Spatial zonations for regional assessment of seismic site effects in the Seoul metropolitan area

Earthquake-induced hazards are profoundly affected by site effects related to the amplification of ground motions, which are strongly influenced by local geologic conditions such as soil thickness or bedrock depth and soil stiffness. In this study, an integrated geographic information system (GIS)-based system for geotechnical data, called the geotechnical information system (GTIS), was developed to establish a regional counterplan against earthquake ground motions in the Seoul metropolitan area. In particular, to reliably predict spatial geotechnical information, a procedural methodology for building the GTIS within a GIS framework was developed and applied to the Seoul area in Korea. To build the GTIS, pre-existing geotechnical data were collected in and around the study area, and then a walk-over site survey was conducted to acquire surface geo-knowledge data. In addition, the representative shear wave velocities for geotechnical layers were derived by statistically analyzing many seismic test data in Korea. The GTIS was used in a practical application to estimate site effects in the study area; seismic zoning maps of geotechnical earthquake parameters, such as the depth to bedrock and the site period, were created and presented as a regional synthetic strategy for earthquake risk assessment. Furthermore, seismic zonation of site classification was also performed to determine the site amplification coefficients for seismic design and seismic performance evaluation at any site and administrative sub-unit in the study area. The methodology and results of the case study of seismic zonations in the Seoul area verified that the GIS-based GTIS can be very useful for the regional estimation of seismic risk and also to support decisions regarding seismic hazard mitigation, particularly in the metropolitan area.     

Non-uniqueness in surface-wave inversion and consequences on seismic site response analyses

Surface-wave tests are based on the solution of an inverse problem for shear-wave velocity profile identification from the experimentally measured dispersion curve. The main criticisms for these testing methodologies are related to the inverse problem solution and arise from the possible equivalence of different shear-wave velocity profiles. In this paper, some implications of solution non-uniqueness for seismic response studies are investigated using both numerical simulations and experimental data. A Monte Carlo approach for the inversion problem has been used to obtain a set of equivalent shear-wave velocity models. This selection is based on a statistical test which takes into account both data uncertainty and model parameterization. This set of solutions (i.e., soil profiles) is then used to evaluate the seismic response with a conventional one-dimensional analysis. It is shown that equivalent profiles with respect to surface-wave testing are equivalent also with respect to site amplification, thus countering the criticism related to inversion uncertainty for the engineering use of surface-wave tests.     

Site response analyses using downhole arrays at various seismic hazard levels of Singapore

Local site conditions can significantly influence the characteristics of seismic ground motions. In this study, site response analyses using one-dimensional linear elastic (LE), equivalent-linear (EQL) and nonlinear (NL) approaches are performed at different seismic hazard levels of Singapore. Two seismic stations, namely, the KAP and BES stations located at soft soil sites, are selected from the national network of Singapore. Firstly, site response estimates using the LE, EQL (SHAKE04) and NL (DEEPSOIL) approaches are compared with the borehole recordings. Results show favorable matches between the predictions and the observations at the KAP site, while under-predictions are observed for all the three site effect approaches at the BES site. Secondly, the applicability of the LE, EQL and NL models is examined at different hazard levels of Singapore. It is found that for the hazard level at a return period of 475 years, the computed maximum strain (γ_max) is 0.06% and then the EQL model can provide accurate site response predictions. However, for the hazard level at a return period of 2475 years, the calculated γ_max is larger than 2%, resulting in notable differences in the predictions of different site response models. This study highlights the importance of site effects in seismic hazard analysis of Singapore.     

An exploration of incorporating site response into PSHA—Part I: Issues related to site response analysis methods

