Two Empirical Double-Corner-Frequency Source Spectra and Their Physical Implications

Author(s):  
Chen Ji ◽  
Ralph J. Archuleta

ABSTRACT We introduce double-corner-frequency (DCF) source spectral models JA19 and JA19_2S, which, in conjunction with a stochastic ground-motion model, can reproduce the mean peak ground acceleration (PGA) and mean peak ground velocity (PGV) of the Next Generation Attenuation-West 2 database for magnitudes 3.3–7.3. Their displacement amplitude spectrum remains constant for frequencies less than fc1, decays as f−1 between fc1 and fc2, and decays as f−2 for frequencies greater than fc2. The model JA19 is self-similar. Its two corner frequencies fc1 and fc2 scale with moment magnitude (M) as (1) log(fc1(M))=1.754−0.5M and (2) log(fc2(M))=3.250−0.5M. We find that relation (1) is consistent with the known self-similar scaling relations of the rupture duration (Td), in which Td=1/(πfc1). Relation (2) may reflect the scaling relation of the average rise time (TR), where TR∼0.8/(fc2). Stochastic simulations of ground motion using JA19 cannot reproduce the sharp change in magnitude dependence of PGA and PGV at M 5.3, suggesting a breakdown of self-similarity. The magnitude dependence of PGA and PGV and this change in slope is well explained by JA19_2S, which results from perturbing the fc1 scaling relationship in JA19. For JA19_2S: log(fc1(M))=1.474−0.415M for M≤5.3; log(fc1(M))=2.375−0.585M for M>5.3. The scaling relation for fc2 is unchanged. When fc1≪fc2, the scaled energy (ratio of radiated energy and seismic moment) scales with M0fc12fc2. The scaled energy of JA19 is 2.2×10−5, independent of magnitude. Because JA19_2S is not self-similar, its scaled energy is 2.2–4.7×10−5, increasing 2.2 times, when magnitude increases from 3.3 to 5.3, and, subsequently decreasing 2.2 times, as magnitude further increases from 5.3 to 7.3. Both agree with the global average (∼3×10−5) reported previously. Using our proposed empirical models, the standard deviation of average static stress drop from seismological studies can be significantly greater than the standard deviation of the stress parameter used to estimate PGA and PGV.

1983 ◽  
Vol 73 (1) ◽  
pp. 83-96 ◽  
Author(s):  
Michel Campillo ◽  
Michel Bouchon

abstract We present a study of the seismic radiation of a physically realistic source model—the circular crack model of Madariaga—at close distance range and for vertically heterogeneous crustal structures. We use this model to represent the source of small strike-slip earthquakes. We show that the characteristics of the radiated seismic spectra, like the corner frequency, are strongly affected by the presence of the free surface and by crustal layering, and that they can be considerably different from the ones of the homogeneous-medium far-field solution. The vertical and radial displacement spectra are the most strongly affected. We use this source model to calculate the decay of peak ground velocity with epicentral distance and source depth for small strike-slip earthquakes in California. For distances between 10 and 80 km, the peak horizontal velocity decay is of the form r−1.25 for a 4-km hypocentral depth and r−1.65 for deeper sources. The predominance of supercritically reflected arrivals beyond epicentral distances of 70 to 80 km produces a sharp change in the rate of decay of the ground motion. For most of the cases considered, the peak ground velocity increases between 80 and 100 km. We also show that the S-wave velocity in the source layer is the lower limit of phase velocities associated with significant ground motion.


2021 ◽  
Author(s):  
Grigorios Lavrentiadis ◽  
Norman A. Abrahamson ◽  
Nicolas M. Kuehn

