A Far-Field Ground-Motion Model for the North Australian Craton from Plate-Margin Earthquakes

ABSTRACT The Australian territory is just over 400 km from an active convergent plate margin with the collision of the Sunda–Banda Arc with the Precambrian and Palaeozoic Australian continental crust. Seismic energy from earthquakes in the northern Australian plate-margin region are channeled efficiently through the low-attenuation North Australian craton (NAC), with moderate-sized (Mw≥5.0) earthquakes in the Banda Sea commonly felt in northern Australia. A far-field ground-motion model (GMM) has been developed for use in seismic hazard studies for sites located within the NAC. The model is applicable for hypocentral distances of approximately 500–1500 km and magnitudes up to Mw 8.0. The GMM provides coefficients for peak ground acceleration, peak ground velocity, and 5%-damped pseudospectral acceleration at 20 oscillator periods from 0.1 to 10 s. A strong hypocentral depth dependence is observed in empirical data, with earthquakes occurring at depths of 100–200 km demonstrating larger amplitudes for short-period ground motions than events with shallower hypocenters. The depth dependence of ground motion diminishes with longer spectral periods, suggesting that the relatively larger ground motions for deeper earthquake hypocenters may be due to more compact ruptures producing higher stress drops at depth. Compared with the mean Next Generation Attenuation-East GMM developed for the central and eastern United States (which is applicable for a similar distance range), the NAC GMM demonstrates significantly higher short-period ground motion for Banda Sea events, transitioning to lower relative accelerations for longer period ground motions.

Relational database used for ground-motion model development in the NGA-sub project

The Next-Generation Attenuation for subduction zone regions project (NGA-Sub) has developed data resources and ground motion models for global subduction zone regions. Here we describe the NGA-Sub database. To optimize the efficiency of data storage, access, and updating, data resources for the NGA-Sub project are organized into a relational database consisting of 20 tables containing data, metadata, and computed quantities (e.g. intensity measures, distances). A database schema relates fields in tables to each other through a series of primary and foreign keys. Model developers and other users mostly interact with the data through a flatfile generated as a time-stamped output of the database. We describe the structure of the relational database, the ground motions compiled for the project, and the means by which the data can be accessed. The database contains 71,340 three-component records from 1880 earthquakes from seven global subduction zone regions: Alaska, Central America and Mexico, Cascadia, Japan, New Zealand, South America, and Taiwan. These data were processed on a component-specific basis to minimize noise effects in the data and remove baseline drifts. Provided ground motion intensity measures include peak acceleration, peak velocity, and 5%-damped pseudo-spectral accelerations for a range of oscillator periods.

Probabilistic Seismic Hazard Analysis Based on Arias Intensity in the North–South Seismic Belt of China

ABSTRACT Arias intensity (IA), as an important seismic parameter, which contains the information of amplitude, frequencies, and duration of ground motion, plays a crucial role in characterizing seismic hazard such as earthquake-induced landslides. In this article, we conducted probabilistic seismic hazard analysis (PSHA) based on IA in China’s north–south seismic belt. We adopted the seismic sources and seismicity parameters used in the fifth generation of the Seismic Ground Motion Parameter Zoning Map of China, and two ground-motion model of IA. The results show that the values of IA are greater than 0.11 m/s in most regions of the north–south seismic belt. The provincial capital cities and most prefecture-level cities in the seismic zone are located in the region with IA-values greater than 0.32 m/s. The values of IA are above 0.54 m/s in the region around the main fault zone. This means that the north–south seismic belt is prone to extremely high-seismic hazard, particularly earthquake-induced landslides. Therefore, it is important to strengthen the evaluation and prevention of earthquake-induced landslides in this area. As we have found significant differences in the values of IA calculated from different ground-motion model, it is necessary to study the ground-motion model of IA for the western geological environment of China. In addition, the PSHA based on IA gives more consideration to the influence of large earthquakes than that based on peak ground acceleration. Therefore, IA plays an important role in seismic design of major engineering projects. The results of this article are of great scientific significance for understanding the seismic hazard of the north–south seismic belt.

Ground motions variability in Israel from 3-D simulations of M 6 and M 7 earthquakes

