Search USGS⌕ Search

SEARCH · Search USGS

Results for “Soil Dynamics and Earthquake Engineering”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

Optimal-adaptive filters for modelling spectral shape, site amplification, and source scaling

Optimal filtering techniques have been used successfully in various areas in science and engineering. They are based on statistical properties of the signal and the noise, and stochastic approximation theory. In addition to filtering, optimal filters can also be used for smoothing, prediction, and system identification. This paper introduces some applications of optimal filtering techniques to earthquake engineering by using the so-called ARMAX models. Three applications are presented: (a) spectral modelling of ground accelerations, (b) site amplification (i.e., the relationship between two records obtained at different sites during an earthquake), and (c) source scaling (i.e., the relationship between two records obtained at a site during two different earthquakes). A numerical example for each application is presented by using recorded ground motions. The results show that the optimal filtering techniques provide elegant solutions to above problems, and can be a useful tool in earthquake engineering.

Soil Dynamics and Earthquake Engineering↗

Processing of strong-motion accelerograms: Needs, options and consequences

Recordings from strong-motion accelerographs are of fundamental importance in earthquake engineering, forming the basis for all characterizations of ground shaking employed for seismic design. The recordings, particularly those from analog instruments, invariably contain noise that can mask and distort the ground-motion signal at both high and low frequencies. For any application of recorded accelerograms in engineering seismology or earthquake engineering, it is important to identify the presence of this noise in the digitized time-history and its influence on the parameters that are to be derived from the records. If the parameters of interest are affected by noise then appropriate processing needs to be applied to the records, although it must be accepted from the outset that it is generally not possible to recover the actual ground motion over a wide range of frequencies. There are many schemes available for processing strong-motion data and it is important to be aware of the merits and pitfalls associated with each option. Equally important is to appreciate the effects of the procedures on the records in order to avoid errors in the interpretation and use of the results. Options for processing strong-motion accelerograms are presented, discussed and evaluated from the perspective of engineering application. ?? 2004 Elsevier Ltd. All rights reserved.

Soil Dynamics and Earthquake Engineering↗

Examination of differences between three SPT-based seismic soil liquefaction triggering relationships

The preceding companion paper presented the updating of the seismic soil liquefaction triggering relationship of Cetin et al. [1], and compared the resulting updated relationship with the earlier version. In this second paper, a detailed cross-comparison is made between three triggering relationships: (1) Seed et al. [2], as slightly updated by the NCEER Working Group (Youd et al. [3]), (2) Boulanger and Idriss [4], and (3) Cetin et al. [5]. Differences between these three triggering relationships, and the apparent causes of them are examined. Also studied are the impacts of these differences on levels of conservatism with regard to evaluation of liquefaction triggering hazard, and the resulting risks for engineering projects.

Soil Dynamics and Earthquake Engineering↗

GPS in dynamic monitoring of long-period structures

Global Positioning System (GPS) technology with high sampling rates (??? 10 samples per second) allows scientifically justified and economically feasible dynamic measurements of relative displacements of long-period structures-otherwise difficult to measure directly by other means, such as the most commonly used accelerometers that require post-processing including double integration. We describe an experiment whereby the displacement responses of a simulated tall building are measured clearly and accurately in real-time. Such measurements can be used to assess average drift ratios and changes in dynamic characteristics, and therefore can be used by engineers and building owners or managers to assess the building performance during extreme motions caused by earthquakes and strong winds. By establishing threshold displacements or drift ratios and identifying changing dynamic characteristics, procedures can be developed to use such information to secure public safety and/or take steps to improve the performance of the building. Published by Elsevier Science Ltd.

Soil Dynamics and Earthquake Engineering↗

Evaluation of simulated ground motions using probabilistic seismic demand analysis: CyberShake (ver. 15.12) simulations for Ordinary Standard Bridges

