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.     

Geodynamically induced variations in the emission of CO2 gas at San Faustino (Central Apennines,Italy)

The Central Apennines are affected by frequent earthquakes of moderate magnitude that occur mainly within the upper part of the crust at depths of <15 km. A large number of cold gas emissions that are rich in CO₂ are also found in the region. One particular vent with a high rate of degassing was equipped with a sensor to measure flow rates, which were recorded for a number of different periods between 2005 and 2010. Factors that could affect potentially CO₂ flow rates include barometric pressure, atmospheric temperature, precipitation and local seismicity. Our analysis indicates that the periods of anomalous flow rate were related not to the environmental factors but probably to the deformative processes of the crust associated with the local seismicity. Local seismic events as expression of geodynamic processes occurred always before and during these anomalous gas flow periods. This correlation exists only for events that occurred eastwards of the gas emission site close to the Martana fault zone. We herein consider this correlation as indication for a continuous interaction between the field of static strain and the deep fluid pressure. An approximation of the fluid pressure transmission towards the gas emission site gives reasonable values of 1–10 m² sec^?1. To make comparisons with the long‐term effects of the static strain, we also recorded the short‐term effects of the dynamic release of strain induced by the series of strong earthquakes that took place in L’Aquila in 2009. We detected a significant anomalous flow rate that occurred at the same time as this seismic sequence, during which widespread degassing was induced around the focal zone.     

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.     

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.     

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     

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.     

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.     

Spectral Attenuation Relations at Soft Sites Based on Existing Attenuation Relations for Rock Sites

A simple and general method based on well-known random vibration theory is used to compute spectral attenuation relations at soft sites based on existing spectral attenuation relations at rock sites. The method consists of: (1) computation, for given magnitude and distance, of the expected Fourier amplitude spectrum associated with the median rock response spectrum computed with the attenuation relation; (2) inclusion of site effects characterized by a frequency-dependent, linear, or nonlinear transfer function; and (3) computation of the response spectrum at the soft site.     

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.     

EUROCODE 8 DESIGN RESPONSE SPECTRA EVALUATION USING THE K-NET JAPANESE DATABASE

The Eurocode 8 (EC8) currently proposes two standard shapes for the design response spectra. Type 1 spectra are enriched in long period and are suggested for high seismicity regions. Conversely, Type 2 spectra are proposed for low to moderate seismicity areas (like France), and exhibit both a larger amplification at short period, and a much smaller long period contents, with respect to Type 1 spectra. These propositions, however, were constrained using few events mostly recorded on analogical instruments. In the present study, we use the Japanese high quality digital K-net array in order to evaluate the proposed ECS response spectra. Furthermore, all K-net stations have geotechnical characterisation. We first constructed a database of shallow events, depth less than 25 km, to avoid subduction related records. The database spans six years of seismicity from 1996 until 2003. Thus, 591 events were selected with moment magnitude between 4 and 7.3, recorded at 691 stations, giving a total of 6812 two horizontal components accelerograms. Using these records, we computed spectral ground-motion prediction equations and we used them to review the shape of the proposed EC8 spectra. In particular, we studied the plateau-PGA ratio level, the period interval where this plateau is constant, and site amplification effects. The results show surprisingly that the Type 2 rock better envelope the Japanese data. Another interesting observation is that the K-net data corresponding to all soil classes are rich in short periods around 0.1 s. This characteristic has not been observed in other worldwide databases. Normalised empirical predictions show a widening of the plateau as the soil conditions degrade. This suggests that the Type 2 EC8 spectra do not cover enough the long periods for EC3-soil classes C, D and E. Finally, the computed ground-motion prediction equations show that the peak ground acceleration (PGA) is nearly invariant to the soil conditions. Soil effects are mainly seen in the shape and plateau level.     

Seismic Site Response Modeling for Three Missouri River Highway Bridges

