Geology topics

Alan Yong

Publications and source records attributed to Alan Yong.

At least 19 recordsLinked to original sources

USGS2026

Site-specific amplifications in Northwestern Turkiye: A generic approach

In this study, we derive generic site amplification functions applicable to Northwestern Türkiye and follow guidance on site classifications established by the United States National Earthquake Hazards Reduction Program (NEHRP). We employ the one-dimensional (1-D) equivalent linear ground response analysis method and focus on recordings from a select suite of free-field strong motion (SM) stations. These particular SM stations are associated with a comprehensive geotechnical site conditions dataset comprising 76 shear-wave velocity profiles. Despite their varying resolutions and depths, this diversity allows for a more representative analysis across different site conditions, enhancing the robustness of our findings. For analyses, we utilize a robust set of global ground motion records, selected from both within and outside of Türkiye, to serve as strong and weak input bedrock motions. Site amplification is determined by propagating these bedrock records through 76 site-specific Afet ve Acil Durum Yönetimi Başkanlığı (AFAD) station profiles and calculating the spectral acceleration ratio between the ground surface and the input motion. Our analyses reveal substantial disparities between the amplification functions derived from 1-D site response analyses using strong and weak input motions. These disparities are central to our analysis and inform the subsequent comparison with site amplification functions reported in the literature. We compare the generic site amplification functions derived in this study with those reported in the literature for regions with similar geological and seismic conditions to Northwestern Türkiye. For example, site amplification functions for NEHRP Site Classes C and D have been documented for Greece and the United States utilizing frequency-dependent methods such as the quarter-wavelength approximation. Our findings reveal that the generic site amplification functions derived in this study exhibit differences in peak amplitudes and predominant frequencies compared to those in the literature, reflecting regional variations in geotechnical and seismic properties. This study represents an initial effort to formulate generic site amplification functions for Türkiye, with a particular emphasis on the Northwestern region. A thorough literature review indicated no prior studies have proposed such functions specifically for this area.

Northwestern Turkiye↗

USGS2025

Special issue "Effects of surface geology on seismic motion (ESG): General state-of-research"

To understand and predict the behaviors of strong ground motions incurred during devastating earthquakes, studies focused on the “effects of surface geology (ESG) on seismic motions” (following tradition, herein contracted as: “ESG”) have progressed steadily in the last three decades. These ESG studies typically involve in situ measurements of the wavefield and commonly apply both analytical and computational approaches. Concurrently, improvements in ESG related research can be readily attributed to the proliferation of openly accessible strong motion data, as well as advances in computational power. Nevertheless, there remain significant shortfalls in our understanding of the epistemic and aleatory uncertainties associated with the ESG as demonstrated by phenomena from recent deadly earthquake-related site effects. Thus, investigations toward addressing these deficiencies should be underscored in future earthquake related disaster mitigation efforts. This special issue of Earth, Planets, and Space (EPS) is dedicated to the ongoing efforts by the International Association of Seismology and Physics of the Earth’s Interior (IASPEI) and the International Association of Earthquake Engineering (IAEE) that champion activities promoting studies related to the ESG.

Earth, Planets and Space↗

USGS2024

Are the horizontal-to-vertical spectral ratios of earthquakes and microtremors the same?

We consider the similarities and differences between earthquake and microtremor horizontal‐to‐vertical spectral ratios (eHVSR and mHVSR, respectively) using a dataset of 161 sites in southern California. Quantitative comparisons are made in terms of the eHVSR and mHVSR lognormal median curves, as well as the frequencies and amplitudes associated with the fundamental‐ and higher‐mode resonances where present. The results show only 58% of the eHVSR–mHVSR pairs agree in terms of their median curve and only 25% of the eHVSR–mHVSR pairs agree in terms of shared resonances, which increases to 68% if flat HVSRs are considered equivalent. Furthermore, while the shared resonances match very well in terms of frequency (root mean square error, RMSE, <0.11 Hz), the amplitudes of those resonances do not agree (RMSE >1.6). These findings demonstrate that while eHVSR and mHVSR agree at some sites, they are not equivalent at all sites. To investigate if the agreement between eHVSR and mHVSR could be related to features of the microtremor data, earthquake recordings, and/or the site conditions, three machine learning (ML) models at varying levels of interpretability are presented. The ML models—which include multivariate logistic regression, gradient‐boosted trees, and support vector machines—show only partial success at using site‐specific data to predict whether eHVSR and mHVSR will likely agree in terms of their median curve (accuracy of 78%) and number of resonances (accuracy of 84%). Therefore, we conclude that while eHVSR and mHVSR can be quite similar in terms of resonant frequencies at some sites, they are not identical at all sites. Furthermore, preliminary evidence shows that the agreement of eHVSR and mHVSR can be predicted a priori given features of the microtremor measurements, earthquake recordings, and site conditions, although a larger dataset will be necessary for developing a robust predictive model.

