Deterministic 3D Ground-Motion Simulations (0–5 Hz) and Surface Topography Effects of the 30 October 2016 Mw 6.5 Norcia, Italy, Earthquake

Author(s):  
Arben Pitarka ◽  
Aybige Akinci ◽  
Pasquale De Gori ◽  
Mauro Buttinelli

ABSTRACT The Mw 6.5 Norcia, Italy, earthquake occurred on 30 October 2016 and caused extensive damage to buildings in the epicentral area. The earthquake was recorded by a network of strong-motion stations, including 14 stations located within a 5 km distance from the two causative faults. We used a numerical approach for generating seismic waves from two hybrid deterministic and stochastic kinematic fault rupture models propagating through a 3D Earth model derived from seismic tomography and local geology. The broadband simulations were performed in the 0–5 Hz frequency range using a physics-based deterministic approach modeling the earthquake rupture and elastic wave propagation. We used SW4, a finite-difference code that uses a conforming curvilinear mesh, designed to model surface topography with high numerical accuracy. The simulations reproduce the amplitude and duration of observed near-fault ground motions. Our results also suggest that due to the local fault-slip pattern and upward rupture directivity, the spatial pattern of the horizontal near-fault ground motion generated during the earthquake was complex and characterized by several local minima and maxima. Some of these local ground-motion maxima in the near-fault region were not observed because of the sparse station coverage. The simulated peak ground velocity (PGV) is higher than both the recorded PGV and predicted PGV based on empirical models for several areas located above the fault planes. Ground motions calculated with and without surface topography indicate that, on average, the local topography amplifies the ground-motion velocity by 30%. There is correlation between the PGV and local topography, with the PGV being higher at hilltops. In contrast, spatial variations of simulated PGA do not correlate with the surface topography. Simulated ground motions are important for seismic hazard and engineering assessments for areas that lack seismic station coverage and historical recordings from large damaging earthquakes.

2021 ◽  
Vol 111 (5) ◽  
pp. 2441-2462 ◽  
Author(s):  
Rosemary Fayjaloun ◽  
Mayssa Dabaghi ◽  
Cecile Cornou ◽  
Mathieu Causse ◽  
Yang Lu ◽  
...  

ABSTRACT Lebanon is a densely populated country crossed by major faults. Historical seismicity shows the potential of earthquakes with magnitudes >7, but large earthquakes have never been instrumentally recorded in Lebanon. Here, we propose a method to simulate near-fault broadband ground motions for a potential Mw 7 earthquake on the Yammouneh fault (YF)—the largest branch of the Dead Sea Transform fault that bisects Lebanon from north to south. First, we performed the first 3D tomography study of Lebanon using ambient noise correlation, which showed that Lebanon could be approximated by a 1D velocity structure for low-frequency (LF) ground-motion simulation purposes. Second, we generated suites of kinematic rupture models on the YF, accounting for heterogeneity of the rupture process, and uncertainty of the rupture velocity and hypocenter location. The radiated seismic energy was next propagated in the inferred 1D velocity model to obtain suites of LF ground motions (<1 Hz) at four hypothetical near-fault seismic stations. These LF simulations included the main features of near-fault ground motions, such as the impulsive character of ground velocity due to the rupture directivity or fling-step effects (so-called pulse-like ground motions). Third, to obtain broadband ground motions (up to 10 Hz), we proposed a hybrid technique that combined the simulated LF ground motions with high-frequency (HF) stochastic simulations, which were empirically calibrated using a worldwide database of near-fault recordings. Contrary to other hybrid approaches, in which the LF and HF motions are generally computed independently, the characteristics of stochastic HF ground motions were conditioned on those of LF ground motions (namely on the characteristics of the velocity pulse, if it existed, or on the absence of a pulse). The simulated peak ground accelerations were in agreement with the ones reported in the Next Generation Attenuation-West2 (NGA-West2) database for similar magnitude and distances and with three NGA-West2 ground-motion prediction equations.


2019 ◽  
Vol 109 (6) ◽  
pp. 2384-2398 ◽  
Author(s):  
Xiaoyu Chen ◽  
Dongsheng Wang ◽  
Rui Zhang