Abstract A new non-ergodic ground-motion model (GMM) for effective amplitude spectral (EAS) values for California is presented in this study. EAS, which is defined in Goulet et al. (2018), is a smoothed rotation-independent Fourier amplitude spectrum of the two horizontal components of an acceleration time history. The main motivation for developing a non-ergodic EAS GMM, rather than a spectral acceleration GMM, is that the scaling of EAS does not depend on spectral shape, and therefore, the more frequent small magnitude events can be used in the estimation of the non-ergodic terms. The model is developed using the California subset of the NGAWest2 dataset Ancheta et al. (2013). The Bayless and Abrahamson (2019b) (BA18) ergodic EAS GMM was used as backbone to constrain the average source, path, and site scaling. The non-ergodic GMM is formulated as a Bayesian hierarchical model: the non-ergodic source and site terms are modeled as spatially varying coefficients following the approach of Landwehr et al. (2016), and the non-ergodic path effects are captured by the cell-specific anelastic attenuation attenuation following the approach of Dawood and Rodriguez-Marek (2013). Close to stations and past events, the mean values of the non-ergodic terms deviate from zero to capture the systematic effects and their epistemic uncertainty is small. In areas with sparse data, the epistemic uncertainty of the non-ergodic terms is large, as the systematic effects cannot be determined. The non-ergodic total aleatory standard deviation is approximately 30 to 40% smaller than the total aleatory standard deviation of BA18. This reduction in the aleatory variability has a significant impact on hazard calculations at large return periods. The epistemic uncertainty of the ground motion predictions is small in areas close to stations and past event.


2020 ◽  
Vol 224 (3) ◽  
pp. 1753-1779
Author(s):  
Marta Pischiutta ◽  
Aybige Akinci ◽  
Elisa Tinti ◽  
André Herrero

SUMMARY On 24 August 2016 at 01:36 UTC a ML6.0 earthquake struck several villages in central Italy, among which Accumoli, Amatrice and Arquata del Tronto. The earthquake was recorded by about 350 seismic stations, causing 299 fatalities and damage with macroseismic intensities up to 11. The maximum acceleration was observed at Amatrice station (AMT) reaching 916 cm s–2 on E–W component, with epicentral distance of 15 km and Joyner and Boore distance to the fault surface (RJB) of less than a kilometre. Motivated by the high levels of observed ground motion and damage, we generate broad-band seismograms for engineering purposes by adopting a hybrid method. To infer the low frequency seismograms, we considered the kinematic slip model by Tinti et al . The high frequency seismograms were produced using a stochastic finite-fault model approach based on dynamic corner-frequency. Broad-band synthetic time-series were therefore obtained by merging the low and high frequency seismograms. Simulated hybrid ground motions were compared both with the observed ground motions and the ground-motion prediction equations (GMPEs), to explore their performance and to retrieve the region-specific parameters endorsed for the simulations. In the near-fault area we observed that hybrid simulations have a higher capability to detect near source effects and to reproduce the source complexity than the use of GMPEs. Indeed, the general good consistency found between synthetic and observed ground motion (both in the time and frequency domain), suggests that the use of regional-specific source scaling and attenuation parameters together with the source complexity in hybrid simulations improves ground motion estimations. To include the site effect in stochastic simulations at selected stations, we tested the use of amplification curves derived from HVRSs (horizontal-to-vertical response spectra) and from HVSRs (horizontal-to-vertical spectral ratios) rather than the use of generic curves according to NTC18 Italian seismic design code. We generally found a further reduction of residuals between observed and simulated both in terms of time histories and spectra.


1987 ◽  
Vol 58 (2) ◽  
pp. 53-58 ◽  
Author(s):  
Otto W. Nuttli ◽  
David S. Bowling ◽  
J. E. Lawson ◽  
Randall Wheeler

Abstract Strong-motion records from a velocity meter were recorded at an epicentral distance of 4.5 km from the January 12, 1984 eastern Missouri earthquake of mb = 3.0. Peak values of ground velocity, associated with only one or two wave cycles, are: transverse component, 0.18 cm/sec; radial component, 0.16 cm/sec; vertical component, 0.12 cm/sec. The levels of the sustained motion, which extends from the onset of S to about 0.4 sec later, are: transverse component, 0.064 cm/sec; radial component, 0.061 cm/sec; vertical component, 0.061 cm/sec. The data are consistent with a spectral scaling relation assuming either a 3.5 or 4.0 slope of the logarithm of the seismic moment versus the logarithm of the corner frequency, but cannot be used to choose between the two relations.


Author(s):  
Fabio Sabetta ◽  
Antonio Pugliese ◽  
Gabriele Fiorentino ◽  
Giovanni Lanzano ◽  
Lucia Luzi

