Rupture Process and Strong-Motion Generation of the 2014 Iquique, Northern Chile, Earthquake

2016 ◽  
Vol 10 (03) ◽  
pp. 1640008 ◽  
Author(s):  
Wataru Suzuki ◽  
Nelson Pulido ◽  
Shin Aoi

To investigate the rupture process and strong-motion generation of the [Formula: see text] 8.2 Iquique, Northern Chile, earthquake in 2014, we estimated kinematic source models from waveform inversion and back-projection analyses using strong-motion records. A slip model derived from the waveform inversion using the low-frequency (0.02–0.125[Formula: see text]Hz) velocities is characterized by a large slip area localized 50[Formula: see text]km south of the epicenter with a peak slip of 10[Formula: see text]m, and a deeper slip area with a peak slip above 2[Formula: see text]m located below the coast. The main rupture of these areas started 25[Formula: see text]s after the initial break generating two pronounced phases observed in most of the records. The landward slip area ruptured for about 10[Formula: see text]s generating the first impulsive phase, while the offshore largest slip area ruptured for 20[Formula: see text]s creating a longer duration phase observed later. Results from a back-projection analysis based on stacking of envelopes of 5–10-Hz accelerations indicate that the high-frequency radiation propagated down-dip towards the coast, reaching its maximum value from 25[Formula: see text]s to 40[Formula: see text]s, far away from the shallow main slip area obtained from low-frequency waveform inversion. Our results suggest a clear depth dependence of the seismic wave radiation during the Iquique earthquake.

1997 ◽  
Vol 87 (4) ◽  
pp. 904-917 ◽  
Author(s):  
Yasumaro Kakehi ◽  
Kojiro Irikura

Abstract We investigate the process of high-frequency (1 to 10 Hz) radiation on the fault plane of the 1993 Hokkaido-Nansei-Oki, Japan, earthquake (MW = 7.5) from the envelope inversion of strong-motion acceleration seismograms. For the analysis, empirical Green's functions are used because theoretical approach is not available for such high frequencies. The source is modeled with two fault planes with different strike angles. The rupture process of this earthquake is very complex in terms of high-frequency wave generation. The rupture, which started on the northern fault plane, had a delay of about 10 sec or propagated very slowly between the northern and southern fault planes. High-frequency wave radiation is large at the northern and southern edges of the source region. Deceleration of rupture is also observed there. This is interpreted to be associated with stopping of rupture. Another high-frequency wave radiation area is found at the center of the northern fault plane, where discontinuity in the depth distribution of aftershocks suggests an existence of a barrier. The areas of high- and low-frequency wave radiation are not correlated. This is considered to result from the complexity of rupture process. We cannot distinguish between westward and eastward dip of the southern fault plane because of one-sided station distribution.


2019 ◽  
Vol 71 (1) ◽  
Author(s):  
Hisahiko Kubo ◽  
Asako Iwaki ◽  
Wataru Suzuki ◽  
Shin Aoi ◽  
Haruko Sekiguchi

Abstract In this study, we investigate the source rupture process of the 2018 Hokkaido Eastern Iburi earthquake in Japan (MJMA 6.7) and how the ground motion can be reproduced using available source and velocity models. First, we conduct a multiple-time-window kinematic waveform inversion using strong-motion waveforms, which indicates that a large-slip area located at a depth of 25–30 km in the up-dip direction from the hypocenter was caused by a rupture propagating upward 6–12 s after its initiation. Moreover, the high-seismicity area of aftershocks did not overlap with the large-slip area. Subsequently, using the obtained source model and a three-dimensional velocity structure model, we conduct a forward long-period (< 0.5 Hz) ground-motion simulation. The simulation was able to reproduce the overall ground-motion characteristics in the sedimentary layers of the Ishikari Lowland.


Author(s):  
Deyu Yin ◽  
Yun Dong ◽  
Qifang Liu ◽  
Jingke Wu ◽  
Huasheng Sun ◽  
...  

ABSTRACT We estimated the areas exhibiting high-frequency (1∼10  Hz) wave radiation on the fault plane of the 2008 Wenchuan earthquake, by applying envelope inversion to strong-motion acceleration records. The corrected records of two small earthquakes are adopted as the empirical Green’s functions. Considering the change in the rupture pattern of the Wenchuan earthquake from southwest to northeast, the records of small earthquakes dominated by thrust and strike-slip are utilized as the empirical Green’s function for the southwestern and northeastern fault sections, respectively. The results are as follows: (1) According to the high-frequency wave radiation, the rupture process is complex. High-frequency waves radiated strongly in six areas: around the initial rupture point, along the north and south edges of the fault plane, near the area of intersection with the cross-cutting Xiaoyudong fault, south of Nanba, and near the area of Qingchuan. In total, these areas can be divided into three cases. In the first situation, high-frequency waves radiated strongly around the initial rupture area, which may be associated with the initiation of rupture and a high stress drop. The second location is near the periphery of the fault, which is associated with the termination of rupture. The third condition comprises high-frequency waves near the intersection with the cross-cutting Xiaoyudong fault. This area as a geometric barrier, and the surface rupture is observed. (2) The distribution patterns of the high- and low-frequency radiation intensity differ on the fault plane. From the hypocenter to the point of intersection with the Xiaoyudong fault, the high-frequency wave is located around the area with large slip value. In other areas, the distribution of the high- and low-frequency radiation is no obvious relationship. This different characteristic indicates the complexity of the rupture process.


