Site Response Models Research Articles

Site effect assessment using microtremor measurement, equivalent linear method, and artificial neural network (case study: Babol, Iran)

Site effect assessment is an important procedure for a reliable site-specific hazard assessment. Fundamental frequency is a very important factor that must be considered in a construction site for examining the potential damage resulting from earthquake. In the two last decades, the microtremor H/V spectral ratio method has been widely used for site effect studies. Microtremor measurement is fast, applicable, and cost-effective. Microtremor measurement was undertaken at 60 stations in the Babol, north of Iran, during 2011 and 2012. Regarding Nakamura’s method, H/V spectral ratios, fundamental frequency, and amplification factor were calculated for all microtremor stations. The obtained results were controlled with SESAME guidelines. In order to assess the accuracy of the microtremor measurements and its application in site effect evaluation, a preliminary site response modeling was carried out using the equivalent linear methods at some stations. More than 18 boreholes in the study area were modeled as multilayer column, overlaid on bedrock. Also, the artificial neural networks (ANN) with different inputs including thickness, type of soil, unit weight, shear wave velocity, and max shear modulus of soil layer were trained. Then, the trained ANN predicts the fundamental frequency, and the output results were compared with microtremor measurement. The results showed that microtremor measurement approach provides an acceptable means of site effect evaluation in the study area.

VS30 and Spectral Response from Collocated Shallow, Active-, and Passive-Source VS Data at 27 Sites in Puerto Rico

Shear-wave velocity (VS) and time-averaged shear-wavevelocity to 30 m depth (VS30) are the key parameters used in seismic site response modeling and earth- quake engineering design. Where VS data are limited, available data are often used to develop and refine map-based proxy models of VS30 for predicting ground-motion intensities. In this paper, we present shallow VS data from 27 sites in Puerto Rico. These data were acquired using a multimethod acquisition approach consisting of noninvasive, collocated, active-source body-wave (refraction/reflection), active- source surface wave at nine sites, and passive-source surface-wave refraction micro- tremor (ReMi) techniques. VS-versus-depth models are constructed and used to calculate spectral response plots for each site. Factors affecting method reliability are analyzed with respect to site-specific differences in bedrock VS and spectral response. At many but not all sites, body- and surface-wave methods generally determine similar depths to bedrock, and it is the difference in bedrock VS that influences site ampli- fication. The predicted resonant frequencies for the majority of the sites are observed to be within a relatively narrow bandwidth of 1-3.5 Hz. For a first-order comparison of peak frequency position, predictive spectral response plots from eight sites are plot- ted along with seismograph instrument spectra derived from the time series of the 16 May 2010 Puerto Rico earthquake. We show how a multimethod acquisition approach using collocated arrays compliments and corroborates VS results, thus adding confi- dence that reliable site characterization information has been obtained.

Lateral Response of Buried Pipeline through Three-Soil Site to Random Ground Motion

Using the pseudo-excitation method, a simple model was built to study the seismic response of buried pipeline laid through three-soil site for vertical incidence of seismic waves. Seismic ground motions are considered as stationary stochastic processes, and the cross-PSD function was established based on auto-PSD model, coherence model, and site response model. The results show that the peak bending stress depends on the properties of adjacent soils and the width of middle-soil, and appears in the pipeline near the soil interfaces.

Earthquake site effect modeling in the Granada basin using a 3-D indirect boundary element method

In this work, we develop the indirect boundary element method (IBEM) to simulate the seismic site response in a realistic, large-scale 3-D sedimentary basin. Most previous applications of boundary element method have used full-space Green’s functions for wave propagation between element points. We use half-space Green’s functions, which include the seismic wavefield interactions at the free surface and require only the boundary elements of the basin interface. In this way, the size of the matrix equation for solution in the IBEM can be reduced to approximately a quarter of that using full-space Green’s functions. The site response modeling of the Granada basin in southern Spain using the IBEM shows that the basin-induced scattering waves were identified as propagating back and forth inside the basin. The scattered waves also generate surface waves that are weakly propagated outside of the basin. The wave propagation inside and outside of the basin shows different patterns. We observe that the scattered wave is locally amplified, and its propagation direction deviates from that of the incident waves propagation direction. Therefore, the computed seismic response in the basin could provide us with good estimates of the seismic motion.

Site Effects in Structural Response Predictions of Inelastic SDOF Oscillators

