Earthquake dynamics on circular faults: a review 1970–2015

We review the development of earthquake dynamics taken from the point of view of the origin of seismic radiation instead of the detailed study of rupture propagation on complex surfaces. Many features of seismic radiation can be explained by simple models that serve as elementary canonical problems. Some of these properties are very well known like the fact that at low frequencies, the seismic spectrum is proportional to the seismic moment. At high frequencies, on the other hand, radiation is generated by the motion of the rupture front, in particular stopping phases and geometrical obstacles (barriers). A rupture front moving at constant speed does not produce far-field radiation. For many practical applications, for determining source size and stress drop, for example, it is not necessary to determine geometrical details of the rupture. For such cases, a simple circular crack model is quite sufficient. An improvement on this method is to do dynamic inversion on simple, elliptical-shaped sources and letting the rupture start arbitrarily from a point on the fault. This problem can be solved nowadays with finite differences and a variety of inversion techniques.     

Variability of Antarctic Intermediate Water south of Australia and its relationship with frontal waves

A streamfunction EOF method is applied to a time series of hydrographic sections in the Southern Ocean south of Australia to study water mass variations. Results show that there are large thermohaline variations north of the Subantarctic Front (SAF) at 300–1500 dbar level, indicating upwelling and downwelling of the Antarctic Intermediate Water (AAIW) along isopycnal surfaces. Based on the latest altimeter product, Absolute Dynamic Topography, a mechanism due to frontal wave propagation is proposed to explain this phenomenon, and an index for frontal waves is defined. When the frontal wave is in positive (negative) phase, the SAF flows northeastward (southeastward) with the fresh AAIW downwelling (upwelling). Such mesoscale processes greatly enhance cross-frontal exchanges of water masses. Spectral analysis shows that frontal waves in the Southern Ocean south of Australia are dominated by a period of about 130 days with a phase speed of 4 cm/s and a wavelength of 450 km.     

Three-dimensional admittance analysis of lithospheric elastic thickness over the Louisville Ridge

Using bathymetry and altimetric gravity anomalies, a 1° 9 1° lithospheric effective elastic thickness(Te) model over the Louisville Ridge and its adjacent regions is calculated using the moving window admittance technique. For comparison, three bathymetry models are used: general bathymetric charts of the oceans, SIO V15.1,and BAT_VGG. The results show that BAT_VGG is more suitable for calculating T e than the other two models. T e along the Louisville Ridge was re-evaluated. The southeast of the ridge has a medium Te of 10–20 km, while Te increases dramatically seaward of the Tonga-Kermadec trench as a result of the collision of the Pacific and IndoAustralian plates.     

High-resolution crustal structure of the Yinchuan basin revealed by deep seismic reflection profiling: implications for deep processes of basin

The Yinchuan basin, located on the western margin of the Ordos block, has the characteristics of an active continental rift. A NW-striking deep seismic reflection profile across the center of Yinchuan basin precisely revealed the fine structure of the crust. The images showed that the crust in the Yinchuan basin was characterized by vertical stratifications along a detachment located at a two-way travel time(TWT) of 8.0 s.The most outstanding feature of this seismic profile was the almost flat Mohorovicˇic′ discontinuity(Moho) and a high-reflection zone in the lower crust. This sub-horizontal Moho conflicts with the general assumption of an uplifted Moho under sedimentary basins and continental rifts, and may indicate the action of different processes at depth during the evolution of sedimentary basins or rifts.We present a possible interpretation of these deep processes and the sub-horizontal Moho. The high-reflection zone, which consists of sheets of high-density, mantlederived materials, may have compensated for crustal thinning in the Yinchuan basin, leading to the formation of a sub-horizontal Moho. These high-density materials may have been emplaced by underplating with mantlesourced magma.     

Multiscale anisotropic diffusion for ringing artifact suppression in geophysical deconvolution data

Ringing artifact degradations always appear in the deconvolution of geophysical data. To address this problem, we propose a postprocessing approach to suppress ringing artifacts that uses a novel anisotropic diffusion based on a stationary wavelet transform (SWT) algorithm. In this paper, we discuss the ringing artifact suppression problem and analyze the characteristics of the deconvolution ringing artifact. The deconvolution data containing ringing artifacts are decomposed into different SWT subbands for analysis, and a new multiscale adaptive anisotropic filter is developed to suppress these degradations. Finally, we demonstrate the performance of the proposed method and describe the experiments in detail.     

Dispersion function of Rayleigh waves in porous layered half-space system

