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 ASSESSMENT OF UNITED ARAB EMIRATES AND ITS SURROUNDINGS

This paper presents the seismic hazard assessment and seismic zoning of the United Arab Emirates (UAE) and its surroundings based on the probabilistic approach. The area that has been studied lies between 50°E-60°E and 20°N-30°N and spans several Gulf countries. First, the tectonics of the area and its surroundings is reviewed. An updated catalogue, containing both historical and instrumental events is used. Seismic source regions are modelled and relationships between earthquake magnitude and earthquake frequency is established. A modified attenuation relation for Zagros region is adopted. Seismic hazard assessment is then carried out for 20 km interval grid points. Seismic hazard maps of the studied area based on probable Peak Ground Acceleration (PGA) for 10% probability of exceedance for time-spans of 50, 100 and 200 years are shown. A seismic zone map is also shown for a 475-year return period. Although the results of the seismic hazard assessment indicated that UAE has moderate to low seismic hazard levels, nevertheless high seismic activities in the northern part of UAE warrant attention. The northern Emirates region is the most seismically active part of UAE. The PGA on bedrock in this region ranges between 0.22 g for a return period of 475 years to 0.38 g for a return period of 1900 years. This magnitude of PGA, together with amplification from local site effect, can cause structural damage to key structures and lifeline systems.     

TUNING OF SEISMIC HAZARD ESTIMATES BY MEANS OF OBSERVED SITE INTENSITIES

Probabilistic seismic hazard is usually assessed by means of computer codes using seis-mogenic sources, parametric catalogues, seismicity rates and attenuation relationships. All these ingredients are conditioned by expert judgement that influences the final results. Even the attenuation relationships, though strictly based on experimental data, are considered a weak point due to the difficulty of modelling the interaction between seismic energy radiation and site response, and because earthquakes do not usually repeat themselves according to one theoretical model. Recently, methods making wide use of site intensity data have been developed in regions such as Italy, where the observed seismic history at selected sites is quite exhaustive. We analyzed these observations to assess seismic hazard at about 600 sites. We used a probabilistic counting technique, integrating the observations (when necessary), with computed shakings obtained from a logistic-type attenuation model. The results were then compared with the estimates provided by the recent seismic hazard map of Italy, compiled according to the traditional probabilistic seismic hazard approach. The match shows significant differences for some sites. A tentative explanation which seems to point to three alternatives is provided: (1) The mismatch between the two methodologies might appear because the stationary assumption has a poor fit with reality (at least in certain areas); (2) Some sites show a response that is systematically different from the average values predicted using attenuation relationships; (3) The definition of seismogenic zones leads to a bias in the seismicity rate estimate.     

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.     

EARTHQUAKE DAMAGE SCENARIOS AND SEISMIC HAZARD OF MESSINA,NORTH-EASTERN SICILY (ITALY) AS INFERRED FROM HISTORICAL DATA

A study aimed at evaluating earthquake damage scenarios and seismic hazard of Messina using historical data, is presented. The analysis of coeval reports allowed us to reconstruct the seismic history of the city and to obtain a homogeneous earthquake site catalogue based on intensity assessed by the European Macroseismic Scale 1998. In the last 1200 years Messina was destroyed once (1908, intensity X-XI EMS) and suffered effects estimated between intensities VII and IX EMS many times (e.g. 853, 1169, 1494, 1509, 1599, 1693, 1783, 1894, 1909). Destruction or severe damage which affected the city are mainly related to earthquakes occurring in the Messina Straits and Southern Calabria, while slighter, moderate effects are usually due to shocks taking place in the seismogenic sources of SE Sicily, Gulf of Patti and Northern Calabria. The damage scenarios of the most relevant events, delineated using coeval urban plans of the city, showed that damage distribution is strongly conditioned by the different soil response. A probabilistic seismic hazard assessment was obtained by using site observed intensities: The expected intensity in a time span of 50 years (i.e. maximum intensity characterised by at least 10% exceedance probability in 50 years) is DC EMS; the expected intensity in a time span of 300 years (10% exceedance probability in 300 years) is X EMS.     

HINTS ABOUT SITE AMPLIFICATION EFFECTS COMPARING MACROSEISMIC HAZARD ESTIMATE WITH MICROTREMOR MEASUREMENTS: THE AGRI VALLEY (ITALY) EXAMPLE

We approach from a new standpoint the problem of estimating seismic hazard for some towns and villages located in Val d'Agri area (Southern Italy) that in the past have been affected by several seismic events. The estimates are carried out using a method that is based on the analysis of site seismic history extracted from macroseismic catalogues. To study the influence of site effects two different procedures have been performed: in the first, seismic hazard estimates have been deduced from epicentral data only, in the second, intensity data actually observed at the site are also considered. The difference between the two estimates can be correlated with local variations of seismic response due to local geological features which are responsible for possible cases of amplification. In order to validate the presence of such correlation, seismic hazard estimates have been compared to site amplification measurements obtained by using the HVSR (Horizontal to Vertical Spectral Ratio) technique. Our findings reveal a good correlation between seismic hazard enhancements and the presence of site amplification effects. The application of this kind of analysis to the Val d'Agri area has pointed out that the joint estimates of site seismic hazard enhancement and HVSR measurements could be a helpful tool to identify problems related to seismic microzonation.     

