An open-source site database of strong-motion stations in Japan: K-NET and KiK-net (v1.0.0)

2021 ◽  
pp. 875529302098802
Author(s):  
Chuanbin Zhu ◽  
Graeme Weatherill ◽  
Fabrice Cotton ◽  
Marco Pilz ◽  
Dong Youp Kwak ◽  
...  

This article describes an open-source site database for a total number of 1742 earthquake recording sites in the K-NET (Kyoshin network) and KiK-net (Kiban Kyoshin network) networks in Japan. This database contains site characterization parameters directly derived from available velocity profiles, including average wave velocities, bedrock depths, and velocity contrast. Meanwhile, it also consists of earthquake horizontal-to-vertical spectral ratio (HVSR) and peak parameters, for example, peak frequency, amplitude, width, and prominence. In addition, the site database also comprises topographic and geological proxies inferred from regional models or maps. Each parameter is derived in a consistent manner for all sites. This site database can benefit the application of machine learning techniques in studies on site amplification. Besides, it can facilitate, for instances, the search of the optimal site parameter(s) for the prediction of site amplification, the development and testing of ground-motion models or methodologies, as well as investigations on spatial or regional variability in site response. All resources (the site database, earthquake HVSR data at all sites, and the MATLAB script for peak identification) can be freely accessed via: https://doi.org/10.5880/GFZ.2.1.2020.006

2021 ◽  
Author(s):  
Olga-Joan Ktenidou ◽  
Faidra Gkika ◽  
Erion-Vasilis Pikoulis ◽  
Christos Evangelidis

<p>Although it is nowadays desirable and even typical to characterise site conditions in detail at modern recording stations, this is not yet a general rule in Greece, due to the large number and geographical dispersion of stations. Indeed, most of them are still characterised merely through geological descriptions or proxy-based parameters, rather than through in-situ measurements. Considering: 1. the progress made in recent years with sophisticated ground motion models and the need to define region-specific rock conditions based on data, 2. the move towards large open-access strong-motion databases that require detailed site metadata, and 3. that Greek-provenance recordings represent a significant portion of European seismic data, there are many reasons to improve our understanding of site response at these stations. Moreover, it has been shown recently in several regions that even sites considered as rock can exhibit amplification and ground motion variability, which has given rise to more scientific research into the definition of reference sites. For Greece, in-situ-characterisation campaigns for the entire network would impose unattainable time/budget constraints; so, instead, we implement alternative empirical approaches using the recordings themselves, such as the horizontal-to-vertical spectral ratio technique and its variability. We present examples of 'well-behaved', typical rock sites, and others whose response diverges from what is assumed for their class.</p><p> </p>


Author(s):  
Paolo Bergamo ◽  
Conny Hammer ◽  
Donat Fäh

ABSTRACT We address the relation between seismic local amplification and topographical and geological indicators describing the site morphology. We focus on parameters that can be derived from layers of diffuse information (e.g., digital elevation models, geological maps) and do not require in situ surveys; we term these parameters as “indirect” proxies, as opposed to “direct” indicators (e.g., f0, VS30) derived from field measurements. We first compiled an extensive database of indirect parameters covering 142 and 637 instrumented sites in Switzerland and Japan, respectively; we collected topographical indicators at various spatial extents and focused on shared features in the geological descriptions of the two countries. We paired this proxy database with a companion dataset of site amplification factors at 10 frequencies within 0.5–20 Hz, empirically measured at the same Swiss and Japanese stations. We then assessed the robustness of the correlation between individual site-condition indicators and local response by means of statistical analyses; we also compared the proxy-site amplification relations at Swiss versus Japanese sites. Finally, we tested the prediction of site amplification by feeding ensembles of indirect parameters to a neural network (NN) structure. The main results are: (1) indirect indicators show higher correlation with site amplification in the low-frequency range (0.5–3.33 Hz); (2) topographical parameters primarily relate to local response not because of topographical amplification effects but because topographical features correspond to the properties of the subsurface, hence to stratigraphic amplification; (3) large-scale topographical indicators relate to low-frequency response, smaller-scale to higher-frequency response; (4) site amplification versus indirect proxy relations show a more marked regional variability when compared with direct indicators; and (5) the NN-based prediction of site response is the best achieved in the 1.67–5 Hz band, with both geological and topographical proxies provided as input; topographical indicators alone perform better than geological parameters.


2020 ◽  
Vol 36 (1) ◽  
pp. 87-110 ◽  
Author(s):  
Youssef M. A. Hashash ◽  
Okan Ilhan ◽  
Behzad Hassani ◽  
Gail M. Atkinson ◽  
Joseph Harmon ◽  
...  