There is a need for benchmarking and validating simulated ground motions in order for them to be utilized by the engineering community. Such validation may be geared towards a specific ground motion simulation method, a target engineering application, and a specific location; the validation presented herein focuses on a bridge engineering application in southern California. Catalogs of simulated ground motions representing a 200,000-year forecast are selected from the Southern California Earthquake Center CyberShake version 15.12 database for five sites in Southern California (~20,000 unscaled ground motions per site). They are used in Non-Linear Time History Analysis (NLTHA) of four Ordinary Standard Bridge structures. For each site, these data are used to obtain simulation-based Engineering Demand Parameter (EDP) hazard curves. These are compared against EDP hazard curves that are constructed using conventional methods based on empirical models, i.e., using recorded ground motions through Incremental Dynamic Analysis and integration over the Intensity Measure (IM) hazard curve. The two sets of simulation-based and conventional EDP hazard curves are compared at various return periods. To further account for the differences between simulated and recorded ground motions, direct comparisons are also made between IM hazard curves for simulated and recorded catalogs, as well as the EDP versus IM data obtained from NLTHA of the bridges. We observe that CyberShake simulates motions that yield similar EDP values compared to empirical data for shorter return periods. For longer return periods, however, EDPs from the simulation-based analysis tend to be lower than the EDPs obtained from utilizing recorded ground motions for short-period bridges, while the opposite is the case for long-period bridges. It is recommended that validation efforts go beyond IM levels and also include comparisons of the relations between IMs and EDPs. Finally, site-specific relations are proposed that correlate the ratio between the two types of EDPs (simulation-based and conventional) with the hazard level, shallow site condition, and site basin depth.

Soil Dynamics and Earthquake Engineering Journal↗

Multichannel analysis of surface wave method with the autojuggie

The shear (S)-wave velocity of near-surface materials and its effect on seismic-wave propagation are of fundamental interest in many engineering, environmental, and groundwater studies. The multichannel analysis of surface wave (MASW) method provides a robust, efficient, and accurate tool to observe near-surface S-wave velocity. A recently developed device used to place large numbers of closely spaced geophones simultaneously and automatically (the 'autojuggie') is shown here to be applicable to the collection of MASW data. In order to demonstrate the use of the autojuggie in the MASW method, we compared high-frequency surface-wave data acquired from conventionally planted geophones (control line) to data collected in parallel with the automatically planted geophones attached to steel bars (test line). The results demonstrate that the autojuggie can be applied in the MASW method. Implementation of the autojuggie in very shallow MASW surveys could drastically reduce the time required and costs incurred in such surveys. ?? 2003 Elsevier Science Ltd. All rights reserved.

Soil Dynamics and Earthquake Engineering↗

Response of pendulums to complex input ground motion

Dynamic response of most seismological instruments and many engineering structures to ground shaking can be represented via response of a pendulum (single-degree-of-freedom oscillator). In most studies, pendulum response is simplified by considering the input from uni-axial translational motion alone. Complete ground motion however, includes not only translational components but also rotations (tilt and torsion). In this paper, complete equations of motion for three following types of pendulum are described: (i) conventional (mass-on-rod), (ii) mass-on-spring type, and (iii) inverted (astatic), then their response sensitivities to each component of complex ground motion are examined. The results of this study show that a horizontal pendulum similar to an accelerometer used in strong motion measurements is practically sensitive to translational motion and tilt only, while inverted pendulum commonly utilized to idealize multi-degree-of-freedom systems is sensitive not only to translational components, but also to angular accelerations and tilt. For better understanding of the inverted pendulum's dynamic behavior under complex ground excitation, relative contribution of each component of motion on response variants is carefully isolated. The systematically applied loading protocols indicate that vertical component of motion may create time-dependent variations on pendulum's oscillation period; yet most dramatic impact on response is produced by the tilting (rocking) component. ?? 2007 Elsevier Ltd. All rights reserved.

Soil Dynamics and Earthquake Engineering↗

Improved computational methods for probabilistic liquefaction hazard analysis

Current procedures for analysis of and design against liquefaction hazards focus primarily on the use of probabilistic ground motions at a single ground-shaking hazard level, with the cyclic loading represented by a peak ground acceleration ( PGA ) corresponding to a target return period and a single representative moment magnitude M w . These parameters are typically used in conjunction with deterministic simplified procedures for estimating liquefaction triggering and consequences that largely neglect the considerable uncertainties inherent to liquefaction problems. Motivated by these limitations and the resulting inconsistencies in liquefaction design levels, early methods for probabilistic liquefaction hazard analysis (PLHA) were proposed that incorporate the full ground motion hazard space, integrated with probabilistic liquefaction triggering models. Such methods provide liquefaction factor of safety ( FS L ) hazard curves for standard penetration test (SPT) data. Recognizing the increased use of higher-resolution cone penetrometer test (CPT) data in engineering analysis and design and the potential computational challenges it presents, an expanded suite of probabilistic triggering models, and wider availability of more detailed seismic hazard data, an improved PLHA computational methodology is presented in this study. The methods described utilize the U.S. Geological Survey National Seismic Hazard Model directly through web services to obtain PGA hazard and disaggregation calculations for any site and average shear wave velocity in the upper 30 m of the site ( V s 30 ) in the conterminous United States, to reconstruct the PGA hazard space for use in PLHA calculations, and to employ array calculations for efficient liquefaction hazard curve estimates for the thousands of CPT measurements in a given profile. The framework presented here is modular in nature, and discussion on the use of alternative models and extension to probabilistic liquefaction consequence evaluation is presented. The importance of appropriate representation of probabilistic liquefaction model uncertainties is also highlighted, along with the impacts of different levels of uncertainty on the PLHA calculation. Finally, a potential roadmap for incorporation of the PLHA framework in seismic provisions is presented, with an illustration of how it can address current limitations and impacts in liquefaction hazard analysis and design.