SENSITIVITY OF PARAMETERS FOR PROBABILISTIC SEISMIC HAZARD ANALYSIS USING A LOGIC TREE APPROACH

The sensitivity of different parameters used in probabilistic seismic hazard calculation is investigated by different logic tree runs with alternative magnitude sets, source zone models and attenuation relations, and with different sets of values for the seismicity parameters and the <r-value. Also the influence from the different parameters on the hazard uncertainty, represented by fractiles, is investigated. The calculations are made for peak ground acceleration at a site near Aachen in the Lower Rhine Embayment. The model where the site is located in a larger source zone gives lower hazard values. This is typical for the case where the seismicity near the site is high relative to its surrounding. The hazard curves for the different attenuation functions are similar, an effect of the similarity of the functions themselves. However, a large sensitivity of this parameter is indicated for small mean return periods. An increased α-value implies a moderate increase of the hazard at long mean return periods. The hazard is increasing for decreased focal depth, decreased β-value and increased maximum expected magnitude, respectively. However, the effects are noteworthy only at low hazard levels for variations in the focal depth and to some extent in the maximum expected magnitude. Finally, decreasing the minimum magnitude thought to be of engineering relevance causes a drastic increase of the hazard at small mean return periods.     

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.     

PSHA for Broad-Band Strong Ground-Motion Hazards in the Saronikos Gulf,Greece, from Potential Earthquake with Synthetic Ground Motions

We perform a probabilistic seismic hazard analysis (PSHA) for broad-band strong ground motion within the Saronikos Gulf region, Greece, from potential earthquakes along the 30 km long Aigina fault in the northern part of the Gulf. We perform the PSHA utilizing empirical Green’s functions (EGFs) merged with synthetic Green’s functions (SGFs) along with models of finite rupture in place of standard “attenuation relations.” Our approach considers all significant magnitudes for PSHA and full broadband ground motion simulations. Calculations are source and site specific, and could reduce uncertainties in estimating standard engineering parameters. We use a range of rupture scenarios for all significant magnitude earthquakes along the fault. The hazard calculation is for frequencies 0.0– 15.0 Hz. Recordings of small earthquakes from an onshore/offshore local seismic array were used as EGFs for frequencies of 1.5–15.0 Hz, the finite difference code E3D was utilized to synthesize SGFs for frequencies 0.0–1.5 Hz, and an algorithm for merging the EGFs with SGFs was developed. The full-waveform calculations are important for non-linear dynamic analysis of structures in the coastal zone and potential hazard to long period structures. Results of proposed PSHA identify 2%, 10%, and 50% hazard at the selected sites of Saronikos Gulf.

Finally, we compare our PSHA results to those obtained by standard practice which involves prediction equations (GMPEs) recently developed in the Next Generation Attenuation (NGA) project and empirical predictive attenuation relations proposed for Greece. We believe that differences with the NGA results are due to site- and source-specific information utilized in this study, and incorporation of this information may significantly reduce the uncertainty in seismic hazard calculations.     

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.     

Predicting Probabilistic-Based Strong Ground Motion Time Series for Citadel of Arg-E-Bam (South-East of Iran)

The main objective of this article is to present a probabilistic-based strong motion compatible with the source-path and site soil condition given the probability of exceedence for citadel of Arg-e-Bam site bed rock (South-East of Iran). A Fourier amplitude spectral attenuation relation for bed rock beneath the site is proposed which permits the estimation of time-histories through a probabilistic seismic hazard analysis procedure. Due to lack of data, the two well-known simulation techniques, point source and finite fault models have been used for generating hundreds of strong motion as input data. Tens of model parameter values such as stress-drop nucleation points were used, in each specified magnitude-distance, to reduce the uncertainty effects inherently existing in seismological/geological parameters. The proposed attenuation relation is validated by comparing the estimated strong motion, in the form of Fourier amplitude spectral, using the proposed attenuation relation with those of recorded ground motion data at three stations far away from the assumed source so that the results would not be influenced by the near source problems such as directivity and fling step. The results of proposed technique is assessed by comparing the estimated response spectra, with 10% probability of exceedence and 5% damping ratio, with those of traditional uniform hazard spectra. The proposed technique is supposed to be used in retrofitting procedure of international historical adobe structures in Arg-e-Bam site, which have been damaged during the destructive Bam earthquake 2003, Iran     

MODELLING OF THE GROUND MOTION AT RUSSE SITE (NE BULGARIA) DUE TO THE VRANCEA EARTHQUAKES