This article evaluates linear simulation-based and empirical site amplification models including site natural period dependency parameters to account for the distinctive amplification behavior near site fundamental frequencies resulting from the sharp impedance contrast between soil and underlying hard bedrock in central and eastern North America (CENA). The simulation-based amplification models are developed using 581,685 frequency-domain linear analyses generated from a parametric study and include VS30-scaling and site natural period ( Tnat) parameters. The empirical models are derived from residuals analyses of ground-motion models for two reference conditions: B/C boundary ( VS30 = 760 m/s) and CENA hard-rock condition ( VS = 3000 m/s). The simulation-based and empirical models are compared for 8 site profiles in CENA to measured horizontal-to-vertical (H/V) component response spectral (RS) ratios, the mean of linear simulations for similar sites, and one-dimensional (1D) linear site response analysis for four of these sites. Comparisons between observed and estimated site amplification behaviors highlight model dependency on Tnat in CENA. Model consistencies and differences related to the distinct linear amplification features near site fundamental frequency are discussed.


2020 ◽  
Vol 110 (6) ◽  
pp. 2892-2911
Author(s):  
Eri Ito ◽  
Kenichi Nakano ◽  
Fumiaki Nagashima ◽  
Hiroshi Kawase

ABSTRACT The main purpose of the site classification or velocity determination at a target site is to obtain or estimate the horizontal site amplification factor (HSAF) at that site during future earthquakes because HSAF would have significant effects on the strong-motion characteristics. We have been investigating various kinds of methods to delineate the S-wave velocity structures and the subsequent HSAF, as precisely as possible. After the advent of the diffuse field concept, we have derived a simple formula based on the equipartitioned energy density observed in the layered half-space for incident body waves. In this study, based on the diffuse field concept, together with the generalized spectral inversion technique (GIT), we propose a method to directly estimate the HSAF of the S-wave portion from the horizontal-to-vertical spectral ratio of earthquakes (eHVSRs). Because the vertical amplification is included in the denominator of eHVSR, it cannot be viewed as HSAF without correction. We used GIT to determine both the HSAF and the vertical site amplification factor (VSAF) simultaneously from strong-motion data observed by the networks in Japan and then deduced the log-averaged vertical amplification correction function (VACF) from VSAFs at a total of 1678 sites in which 10 or more earthquakes have been observed. The VACF without a category has a constant amplitude of about 2 in the frequency range from 1 to 15 Hz. By multiplying eHVSR by VACF, we obtained the simulated HSAF. We verified the effectiveness of this correction method using data from observation sites not used in the aforementioned averaging in the frequency range from 0.12 to 15 Hz.


2019 ◽  
Vol 35 (2) ◽  
pp. 849-881 ◽  
Author(s):  
Grace A. Parker ◽  
Jonathan P. Stewart ◽  
Youssef M. A. Hashash ◽  
Ellen M. Rathje ◽  
Kenneth W. Campbell ◽  
...  

We present empirical linear site amplification models conditioned on time-averaged shear wave velocity in the upper 30 m ( VS30) for central and eastern North America. The models are derived from ground motion data and site condition information from the NGA-East project and are intended for use with reference rock ground motion models. Site amplification is found to scale with VS30 for intermediate to stiff site conditions ( VS30 > 300 m/s) in a weaker manner than for active tectonic regions such as the western United States. For stiff sites ( >800 m/s), we find differences in site amplification for previously glaciated and nonglaciated regions, with nonglaciated sites having lower amplification. The models were developed using a combination of least-squares, mixed effects, and Bayesian techniques; the latter show that accounting for predictor uncertainty does not appreciably affect the median model but decreases model dispersion. Our VS30-scaling models are modular and additive to simulation-based models for the nonlinear components of site response. A limitation of the present models is that they do not account for site-specific resonance effects.


2021 ◽  
Author(s):  
Claudia Mascandola ◽  
Giovanni Lanzano ◽  
Francesca Pacor

<p>The rapid increase of seismic waveforms, due to the increment of seismic stations and continuous real-time streaming to data centres, leads to the need for automatic procedures aimed at supporting data processing and data quality control. In this study, we propose a semi-automatic procedure for the consistency check of large strong-motion datasets, classifying the anomalies observed on the residuals analysis and identifying the possible causes.</p><p>The data collected in the strong-motion databases are usually arranged as parametric tables (called flatfiles), used to disseminate the Intensity Measures (IMs) and the associated metadata of the processed waveforms. This is the current practice for the ITalian ACcelerometric Archive (ITACA, D’Amico et al., 2020) and Engineering Strong Motion (ESM; Lanzano et al. 2019a) databases. The adopted criteria for flatfile compilation are designed to collect IMs and related metadata in a uniform, updated, and traceable way, with the aim of providing datasets useful to develop Ground Motion Models (GMMs) for Probabilistic Seismic Hazard Assessment (PSHA) and engineering applications. Therefore, the consistency check of the flatfiles is a crucial task to improve the quality of the products provided by the waveform services.</p><p>The proposed procedure is based on the residual distributions obtained from ad-hoc ground motion prediction equations for the ordinates of the 5% damped acceleration response spectra. In this study, we focus on the active shallow crust events in ITACA, considering the ITA18 ground motion model (Lanzano et al., 2019b) as a reference for Italy. The total residuals, computed as logarithm difference between observations and predictions, are decomposed in between-event, between-station and event-and-station corrected residuals by applying a mixed-effect regression (Bates et al., 2015). This is the common practice for the (partial) removal of the ergodic assumption in empirical GMMs (e.g., Stafford 2014), where the contribution of the systematic corrective effects of event and station on aleatory variability are identified and shifted to the epistemic uncertainty. Afterward, the proposed procedure is applied to raise a warning in case of anomalous residual values. Warnings are provided when the normalized residuals exceed a certain threshold, in three ranges of periods (i.e., 0.01-0.15 s, 0.15-1 s, 1-5 s). The causes of warnings may be several and may concern the event, the site, the waveform, or a combination of them. Among the possible sources of anomalous trends, the more common are: preliminary or inaccurate event localization or magnitude, wrong soil category assigned based on proxies, misleading tectonic regime assigned to the earthquake, and fault directivity that may cause strong-ground motion amplification in certain directions. Warnings may also raise for peculiarities in the site-response (e.g., large amplifications/de-amplifications at certain frequency-bands) and to the occurrence of near-source effects in the waveforms (see Pacor et al., 2018). Based on the raised warnings, a decision tree classifier is developed to identify the common anomaly sources and to support the consistency check of the semi-automatic procedure.</p><p>This study may help to enhance the waveform services and related products, besides reducing the variability of ground motion models and guiding decisions for site characterization studies and network maintenance.</p>