We study how the inelastic structural response predicted via synthetic seismograms is affected by the selection of site response models in ground motion simulations. We first generate synthetics for multiple scenarios and site conditions in Southern California using attenuation relations, site specific linear, vis-coelastic and nonlinear analyses, and estimate the ground motion variability that results from the soil model selection. We next use bilinear single degree-of-freedom oscillators to demonstrate how this variability propagates to the inelastic structural response predictions. Results show high bias and scatter of the inelastic displacement ratio predicted using the empirical and linear elastic site response models relative to the nonlinear, for periods close to the fundamental period of the site. For the synthetic motions and sites used, we derive empirical correlations between the amount of bias and period range where it manifests, and selected input motion and site parameters.

NEES integrated seismic risk assessment framework (NISRAF)

The paper presents an integrated and extensible framework for assessment of the impact of earthquakes on civil infrastructure systems, particularly buildings and bridges. The framework, referred to as NEES Integrated Seismic Risk Assessment Framework (NISRAF), is developed with a focus on improving the reliability of earthquake assessment results. The components are structural fragility assessment using measured data and hybrid simulation, hazard characterization by free-field site response analysis, and integrated impact assessment. The hazard and fragility components are refined by employing nonlinear site response tools and model updating techniques, respectively. Several of these NISRAF components are tailored to achieve seamless integration and to arrive at an operational system. The novelty of the developed framework is primarily the integration of the various components of earthquake impact assessment, which have not been deployed in such an application before. The framework has been built and demonstrated via applications to a test bed in California. Earthquake impact assessment results using the generated hazard map and fragility curves correlate well with field reports, indicating the efficacy of the approach.

An exploration of incorporating site response into PSHA-part II: Sensitivity of hazard estimates to site response approaches

This paper examines the sensitivity of seismic hazard analyses to various site response analysis procedures. Site effects are incorporated in the hazard calculations using a probabilistic approach and specifically the methodology of Bazzurro and Cornell [1] for the transformation of a generic ground-motion prediction equation to a site-specific one. The paper explores the sensitivity of the median amplification function, its standard deviation and the resulting surface hazard curve, to different methods of site response analysis and model input parameters. The computed site-specific surface hazard curves are also compared with those obtained from a generic soil ground-motion prediction equation. For the two sites investigated, it is shown that the choice of equivalent linear or nonlinear analysis with different constitutive model parameters has a significant impact on the hazard results. The sandy site was seen to be more sensitive to the site response analysis approach employed than the clayey site.

A Study of Local Site Effects on Seismic Design of Buildings

In present work, first the multi-storeyed building was analysing taking medium site condition in zone II. Sa/g and Base shear for that building was calculated from IS 1893:2002 (Part 1).After that site specific ground response analysis is carried out for Amravati and Nagpur region soil sites. For development of site specific response spectrum ProSHAKE software is used ProSHAKE is based on one dimensional geotechnical site response model. Site response study is usually carried out using the geotechnical data or Geophysical data (shear wave velocity). Basic aim of present work is to study local site effects on seismic analysis of multi-storeyed Buildings, for this purpose local soil conditions are considered. In the present work various local sites in Amravati and Nagpur region are studied.

Response of underground pipeline through fault fracture zone to random ground motion

It is assumed that a pipeline is laid through a vertical fault fracture zone, and is excited by seismic ground motion modelled as stationary stochastic process. For horizontal incidence of waves, the cross-PSD (Power Spectral Density) function is developed using wave propagation theory, while for vertical incidence of waves the cross-PSD function is composed by auto-PSD model, coherence model and site response model. As the seismic input, the cross-PSD function is used to calculate the the axial and lateral seismic responses of underground pipeline through the fracture zone. The results show that the incident directions of seismic waves, width and soil property of the fracture zone have great influence on underground pipeline. It is suggested that the flexible joints with appropriate stiffness should be added into the pipeline near the interfaces between the fracture zone and the surrounded media.

Multiscale Site-Response Mapping: A Case Study of Parkfield, California