This paper presents a series of analyses for the evaluation of the ground response of two NEHRP class D sites, subjected to shaking by a large number of strong ground-motion records. The two investigated sites have very distinct profiles, but they are characterised by almost identical V_s30 values. The site response analyses are performed using various methods of analysis and input parameters in order to explore the sensitivity of the ground response estimates and to identify the dominating parameters. Equivalent linear analysis is performed using different sets of dynamic soil properties curves, while nonlinear analysis is performed using different target dynamic soil curves, viscous damping formulations and fitting procedures for the constitutive model parameters. Particular focus is given to the sensitivity of the response when soil sites are subjected to high-intensity shaking, a subject of particular interest when the prediction of surface ground motions with low annual probabilities of exceedance is the target of probabilistic seismic hazard analyses (PSHA). The site response analysis results of this paper are incorporated into the probabilistic framework of Bazzurro and Cornell [1] in our companion paper in order to assess their impact on the final soil surface hazard calculation.     

Evaluation of seismic site factors in the Mississippi Embayment. I. Estimation of dynamic properties

Deep deposits of the Mississippi Embayment, overlying the New Madrid Seismic Zone, present unique challenges for the estimation of local site effects on propagated ground motion. Limited information is available on the dynamic properties of these deposits. This paper develops generalized depth dependent modulus degradation and damping curves specific to the embayment deposits. Depth dependent small strain damping is estimated using weak motion recordings during the Enola earthquake, 2001. Modulus degradation and damping curves are further constrained using limited laboratory test data of embayment soils at low confining pressure. At high confining pressures embayment specific data is unavailable and test data from outside the embayment is used in developing the dynamic properties. The representative modulus degradation and damping curves are used in non-linear and equivalent linear one-dimensional site response analyses. A companion paper describes a large-scale probabilistic seismic hazard analysis study in the Mississippi Embayment that integrates non-linear site effects.     

From regional- to local-scale seismic hazard assessment: examples from Southern Spain

An approach that relates results from a regional seismic hazard assessment study with local-scale site-effect characterizations in an area of low-to-moderate seismic activity such as Andalusia (southern Spain), is presented. Results of a previous probabilistic seismic hazard analysis of Andalusia on rock conditions are disaggregated to infer hazard controlling earthquakes for different target motions. A collection of controlling magnitude-distance pairs and the corresponding site-specific response spectra at main capital cities of the region are obtained. These spectra are first-order approximations to expected seismic actions required in local earthquake risk assessments. In addition, results of independent, local-scale studies developed in Almeria City (SE Andalusia) are used to derive an updated seismic zonation of the city. These include predominant soil period estimates and shear-wave velocity profiles at different locations. If a local seismic risk assessment study or an earthquake-resistant structural design is to be developed, it may be recommended the use of different seismic actions on sites characterized by distinct response to seismic shaking (as derived from the seismic zonation). The seismic action related to worst-case scenarios may be modeled through a hazard-consistent response spectrum, obtained by hazard disaggregation at the spectral acceleration with period matching the prevailing resonant period of the target site or structure.     

Evaluation of site conditions for the Ankara Basin of Turkey based on seismic site characterization of near-surface geologic materials

Regional site conditions relevant for seismic hazard studies can be derived from various geologic, seismologic, and geotechnical sources. In this study, site conditions are derived for the Ankara Basin in Turkey by merging in-situ seismic measurements of dynamic properties, geologic information, and some geotechnical boring information. Field seismic refraction surveys were performed at 259 sites in the project area to classify and characterize Plio–Pleistocene fluvial deposits and Quaternary alluvial and terrace deposits. The shear-wave velocity profiles of the near-surface geologic units are used to characterize site classes according to the International Building Code [International Code Council, ICC. International Building Code. Structural and fire- and life-safety provisions (seismic, wind, accessibility, egress, occupancy and roof codes), 2006. Whittier, CA.] and the Turkish Seismic Code [Ministry of Public Works and Settlement, 1998. Turkish Seismic Code, Specification for Structures to be Built in Disaster Areas, Ankara, Turkey], and to develop a regional model for the average shear-wave velocity in the top 30 m. The resulting maps of site class indicate that the classification system for the Turkish Seismic Code results in a significant portion of the Ankara Basin being classified as Z4, the softest site class. The International Building Code site classification system results in most of the Ankara Basin being classified as D, stiff soil. These differences are caused by the Turkish Seismic Code incorporating information from only the surface layer, while the International Building Code incorporates information from the top 30 m.     