Soil Dynamics and Earthquake Engineering↗

Pseudospectral modeling and dispersion analysis of Rayleigh waves in viscoelastic media

Multichannel Analysis of Surface Waves (MASW) is one of the most widely used techniques in environmental and engineering geophysics to determine shear-wave velocities and dynamic properties, which is based on the elastic layered system theory. Wave propagation in the Earth, however, has been recognized as viscoelastic and the propagation of Rayleigh waves presents substantial differences in viscoelastic media as compared with elastic media. Therefore, it is necessary to carry out numerical simulation and dispersion analysis of Rayleigh waves in viscoelastic media to better understand Rayleigh-wave behaviors in the real world. We apply a pseudospectral method to the calculation of the spatial derivatives using a Chebyshev difference operator in the vertical direction and a Fourier difference operator in the horizontal direction based on the velocity-stress elastodynamic equations and relations of linear viscoelastic solids. This approach stretches the spatial discrete grid to have a minimum grid size near the free surface so that high accuracy and resolution are achieved at the free surface, which allows an effective incorporation of the free surface boundary conditions since the Chebyshev method is nonperiodic. We first use an elastic homogeneous half-space model to demonstrate the accuracy of the pseudospectral method comparing with the analytical solution, and verify the correctness of the numerical modeling results for a viscoelastic half-space comparing the phase velocities of Rayleigh wave between the theoretical values and the dispersive image generated by high-resolution linear Radon transform. We then simulate three types of two-layer models to analyze dispersive-energy characteristics for near-surface applications. Results demonstrate that the phase velocity of Rayleigh waves in viscoelastic media is relatively higher than in elastic media and the fundamental mode increases by 10-16% when the frequency is above 10. Hz due to the velocity dispersion of P and S waves. ?? 2011 Elsevier Ltd.

Soil Dynamics and Earthquake Engineering↗

An empirical method to estimate shear wave velocity of soils in the New Madrid seismic zone

In this study, a set of charts are developed to estimate shear wave velocity of soils in the New Madrid seismic zone (NMSZ), using the standard penetration test (SPT) N values and soil depths. Laboratory dynamic test results of soil samples collected from the NMSZ showed that the shear wave velocity of soils is related to the void ratio and the effective confining pressure applied to the soils. The void ratio of soils can be estimated from the SPT N values and the effective confining pressure depends on the depth of soils. Therefore, the shear wave velocity of soils can be estimated from the SPT N value and the soil depth. To make the methodology practical, two corrections should be made. One is that field SPT N values of soils must be adjusted to an unified SPT N??? value to account the effects of overburden pressure and equipment. The second is that the effect of water table to effective overburden pressure of soils must be considered. To verify the methodology, shear wave velocities of five sites in the NMSZ are estimated and compared with those obtained from field measurements. The comparison shows that our approach and the field tests are consistent with an error of less than of 15%. Thus, the method developed in this study is useful for dynamic study and practical designs in the NMSZ region. Copyright ?? 1996 Elsevier Science Limited.

Soil Dynamics and Earthquake Engineering↗

Effect of tilt on strong motion data processing

In the near-field of an earthquake the effects of the rotational components of ground motion may not be negligible compared to the effects of translational motions. Analyses of the equations of motion of horizontal and vertical pendulums show that horizontal sensors are sensitive not only to translational motion but also to tilts. Ignoring this tilt sensitivity may produce unreliable results, especially in calculations of permanent displacements and long-period calculations. In contrast to horizontal sensors, vertical sensors do not have these limitations, since they are less sensitive to tilts. In general, only six-component systems measuring rotations and accelerations, or three-component systems similar to systems used in inertial navigation assuring purely translational motion of accelerometers can be used to calculate residual displacements. ?? 2004 Elsevier Ltd. All rights reserved.