Rayleigh wave exploration is based on an elastic layered half-space model. If practical formations contain porous layers, these layers need to be simplified as an elastic medium. We studied the effects of this simplification on the results of Rayleigh wave exploration. Using a half-space model with coexisting porous and elastic layers, we derived the dispersion functions of Rayleigh waves in a porous layered half-space system with porous layers at different depths, and the problem of transferring variables to matrices of different orders is solved. To solve the significant digit overflow in the multiplication of transfer matrices, we propose a simple, effective method. Results suggest that dispersion curves differ in a lowfrequency region when a porous layer is at the surface; otherwise, the difference is small.     

Nonlinear inversion of electrical resistivity imaging using pruning Bayesian neural networks

Conventional artificial neural networks used to solve electrical resistivity imaging (ERI) inversion problem suffer from overfitting and local minima. To solve these problems, we propose to use a pruning Bayesian neural network (PBNN) nonlinear inversion method and a sample design method based on the K-medoids clustering algorithm. In the sample design method, the training samples of the neural network are designed according to the prior information provided by the K-medoids clustering results; thus, the training process of the neural network is well guided. The proposed PBNN, based on Bayesian regularization, is used to select the hidden layer structure by assessing the effect of each hidden neuron to the inversion results. Then, the hyperparameter α _k , which is based on the generalized mean, is chosen to guide the pruning process according to the prior distribution of the training samples under the small-sample condition. The proposed algorithm is more efficient than other common adaptive regularization methods in geophysics. The inversion of synthetic data and field data suggests that the proposed method suppresses the noise in the neural network training stage and enhances the generalization. The inversion results with the proposed method are better than those of the BPNN, RBFNN, and RRBFNN inversion methods as well as the conventional least squares inversion.     

Seismic evaluation of rocking structures through performance assessment and fragility analysis

Numerical studies have been conducted for low- and medium-rise rocking structures to investigate their efficiency as earthquake-resisting systems in comparison with conventional structures. Several non-linear time-history analyses have been performed to evaluate seismic performance of selected cases at desired ground shaking levels, based on key parameters such as total and flexural story drifts and residual deformations. The Far-field record set is selected as input ground motions and median peak values of key parameters are taken as best estimates of system response. In addition, in order to evaluate the probability of exceeding relevant damage states, analytical fragility curves have been developed based on the results of the incremental dynamic analysis procedure. Small exceedance probabilities and acceptable margins against collapse, together with minor associated damages in main structural members, can be considered as superior seismic performance for medium-rise rocking systems. Low-rise rocking systems could provide significant performance improvement over their conventional counterparts notwithstanding certain weaknesses in their seismic response.     

Newmark design spectra considering earthquake magnitudes and site categories

Newmark design spectra have been implemented in many building codes, especially in building codes for critical structures. Previous studies show that Newmark design spectra exhibit lower amplitudes at high frequencies and larger amplitudes at low frequencies in comparison with spectra developed by statistical methods. To resolve this problem, this study considers three suites of ground motions recorded at three types of sites. Using these ground motions, influences of the shear-wave velocity, earthquake magnitudes, source-to-site distances on the ratios of ground motion parameters are studied, and spectrum amplification factors are statistically calculated. Spectral bounds for combinations of three site categories and two cases of earthquake magnitudes are estimated. Site design spectrum coefficients for the three site categories considering earthquake magnitudes are established. The problems of Newmark design spectra could be resolved by using the site design spectrum coefficients to modify the spectral values of Newmark design spectra in the acceleration sensitive, velocity sensitive, and displacement sensitive regions.     

Water quality assessment of the Huaihe River segment of Bengbu (China) using multivariate statistical techniques

In the study, multivariate statistical methods including principal component analysis (PCA)/factor analysis (FA) and cluster analysis (CA) were applied to analyze surface water quality data sets obtained from the Huaihe River segment of Bengbu (HRSB) and generated during 2 years (2011–2012) monitoring of 19 parameters at 7 sampling sites. The results of PCA for 7 sampling sites revealed that the first four components of PCA showed 94.89% of the total variance in the data sets of HRSB. The Principal components (Factors) obtained from FA indicated that the parameters for water quality variations were mainly related to heavy metals (Pb, Mn, Zn and Fe) and organic related parameters (COD, PI and DO). The results revealed that the major causes of water quality deterioration were related to inflow of industrial, domestic and agricultural effluents into the Huaihe River. Three significant sampling locations—(sites 2, 3 and 4), (sites 1 and 5) and (sites 6 and 7)—were detected on the basis of similarity of their water quality. Thus, these methods were believed to be valuable to help water resources managers understand complex nature of water quality issues and determine the priorities to improve water quality.