The scale of previously proposed methods for mapping site-response ranges from global coverage down to individual urban regions. Typically, spatial coverage and accuracy are inversely related. We use the densely spaced strong-motion stations in Parkfield, California, to estimate the accuracy of different site-response mapping methods and demonstrate a method for integrating multiple site-response estimates from the site to the global scale. This method is simply a weighted mean of a suite of different estimates, where the weights are the inverse of the variance of the individual estimates. Thus, the dominant site-response model varies in space as a function of the accuracy of the different models. For mapping applications, site-response models should be judged in terms of both spatial coverage and the degree of correlation with observed amplifications. Performance varies with period, but in general the Parkfield data show that: (1) where a velocity profile is available, the square-root-of-impedance (SRI) method outperforms the measured V S 30 (30 m divided by the S -wave travel time to 30 m depth) and (2) where velocity profiles are unavailable, the topographic slope method outperforms surficial geology for short periods, but geology outperforms slope at longer periods. We develop new equations to estimate site response from topographic slope, derived from the Next Generation Attenuation (NGA) database.

A Wavelet-based Seismogram Inversion Algorithm for the In Situ Characterization of Nonlinear Soil Behavior

We present a full waveform inversion algorithm of downhole array seismogram recordings that can be used to estimate the inelastic soil behavior in situ during earthquake ground motion. For this purpose, we first develop a new hysteretic scheme that improves upon existing nonlinear site response models by allowing adjustment of the width and length of the hysteresis loop for a relatively small number of soil parameters. The constitutive law is formulated to approximate the response of saturated cohesive materials, and does not account for volumetric changes due to shear leading to pore pressure development and potential liquefaction. We implement the soil model in the forward operator of the inversion, and evaluate the constitutive parameters that maximize the cross-correlation between site response predictions and observations on ground surface. The objective function is defined in the wavelet domain, which allows equal weight to be assigned across all frequency bands of the non-stationary signal. We evaluate the convergence rate and robustness of the proposed scheme for noise-free and noise-contaminated data, and illustrate good performance of the inversion for signal-to-noise ratios as low as 3. We finally employ the proposed scheme to downhole array data, and show that results compare very well with published data on generic soil conditions and previous geotechnical investigation studies at the array site. By assuming a realistic hysteretic model and estimating the constitutive soil parameters, the proposed inversion accounts for the instantaneous adjustment of soil response to the level and strain and load path during transient loading, and allows results to be used in predictions of nonlinear site effects during future events.

Earth modeling and estimation of the local seismic ground motion due to site geology in complex volcanoclastic areas

Volcanic areas often show complex behaviour as far as seismic waves propagation and seismic motion at surface are concerned. In fact, the finite lateral extent of surface layers such as lava flows, blocks, differential welding and/or zeolitization within pyroclastic deposits, introduces in the propagation of seismic waves effects such as the generation of surface waves at the edge, resonance in lateral direction, diffractions and scattering of energy, which tend to modify the amplitude as well as the duration of the ground motion. The irregular topographic surface, typical of volcanic areas, also strongly influences the seismic site response. Despite this heterogeneity, it is unfortunately a common geophysical and engineering practice to evaluate even in volcanic environments the subsurface velocity field with monodimensional investigation method (i.e. geognostic soundings, refraction survey, down-hole, etc.) prior to the seismic site response computation which in a such cases is obviously also made with 1D algorithms. This approach often leads to highly inaccurate results. In this paper we use a different approach, i.e. a fully 2D P-wave Cturning rayE tomographic survey followed by 2D seismic site response modeling. We report here the results of this approach in three sites located at short distance from Mt. Vesuvius and Campi Flegrei and characterized by overburdens constituted by volcanoclastic deposits with large lateral and vertical variations of their elastic properties. Comparison between 1D and 2D Dynamic Amplification Factor shows in all reported cases entirely different results, both in terms of peak period and spectral contents, as expected from the clear bidimensionality of the geological section. Therefore, these studies suggest evaluating carefully the subsoil geological structures in areas characterized by possible large lateral and vertical variations of the elastic properties in order to reach correct seismic site response curves to be used for engineering projects.

Site Effects for the Sacramento-San Joaquin Delta

Seismic site response and site effects models are presented for levees in the Sacramento-San Joaquin Delta where the subsurface soils include thick deposits of highly organic soils. Sources of uncertainty that contribute to the variation of seismic wave amplification are investigated, including variations in the input ground motions, soil profiles, and dynamic soil properties through Monte Carlo simulations of equivalent-linear site response analyses. Regression models for seismic wave amplification for levees in the Delta are presented that range from a function of peak outcrop acceleration alone to a vector of response spectra ordinates and soil profile parameters. The site effects models were incorporated into a probabilistic seismic hazard analysis for a representative location, and the relative impacts of the various models on the computed hazard are evaluated.