Three highway bridges spanning the Missouri River flood plain were selected for evaluation of seismic site response for moderate size earthquakes emanating from the New Madrid Seismic Zone (NMSZ) in the Midwestern United States. The NMSZ is known to be capable spawning earthquakes larger than magnitude (M) 7.0, four of which occurred in a three-month period between 1811 and 1812, and the M_w 6.0 earthquake of October 1895 centered near Charleston, Missouri. This study evaluated the likely impacts of long period motion of these historic earthquakes on three long-span highway bridges using geotechnical data obtained from recent investigations. Our results suggest site amplification between 6× and 9×, depending on the magnitude and epicentral distance. We believe that threshold magnitude for serious foundation failure and damage to these bridges is between M_w 6.5 and 6.6. Above these magnitudes widespread liquefaction is predicted, which would effect the peak horizontal acceleration and spectral accelerations, causing the ground motions to be different than predicted. Increase in amplification of the response spectra also should be expected where the periods are higher than 1.0 sec. Therefore, M_w 6.5+ earthquakes at ranges 210–260 km could be expected to engender resonant frequency problems for multiple span bridges and tall buildings (10 to 25 stories) in channel corridors containing 20 to 46 m of unconsolidated sediment.     

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.     

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.     

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.     

Site Classification Assessment for Estimating Empirical Attenuation Relationships for Central-Northern Italy Earthquakes

The aim of this article is to investigate the ground motion attenuation of the most industrialized and populated regions of Italy, evaluating the capability of different approaches to estimate site dependent models. The 5.2 local magnitude earthquake on November 24, 2004 shocked the areas of Northern Italy producing damage of about 215 million euros. The data set, including 243 earthquakes of local magnitude up to 5.2, has been collected in the period December 2002–October 2005 by 30 three-component seismic stations managed by Istituto Nazionale di Geofisica e Vulcanologia, Sezione di Milano (INGV-MI). Empirical attenuation relationships have been estimated for horizontal peak ground velocity (PGHV), acceleration (PGHA), displacement (PGHD), and for response spectral acceleration (SA) for periods between 0.1 and 1.5 s. To estimate suitable attenuation models, in particular for sites characterized by thick sedimentary geological formations, a soil discrimination based on EU8 code can lead to wrong evaluations. On the contrary, a classification based on H/V spectral ratios of seismic ambient noise (NHV) allows the models to fit better real and predicted data and to reduce the uncertainties of the process. For each receiver, NHV have been strengthened by additional H/V spectral ratio of earthquake data (EHV), calculated considering different portions of the analysed signals. In order to validate the PGHA attenuation relationship for greater magnitudes, accelerometric records, relative to Central-Northern Italy strong motions occurring in the last 30 years, have been collected and superimposed to our attenuation curves.     

Site Effects at Long Periods from Digital Strong Motion Records of the KiK-net,Japan

The high-quality digital records of the Japanese KiK-net were examined, with the aim of studying the influence of local site conditions on the displacement spectral ordinates at long periods. The attention was limited to those records for which the velocity profiles up to 100–200 m depth were known, and corresponding surface and borehole accelerograms were available. Based on the available records and with the support of 1D numerical simulations, different aspects that may have an influence on the amplification of long period spectral ordinates were studied, including the bedrock velocity profile, the site classification using V_s,30, and the earthquake magnitude.

Small amplification factors at long periods were found, ranging from 1 to 1.3, with median value from 1.05 to 1.1, for Eurocode 8 site classes B and C, respectively. Only for two records on soft soils (at the boundary between Eurocode 8 classes C and D), from small magnitude earthquakes, large amplification factors were obtained, up to about 4. A good correlation was found of the amplification levels with the response spectral ratio D(T₀)/D(10), where D(T₀) and D(10) are displacement spectral ordinates of the input signal at bedrock, at the fundamental period T₀ of the soil profile and at T = 10 s, respectively.     

Effects of Soil Characterization on the Seismic Input

This paper is aimed at determining the effects of the soil characterization on the seismic input to use for seismic assessment. Three different soil profiles have been assembled to represent the stratigraphies found through a proper experimental investigation, carefully described, and alternative seismic site response analyses have been performed. The surface spectra obtained from the seismic site response analysis (SRA) are very different from each other, thus evidencing the importance of carefully describing soil stratigraphy. Furthermore, the comparison among the surface records found for different return periods has shown a limited sensitivity of the SRA to the seismic intensity.     

Tsunamigenic Analysis in and around Makran

The Makran region is only segment east of the Mediterranean and west of Andaman arc in which subduction of oceanic lithosphere is still an ongoing process. We applied statistical models for the possible explanation of earthquakes in and around Makran, and described a method of Fault Plane Solutions to present seismicity picture in the investigated region, and their possible tsunamigenic effects on Gwadar Coast. We also applied the Community Modeling Interface for Tsunami (ComMIT) for the estimation of tsunami hazard along the coast of Gwader due to the Arabian and Eurasian plate boundary earthquakes along the Makran subduction zone.