Abstract Large‐amplitude and long‐period pulses are observed in velocity time histories of near‐fault ground‐motion records. The pulses in these records have significant damage effect on flexible structures due to their long‐period property; therefore, more attention should be paid to the frequency components in the ground motion. Based on the identification of frequency components in the original record, a new method based on the Hilbert–Huang transform (HHT) is proposed here. A ground‐motion record can be decomposed into several intrinsic mode functions (IMFs) that carry different frequency components by the HHT without contamination from any a prior function. Only two fixed parameters, the peak ground velocity (PGV)/peak ground acceleration (PGA) ratio and the energy change of every IMF, are used to classify pulse‐like ground‐motion records. The inherent pulses of these records can also be extracted, based on the selection of IMFs for which PGV/PGA ratios are larger than 0.12 and energy changes that are greater than 0.1. For multipulse cases, all the pulses can be captured after extracting once, and the time course of inherent pulses can also be obtained. Then, pulse periods are calculated based on the solutions of instantaneous frequency of the peak for the extracted pulses. All the periods obtained using the HHT method can be verified by the results obtained from Baker’s wavelet method. The 24 controversial records that are discussed in previous studies are examined here as well. The HHT method is a complete procedure that includes the classification of pulse‐like ground motions, the extraction of velocity pulses, and the solution of pulse periods. It works well for multipulse records, especially because it can provide the exact timing of all the inherent pulses.


2021 ◽  
Author(s):  
Joanna Holmgren ◽  
Maximilian Werner

<p>The United Downs enhanced geothermal system in Cornwall, UK, has induced several microseismic events since flow testing began in August 2020, targeting a granitic intrusion at 5 km depth. As of January 2021, two events exceeding local magnitudes (M<sub>L</sub>) 1.5 have occurred, highlighting the associated seismic hazard and providing initial data for a preliminary assessment of the region’s ground motion response. However, with only one national seismic station publicly available within 90 km of the site, public data are scarce. In an effort to involve the surrounding communities in the geothermal project, United Downs provided Raspberry Shake 1D or 4D (one vertical geophone, with 4D containing an additional three accelerometers) seismographs to nearby schools, increasing the number of publicly available seismic stations to ten within 15 km of the site. In this study, we assess the ground motions recorded by the Raspberry Shake stations and evaluate their utility for probing ground motions models (GMMs) and the effects of the local geology.</p><p>171 earthquakes between M<sub>L</sub> -1.3 to 1.7 originating at United Downs have been recorded to date, with 37 events above M<sub>L</sub> 0.0. Unfortunately, the accelerometer components of the Raspberry Shake instruments contained too high background noise levels to be useable, leaving only the vertical geophone component to be analysed for each of the instruments. We find that while the peak ground velocity (PGV) values are in line with those predicted from the Douglas et al. (2013) geothermal GMM, the area experiences higher peak ground acceleration (PGA) than expected. We also find that the observed PGVs and PGAs match the region’s geological features, consisting of a combination of igneous intrusions and sedimentary sandstones and mudstones. For sparse national seismic networks, Raspberry Shake stations can provide a quick initial evaluation of seismic events and their ground motions before industry releases private data for more detailed analyses.</p>


2021 ◽  
Author(s):  
Aybige Akinci ◽  
Daniele Cheloni ◽  
AHMET ANIL DINDAR

Abstract On 30 October 2020 a MW 7.0 earthquake occurred in the eastern Aegean Sea, between the Greek island of Samos and Turkey’s Aegean coast, causing considerable seismic damage and deaths, especially in the Turkish city of Izmir, approximately 70 km from the epicenter. In this study, we provide a detailed description of the Samos earthquake, starting from the fault rupture to the ground motion characteristics. We first use Interferometric Synthetic Aperture Radar (InSAR) and Global Positioning System (GPS) data to constrain the source mechanisms. Then, we utilize this information to analyze the ground motion characteristics of the mainshock in terms of peak ground acceleration (PGA), peak ground velocity (PGV), and spectral pseudo-accelerations. Modelling of geodetic data shows that the Samos earthquake ruptured a NNE-dipping normal fault located offshore north of Samos, with up to 2.5-3 m of slip and an estimated geodetic moment of 3.3 ⨯ 1019 Nm (MW 7.0). Although low PGA were induced by the earthquake, the ground shaking was strongly amplified in Izmir throughout the alluvial sediments. Structural damage observed in Izmir reveals the potential of seismic risk due to the local site effects. To better understand the earthquake characteristics, we generated and compared stochastic strong ground motions with the observed ground motion parameters as well as the ground motion prediction equations (GMPEs), exploring also the efficacy of the region-specific parameters which may be improved to better predict the expected ground shaking from future large earthquakes in the region.