Seismic vulnerability estimation of the building considering seismic environment and local site condition

Introduction The estimation of damage probability distribution among different damage states of rein-forced concrete buildings is a key component of earthquake loss estimation for modern city or a group of cities. With the development of city, the reinforced concrete buildings are major compo-nent parts of modern cities. Vulnerability estimates for these kinds of buildings are of importance to those responsible for civil protection, relief, and emergency services to enable adequate contin-genc…     

Introduction of conditional mean spectrum and conditional spectrum in the practice of seismic safety evaluation in China

This paper aims at implementing and introducing the use of conditional mean spectrum (CMS) and conditional spectrum (CS) as the main input parameters in the practice of seismic safety evaluation (SSE) in China, instead of the currently used uniform hazard spectrum (UHS). For this purpose, a procedure for M-R-epsilon seismic hazard deaggregation in China was first developed. For illustration purposes, two different typical sites in China, with one to two dominant seismic zones, were considered as examples to carry out seismic hazard deaggregation and illustrate the construction of CMS/CS. Two types of correlation coefficients were used to generate CMS and the results were compared over a vibration period range of interest. Ground motion records were selected from the NSMONS (2007–2015) and PEER NGA-West2 databases to correspond to the target CMS and CS. Hazard consistency of the spectral accelerations of the selected ground motion records was evaluated and validated by computing the annual exceedance probability rate of the response spectra and comparing the results to the hazard curve corresponding to each site of concern at different periods. The tools developed in this work and their illustrative application to specific case studies in China are a first step towards the adoption of CMS and CS into the practice of seismic safety evaluation in this country.     

Probabilistic seismic hazard assessment for South-Eastern France

The accurate evaluation and appropriate treatment of uncertainties is of primary importance in modern probabilistic seismic hazard assessment (PSHA). One of the objectives of the SIGMA project was to establish a framework to improve knowledge and data on two target regions characterized by low-to-moderate seismic activity. In this paper, for South-Eastern France, we present the final PSHA performed within the SIGMA project. A new earthquake catalogue for France covering instrumental and historical periods was used for the calculation of the magnitude-frequency distributions. The hazard model incorporates area sources, smoothed seismicity and a 3D faults model. A set of recently developed ground motion prediction equations (GMPEs) from global and regional data, evaluated as adequately representing the ground motion characteristics in the region, was used to calculate the hazard. The magnitude-frequency distributions, maximum magnitude, faults slip rate and style-of-faulting are considered as additional source of epistemic uncertainties. The hazard results for generic rock condition (Vs30 = 800 m/s) are displayed for 20 sites in terms of uniform hazard spectra at two return periods (475 years and 10,000 years). The contributions of the epistemic uncertainties in the ground motion characterizations and in the seismic source characterization to the total hazard uncertainties are analyzed. Finally, we compare the results with existing models developed at national scale in the framework of the first generation of models supporting the Eurocode 8 enforcement, (MEDD 2002 and AFPS06) and at the European scale (within the SHARE project), highlighting significant discrepancies at short return periods.     

Validation through HVSR measurements of a method for the quick detection of site amplification effects from intensity data: an application to a seismic area in Northern Italy

An instrumental validation is attempted of an innovative approach devoted to the quick individuation, from macroseismic data, of site amplification phenomena able to significantly modify seismic hazard levels expected on the basis of average propagation effects only. According to this methodology, two evaluations of hazard are performed at each investigated locality: the former, obtained by epicentral intensity data ‘reduced’ at the site through a probabilistic attenuation function and, the latter, computed by integrating such data with seismic effects actually observed at the site during past earthquakes. The comparison, for each locality, between these two hazard estimates allow to orientate the identification of those sites where local amplifications of earthquake ground motion could be significant. In order to check such methodology, indications obtained in this way from macroseismic data are compared with the estimates of transfer functions performed through the HVSR technique applied to microtremors. Results concerning municipalities located in a seismic area of Northern Italy indicate a good agreement between macroseismic and instrumental estimates.     

