Seismic radiation from dynamic coalescence, and the reconstruction of dynamic source parameters on a planar fault

The dynamic coalescence of two mode II cracks on a planar fault is simulated here using the elastodynamic boundary integral equation method. We focus on the complexity of the resultant slip rate and seismic radiation in the crack coalescence model (CCM) and on the reconstruction of a single crack model (SCM) that can reproduce the CCM waveforms from heterogeneous source parameters rather than coalescence. Simulation results reveal that localized higher slip rates are generated by coalescence as a result of stress interaction between the approaching crack tips. The synthesized seismic radiation exhibits a distinct coalescence phase that has striking similarities to stopping phases in the radiation and propagation properties. The corresponding SCM yields a singular increase in the stress drop distribution, which is accompanied by a sudden decrease in it across the point of coalescence in the CCM. This implies that the generation of high-frequency radiation is more efficient from coalescence than from stopping, although both phenomena exhibit the same strong  ω⁻²  -type displacement spectra.     

Retrieving earthquake signature in grace gravity solutions

The GRACE satellites have been orbiting the Earth since 2002, monitoring the time variable gravity field. Some of the observed fluctuations are due to geodynamic causes, but they are often hidden in the complex signal, composed of hydrology, ocean, atmosphere, and geodynamics, the signal of geodynamic origin being usually the smallest. In addition, dealiasing residuals and noise make the separation of the signal from the different causes more difficult. We proposed a method based on the Empirical Orthogonal Function decomposition to extract the signal of physical origin, under the hypothesis that the physical signal is spatially more consistent than the noise and aliasing incomplete correction. We used synthetic geoid variations associated with earthquakes located at nearly 2000 positions at the Earth surface, based on several examples of large actual subduction events. We show that, with the present day accuracy, we can retrieve the geoid variations associated with more than 98 per cent of the earthquakes of magnitude 9 or above, around 60 per cent for magnitude 8.8, 40 per cent for magnitude 8.6 and 33 per cent for magnitude 8.3. Some events, with the right properties and location, can be detected with magnitude as low as 8. We then applied the method to the GRACE solutions, and retrieved the Hokkaido event (2003) and the Sumatra event (2004), which is in agreement with the retrieval rates mentioned here above.     

2015年4月25日尼泊尔廓尔喀M_S8.1级地震强地面运动

2015年4月25日尼泊尔廓尔喀地区发生MS8.1级地震,本文分析了加德满都强震台站的强地震动记录特征,其水平向地震动表现为明显的脉冲地震动,脉冲周期约为6.0s,反映了近场地震动的方向性效应;处于深厚沉积层上的加德满都谷地对地震动有一定的放大作用,且主震的场地峰值频率向低频段偏移,出现明显的场地非线性反应;采用随机有限断层方法模拟的空间地震动分布与宏观地震烈度图符合较好,可以为分析宏观震害提供参考;通过与尼泊尔设计反应谱对比,得出在短周期与长周期段,加德满都台站的水平向地震动的反应谱远高于设计反应谱,说明即使严格按照尼泊尔建筑抗震规范设计施工的建筑结构也难以抵御此次地震的破坏。     

Multifractal characteristic analysis of near-fault earthquake ground motions

This study aims to reveal the multi-scaling behavior and quantify the irregularity of near-fault earthquake ground motions from a new perspective of multifractal theory. Based on multifractal detrended fluctuation analysis, the multifractal characteristic parameters of acceleration time series for typical near-fault ground motions are calculated, and their correlations with two period parameters (i.e., mean period T_m and characteristic period T_c) and box-counting fractal dimensions are analyzed. Numerical results of strong nonlinear dependence of generalized Hurst exponents h(q) upon the fluctuation orders q indicate that near-fault ground motions present the multifractal properties and long-range correlation obviously. Furthermore, the scaling exponent h(2) of near-fault records has a strong correlation with their periods T_m and T_c, and strongly negative correlation with their box dimension. Moreover, h(2) can be regarded as a measure of frequency content and irregularity degree of strong earthquake ground motions. Finally, it is pointed out that the long-range correlation of small and large fluctuation is the major source of multifractality of near-fault ground motions.     