2019 ◽  
Vol 36 (1) ◽  
pp. 111-137 ◽  
Author(s):  
Boqin Xu ◽  
Ellen M Rathje ◽  
Youssef Hashash ◽  
Jonathan Stewart ◽  
Kenneth Campbell ◽  
...  

Small-strain damping profiles developed from geotechnical laboratory testing have been observed to be smaller than the damping inferred from the observed site amplification from downhole array recordings. This study investigates the high-frequency spectral decay parameter ( κ0) of earthquake motions from soil sites and evaluates the use of κ0 to constrain the small-strain damping profile for one-dimensional site response analysis. Using data from 51 sites from the Kiban-Kyoshin strong motion network (KiK-net) array in Japan and six sites from California, a relationship was developed between κ0 at the surface and both the 30-m time-averaged shear wave velocity ( V s30) and the depth to the 2.5 km/s shear wave velocity horizon ( Z2.5). This relationship demonstrates that κ0 increases with decreasing V s30 and increasing Z2.5. An approach is developed that uses this relationship to establish a target κ0 from which to constrain the small-strain damping profile used in one-dimensional site response analysis. This approach to develop κ0-consistent damping profiles for site response analysis is demonstrated through a recent site amplification study of Central and Eastern North America for the NGA-East project.


2006 ◽  
Vol 22 (1) ◽  
pp. 169-188 ◽  
Author(s):  
Sheri Molnar ◽  
John F. Cassidy

The applicability of the microtremor spectral ratio method is examined by comparing microtremor and weak-motion earthquake site responses at seven permanent strong-motion sites in Victoria, British Columbia. For each site, a weak-motion earthquake standard spectral ratio (bedrock reference), the average horizontal-to-vertical spectral ratio of up to five weak-motion earthquakes, and the average microtremor (Nakamura method) spectral ratio are compared. The geologic setting of Victoria is ideal for site response studies with a near-surface high impedance contrast between thin geologic layers of Victoria clay (about 11 m maximum in this study) and Pleistocene till or bedrock. Regardless of excitation source (weak-motion earthquakes or microtremors) and spectral ratio method, similar peak amplitudes and fundamental frequencies were found. Thicker material (>10 m) sites displayed higher peak amplitudes (up to six times amplification) at frequencies of 2–5 Hz compared to sites with a thin lens of material (<3 m) over bedrock that showed peak amplitudes at frequencies of >8 Hz.


Author(s):  
Ruibin Hou ◽  
John X. Zhao

ABSTRACT This article presents a nonlinear site amplification model for ground-motion prediction equations (GMPEs), using site period as site-effect proxy based on the measured shear-wave velocity profiles of selected KiK-net and K-NET sites in Japan. This model was derived using 1D equivalent-linear site-response analysis for a total of 516 measured soil-site shear-wave velocity profiles subjected to a total of 912 components of rock-site records. The modulus reduction and damping curves for each soil layer were assigned based on the soil-type description for a particular layer. The site period and site impedance ratio affect both the linear and nonlinear parts of this study, and were used as the site parameters in the 1D amplification model. A large impedance ratio enhances the amplification ratios when the site responds elastically and enhances the nonlinear response when the site develops a significant nonlinear response. The effects of moment magnitude and source distance on the linear part of the 1D amplification model were also incorporated in the model. To implement the 1D amplification model into GMPEs, a model adjustment is required to match the GMPE amplification ratio at weak motion and to retain the nonlinear amplification ratio at the strong motion of the 1D model. The two-step adjustment method by Zhao, Hu, et al. (2015) was adopted in this study with significant modifications. It is not possible to obtain a credible second-step adjustment parameter using the GMPEs dataset only. We proposed three methods for calculating the scale factors. Method 1 is a constant angle in a 30°–60° range for all spectral periods; method 2 was based on the GMPE dataset and 1-D model parameters; and method 3 was based on the strong-motion records used for the 1D site modeling. A simple second-step adjustment factor leads to smoothing amplification ratios and soil-site spectrum.


Sign in / Sign up

Export Citation Format

Share Document