Sensitivity of site response analysis on the number of ground motion records and implications for PSHA

This paper investigates issues related to the number of ground-motion records required for the performance of site response analysis and the inclusion of the site-specific amplification function within probabilistic seismic hazard calculations (PSHA). It explores the minimum number of records required for a robust estimation of the median and standard deviation of the site amplification function, as well as the impact of the selected ground-motion suites on the results of PSHA. Site response analyses are performed using both equivalent linear and nonlinear methodologies. Although the median amplification was observed to be relatively easy to capture, the standard deviation was seen to fluctuate considerably, especially when suites of few records were used. It was observed that in the case of the nonlinear site response analysis 10 records provide relatively stable estimates of the hazard curves for the majority of periods, while in the case of the equivalent linear analysis 20 records or more are required to achieve a similar level of accuracy.     

The reference spectral noise ratio method to evaluate the seismic response of a site

The seismic response of several sites in the Durban area of South Africa has been investigated using a refined version of the spectral ratio (H/V) generated by cultural seismic noise. Two samples of noise separated by an interval of several minutes for the vertical and the two orthogonal horizontal components were taken at each site. Though the two samples differed noticeably in the time as well as in the frequency domain, when the components are considered separately, the ratios of the horizontal to the vertical spectra were very similar for most of the pairs. This similitude was obvious not only in the frequencies of the peaks, but also in the range of amplifications involved. These observations have been used to introduce a refinement to the Nakamura Method. This refinement consists in using as a reference the H/V spectral ratio of a well-known, hard rock site. All the other soft sites of interest are then compared with this standard reference site. This comparison provides valuable information on the frequencies and amplification factors introduced by the critical soft sediments. This refined method will be referred here as the ‘Reference spectral noise ratio’.This method, is compared to a similar one devised to investigate the response of horizontal layers to longitudinal seismic waves, using the ratio of the spectra generated by long- period waves.     

Dynamic properties of the Quaternary sedimentary rocks and their influence on seismic site effects. Case study in Bucharest City,Romania

Bucharest is one of the cities most affected by earthquakes in Europe. Situated at 150–170 km distance from Vrancea epicentral zone, Bucharest had suffered many damages due to high energy Vrancea intermediate-depth earthquakes. For example, the 4 March 1977 event produced the collapse of 32 buildings with 8–12 levels, while more than 150 old buildings with 6–9 levels were seriously damaged. The studies done after this earthquake had shown the importance of the surface geological structure upon ground motion parameters. New seismic measurements are performed in Bucharest area aiming at defining better elastic and dynamic properties of the shallow sedimentary rocks. Down-hole seismic measurements were performed in a number of 10 cased boreholes drilled in the Bucharest City area. Processing and interpretation of the data lead to the conclusion that shallow sedimentary rocks can be considered weak in the area, down to 150–200 m depth. Seismic wave velocity values and bulk density values presented in the paper associated with local geology are useful primary data in the seismic microzonation of Bucharest City. They are used as 1D models to derive transfer functions and response spectra for the stack of sedimentary rocks in several parts of Bucharest area, leading to a better knowledge of the local site amplification and associated frequency spectra. In a recent study the H/V spectral ratio using Nakamura’s method was applied on the seismic noise measurements in 22 sites in Bucharest City in order to derive the fundamental period associated with these sites. The values confirm the previous results, showing a dominant resonance in the period range of 1.25–1.75 s. The fundamental periods obtained with Nakamura’s method are in good agreement with those computed on the basis of geological and geotechnical data in boreholes, which show an increase of the fundamental period in the Bucharest area from south to north, in the same direction as the increase of the thickness of the Quaternary deposits above the Fratesti layer which is considered the bedrock in the area.     

Effect of site amplification on inelastic seismic response