AbstractThis work presents an up-to-date model for the simulation of non-stationary ground motions, including several novelties compared to the original study of Sabetta and Pugliese (Bull Seism Soc Am 86:337–352, 1996). The selection of the input motion in the framework of earthquake engineering has become progressively more important with the growing use of nonlinear dynamic analyses. Regardless of the increasing availability of large strong motion databases, ground motion records are not always available for a given earthquake scenario and site condition, requiring the adoption of simulated time series. Among the different techniques for the generation of ground motion records, we focused on the methods based on stochastic simulations, considering the time- frequency decomposition of the seismic ground motion. We updated the non-stationary stochastic model initially developed in Sabetta and Pugliese (Bull Seism Soc Am 86:337–352, 1996) and later modified by Pousse et al. (Bull Seism Soc Am 96:2103–2117, 2006) and Laurendeau et al. (Nonstationary stochastic simulation of strong ground-motion time histories: application to the Japanese database. 15 WCEE Lisbon, 2012). The model is based on the S-transform that implicitly considers both the amplitude and frequency modulation. The four model parameters required for the simulation are: Arias intensity, significant duration, central frequency, and frequency bandwidth. They were obtained from an empirical ground motion model calibrated using the accelerometric records included in the updated Italian strong-motion database ITACA. The simulated accelerograms show a good match with the ground motion model prediction of several amplitude and frequency measures, such as Arias intensity, peak acceleration, peak velocity, Fourier spectra, and response spectra.


2021 ◽  
Vol 37 (1_suppl) ◽  
pp. 1420-1439
Author(s):  
Albert R Kottke ◽  
Norman A Abrahamson ◽  
David M Boore ◽  
Yousef Bozorgnia ◽  
Christine A Goulet ◽  
...  

Traditional ground-motion models (GMMs) are used to compute pseudo-spectral acceleration (PSA) from future earthquakes and are generally developed by regression of PSA using a physics-based functional form. PSA is a relatively simple metric that correlates well with the response of several engineering systems and is a metric commonly used in engineering evaluations; however, characteristics of the PSA calculation make application of scaling factors dependent on the frequency content of the input motion, complicating the development and adaptability of GMMs. By comparison, Fourier amplitude spectrum (FAS) represents ground-motion amplitudes that are completely independent from the amplitudes at other frequencies, making them an attractive alternative for GMM development. Random vibration theory (RVT) predicts the peak response of motion in the time domain based on the FAS and a duration, and thus can be used to relate FAS to PSA. Using RVT to compute the expected peak response in the time domain for given FAS therefore presents a significant advantage that is gaining traction in the GMM field. This article provides recommended RVT procedures relevant to GMM development, which were developed for the Next Generation Attenuation (NGA)-East project. In addition, an orientation-independent FAS metric—called the effective amplitude spectrum (EAS)—is developed for use in conjunction with RVT to preserve the mean power of the corresponding two horizontal components considered in traditional PSA-based modeling (i.e., RotD50). The EAS uses a standardized smoothing approach to provide a practical representation of the FAS for ground-motion modeling, while minimizing the impact on the four RVT properties ( zeroth moment, [Formula: see text]; bandwidth parameter, [Formula: see text]; frequency of zero crossings, [Formula: see text]; and frequency of extrema, [Formula: see text]). Although the recommendations were originally developed for NGA-East, they and the methodology they are based on can be adapted to become portable to other GMM and engineering problems requiring the computation of PSA from FAS.


1983 ◽  
Vol 73 (6A) ◽  
pp. 1499-1511
Author(s):  
Paul Silver

Abstract A method is proposed for retrieving source-extent parameters from far-field body-wave data. At low frequency, the normalized P- or S-wave displacement amplitude spectrum can be approximated by |Ω^(r^,ω)| = 1 − τ2(r^)ω2/2 where r^ specifies a point on the focal sphere. For planar dislocation sources, τ2(r^) is linearly related to statistical measures of source dimension, source duration, and directivity. τ2(r^) can be measured as the curvature of |Ω^(r^,ω)| at ω = 0 or the variance of the pulse Ω^(r^,t). The quantity ωc=2τ−1(r^) is contrasted with the traditional corner frequency ω0, defined as the frequency at the intersection of the low- and high-frequency trends of |Ω^(r^,ω)|. For dislocation models without directivity, ωc(P) ≧ ωc(S) for any r^. A mean corner frequency defined by averaging τ2(r^) over the focal sphere, ω¯c=2<τ2(r^)>−1/2, satisfies ωc(P) > ωc(S) for any dislocation source. This behavior is not shared by ω0. It is shown that ω0 is most sensitive to critical times in the rupture history of the source, whereas ωc is determined by the basic parameters of source extent. Evidence is presented that ωc is the corner frequency measured on actual seismograms. Thus, the commonly observed corner frequency shift (P-wave corner greater than the S-wave corner), now viewed as a shift in ωc is simply a result of spatial finiteness and is expected to be a property of any dislocation source. As a result, the shift cannot be used as a criterion for rejecting particular dislocation models.