California↗

USGS2024

Application of non-stationary shear-wave velocity randomization approach to predict 1D seismic site response and its variability at two downhole array recordings

Accounting for uncertainties in seismic site response is crucial to improving the performance of one-dimensional (1D) ground response analyses (GRAs) at downhole array recording sites. In addition to site effects, uncertainties in 1D-GRAs can also be contributed from the seismic source and/or path. Though often representing not more than one percent of the distance (path) from the source, site conditions are known to have an enormous influence on ground shaking. In this study, we focus on the site shear-wave velocity ( V S ) structure, which is the main ingredient for estimating the variability of site response. As such, V S can manifest aleatory uncertainties related to the effects of small-scale spatial heterogeneities within the near surface, thus V S can substantially modify ground shaking during earthquakes. We apply a novel V S randomization approach to propagate the small-scale heterogeneities of V S to estimate seismic site response within a non-stationary probabilistic framework. The randomization approach generates samples of V S profiles that are used to perform several 1D-GRAs and obtain an averaged site response and related variability. The proposed method is implemented on data recorded at two downhole array sites with different subsurface soil conditions: a soft soil site on Treasure Island (California, United States of America) and a rock outcrop site in Cadarache (South-East France). We show that synthetic surface-to-borehole transfer functions from 1D-GRAs provide an acceptable fit to the empirical transfer functions from low-motion earthquake records and succeed in reproducing most of the site-specific seismic response variability. The remaining mismatch between transfer functions is likely due to insufficient precision on the seismic bedrock and the impedance contrast. The variability in site response is discussed with emphasis on the role of V S small-scale heterogeneities, attenuation, and input motion incidence angle in ground motion variability for the site and soil conditions at both locations.

California↗

USGS2024

Developing international standards and guidelines for disseminating and validating simulated ground motions

We are leading an effort to develop international standards and guidelines for curating, disseminating, and validating simulated ground-motion data. This effort is organized as a working group within the Consortium of Organizations for Strong Motion Observation Systems (COSMOS). In 2022 we held online workshops on 7 and 8 June and 20 October. The first workshop focused on curating and disseminating simulated ground-motion data. The second workshop focused on validating simulated ground-motions for engineering applications. About 100 people participated in each of the workshops with strong representation from North America, Europe, and western Asia. In the coming year, we intend to form a technical committee to draft the international guidelines and standards while continuing to engage stakeholders from across the globe through online and in-person meetings. The key points from the workshops include: (1) Numerous groups are generating simulated earthquake ground motions and making them openly available; however there is very little coordination among groups to provide consistent interfaces for searching and retrieving data; (2) Standardizing interfaces for metadata and data access should consider agile approaches that can adapt to changing capabilities and user needs while building upon existing efforts; (3) A distributed architecture is preferred by most participants to allow institutions to host and manage their own data while broadcasting their holdings to a combined catalog; (4) Validation of ground-motion simulations applies to the entire workflow for simulating earthquake ground-motions, including the rupture model, seismic velocity model, and seismic wave propagation software; (5) Metrics to evaluate the validation are application dependent; more research is needed to tie structural response characteristics to ground-motion characteristics; and (6) Validation results should provide a clear, transparent, and quantitative assessment of the simulated ground motions.

Conference Paper↗

USGS2022

Introduction to the special issue of the Consortium of Organizations for Strong Motion Observation Systems (COSMOS) international guidelines for applying noninvasive geophysical techniques to characterize seismic site conditions