In Israel, due to low seismicity rates and sparse seismic network, the temporal and spatial coverage of ground motion data is insufficient to estimate the variability of moderate-strong (M > 6) ground motions required to construct a local ground motion model (GMM). To fill this data gap and to study the ground motions variability of M > 6 events, we performed a series of 3-D numerical simulations of M 6 and M 7 earthquakes. Based on the results of the simulations, we developed a statistical attenuation model (AM) and studied the residuals between simulated and AM PGVs and the single station variability. We also compared the simulated ground motions with a global GMM in terms of peak ground velocity (PGV) and significant duration (Ds 595). Our results suggest that the AM was unable to fully capture the simulated ground motions variability, mainly due to the incorporation of super-shear rupture and effects of local sedimentary structures. We also show that an imported GMM considerably deviates from simulated ground motions. This work sets the basis for future development of a comprehensive GMM for Israel, accounting for local sources, path, and site effects.

Estimating Ground-Motion Models via Bayesian Inference using Stan

We provide a simple introduction to the estimation of ground-motion models via Bayesian inference and the probabilistic programming language Stan.We show one ca implement a simple ground-motion model in Stan, and how can run the program from the computer environment R.We also show how one can access the results, and plot summaries of estimated parameters.A large number of different Stan models for the development https://github.com/pstafford/StanGMMTutorial.

A Ground-Motion Model for GNSS Peak Ground Displacement

ABSTRACT We present an updated ground-motion model (GMM) for Mw 6–9 earthquakes using Global Navigation Satellite Systems (GNSS) observations of the peak ground displacement (PGD). Earthquake GMMs inform a range of Earth science and engineering applications, including source characterization, seismic hazard evaluations, loss estimates, and seismic design standards. A typical GMM is characterized by simplified metrics describing the earthquake source (magnitude), observation distance, and site terms. Most often, GMMs are derived from broadband seismometer and accelerometer observations, yet during strong shaking, these traditional seismic instruments are affected by baseline offsets, leading to inaccurate recordings of low-frequency ground motions such as displacement. The incorporation of geodetic data sources, particularly for characterizing the unsaturated ground displacement of large-magnitude events, has proven valuable as a complement to traditional seismic approaches and led to the development of an initial point-source GMM based on PGD estimated from high-rate GNSS data. Here, we improve the existing GMM to more effectively account for fault finiteness, slip heterogeneity, and observation distance. We evaluate the limitations of the currently available GNSS earthquake data set to calibrate the GMM. In particular, the observed earthquake data set is lacking in observations within 100 km of large-magnitude events (Mw>8), inhibiting evaluation of fault dimensions for earthquakes too large to be represented as point sources in the near field. To that end, we separately consider previously validated synthetic GNSS waveforms within 10–1000 km of Mw 7.8–9.3 Cascadia subduction zone scenario ruptures. The synthetic data highlight the importance of fault distance rather than point-source metrics and improve our preparedness for large-magnitude earthquakes with spatiotemporal qualities unlike those in our existing data set.

A Non-Ergodic Ground-Motion Model of Fourier Amplitude Spectra for France

We used an ergodic ground-motion model (GMM) of California of Bayless and Abrahamson (Bull Seismol Soc Am 109(5):2088–2105, 2019) as a backbone model and incorporated the varying-coefficient model (VCM), with a modification for anisotropic path effects, to develop a new non-ergodic GMM for France based on the French RESIF data set (1996-2016). Most of the earthquakes in this database have small-to-moderate magnitudes (M2.0 – M5.2). We developed the GMM for the smoothed effective amplitude spectrum (EAS) rather than for elastic spectral acceleration because it allows the use of small magnitude data to constrain linear effects of the path and site without the complication of the scaling being affected by differences in the response spectral shape. For the VCM, the coefficients of GMM can vary by geographical location and they are estimated using Gaussian-process regression. There is a separate set of coefficients for each source and site coordinate, including both the mean coefficients and the epistemic uncertainty in the coefficients. We further modify the anelastic attenuation term of a GMM by the cell-specific approach of Kuehn et al. (Bull Seismol Soc Am 109 (2): 575–585, 2019) to allow for azimuth-dependent attenuation for each source which reduces the standard deviation of the residuals at long distances. As an example, we compute the 5Hz seismic hazard for two sites using the non-ergodic EAS GMM. At the 1 10-4 annual frequency of exceedance hazard level, there can be a large difference between the ergodic hazard and the non-ergodic hazard if the site is close to the available data. The combination of the non-ergodic median ground motion and the reduced aleatory variability can have large implications for seismic-hazard estimation for long return periods. For some sites, the estimated hazard will increase and for other sites the estimated hazard will decrease compared to the traditional ergodic GMM approach. Due to the skewed distribution of the epistemic uncertainty of the hazard, more of the sites will see a decrease in the mean hazard mean hazard at the 1 10-4 hazard level than will see an increase as a result of using the non-ergodic GMM.