Learning of Dynamic Soil Behavior from Downhole Arrays

An increasing number of downhole arrays are deployed to measure motions at the ground surface and within the soil profile. Measurements from these arrays provide an opportunity to improve site response models and to better understand underlying dynamic soil behavior. Parametric inverse analysis approaches have been used to identify constitutive model parameters to achieve a better match with field observations. However, they are limited by the selected material model. Nonparametric inverse analysis approaches identify averaged soil behavior between measurement locations. A novel inverse analysis framework, self-learning simulations (SelfSim), is employed to reproduce the measured downhole array response while extracting the underlying soil behavior of individual soil layers unconstrained by prior assumptions of soil behavior. SelfSim is successfully applied to recordings from Lotung and La Cienega. The extracted soil behavior from few events can be used to reliably predict the measured response for other events. The field extracted soil behavior shows dependencies of shear modulus and damping on cyclic shear strain level, number of loading cycles, and strain rate that are similar qualitatively to those reported from laboratory studies but differ quantitatively.

Site response zones and short-period earthquake ground motion projections for the Las Vegas Basin

A deterministic seismic hazard analysis was conducted to address the effect of local soil conditions on earthquake-induced strong ground motion in the Las Vegas Basin, Nevada (US). Using a large geological and geotechnical database, two response units were defined: a fine-grained unit, predominantly clay; and a coarse-grained unit, predominantly gravel. A moderate number of high-quality shallow shear wave velocity measurements were collected from which characteristic shear wave velocity profiles were developed for each response unit. An equivalent-linear one-dimensional site response model was used. The model was calibrated using a basin-wide, small-strain ground motion database. Calibration tests showed that ground motion projections become increasingly conservative with increasing ground-motion amplitude. Projections were overconservative for the coarse-grained response unit, likely due to the sparseness of the velocity database. For the earthquake response analyses, historical ground motions were used to model characteristic ‘bedrock’ motion for earthquakes on 10 faults judged to be critical. Response spectral envelopes were generated for each unit through Monte-Carlo simulations. For the fine-grained response unit, 95th percentile peak ground acceleration, peak spectral acceleration and predominant period were 310 cm/s2, 1100 cm/s2, and 0.29 s, respectively. With respect to codified design spectra, projections are lower at short periods and higher at long periods. Projections of peak spectral accelerations for the coarse-grained response unit, were more than double that of codified spectra; however, they are believed to be overconservative. Near-fault effects and basin-edge effects, though potentially important, were not considered in these analyses.

Quantifying Nonlinearity Susceptibility via Site-Response Modeling Uncertainty at Three Sites in the Los Angeles Basin

The effects of near-surface soil stratigraphy on the amplitude and frequency content of ground motion are accounted for in most modern U.S. seismic design codes for building structures as a function of the soil conditions prevailing in the area of interest. Nonetheless, currently employed site-classification criteria do not adequately describe the nonlinearity susceptibility of soil formations, which prohibits the development of standardized procedures for the computationally efficient integration of nonlinear ground response analyses in broadband ground-motion simulations. In turn, the lack of a unified methodology for nonlinear site-response analyses affects both the prediction accuracy of site-specific ground-motion intensity measures and the evaluation of site-amplification factors when broadband simulations are used for the development of hybrid attenuation relations. In this article, we introduce a set of criteria for quantification of the nonlinearity susceptibility of soil profiles based on the site conditions and incident ground-motion characteristics, and we implement them to identify the least complex ground response prediction methodology required for the simulation of nonlinear site effects at three sites in the Los Angeles basin. The criteria are developed on the basis of a comprehensive nonlinear site-response modeling uncertainty analysis, which includes both detailed soil profile descriptions and statistical adequacy of ground-motion time histories. Approximate and incremental nonlinear models are implemented, and the limited site-response observations are initially compared to the ensemble site-response estimates. A suite of synthetic ground motions for rupture scenarios of weak, medium, and large magnitude events (M 3.5–7.5) is next generated, parametric studies are conducted for each fixed magnitude scenario by varying the source-to-site distance, and the variability introduced in ground-motion predictions is quantified for each nonlinear site-response methodology. A frequency index is developed to describe the frequency content of incident ground motion relative to the resonant frequencies of the soil profile, and this index is used in conjunction with the rock-outcrop acceleration peak amplitude (PGA_(RO)) to identify the site conditions and ground-motion characteristics where incremental nonlinear analyses should be employed in lieu of approximate methodologies. We show that the proposed intensity-frequency representation of ground motion may be implemented to describe the nonlinearity susceptibility of soil formations in broadband simulations by accounting both for the magnitude-distance-orientation characteristics of seismic motion and the profile stiffness characteristics. The synthetic ground-motion predictions are next used for the development of site-amplification factors for the alternative site-response methodologies, and the results are compared to published site factors of attenuation relations. For the site conditions investigated, currently established amplification factors compare well with synthetic simulations for class C and D site conditions, while long-period amplification factors are overestimated by a factor of 1.5 at the class E site, where site-specific nonlinear analyses should be employed for levels of PGA_(RO)>0.2g.

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