Soil Dynamics and Earthquake Engineering↗

Assessment of seismic risk in Tashkent, Uzbekistan and Bishkek, Kyrgyz Republic

The impact of earthquakes in urban centers prone to disastrous earthquakes necessitates the analysis of associated risk for rational formulation of contingency plans and mitigation strategies. In urban centers the seismic risk is best quantified and portrayed through the preparation of 'Earthquake damage and Loss Scenarios'. The components of such scenarios are the assessment of the hazard, inventories and the vulnerabilities of elements at risk. For the development of earthquake risk scenario in Tashkent-Uzbekistan and Bishkek-Kyrgyzstan an approach based on spectral displacements is utilized. This paper will present the important features of a comprehensive study, highlight the methodology, discuss the results and provide insights to the future developments. ?? 2005 Elsevier Ltd. All rights reserved.

Soil Dynamics and Earthquake Engineering↗

Time-domain representation of frequency-dependent foundation impedance functions

Foundation impedance functions provide a simple means to account for soil-structure interaction (SSI) when studying seismic response of structures. Impedance functions represent the dynamic stiffness of the soil media surrounding the foundation. The fact that impedance functions are frequency dependent makes it difficult to incorporate SSI in standard time-history analysis software. This paper introduces a simple method to convert frequency-dependent impedance functions into time-domain filters. The method is based on the least-squares approximation of impedance functions by ratios of two complex polynomials. Such ratios are equivalent, in the time-domain, to discrete-time recursive filters, which are simple finite-difference equations giving the relationship between foundation forces and displacements. These filters can easily be incorporated into standard time-history analysis programs. Three examples are presented to show the applications of the method.

Soil Dynamics and Earthquake Engineering↗

Simple equations guide high-frequency surface-wave investigation techniques

We discuss five useful equations related to high-frequency surface-wave techniques and their implications in practice. These equations are theoretical results from published literature regarding source selection, data-acquisition parameters, resolution of a dispersion curve image in the frequency-velocity domain, and the cut-off frequency of high modes. The first equation suggests Rayleigh waves appear in the shortest offset when a source is located on the ground surface, which supports our observations that surface impact sources are the best source for surface-wave techniques. The second and third equations, based on the layered earth model, reveal a relationship between the optimal nearest offset in Rayleigh-wave data acquisition and seismic setting - the observed maximum and minimum phase velocities, and the maximum wavelength. Comparison among data acquired with different offsets at one test site confirms the better data were acquired with the suggested optimal nearest offset. The fourth equation illustrates that resolution of a dispersion curve image at a given frequency is directly proportional to the product of a length of a geophone array and the frequency. We used real-world data to verify the fourth equation. The last equation shows that the cut-off frequency of high modes of Love waves for a two-layer model is determined by shear-wave velocities and the thickness of the top layer. We applied this equation to Rayleigh waves and multi-layer models with the average velocity and obtained encouraging results. This equation not only endows with a criterion to distinguish high modes from numerical artifacts but also provides a straightforward means to resolve the depth to the half space of a layered earth model. ?? 2005 Elsevier Ltd. All rights reserved.

Soil Dynamics and Earthquake Engineering↗

Spatial correlation of shear-wave velocity in the San Francisco Bay Area sediments

Ground motions recorded within sedimentary basins are variable over short distances. One important cause of the variability is that local soil properties are variable at all scales. Regional hazard maps developed for predicting site effects are generally derived from maps of surficial geology; however, recent studies have shown that mapped geologic units do not correlate well with the average shear-wave velocity of the upper 30 m, V s (30). We model the horizontal variability of near-surface soil shear-wave velocity in the San Francisco Bay Area to estimate values in unsampled locations in order to account for site effects in a continuous manner. Previous geostatistical studies of soil properties have shown horizontal correlations at the scale of meters to tens of meters while the vertical correlations are on the order of centimeters. In this paper we analyze shear-wave velocity data over regional distances and find that surface shear-wave velocity is correlated at horizontal distances up to 4 km based on data from seismic cone penetration tests and the spectral analysis of surface waves. We propose a method to map site effects by using geostatistical methods based on the shear-wave velocity correlation structure within a sedimentary basin. If used in conjunction with densely spaced shear-wave velocity profiles in regions of high seismic risk, geostatistical methods can produce reliable continuous maps of site effects.

