Prediction of the severity of ground failure in Quaternary deposits is a critical component of hazard studies. Model development in our project is focused on design and application of methods for quantitative assessment of ground deformation potential.
Modeling Shear Wave Velocity Structure Using Spectral Analysis of Surface Waves (SASW) techniques

Collaborative with Woods Hole Coastal and Marine Science Center
The SASW technique is ideal for the assessment of ground failure in the coastal and offshore environment principally because of the portability of the apparatus and the non-invasive nature of the technique. For onshore deployments we use a spectrum analyzer and commercially available vertical 1-Hz sensors. For offshore investigation, we use custom designed ocean bottom seismometers (OBS) built into a ship deployable seafloor SASW system. The OBS units are arrayed in a linear multi-station configuration and mounted onto a linear support frame. The surface wave source for onshore and offshore work is an electromechanical shaker. For offshore deployment, the shaker is housed in an air filled vessel and deployed at the end of the support frame. We routinely use the SASW method to develop profiles of shear wave velocity in the upper 30+ meters of the earth. These profiles can be used as part of a liquefaction hazard assessment and as the basis for predicting response of sites to incoming earthquake ground motions.
Probabilistic Liquefaction Triggering Assessment
Collaborative research with U.C. Berkeley; PEER; Kobe University; and the China Seismological Bureau
The probabilistic models for liquefaction assessment are being advanced along two tracks. First, global data sets are being assembled to identify the ground conditions and earthquake stresses at the majority of the world's documented ground failure sites. To do this we have cataloged existing data sets, and have revisited and tested over 300 sites in China, Japan, Taiwan and the USA. The second track involves modeling these data through a Bayesian updating and structural reliability engine used to assess the probability of failure in a component or system. Here, the system is a geotechnical model of soil resistance to liquefaction occurrence, and the loads are associated with earthquake motion.
At first step, we are gathering global databases for field measurements of liquefaction resistance by cataloging existing data and collecting new data for penetration-based, and shear wave velocity-based, liquefaction evaluation measurements. To acquire new Vs field profiles, we use a variety of active and passive surface wave methods and analysis techniques at liquefaction evaluation sites already documented by drill logs
To correlate these databases with likelihood of initiation of seismic-soil liquefaction, we utilize high-order probabilistic tools (Bayesian updating) developed for structural engineering reliability. A multi-parameter limit-state function for liquefaction triggering is modeled and evaluated based on the means, distributions and uncertainties of each model-parameter. Each case history is then sub-divided into 'quality'-ranking categories based on the conjugate-uncertainties of the estimated earthquake induced stress and the measured field-based liquefaction resistance of the ground. A low-pass cut-off in the coefficient of variation is used filter-out poorly constrained sites. Finally for the probabilistic analysis, the Bayesian updating procedure is used to iteratively compute coefficients for the limit-state function that minimize model error. The intended outcome of this effort is a new evaluation of liquefaction-triggering boundaries in light of global data sets and modern limit-state probabilistic tools.
Probabilistic Assessment of Dynamic Displacements
This project is focused on developing probabilistic models for multidirectional seismic shear displacements in soil for liquefiable and non-liquefiable deposits. Bayesian methods, described above, and Newmark-type models for computing seismic slope displacements are used to estimate the amplitude of ground failures during earthquakes. Learn more: “Seismic displacement of gently-sloping coastal and marine sediment under multidirectional earthquake loading”
This computer modeling effort is part of the USGS Pacific Coastal and Marine Science Center Global Geoengineering Research.
Global Geoengineering Research
Global Geoengineering Research: Methods of study
Below are publications associated with this project.
Examination of differences between three SPT-based seismic soil liquefaction triggering relationships
Seismic velocity site characterization of 10 Arizona strong-motion recording stations by spectral analysis of surface wave dispersion
Seismic displacement of gently-sloping coastal and marine sediment under multidirectional earthquake loading
Summary of SPT based field case history data of CETIN (2016) database
A comparative study of surface waves inversion techniques at strong motion recording sites in Greece
Shear Wave Velocity and Site Amplification Factors for 25 Strong-Motion Instrument Stations Affected by the M5.8 Mineral, Virginia, Earthquake of August 23, 2011
Shear-wave velocity and site-amplification factors for 50 Australian sites determined by the spectral analysis of surface waves method
Surface wave site characterization at 27 locations near Boston, Massachusetts, including 2 strong-motion stations
Shear-wave velocity-based probabilistic and deterministic assessment of seismic soil liquefaction potential
- Overview
Prediction of the severity of ground failure in Quaternary deposits is a critical component of hazard studies. Model development in our project is focused on design and application of methods for quantitative assessment of ground deformation potential.
Modeling Shear Wave Velocity Structure Using Spectral Analysis of Surface Waves (SASW) techniques
Sources/Usage: Public Domain. Visit Media to see details.Site ID, location, and average shear wave velocity (left) for surface wave test site 997-Y16A; average dispersion curve in gray, and individual empirical dispersion curves in black (center plot); shear wave velocity profile computed by three inversion methods (right plot). Learn more: USGS OFR 2016-1208 Collaborative with Woods Hole Coastal and Marine Science Center
The SASW technique is ideal for the assessment of ground failure in the coastal and offshore environment principally because of the portability of the apparatus and the non-invasive nature of the technique. For onshore deployments we use a spectrum analyzer and commercially available vertical 1-Hz sensors. For offshore investigation, we use custom designed ocean bottom seismometers (OBS) built into a ship deployable seafloor SASW system. The OBS units are arrayed in a linear multi-station configuration and mounted onto a linear support frame. The surface wave source for onshore and offshore work is an electromechanical shaker. For offshore deployment, the shaker is housed in an air filled vessel and deployed at the end of the support frame. We routinely use the SASW method to develop profiles of shear wave velocity in the upper 30+ meters of the earth. These profiles can be used as part of a liquefaction hazard assessment and as the basis for predicting response of sites to incoming earthquake ground motions.
Probabilistic Liquefaction Triggering Assessment
Plots of cyclic stress ratio (CSR) and effective stress normalized shear wave velocity Vs1, used just to demonstrate types of data modeling. Collaborative research with U.C. Berkeley; PEER; Kobe University; and the China Seismological Bureau
The probabilistic models for liquefaction assessment are being advanced along two tracks. First, global data sets are being assembled to identify the ground conditions and earthquake stresses at the majority of the world's documented ground failure sites. To do this we have cataloged existing data sets, and have revisited and tested over 300 sites in China, Japan, Taiwan and the USA. The second track involves modeling these data through a Bayesian updating and structural reliability engine used to assess the probability of failure in a component or system. Here, the system is a geotechnical model of soil resistance to liquefaction occurrence, and the loads are associated with earthquake motion.
At first step, we are gathering global databases for field measurements of liquefaction resistance by cataloging existing data and collecting new data for penetration-based, and shear wave velocity-based, liquefaction evaluation measurements. To acquire new Vs field profiles, we use a variety of active and passive surface wave methods and analysis techniques at liquefaction evaluation sites already documented by drill logs
To correlate these databases with likelihood of initiation of seismic-soil liquefaction, we utilize high-order probabilistic tools (Bayesian updating) developed for structural engineering reliability. A multi-parameter limit-state function for liquefaction triggering is modeled and evaluated based on the means, distributions and uncertainties of each model-parameter. Each case history is then sub-divided into 'quality'-ranking categories based on the conjugate-uncertainties of the estimated earthquake induced stress and the measured field-based liquefaction resistance of the ground. A low-pass cut-off in the coefficient of variation is used filter-out poorly constrained sites. Finally for the probabilistic analysis, the Bayesian updating procedure is used to iteratively compute coefficients for the limit-state function that minimize model error. The intended outcome of this effort is a new evaluation of liquefaction-triggering boundaries in light of global data sets and modern limit-state probabilistic tools.
Probabilistic Assessment of Dynamic Displacements
This project is focused on developing probabilistic models for multidirectional seismic shear displacements in soil for liquefiable and non-liquefiable deposits. Bayesian methods, described above, and Newmark-type models for computing seismic slope displacements are used to estimate the amplitude of ground failures during earthquakes. Learn more: “Seismic displacement of gently-sloping coastal and marine sediment under multidirectional earthquake loading”
- Science
This computer modeling effort is part of the USGS Pacific Coastal and Marine Science Center Global Geoengineering Research.
Global Geoengineering Research
The Pacific Coastal and Marine Science Center geotechnical group investigates the causes of ground deformation and ground failure as a result of earthquakes, storms, and wave actionGlobal Geoengineering Research: Methods of study
Field methods of study in global geoengineering research - Publications
Below are publications associated with this project.
Examination of differences between three SPT-based seismic soil liquefaction triggering relationships
The preceding companion paper presented the updating of the seismic soil liquefactiontriggering 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) BoulangSeismic velocity site characterization of 10 Arizona strong-motion recording stations by spectral analysis of surface wave dispersion
Vertical one-dimensional shear wave velocity (VS) profiles are presented for strong-motion sites in Arizona for a suite of stations surrounding the Palo Verde Nuclear Generating Station. The purpose of the study is to determine the detailed site velocity profile, the average velocity in the upper 30 meters of the profile (VS30), the average velocity for the entire profile (VSZ), and the National ESeismic displacement of gently-sloping coastal and marine sediment under multidirectional earthquake loading
Gentle sediment-laden slopes are typical of the onshore coastal zone and offshore continental shelf and slope. Coastal sediment are commonly young weakly consolidated materials that are well stratified, have low strength, and can mobilize shear displacements at low levels of stress. Seismically-driven plastic displacements of these sediment pose a hazard to coastal cities, buried onshore utilitiesSummary of SPT based field case history data of CETIN (2016) database
This report provides documentation of the Cetin et al. (2016) field performance case histories, probabilistic maximum likelihood assessment and the sources of differences between the liquefaction triggering resistance estimations (CRR values) of the widely used liquefaction triggering relationships of Seed et al. (1985), Cetin et al. (2004, 2016) and Boulanger and Idriss (2012). Cetin et al. (2016A comparative study of surface waves inversion techniques at strong motion recording sites in Greece
Surface wave method was used for the estimation of Vs vs depth profile at 10 strong motion stations in Greece. The dispersion data were obtained by SASW method, utilizing a pair of electromechanical harmonic-wave source (shakers) or a random source (drop weight). In this study, three inversion techniques were used a) a recently proposed Simplified Inversion Method (SIM), b) an inversion techniqueShear Wave Velocity and Site Amplification Factors for 25 Strong-Motion Instrument Stations Affected by the M5.8 Mineral, Virginia, Earthquake of August 23, 2011
Vertical one-dimensional shear wave velocity (Vs) profiles are presented for 25 strong-motion instrument sites along the Mid-Atlantic eastern seaboard, Piedmont region, and Appalachian region, which surround the epicenter of the M5.8 Mineral, Virginia, Earthquake of August 23, 2011. Testing was performed at sites in Pennsylvania, Maryland, West Virginia, Virginia, the District of Columbia, North CShear-wave velocity and site-amplification factors for 50 Australian sites determined by the spectral analysis of surface waves method
One-dimensional shear-wave velocity (VS ) profiles are presented at 50 strong motion sites in New South Wales and Victoria, Australia. The VS profiles are estimated with the spectral analysis of surface waves (SASW) method. The SASW method is a noninvasive method that indirectly estimates the VS at depth from variations in the Rayleigh wave phase velocity at the surface.Surface wave site characterization at 27 locations near Boston, Massachusetts, including 2 strong-motion stations
The geotechnical properties of the soils in and around Boston, Massachusetts, have been extensively studied. This is partly due to the importance of the Boston Blue Clay and the extent of landfill in the Boston area. Although New England is not a region that is typically associated with seismic hazards, there have been several historical earthquakes that have caused significant ground shaking (forShear-wave velocity-based probabilistic and deterministic assessment of seismic soil liquefaction potential
Shear-wave velocity (Vs) offers a means to determine the seismic resistance of soil to liquefaction by a fundamental soil property. This paper presents the results of an 11-year international project to gather new Vs site data and develop probabilistic correlations for seismic soil liquefaction occurrence. Toward that objective, shear-wave velocity test sites were identified, and measurements made