Quantifying sensitivity of local site response models to statistical variations in soil properties

We perform a combined stochastic-deterministic analysis of local site response using two computer codes, an equivalent linear analysis program SHAKE and a fully nonlinear finite element code SPECTRA. Our goal is to compare the relative sensitivity of the two codes to statistical variations in soil properties. For the case studies, we re-analyze two ground motion records in Lotung, Taiwan, and one ground motion record in Gilroy, California, utilizing the recorded ground motions at the site deterministically as input into the two codes while treating the uncertain soil parameters as random variables. We then obtain empirical cumulative distribution functions of the Arias intensity and acceleration spectrum intensity, two measures of cumulative damage, to compare the relative sensitivity of the two codes to variations in model parameters. We show that the two codes exhibit comparable sensitivities to statistical parameter variations, indicating that even in the presence of fluctuations in the soil parameter values it is possible to pursue a fully nonlinear site response analysis with SPECTRA and benefit from its superior accuracy.

Site-Response Models for Charleston, South Carolina, and Vicinity Developed from Shallow Geotechnical Investigations

The study models the response of near-surface materials in Charleston, South Carolina, and the adjacent area. Geotechnical investigations at 281 locations were made available by local engineering firms. The data used for dynamic site- response analysis were derived from shear-wave velocity measurements at 52 loca- tions. Site response was quantified as the ratio of surface motion to hypothetical hard-rock basement outcrop motion. Scenario earthquake motions were developed with the stochastic model. Acceleration response ratios for 5% critical oscillator damping were computed for 12 frequencies ranging from 0.1 to 30 Hz and for peak ground acceleration. Two features determine the general nature of site response in the study area: the impedance contrast between Mesozoic basement and Cretaceous sediments, and the shallow impedance contrast between Quaternary and Tertiary sediments. Average S- wave velocities in the Quaternary are relatively uniform and range from 150 to 250 m/sec. They are not strongly correlated with surface geology. The velocities of the immediately underlying Tertiary sediments range from 300 to 500 m/sec. Because of the uniformity of velocity in the Quaternary, depth to the Quaternary-Tertiary contact appears to be the most important variable leading to differences in calculated site response. This surface is irregular, and varies in depth from near surface at inland sites to approximately 30 m at sites near the coast. As a consequence, estimated site response in the frequency band 1-10 Hz varies by as much as a factor of 3. Site response at frequencies less than 1 Hz is dominated by the first few resonant har- monics of the entire sedimentary section, with fundamental frequency near 0.2 Hz.

Definition of subsoil structure and preliminary ground response in Aigion city (Greece) using microtremor and earthquakes

An extensive campaign—including detailed geologic and geotechnical surveys both existing and news as well as noise measurements—was conducted along a cross-section in order to define both geometry and soil properties (mainly the shear wave velocity) of the main formations in Aigion city. Aigion city is located in the Gulf of Corinth, Greece, a highly seismic region of the Aegean Sea. The main objective of the accurate 2D soil model is its use in site response modeling and in the interpretation of observations from a vertical down-hole accelerograph array. This model revealed a complex geologic structure with a multi-faulted shear zone related to the Aigion fault. The defined subsurface structure offered the possibility for its correlation with estimated site effects, in terms of spectral ratios. Two different data sets, earthquakes recorded at down-hole accelerograph network and noise measurements at 17 sites, were used. To translate the empirical transfer functions with the geologic structure, the 1D estimates were also computed. All these results are consistent, indicating a satisfactory correlation between the soil model and preliminary site response.