An effective earthquake (Mw 7.9) struck Alaska on 3 November, 2002. This earthquake ruptured 340 km along Susitna Glacier, Denali and Totschunda faults in central Alaska. The peak ground acceleration (PGA) was recorded about 0.32 g at station PS10, which was located 3 km from the fault rupture. The PGA would have recorded a high value, if more instruments had been installed in the region. A numerical study has been conducted to find out the possible ground motion record that could occur at maximum horizontal slip during the Denali earthquake. The current study overcomes the limitation of number of elements to model the Denali fault. These numerical results are compared with observed ground motions. It is observed that the ground motions obtained through numerical analysis are in good agreement with observed ground motions. From numerical results, it is observed that the possible expected PGA is 0.62 g at maximum horizontal slip of Denali fault.


Solid Earth ◽  
2019 ◽  
Vol 10 (3) ◽  
pp. 931-949 ◽  
Author(s):  
Peter Klin ◽  
Giovanna Laurenzano ◽  
Maria Adelaide Romano ◽  
Enrico Priolo ◽  
Luca Martelli

Abstract. During the 2012 seismic sequence of the Emilia region (northern Italy), the earthquake ground motion in the epicentral area featured longer duration and higher velocity than those estimated by empirical-based prediction equations typically adopted in Italy. In order to explain these anomalies, we (1) build up a structural and geophysical 3-D digital model of the crustal sector involved in the sequence, (2) reproduce the earthquake ground motion at some seismological stations through physics-based numerical simulations and (3) compare the observed recordings with the simulated ones. In this way, we investigate how the earthquake ground motion in the epicentral area is influenced by local stratigraphy and geological structure buried under the Po Plain alluvium. Our study area covers approximately 5000 km2 and extends from the right Po River bank to the Northern Apennine morphological margin in the N–S direction, and between the two chief towns of Reggio Emilia and Ferrara in the W–E direction, involving a crustal volume of 20 km thickness. We set up the 3-D model by using already-published geological and geophysical data, with details corresponding to a map at scale of 1:250 000. The model depicts the stratigraphic and tectonic relationships of the main geological formations, the known faults and the spatial pattern of the seismic properties. Being a digital vector structure, the 3-D model can be easily modified or refined locally for future improvements or applications. We exploit high-performance computing to perform numerical simulations of the seismic wave propagation in the frequency range up to 2 Hz. In order to get rid of the finite source effects and validate the model response, we choose to reproduce the ground motion related to two moderate-size aftershocks of the 2012 Emilia sequence that were recorded by a large number of stations. The obtained solutions compare very well to the recordings available at about 30 stations in terms of peak ground velocity and signal duration. Snapshots of the simulated wavefield allow us to attribute the exceptional length of the observed ground motion to surface wave overtones that are excited in the alluvial basin by the buried ridge of the Mirandola anticline. Physics-based simulations using realistic 3-D geomodels show eventually to be effective for assessing the local seismic response and the seismic hazard in geologically complex areas.


2019 ◽  
Vol 2019 ◽  
pp. 1-14 ◽  
Author(s):  
Luqi Xie ◽  
Jing Wu ◽  
Qing Huang ◽  
Chao Tong