2021 ◽  
Author(s):  
Nesrin Yenihayat ◽  
Eser Çaktı ◽  
Karin Şeşetyan

<p>One of the major earthquakes that resulted in intense damages in Istanbul and its neighborhoods took place on 10 July 1894. The 1894 earthquake resulted in 474 losses of life and 482 injuries. Around 21,000 dwellings were damaged, which is a number that corresponds to 1/7 of the total dwellings of the city at that time. Without any doubt, the exact loss of life was higher. Because of the censorship, the exact loss numbers remained unknown. There is still no consensus about its magnitude, epicentral location, and rupture of length. Even though the hardness of studying with historical records due to their uncertainties and discrepancies, researchers should enlighten the source parameters of the historical earthquakes to minimize the effect of future disasters especially for the cities located close to the most active fault lines as Istanbul. The main target of this study is to enlighten possible source properties of the 1894 earthquake with the help of observed damage distribution and stochastic ground motion simulations. In this paper, stochastic based ground motion scenarios will be performed for the 10 July 1894 Istanbul earthquake, using a finite fault simulation approach with a dynamic corner frequency and the results will be compared with our intensity map obtained from observed damage distributions. To do this, in the first step, obtained damage information from various sources has been presented, evaluated, and interpreted. Secondly, we prepared an intensity map associated with the 1894 earthquake based on macro-seismic information, and damage analysis and classification. For generating ground motions with a stochastic finite fault simulation approach, the EXSIM 2012 software has been used. Using EXSIM, several scenarios are modeled with different source, path, and site parameters. Initial source properties have been obtained from findings of our previous study on the simulation of the 26 September 2019 Silivri (Istanbul) earthquake with Mw 5.8. With the comparison of spatial distributions of the ground motion intensity parameters to the obtained damage and intensity maps, we estimate the optimum location and source parameters of the 1894 Earthquake.</p>


Author(s):  
Chen Ji ◽  
Ralph J. Archuleta

Abstract We investigate the relation between the kinematic double-corner-frequency source spectral model JA19_2S (Ji and Archuleta, 2020) and static fault geometry scaling relations proposed by Leonard (2010). We find that the nonself-similar low-corner-frequency scaling relation of JA19_2S model can be explained using the fault length scaling relation of Leonard’s model combined with an average rupture velocity ∼70% of shear-wave speed for earthquakes 5.3 < M< 6.9. Earthquakes consistent with both models have magnitude-independent average static stress drop and average dynamic stress drop around 3 MPa. Their scaled energy e˜ is not a constant. The decrease of e˜ with magnitude can be fully explained by the magnitude dependence of the fault aspect ratio. The high-frequency source radiation is generally controlled by seismic moment, static stress drop, and dynamic stress drop but is further modulated by the fault aspect ratio and the relative location of the hypocenter. Based on these two models, the commonly quoted average rupture velocity of 70%–80% of shear-wave speed implies predominantly unilateral rupture.


2017 ◽  
Vol 33 (3) ◽  
pp. 875-894 ◽  
Author(s):  
Tadahiro Kishida ◽  
Danilo Di Giacinto ◽  
Giuseppe Iaccarino

Numerous time series for small-to-moderate-magnitude (SMM) earthquakes have been recorded in many regions. A uniformly-processed ground-motion database is essential in the development of regional ground-motion models. An automated processing protocol is useful in developing the database for these earthquakes especially when the number of recordings is substantial. This study compares a manual and an automated ground-motion processing methods using SMM earthquakes. The manual method was developed by the Pacific Earthquake Engineering Research Center to build the database of time series and associated ground-motion parameters. The automated protocol was developed to build a database of pseudo-spectral acceleration for the Kiban-Kyoshin network recordings. Two significant differences were observed when the two methods were applied to identical acceleration time series. First, the two methods differed in the criteria for the acceptance or rejection of the time series in the database. Second, they differed in the high-pass corner frequency used to filter noise from the acceleration time series. The influences of these differences were investigated on ground-motion parameters to elucidate the quality of ground-motion database for SMM earthquakes.


Sign in / Sign up

Export Citation Format

Share Document