The available models for eff ective periods of site and structure are reviewed in context of frequency tuning in the inelastic seismic response of soil-structure system. The eff ect of seismic intensity and ductility demand, on the eff ective periods, is investigated, and inelastic site amplifi cation is shown to be strongly correlated to the normalized eff ective period. Two non-dimensional parameters, analogous to the conventional site amplifi cation factors in codes, are defi ned to quantify the inelastic site amplifi cation. It is shown that the inelastic site amplifi cation factor (i.e. ratio of constant ductility spectral ordinates at soil site to those at rock outcrop) is able to represent the site eff ects more clearly, as compared to the inelastic site amplifi cation ratio (i.e. ratio of inelastic spectral ordinates at soil site to the corresponding elastic spectral ordinates at rock outcrop). Further, the peak in the amplifi cation factor corresponding to the eff ective site period diminishes rapidly with increasing ductility demand.     

Regional and local seismic hazard assessment

Recent earthquakes such as the M_JMA 7.2 Hyogo-ken Nambu earthquake and the M 7.4 Kocaeli earthquake demonstrate once again the need to include detailed soil investigation into hazard evaluation, that is the need of microzonation. Seismic hazard assessment evaluated at a regional scale generally does not consider soil effects but only in a limited way using an attenuation law that can be ‘soft soil’ or ‘rock’. However, the relevant role of seismic hazard in the assessment of seismic coefficients for the definition of the actions in seismic codes must be properly considered. That is to say, the level of protection of buildings is proportional to a definite level of hazard (generally considered to be the ground motion with 10% probability of exceedence in 50 years). When a microzonation is performed, this criterion cannot be ignored, therefore, a clear linkage must be established between hazard (regional scale) and microzonation. The crucial point is represented by the reference motion (or input motion) to be used for site effects analysis, that must be compatible with the regional seismic hazard. In this paper, three different approaches for reference motion evaluation are analysed: probabilistic; stochastic; and deterministic. Through the case history of Fabriano microzonation the three approaches are compared. It is shown that each approach presents advantages and disadvantages with respect to the others. For example, the probabilistic approach (the reference motion is directly derived from the expected response spectra for a given return period) is linked with hazard, but produces an overestimation in short periods range, while the deterministic approach correctly simulates the wave propagation, but it ends with a kind of conditional probability. Until now, clear criteria to choose the right approach do not appear to exist and the expert experience is of fundamental importance.     

Shear wave velocities of Mississippi embayment soils from low frequency surface wave measurements

Deep unconsolidated sediments in the Mississippi embayment will influence ground motions from earthquakes in the New Madrid seismic zone. Shear wave velocity profiles of these sediments are important input parameters for modeling wave propagation and site response in this region. Low-frequency, active-source surface wave velocity measurements were performed to develop small-strain shear wave velocity (V_S) profiles at eleven deep soil sites in the Mississippi embayment, from north of New Madrid, Missouri to Memphis, Tennessee. A servo-hydraulic, low-frequency source was used to excite surface wave energy to wavelengths of 600 m, resulting in V_S profiles to depths of over 200 m. The average V_S profile calculated from the eleven sites is in good agreement with common reference V_S profiles that have been used in seismic hazard studies of this region. The variability in V_S profiles is shown to be associated with changes in formation depth and thickness from site-to-site. Using lithologic information at each site, average formation velocities were developed and compared to previous studies. We found average V_S values of about 193 m/s for alluvial deposits, 400 m/s for the Upper Claiborne formations, and 685 m/s for the Memphis Sand formation.     

Conditional spectrum‐based ground motion selection. Part I: Hazard consistency for risk‐based assessments