California↗

On using surface-source downhole-receiver logging to determine seismic slownesses

We present a method to solve for slowness models from surface-source downhole-receiver seismic travel-times. The method estimates the slownesses in a single inversion of the travel-times from all receiver depths and accounts for refractions at layer boundaries. The number and location of layer interfaces in the model can be selected based on lithologic changes or linear trends in the travel-time data. The interfaces based on linear trends in the data can be picked manually or by an automated algorithm. We illustrate the method with example sites for which geologic descriptions of the subsurface materials and independent slowness measurements are available. At each site we present slowness models that result from different interpretations of the data. The examples were carefully selected to address the reliability of interface-selection and the ability of the inversion to identify thin layers, large slowness contrasts, and slowness gradients. Additionally, we compare the models in terms of ground-motion amplification. These plots illustrate the sensitivity of site amplifications to the uncertainties in the slowness model. We show that one-dimensional site amplifications are insensitive to thin layers in the slowness models; although slowness is variable over short ranges of depth, this variability has little affect on ground-motion amplification at frequencies up to 5 Hz.

Soil Dynamics and Earthquake Engineering↗

Research on the middle-of-receiver-spread assumption of the MASW method

The multichannel analysis of surface wave (MASW) method has been effectively used to determine near-surface shear- (S-) wave velocity. Estimating the S-wave velocity profile from Rayleigh-wave measurements is straightforward. A three-step process is required to obtain S-wave velocity profiles: acquisition of a multiple number of multichannel records along a linear survey line by use of the roll-along mode, extraction of dispersion curves of Rayleigh waves, and inversion of dispersion curves for an S-wave velocity profile for each shot gather. A pseudo-2D S-wave velocity section can be generated by aligning 1D S-wave velocity models. In this process, it is very important to understand where the inverted 1D S-wave velocity profile should be located: the midpoint of each spread (a middle-of-receiver-spread assumption) or somewhere between the source and the last receiver. In other words, the extracted dispersion curve is determined by the geophysical structure within the geophone spread or strongly affected by the source geophysical structure. In this paper, dispersion curves of synthetic datasets and a real-world example are calculated by fixing the receiver spread and changing the source location. Results demonstrate that the dispersion curves are mainly determined by structures within a receiver spread. ?? 2008 Elsevier Ltd. All rights reserved.

Soil Dynamics and Earthquake Engineering↗

Low-velocity zone and topography as a source of site amplification effect on Tarzana hill, California

Tarzana station is located in the foothills of the Santa Monica Mountains in California near the crest of a low (<20 m) natural hill with gentle slopes. The hill is about 500 m in length by 130 m in width and is formed of extremely weathered shale at the surface to fresh at depth. Average S-wave is about 250 m/s in the top 17-18 m, and S- and P-wave velocities significantly increase below this depth. According to the NEHRP classification based on VS30???300 m/s it is a site class D. Strong-motion instrumentation at Tarzana consisted of an accelerograph at the top of the hill, a downhole instrument at 60 m depth, and an accelerograph at the base of the hill. More than 20 earthquakes were recorded by at least three instruments at Tarzana from 1998 till 2003. Comparisons of recordings and Fourier spectra indicate strong directional resonance in a direction perpendicular to the strike of the hill. The dominant peaks in ground motion amplification on the top of the hill relative to the base are at frequencies ???3.6 and 8-9 Hz for the horizontal components. Our hypothesis is that the hill acts like a wave trap. This results in an amplification at predominant frequencies f=V/4 h (h is layer's thickness) at f???3.6 Hz for S-waves (using average VS17=246 m/s and h=17 m) and f???7.9 Hz for P-waves (using average VP17=535 m/s and h=17 m). As was shown by Bouchon and Barker [Seismic response of a hill: the example of Tarzana, California. Bull Seism Soc Am 1996;86(1A):66-72], topography of this hill amplifies and polarizes ground motion in the frequency range of 3-5 Hz. Hill acts as a magnifying polarizing glass: It polarizes ground motion in the direction perpendicular to the strike of the hill and also amplifies ground motions that had been also amplified by a low-velocity layer.

Soil Dynamics and Earthquake Engineering↗