The analysis of the ductility and cumulative plastic deformation (CPD) demand of a high-performance buckling-restrained brace (HPBRB) under a strong earthquake and its aftershocks is conducted in this paper. A combination of three continuous excitations with the same ground motion is used to simulate the affection of a strong earthquake and its aftershocks. A six-story HPBRB frame (HPBRBF) is taken as an example to conduct the incremental dynamic analysis (IDA). The seismic responses of the HPBRBF under one, two, and three constant continuous ground motions are compared. The IDA result indicates that the ductility and CPD demand of the BRBs under the three constant continuous ground motions are significantly larger than that excited by only one. Probabilistic seismic demand analysis (PSDA) is performed using seven near-fault ground motions and seven far-fault ground motions to consider the indeterminacy of ground motion. The probabilistic seismic demand curves (PSDCs) for the ductility and CPD demand for the HPBRB under the strong earthquake and its aftershocks are obtained in combining the probabilistic seismic hazard analysis. The results indicate that the AISC threshold value of the CPD with 200 is excessively low for a HPBRBF which suffers the continuous strong aftershocks with near-fault excitations, and a stricter threshold value should be suggested to ensure the ductility and plastic deformation capacity demand of the HPBRB.


2019 ◽  
Vol 109 (4) ◽  
pp. 1265-1281 ◽  
Author(s):  
Arthur J. Rodgers ◽  
Arben Pitarka ◽  
David B. McCallen

Abstract We investigated the effects of fault geometry and assumed minimum shear wavespeed (VSmin) on 3D ground-motion simulations (0–2.5 Hz) in general, using a moment magnitude (Mw) 6.5 earthquake on the Hayward fault (HF). Simulations of large earthquakes on the northeast-dipping HF using the U.S. Geological Survey (USGS) 3D seismic model have shown intensity asymmetry with stronger shaking for the Great Valley Sequence east of the HF (hanging wall) relative to the Franciscan Complex to the west (footwall). We performed simulations with three fault geometries in both plane-layered (1D) and 3D models. Results show that the nonvertical fault geometries result in larger motions on the hanging wall relative to the vertical fault for the same Earth model with up to 50% amplifications in single-component peak ground velocity (PGV) within 10 km of the rupture. Near-fault motions on the footwall are reduced for the nonvertical faults, but less than they are increased on the hanging wall. Simulations assuming VSmin values of 500 and 250  m/s reveal that PGVs are on average 25% higher west of the HF when using the lower VSmin, with some locations amplified by a factor of 3. Increasing frequency content from 2.5 to 5 Hz increases PGV values. Spectral ratios of these two VSmin cases show average amplifications of 2–4 (0.5–1.5 Hz) for the lower VSmin west of the fault. Large differences (up to 2×) in PGV across the HF from previous studies persist even for the case with a vertical fault or VSmin of 250  m/s. We conclude that assuming a VSmin of 500  m/s underestimates intensities west of the HF for frequencies above 0.5 Hz, and that low upper crustal (depth <10  km) shear wavespeeds defined in the 3D model contribute most to higher intensities east of the HF.


2020 ◽  
Vol 110 (6) ◽  
pp. 2828-2842
Author(s):  
Esra Zengin ◽  
Norman Abrahamson

ABSTRACT The velocity pulse in near-fault ground motions has been used as a key characteristic of damaging ground motions. Characterization of the velocity pulse involves three parameters: presence of the pulse, period of the pulse, and amplitude of the pulse. The basic concept behind the velocity pulse is that a large amount of seismic energy is packed into a short time, leading to larger demands on the structure. An intensity measure for near-fault ground motions, which is a direct measure of the amount of energy arriving in short time, called instantaneous power (IP (T1)), is defined as the maximum power of the bandpass-filtered velocity time series measured over a time interval of 0.5T1, in which T1 is the fundamental period of the structure. The records are bandpass filtered in the period band (0.2T1−3T1) to remove the frequencies that are not expected to excite the structure. Zengin and Abrahamson (2020) showed that the drift is better correlated with the IP (T1) than with the velocity pulse parameters for records scaled to the same spectral acceleration at T1. A conditional ground-motion model (GMM) for the IP is developed based on the 5%-damped spectral acceleration at T1, the earthquake magnitude, and the rupture distance. This conditional GMM can be used for record selection for near-fault ground motions that captures the key features of velocity pulses and can lead to a better representation of the median and variability of the maximum interstory drift. The conditional GMM can also be used in a vector hazard analysis for spectral acceleration (T1) and IP (T1) that can be used for more accurate estimation of drift hazard and seismic risk.


Sign in / Sign up

Export Citation Format

Share Document