Knowledge about local seismic site conditions provides critical information to account for site effects that are commonly observed in strong motion recordings. Certainly, other wave propagation effects can influence these observations, which are attributable to variations in material properties of the paths traveled by the waves, as well as the characteristics of the seismic source. However, local geologic conditions, particularly, when under shear-wave excitation, are known to have a strong influence on the behavior of ground shaking in the frequency range that is expected to directly affect the built environment. Thus, shear waves traveling in the shallow subsurface—defined here as tens to hundreds of meters beneath the ground surface—are the main foci for application and research in the earthquake engineering community.

Journal of Seismology↗

USGS2022

Flexible multimethod approach for seismic site characterization

We describe the flexible multimethod seismic site characterization technique for obtaining shear-wave velocity ( V S ) profiles and derivative information, such as the time-averaged V S of the upper 30 m ( V S30 ). Simply stated, the multimethod approach relies on the application of multiple independent noninvasive site characterization acquisition and analysis techniques utilized in a flexible field-based approach through a decision process based on primary factors such as the investigator’s available equipment, field logistics (e.g., available array deployment space, urban versus rural), geologic setting, and funding among other primary factors. In a multimethod approach, sites ideally should be characterized using both active and passive noninvasive (i.e., no drilling and only minimal disturbance to the ground surface) methods. Almost without exception, we recommend the use of at least one active-source technique for analyzing surface waves, which in the current state of practice would include one or more of the following: spectral analysis of surface waves (SASW; commonly Rayleigh waves) and multichannel analysis of surface waves (MASW; Rayleigh and/or Love waves). In addition, passive-source microtremor array methods (MAMs) using linear (one-dimensional) and two-dimensional arrays may be appropriate or even required for characterizing a given site. Similarly, the microtremor horizontal-to-vertical spectral ratio (mHVSR) method can be valuable for identifying shallow rock sites. Finally, P- and SH-wave seismic refraction methods warrant consideration, especially at rock and shallow-rock sites.

Journal of Seismology↗

USGS2022

A review of the microtremor horizontal-to-vertical spectral ratio (MHVSR) method

The single-station microtremor horizontal-to-vertical spectral ratio (MHVSR) method was initially proposed to retrieve the site amplification function and its resonance frequencies produced by unconsolidated sediments overlying high-velocity bedrock. Presently, MHVSR measurements are predominantly conducted to obtain an estimate of the fundamental site frequency at sites where a strong subsurface impedance contrast exists. Of the earthquake site characterization methods presented in this special issue, the MHVSR method is the furthest behind in terms of consensus towards standardized guidelines and commercial use. The greatest challenges to an international standardization of MHVSR acquisition and analysis are (1) the what — the underlying composition of the microtremor wavefield is site-dependent, and thus, the appropriate theoretical (forward) model for inversion is still debated; and (2) the how — many factors and options are involved in the data acquisition, processing, and interpretation stages. This paper reviews briefly a historical development of the MHVSR technique and the physical basis of an MHVSR (the what ). We then summarize recommendations for MHVSR acquisition and analysis (the how ). Specific sections address MHVSR interpretation and uncertainty assessment.

Journal of Seismology↗

USGS2022

An assessment of uncertainties in VS profiles obtained from microtremor observations in the phased 2018 COSMOS blind trials

Site response is a critical consideration when assessing earthquake hazards. Site characterization is key to understanding site effects as influenced by seismic site conditions of the local geology. Thus, a number of geophysical site characterization methods were developed to meet the demand for accurate and cost-effective results. As a consequence, a number of studies have been administered periodically as blind trials to evaluate the state-of-practice on-site characterization. We present results from the Consortium of Organizations for Strong Motion Observation Systems (COSMOS) blind trials, which used data recorded from surface-based microtremor array methods (MAM) at four sites where geomorphic conditions vary from deep alluvial basins to an alpine valley. Thirty-four invited analysts participated. Data were incrementally released to 17 available analysts who participated in all four phases: (1) two-station arrays, (2) sparse triangular arrays, (3) complex nested triangular or circular arrays, and (4) all available geological control site information including drill hole data. Another set of 17 analysts provided results from two sites and two phases only. Although data from one site consisted of recordings from three-component sensors, the other three sites consisted of data recorded only by vertical-component sensors. The sites cover a range of noise source distributions, ranging from one site with a highly directional microtremor wave field to others with omni-directional (azimuthally distributed) wave fields. We review results from different processing techniques (e.g., beam-forming, spatial autocorrelation, cross-correlation, or seismic interferometry) applied by the analysts and compare the effectiveness between the differing wave field distributions. We define the M index as a quality index based on estimates of the time-averaged shear-wave velocity of the upper 10 ( V S10 ), 30 ( V S30 ), 100 ( V S100 ), and 300 ( V S300 ) meters and show its usefulness in quantitative comparisons of V S profiles from multiple analysts. Our findings are expected to aid in building an evidence-based consensus on preferred cost-effective arrays and processing methodology for future studies of seismic site effects.