The conditional spectrum (CS, with mean and variability) is a target response spectrum that links nonlinear dynamic analysis back to probabilistic seismic hazard analysis for ground motion selection. The CS is computed on the basis of a specified conditioning period, whereas structures under consideration may be sensitive to response spectral amplitudes at multiple periods of excitation. Questions remain regarding the appropriate choice of conditioning period when utilizing the CS as the target spectrum. This paper focuses on risk‐based assessments, which estimate the annual rate of exceeding a specified structural response amplitude. Seismic hazard analysis, ground motion selection, and nonlinear dynamic analysis are performed, using the conditional spectra with varying conditioning periods, to assess the performance of a 20‐story reinforced concrete frame structure. It is shown here that risk‐based assessments are relatively insensitive to the choice of conditioning period when the ground motions are carefully selected to ensure hazard consistency. This observed insensitivity to the conditioning period comes from the fact that, when CS‐based ground motion selection is used, the distributions of response spectra of the selected ground motions are consistent with the site ground motion hazard curves at all relevant periods; this consistency with the site hazard curves is independent of the conditioning period. The importance of an exact CS (which incorporates multiple causal earthquakes and ground motion prediction models) to achieve the appropriate spectral variability at periods away from the conditioning period is also highlighted. The findings of this paper are expected theoretically but have not been empirically demonstrated previously. Copyright © 2013 John Wiley & Sons, Ltd.     

Site Response in Las Vegas Valley, Nevada from NTS Explosions and Earthquake Data

We report site response in Las Vegas Valley (LVV) from historical recordings of Nevada Test Site (NTS) nuclear explosions and earthquake recordings from permanent and temporary seismic stations. Our data set significantly improves the spatial coverage of LVV over previous studies, especially in the northern, deeper parts of the basin. Site response at stations in LVV was measured for frequencies in the range 0.2–5.0 Hz using Standard Spectral Ratios (SSR) and Horizontal-Vertical Spectral Ratios (HVR). For the SSR measurements we used a reference site (approximately NEHRP B ``rock' classification) located on Frenchman Mountain outside the basin. Site response at sedimentary sites is variable in LVV with average amplifications approaching a factor of 10 at some frequencies. We observed peaks in the site response curves at frequencies clustered near 0.6, 1.2 and 2.0 Hz, with some sites showing additional lower amplitude peaks at higher frequencies. The spatial pattern of site response is strongly correlated with the reported depth to basement for frequencies between 0.2 and 3.0 Hz, although the frequency of peak amplification does not show a similar correlation. For a few sites where we have geotechnical shear velocities, the amplification shows a correlation with the average upper 30-meter shear velocities, V₃₀. We performed two-dimensional finite difference simulations and reproduced the observed peak site amplifications at 0.6 and 1.2 Hz with a low velocity near-surface layer with shear velocities 600–750 m/s and a thickness of 100–200 m. These modeling results indicate that the amplitude and frequencies of site response peaks in LVV are strongly controlled by shallow velocity structure.     

Probabilistic seismic hazard assessment in SE-Spain based on macroseismic site histories

In many countries such as Spain earthquake databases still mainly comprise macroseismic data from felt effects. The full exploit of this information is of basic importance for seismic risk assessment and emergency planning, given the strict link between macroseismic intensity and damage. A probabilistic procedure specifically developed to handle macroseismic data, mostly relying on site information and seismogenic-source free, has been applied to evaluate seismic hazard in SE-Spain (Alicante-Murcia region). Present seismicity is moderate-low with largest magnitudes slightly over Mw5.0. The historical record includes very destructive earthquakes, maximum EMS98 intensities reaching IX–X and X in the nineteenth century (e.g., Torrevieja 1829 earthquake). Very recently, two events in the area on 11 May 2011 (Mw4.5, Mw5.2) killed nine people, injured 300, and produced important damage in the city of Lorca. Regional hazard maps for the area together with specific hazard curves at selected localities are obtained. Results are compared with the maximum observed intensities in the period 1300–2012, and with the values in the seismic hazard map from the Spanish Building Code in force. In general, the maximum felt intensity values are closer to the hazard values calculated for 2 % probability of exceedance in 50 years, using felt and expected intensity. The intensity-based probabilistic hazard maps obtained through the applied approach reduce the inherent smoothing of those based on standard probabilistic seismic hazard assessment approaches for the region, allowing identifying possible over- or sub-estimates of site hazard values, providing very valuable information for risk reduction strategies or for future updates of the building code hazard maps.