Deep structure of the Nojima Fault, southwest Japan, estimated from borehole observations of fault-zone trapped waves

To estimate the deep structure of the southern part of the Nojima Fault, southwest Japan without the influence of near-surface structures, we analyzed the Love-wave-type fault-zone trapped waves (LTWs) recorded by a borehole seismometer at 1800 m depth. We examined the polarization, dispersion, and dominant frequency of the wavetrain following the direct S-wave in each seismogram to identify the LTW. We selected eight candidates for typical LTWs from 462 records. Because the duration of the LTW increases with hypocentral distance, we infer that the low velocity fault-zone of the Nojima Fault continues towards the seismogenic depth. In addition, since the duration of the LTW increases nonlinearly with hypocentral distance, we infer that the S-wave velocity of the fault-zone increases with depth. The location of events showing the LTW indicates that the fault-zone dips to the southeast at 75° and continues to a depth of approximately 10 km. We assumed a uniform low-velocity waveguide to estimate the average structure of the fault-zone. We estimated the average width, S-wave velocity, and Q_s of the fault-zone by comparing an analytical solution of the LTW with measured data. The average width, S-wave velocity, and Q_s of the fault-zone are 150 to 290 m, 2.5 to 3.2 km/s, and 40 to 90, respectively. Hence the fault-zone structure with a larger width and smaller velocity reduction than the fault-zone model estimated by previous surface observation is more suitable to represent the average fault-zone structure of the Nojima fault. The present study also indicated that the shallow layers and/or a shallow fault-zone structure drastically changes the characteristics of the LTW recorded at the surface, and therefore cause a discrepancy in the fault-zone model between the borehole observation and surface observation.     

冲击凿岩钎杆的瞬态动力学分析

以ANSYS为平台，按照真实工作状态对凿岩钎杆进行了瞬态动力学分析，并对不同活塞形状下分析结果数据进行了比较，为冲击凿岩系统设计提供了理论依据。     

Modeling the seismic response of basin-edge structure by a sequential waveform method

A sequential waveform method is developed to simulate the seismic response of basin-edge structure excited by a plane incident P-wave. The full procedure involves: (a) a previous parameterization of the investigated model using the seismic wave velocities and depths of the sedimentary stratifications; (b) an input motion determined from the records at stations installed on hard rock; (c) forward computation of the P-SV elastic wave field by means of a two-dimensional finite difference (FD) method; (d) the optimization of the model vector using simulated annealing technique and comparing the simulated seismic response of the tested structure with the observed wave field; (e) the correction of the initial model by trial-and-error by testing the differences between synthetics and observed data, and (f) the final solution obtained by iteration using the conjugate gradient algorithm. The search of an optimal basin-edge model has been parallel processed by varying the shapes and velocities of strata on the basis of the fitting of relative timing, amplitude and phase between the output and the observed data. The input motion and sensitivity have been checked and the validity of the method has been demonstrated by numeric analysis. Using the teleseismic records generated by 7 earthquakes recorded at 26 broadband seismic stations, we have studied the seismic velocity structure of the southern edge of the Jiyang depression located in the Bohai Bay basin, northern China. Two cross sections show an agreement between the velocity results and the geological sections available in the region. In addition, we obtain evidence of three hidden faults under the sections and features that suggest major extensions at the Paleogene.     

Spherical Slepian functions and the polar gap in geodesy

The estimation of potential fields such as the gravitational or magnetic potential at the surface of a spherical planet from noisy observations taken at an altitude over an incomplete portion of the globe is a classic example of an ill-posed inverse problem. We show that this potential-field estimation problem has deep-seated connections to Slepian's spatiospectral localization problem which seeks bandlimited spherical functions whose energy is optimally concentrated in some closed portion of the unit sphere. This allows us to formulate an alternative solution to the traditional damped least-squares spherical harmonic approach in geodesy, whereby the source field is now expanded in a truncated Slepian function basis set. We discuss the relative performance of both methods with regard to standard statistical measures such as bias, variance and mean squared error, and pay special attention to the algorithmic efficiency of computing the Slepian functions on the region complementary to the axisymmetric polar gap characteristic of satellite surveys. The ease, speed, and accuracy of our method make the use of spherical Slepian functions in earth and planetary geodesy practical.     

Traction image method for irregular free surface boundaries in finite difference seismic wave simulation

In this study, we propose a new numerical method, named as Traction Image method, to accurately and efficiently implement the traction-free boundary conditions in finite difference simulation in the presence of surface topography. In this algorithm, the computational domain is discretized by boundary-conforming grids, in which the irregular surface is transformed into a 'flat' surface in computational space. Thus, the artefact of staircase approximation to arbitrarily irregular surface can be avoided. Such boundary-conforming gridding is equivalent to a curvilinear coordinate system, in which the first-order partial differential velocity-stress equations are numerically updated by an optimized high-order non-staggered finite difference scheme, that is, DRP/opt MacCormack scheme. To satisfy the free surface boundary conditions, we extend the Stress Image method for planar surface to Traction Image method for arbitrarily irregular surface by antisymmetrically setting the values of normal traction on the grid points above the free surface. This Traction Image method can be efficiently implemented. To validate this new method, we perform numerical tests to several complex models by comparing our results with those computed by other independent accurate methods. Although some of the testing examples have extremely sloped topography, all tested results show an excellent agreement between our results and those from the reference solutions, confirming the validity of our method for modelling seismic waves in the heterogeneous media with arbitrary shape topography. Numerical tests also demonstrate the efficiency of this method. We find about 10 grid points per shortest wavelength is enough to maintain the global accuracy of the simulation. Although the current study is for 2-D P-SV problem, it can be easily extended to 3-D problem.