Seismological Research Letters↗

USGS2022

COSMOS Ground-Motion Simulation Working Group workshops #1 and #2

These 2 workshops were held in response to interest generated from sessions on the use of simulated earthquake ground motions at the 2020 and 2021 Consortium of Organizations for Strong Motion Observation Systems (COSMOS) Technical Sessions. The discussions at the Technical Sessions highlighted desires to promote the use of simulated earthquake ground motions for engineering applications and the need to coordinate efforts to validate and disseminate the data. Our first workshop focused on curating and disseminating simulated earthquake ground-motion data. Our second workshop focused on validating simulated ground-motions for engineering applications. Both workshops included a few invited presentations from providers (ground-motion simulators) and engineering user perspectives followed by open discussions. About 100 people participated in each of the workshops.

Report↗

USGS2021

Relational database for horizontal‐to‐vertical spectral ratios

Frequency‐dependent horizontal‐to‐vertical spectral ratios (HVSRs) of Fourier amplitudes from three‐component recordings can provide useful information for site response modeling. However, such information is not incorporated into most ground‐motion models, including those from Next‐Generation Attenuation projects, which instead use the time‐averaged shear‐wave velocity ( ⁠ V S "> V S VS ⁠ ) in the upper 30 m of the site and sediment depth terms. To facilitate utilization of HVSR, we developed a publicly accessible relational database. This database is adapted from a similar repository for V S "> V S VS data and provides microtremor‐based HVSR data (mHVSR) and supporting metadata, but not parameters derived from the data. Users can interact with the data directly within a web portal that contains a graphical user interface (GUI) or through external tools that perform cloud‐based computations. Within the database GUI, the median horizontal‐component mHVSR can be plotted against frequency, with the mean and mean ± one standard deviation (representing variability across time windows) provided. Using external interactive tools (provided as a Jupyter Notebook and an R script), users can replot mHVSR (as in the database) or create polar plots. These tools can also derive parameters of potential interest for modeling purposes, including a binary variable indicating whether an mHVSR plot contains peaks, as well as the fitted properties of those peaks (frequencies, amplitudes, and widths). Metadata are also accessible, which includes site location, details about the instruments used to make the measurements, and data processing information related to windowing, antitrigger routines, and filtering.

Seismological Research Letters↗

USGS2021

A probabilistic framework to model distributions of VS30

The time‐averaged shear‐wave velocity in the upper 30 m depth from the ground surface, or V S 30 "> V S 30 ⁠ , is often used as a predictor to describe local site effects in ground‐motion models. Although V S 30 "> V S 30 is typically determined from in situ measurements, it is not always feasible to obtain such measurements due to project restrictions or site accessibility. This motivates the development and use of proxy‐based V S 30 "> V S 30 predictions that leverage more readily available secondary information such as surface geology, topographic slope, or geomorphic terrain classes to estimate the mean V S 30 "> V S 30 and associated uncertainty. Traditionally, empirical distributions of V S 30 "> V S 30 have been observed to have long right tails, leading to high levels of associated uncertainty. In this study, we present a physical framework that is grounded in fundamental principles of geostatistics and probability to explain the uncertainty and skewness associated with V S 30 "> V S 30 measurements. Specifically, by invoking Lyapunov’s central limit theorem, we hypothesize that the distribution of V S 30 "> V S 30 can be theoretically approximated by a reciprocal–normal distribution. We show that a non‐normal and skewed distribution of V S 30 "> V S 30 is to be expected and is not a sign of measurement error or sampling bias, although sampling bias can exaggerate such skewness. Our framework also enables us to propose the mode as a characteristic value of V S 30 "> V S 30 measurements, as opposed to the mean or median, which can overestimate the most probable value.

Bulletin of the Seismological Society of America↗

USGS2021

Horizontal-to-vertical spectral ratios from California sites: Open-source database and data interpretation to establish site parameters

Frequency-dependent horizontal-to-vertical spectral ratios (HVSR) of Fourier amplitudes from three-component recordings can provide information on one or more site resonant frequencies and relative levels of amplification at those frequencies. Such information is potentially useful for predicting site amplification but is not present in site databases that have been developed over the last 15–20 years for the Next-Generation Attenuation (NGA) projects, which instead use the time-averaged shear-wave velocity (VS) in the upper 30 m of the site (VS30) as the primary site parameter and are supplemented with basin depth terms where available. As a consequence, HVSR parameters are also not used in NGA ground motion models. In order for HVSR-based parameters to be used in future versions of site databases, a publicly accessible repository of this information is needed. We adapt a relational database developed to archive and disseminate VS data to also include HVSR. The database provides relevant microtremor-based HVSR data (mHVSR) and supporting metadata. We consider the most relevant data to be the frequency-dependent mHVSR, where the horizontal is taken as the median component and also as a function of horizontal azimuth (referred to as polar plots). Relevant metadata includes site location information, details about the equipment used to make the measurements, and processing details related to windowing, anti-trigger routines, and filtering. We describe the database schema developed to organize and present this information.

California↗

USGS2020

VS30 and Dominant Site Frequency (⁠fd⁠) as Provisional Station ML Corrections (⁠dML⁠) in California

New seismic stations added to a regional seismic network cannot be used to calculate local magnitude ( ⁠ M L "> M L ⁠ ) until a revised regionwide amplitude decay function is developed. Each station must record a minimum number of local and regional earthquakes that meet specific amplitude requirements prior to recalibration of the amplitude decay function. Station component adjustments ( ⁠ d M L "> d M L ⁠ ; Uhrhammer et al. , 2011 ) are then calculated after inverting for a new regional amplitude decay function, constrained by the sum of d M L "> d M L ⁠ for long‐running stations. Therefore, there can be significant delay between when a new station starts contributing real‐time waveform packets and when data can be included in magnitude determinations. We propose the use of known estimates of seismic site conditions such as the time‐averaged shear‐wave velocity ( ⁠ V S "> V S ⁠ ) of the upper 30 m ( ⁠ V S 30 "> V S 30 ⁠ ) and the site dominant frequency ( ⁠ f d "> f d ⁠ ) to calculate d M L "> d M L ⁠ ⁠ . Previously established d M L "> d M L ⁠ ⁠ , measured V S 30 "> V S 30 ⁠ ⁠ , and f d "> f d data are available for between 126 and 458 horizontal components (east–west and north–south) at 81 seismic stations in the California Integrated Seismic Network; d M L "> d M L ⁠ d M L "> data range from −1.10 to 0.39, V S 30 "> V S 30 ⁠ values range from 202 to 1464    m / s "> 1464 m / s ⁠ , and 440 f d "> f d values are compiled from earthquake and microseismic records that range from 0.13 to 21 Hz. We find V S 30 "> V S 30 ⁠ and d M L "> d M L ⁠ exhibit a positive coefficient of determination ( ⁠ R = 0.59 "> R = 0.59 ⁠ ), indicating that as V S 30 "> V S 30 ⁠ increases, d M L "> d M L ⁠ increases. This implies that greater site amplification (lower V S 30 "> V S 30 ⁠ ⁠ ) results in smaller d M L "> d M L ⁠ ⁠ . f d "> f d and d M L "> d M L ⁠ also generally exhibit a positive correlation ( ⁠ R 2 < 0.56 "> R 2 < 0.56 ⁠ ), which implies lower d M L "> d M L values are related to site resonance at depth‐dependent frequencies. Using the developed relationships, V S 30 "> V S 30 ⁠ or f d "> f d measurements can be used to establish a provisional d M L "> d M L ⁠ for newly established stations. This procedure allows new stations to contribute to regional network M L "> M L determinations immediately without the need to wait until a minimum set of earthquake data has been recorded.

California↗

USGS2020

The U.S. Geological Survey’s Rapid Seismic Array Deployment for the 2019 Ridgecrest Earthquake Sequence

Rapid seismic deployments following large earthquakes capture ephemeral near‐field recordings of aftershocks and ambient noise that can provide valuable data for seismological studies. The U.S. Geological Survey installed 19 temporary seismic stations following the 4 July 2019 M w 6.4 and 6 July 2019 (UTC) M w 7.1 earthquakes near the city of Ridgecrest, California. The stations record the aftershock sequence beginning two days after the mainshock and are expected to remain in the field through approximately January 2020. The deployment augments the permanent seismic network in the area to improve azimuthal coverage and provide additional near‐field observations. This article summarizes the motivation and goals of the deployment; details of station installation, instrumentation, and configurations; and initial data quality and observations from the network. We expect these data to be useful for a range of studies including detailing near‐field variability in strong ground motions, determining stress drops and rupture directivity of small events, imaging the fault zone, documenting the evolution of crustal properties within and outside of the fault zone, and others.

California↗

USGS2019

Length-scale-dependent relationships between VS30 and topographic slopes in southern California

In engineering seismology, the time‐averaged shear‐wave velocity ( ⁠ V S ⁠ ) of the upper 30 m of the crust ( ⁠V S30⁠ ) is the primary parameter used in ground‐motion models to predict seismic site effects. ⁠V S30⁠ is typically derived from in situ recordings of V S ⁠ , although proxy‐based approaches (using geologic and/or geomorphometric classifications) are provisionally adopted when measurement‐based ⁠V S30⁠ are sparse or not readily available. Despite the acceptance of proxy approaches, there are limited studies that examine the empirical relationships between ⁠V S30⁠ and topographic attributes measured from various length scales and different resolutions of the digital elevation model. In this study, we examine the relationships between compiled ⁠V S30⁠ measurements from 218 sites in southern California and topographic metrics of slope and relief measured over various length scales. We find that the correlations between topographic metrics and ⁠V S30⁠ are weak but statistically significant. The correlations are improved when topographic slopes and relief are measured over length scales longer than typical hillslopes and ⁠V S30⁠ sites are separated by different geologic groups. This is likely because ⁠V S30⁠ ⁠ , especially on the rock sites, is better reflected in topographic metrics that capture large‐scale topographic relief, as well as landscape positions such as hilltops and valley bottoms.

California↗

USGS2019

Comparison of site dominant frequency from earthquake and microseismic data in California

An important predictive variable for site amplification is the site dominant frequency ( ⁠ƒ d ⁠ ). At seismic monitoring stations, ƒ d can be calculated from the peak of the horizontal‐to‐vertical spectral ratios (HVSRs) obtained from earthquake recordings (eHVSR). For other sites, ƒ d can be estimated from microseismic (mHVSR) observations. We compare the ƒ d values derived from eHVSR (5% damped response spectra from the Next Generation Attenuation‐West2 [NGA‐West2] database; Ancheta et al. , 2014 ) with those derived from mHVSR (Fourier spectra from Yong et al. , 2013 ) for seismic stations in California. We show that the logarithm of eHVSR ƒ d scales linearly with the logarithm of mHVSR ƒ d ⁠ , with a standard deviation of 0.14log 10 units for mHVSR ƒ d larger than 0.2 Hz. The relationship holds for microseismic surveys at distances up to 300 m away from the seismic stations. The results of this study have beneficial implications for the characterization of site response in modern ground‐motion models as well as in building codes.

California↗

USGS2018

Development of a United States community shear wave velocity profile database

We present the details of a multi-institutional effort to develop an open-access shear-wave velocity ( V S ) profile database (PDB), which will include a public repository for V S profile data in the United States. V S profiles are an essential resource for ground motion modeling and other applications. The minimum requirements for a site to be included in the database are in situ geophysical V S measurements and location metadata (geodetic coordinates and elevation). Other information is included as available, including geotechnical logs, penetration resistance, laboratory test data, ground water elevation, and P-wave velocity profiles. The project is currently at the stage of data collection (over 4500 V S profiles) and prototype data model development. The database will be presented as an online map-based interface with downloadable V S profile and metadata information. This paper serves as a progress report to the geotechnical and earthquake engineering communities, as we seek community engagement and support.

Conference Paper↗