An approach, capable of synthesising strong ground motion from a basic understanding of fault mechanism and of seismic wave propagation in the Earth, is applied to model the seismic input at a set of 25 sites along a chosen profile at Russe, NE Bulgaria, due to two intermediate-depth Vrancea events (August 30, 1986, M _ω=7.2, and May 30, 1990, M _ω=6.9). Accordingly to our results, once a strong ground motion parameter has been selected to characterise the ground motion, it is necessary to investigate the relationships between its values and the features of the earthquake source, the path to the site and the nature of the site. Therefore, a proper seismic hazard assessment requires an appro-priate parametric study to define the different ground shaking scenarios corresponding to the relevant seismogenic zones affecting the given site. Site response assessment is provided simultaneously in frequency and space domains, and thus the applied procedure differs from the traditional engineering approach that discusses the site as a single point. The applied procedure can be efficiently used to estimate the ground motion for different purposes like microzonation, urban planning, retrofitting or insurance of the built environment.     

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.     

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 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.     

Considering Soil-Nonlinearity Effect on Seismic Performance Evaluation of Structures based on Pushover Analysis

According to the most of current seismic codes, nonlinear soil behavior is commonly ignored in seismic evaluation procedure of the structures. To contribute on this matter, a pushover analysis method incorporating the probabilistic seismic hazard analysis (PSHA) is proposed to evaluate the effect of nonlinear soil response on seismic performance of a structure. The PSHA outcomes considering soil nonlinearity effect is involved in the analysis procedures by modifying the site-specific response spectrum. Results showed that incorporation of nonlinear soil behavior leads to an increase in displacement demand of structures which should accurately be considered in seismic design/assessment procedure. Results of implemented procedure are confirmed with the estimated displacement demand including soil-structure interaction (SSI).     

Dome of the Basilica of Santa Maria Degli Angeli in Assisi: Static and Dynamic Assessment

This study presents the results of the static and dynamic assessment of the dome of the Basilica of Santa Maria degli Angeli in Assisi, designed by Galeazzo Alessi. The first section is devoted to an overview of the masonry domes designed by the Italian architect and focuses on the structural solutions adopted in the several cases described to better understand Alessi’s designing skills. In the second part, the drum-dome system is analyzed in order to attain a structural assessment. The static assessment is performed by means of limit analysis and finite element model approaches with non-linear mechanical behavior. The obtained results are consistent with the detected crack pattern and confirm the suitability of the reinforcement steel rings applied to the drum. The seismic assessment has been performed by response spectrum analysis. Due to the lack of specific information, a probabilistic approach for the material mechanical properties was used. The results obtained highlight an adequate seismic response of the structure that can be attributed to the dynamic properties of the slender drum-dome system. This finding justifies the good performance of the structure during the seismic events of 1832 and 1997.     

Extreme weather events in northeastern New Brunswick (Canada) for the period 1950–2012: Comparison of newspaper archive and weather station data

In the context of global warming, changes in extreme weather and climate events are expected, particularly those associated with changes in temperature and precipitation regimes and those that will affect coastal areas. The main objectives of this study were to establish the number of extreme events that have occurred in northeastern New Brunswick, Canada in recent history, and to determine whether their occurrence has increased. By using archived regional newspapers and data from three meteorological stations in a national network, the frequency of extreme events in the study area was established for the time period 1950–2012. Of the 282 extreme weather events recorded in the newspaper archives, 70% were also identified in the meteorological time series analysis. The discrepancy might be explained by the synergistic effect of co‐occurring non‐extreme events, and increased vulnerability over time, resulting from more people and infrastructure being located in coastal hazard zones. The Mann Kendall and Pettitt statistical tests were used to identify trends and the presence of break points in the weather data time series. Results indicate a statistically significant increase in average temperatures and in the number of extreme events, such as extreme hot days, as well as an increase in total annual and extreme precipitation. A significant decrease in the number of frost‐free days and extreme cold days was also found, in addition to a decline in the number of dry days.     

Critical Issues in Earthquake Scenario Loss Modeling

The assessment of human or economic losses due to single events (scenario) may effectively support decision makers in the development of important risk mitigation actions. The study presented herein sheds light on several problems and limitations in the current practice of scenario loss modeling, such as: the number of simulations required to achieve convergence; epistemic and aleatory uncertainty in the ground motion prediction and vulnerability models; and consideration of the earthquake rupture geometry. These issues are investigated using the Metropolitan Area of Lisbon (Portugal), and it has been observed that distinct assumptions in the loss modeling can lead to considerably different results. The findings of this study are also pertinent for probabilistic seismic risk analyses in which a large number of stochastically generated events are employed to assess probabilistic losses.     

A Ground Motion Record Selection Approach Based on Multiobjective Optimization

The evaluation of the seismic safety and reliability of buildings and building contents within a probabilistic framework often requires response history analyses using site-specific ground motion records. The ground motion selection method proposed in this paper addresses this issue by a stochastic search procedure in which record sets are selected such that first- and second-order statistics (median and dispersion) satisfy predefined ground motion spectrum targets over a wide period range. Once a ground motion record set is selected, it can be used for seismic assessment of a broad class of buildings within the target period range at the given location.