Nonlinear stochastic seismic analysis of buried pipeline systems

In this paper, a nonlinear stochastic seismic analysis program for buried pipeline systems is developed on the basis of a probability density evolution method (PDEM). A finite element model of buried pipeline systems subjected to seismic wave propagation is established. The pipelines in this model are simulated by 2D beam elements. The soil surrounding the pipelines is simulated by nonlinear distributed springs and linear distributed springs along the axial and horizontal directions, respectively. The joints between the segmented pipes are simulated by nonlinear concentrated springs. Thereafter, by considering the basic random variables of ground motion and soil, the PDEM is employed to capture the stochastic seismic responses of pipeline systems. Meanwhile, a physically based method is employed to simulate the random ground motion field for the area where the pipeline systems are located. Finally, a numerical example is investigated to validate the proposed program.     

Study on elastic response of structures to near-fault ground motions through record decomposition

Accelerograms recorded near active faults have some important characteristics that make them different from those recorded in far-fault regions. High-frequency components in acceleration records and long-period velocity pulses are among notable specifications of such ground motions. In this paper, a moving average filtering with appropriate cut-off frequency has been used to decompose the near-fault ground motions into two components having different frequency contents: first, Pulse-Type Record (PTR) that possesses long-period pulses; second, the relatively high-frequency BackGround Record (BGR), which does not include large velocity pulses. Comparing the results with those extracted through wavelet analysis shows that moving average filter is an appropriate and efficient tool for near-fault records decomposition. The method is applied to decompose a suite of 91 selected near-fault records and the elastic response of structures is examined through their response to the decomposed parts. The results emphasizes that in contrast with ordinary far-fault earthquake records, response spectra of near-fault ground motions typically have two distinct local peaks, which are representatives of the high- and low-frequency components, i.e., BGR and PTR, respectively. Moreover, a threshold period is identified below which the response of structures is dominated by BGR while PTR controls the response of structures with periods longer than this period.     

Approximation approach to the SRM based on root decomposition in the simulation of spatially varying ground motions

The spectral representation method(SRM) is widely used to simulate spatially varying ground motions.This study focuses on the approximation approach to the SRM based on root decomposition,which can improve the efficiency of the simulation.The accuracy of the approximation approach may be affected by three factors: matrix for decomposition,distribution of frequency interpolation nodes and elements for interpolation.The influence of these factors on the accuracy of this approach is examined and the following conclusions are drawn.The SRM based on the root decomposition of the lagged coherency matrix exhibits greater accuracy than the SRM based on the root decomposition of the cross spectral matrix.The equal energy distribution of frequency interpolation nodes proposed in this study is more effective than the counter pith with an equal spacing.Elements for interpolation do not have much of an effect on the accuracy,so interpolation of the elements of the decomposed matrix is recommended because it is less complicated from a computational efficiency perspective.     

Dominant pulse simulation of near fault ground motions

In this study, a new mathematical model is developed composed of two parts, including harmonic and polynomial expressions for simulating the dominant velocity pulse of near fault ground motions. Based on a proposed velocity function, the corresponding expressions for the ground acceleration and displacement time histories are also derived. The proposed model is then fitted using some selected pulse-like near fault ground motions in the Next Generation Attenuation (NGA) project library. The new model is not only simple in form but also simulates the long-period portion of actual velocity near fault records with a high level of precision. It is shown that the proposed model-based elastic response spectra are compatible with the near fault records in the neighborhood of the prevailing frequency of the pulse. The results indicate that the proposed model adequately simulates the components of the time histories. Finally, the energy of the proposed pulse was compared with the energy of the actual record to confirm the compatibility.