2019 ◽  
Vol 71 (1) ◽  
Author(s):  
Kimiyuki Asano ◽  
Tomotaka Iwata

Abstract The source rupture process of the 2018 Hokkaido Eastern Iburi earthquake (MJMA 6.7) was analyzed by a kinematic waveform inversion method using strong-motion data in 0.04–0.5 Hz. This earthquake occurred close to the Hidaka Collision Zone and the Ishikari depression, where the crustal structure is rather complex. Thus, we used a three-dimensional velocity structure model to compute the theoretical Green’s functions by the finite difference method. A source fault model with strike-angle variation was set based on the spatial distribution of the early aftershocks. The strong-motion stations used for the source inversion were selected based on the result of forward ground motion simulation of a moderate aftershock. The slip in the first 5 s was relatively small, but an area of significant slip with peak slip of 1.7 m was found in the depth range from 22 to 32 km. The rupture propagated upward mainly in the southwest direction. Based on the regional crustal structure and the configuration of the Moho discontinuity, the large-slip area was thought to be located in the lower crust, and its rupture did not reach the upper part of the continental crust. Most of the early aftershocks occurred around the large-slip area. The later aftershocks at the depth shallower than 20 km occurred outside the causative source fault of the mainshock. Three-dimensional ground motion simulation demonstrated that the heterogeneous source process and the three-dimensional basin and crustal velocity structure brought a large velocity pulse to an area to the southwest of the source fault, where the largest PGV was observed during the mainshock. The spatial distribution of the simulated PGV resembled the observed PGV distribution except some sites located inside the Ishikari depression where thick Quaternary soft low-velocity sediments exist at the top of the basin.


Geophysics ◽  
1992 ◽  
Vol 57 (6) ◽  
pp. 854-859 ◽  
Author(s):  
Xiao Ming Tang

A new technique for measuring elastic wave attenuation in the frequency range of 10–150 kHz consists of measuring low‐frequency waveforms using two cylindrical bars of the same material but of different lengths. The attenuation is obtained through two steps. In the first, the waveform measured within the shorter bar is propagated to the length of the longer bar, and the distortion of the waveform due to the dispersion effect of the cylindrical waveguide is compensated. The second step is the inversion for the attenuation or Q of the bar material by minimizing the difference between the waveform propagated from the shorter bar and the waveform measured within the longer bar. The waveform inversion is performed in the time domain, and the waveforms can be appropriately truncated to avoid multiple reflections due to the finite size of the (shorter) sample, allowing attenuation to be measured at long wavelengths or low frequencies. The frequency range in which this technique operates fills the gap between the resonant bar measurement (∼10 kHz) and ultrasonic measurement (∼100–1000 kHz). By using the technique, attenuation values in a PVC (a highly attenuative) material and in Sierra White granite were measured in the frequency range of 40–140 kHz. The obtained attenuation values for the two materials are found to be reliable and consistent.


2019 ◽  
Vol 16 (6) ◽  
pp. 1017-1031 ◽  
Author(s):  
Yong Hu ◽  
Liguo Han ◽  
Rushan Wu ◽  
Yongzhong Xu

Abstract Full Waveform Inversion (FWI) is based on the least squares algorithm to minimize the difference between the synthetic and observed data, which is a promising technique for high-resolution velocity inversion. However, the FWI method is characterized by strong model dependence, because the ultra-low-frequency components in the field seismic data are usually not available. In this work, to reduce the model dependence of the FWI method, we introduce a Weighted Local Correlation-phase based FWI method (WLCFWI), which emphasizes the correlation phase between the synthetic and observed data in the time-frequency domain. The local correlation-phase misfit function combines the advantages of phase and normalized correlation function, and has an enormous potential for reducing the model dependence and improving FWI results. Besides, in the correlation-phase misfit function, the amplitude information is treated as a weighting factor, which emphasizes the phase similarity between synthetic and observed data. Numerical examples and the analysis of the misfit function show that the WLCFWI method has a strong ability to reduce model dependence, even if the seismic data are devoid of low-frequency components and contain strong Gaussian noise.


Sign in / Sign up

Export Citation Format

Share Document