Study on S wave velocity structure beneath part stations in Shanxi Province

Based on S wave records of deep teleseisms on Digital Seismic Network of Shanxi Province, shear wave velocity structures beneath 6 stations were obtained by means of S wave waveform fitting. The result shows that the crust is thick in the studied region, reaching 40 km in thickness under 4 stations. The crust all alternatives high velocity layer with low velocity one. There appear varied velocity structures for different stations, and the stations around the same tectonic region exhibit similar structure characteristics. Combined with dominant depth distribution of many small-moderate earthquakes, the correlation between seismogenic layers and crustal structures of high and low velocity layers has been discussed.

Study on S wave splitting in Dayao earth-quake sequence with M=6.2 and M=6.1 in Yunnan in 2003

The polarization direction of fast wave and the delay time between fast and slow wave were measured for two earthquake sequences occurred continuously on 21 July （M=6.2） and 16 October （M=6.1） in Dayao, Yunnan in 2003 using cross-correlation coefficient method, after determining the high-resolution hypocentral locations of the earthquake sequences using the double-difference earthquake location algorithm. The results indicated that ① The phenomena of S wave splitting are obvious in the two earthquake sequences, and the average polarization directions of fast wave in most stations are almost consistent with regional maximum horizontal compressive stress direction except the station Santai. There are bimodal fast directions in the polarization directions at station Santai and the mean polarization direction is N80°E, indicating an inconsistent phenomenon referred to regional maximum horizontal compressive stress direction. ② There is no apparent relation between delay time and focal depth in the sequences, but the polarization direction show different character in different delay time range. ③ The comparison of S wave splitting results in the two earthquake sequences show that the polarization direction in M=6.2 earthquake sequence is more scattered and its average fast direction is 20° larger than that of M=6.1 sequence, and the delay times between two sequences show a little difference. ④ The spatial variation in S wave splitting polarization direction may be due to the stress disturbance imposed by the M=6.2 and the M=6.1 mainshocks on regional background stress field.

S wave inelastic attenuation and site effect in the northern Tianshan area

In this paper, according to the relationship of observation spectrum with source spectrum, travel-path attenuation factor and site effect, we use the digital seismic data of moderate-small earthquakes to invert the soft rock site effect and S wave inelastic attenuation under the constraint of site effect, and determine the parameters of path and site. The parameters obtained by this method exclude the topographic effect of seismic stations, and truly reflect the seismic effect of free soft rock sites and the S wave inelastic attenuation in the studied region. The average soft rock site effect is about 1.5 times in the frequency domain of 2-4 Hz and the S wave quality factor is obtained as Qs=278f0.346, which can be directly applied to the stochastic modeling of ground motions in the studied region.

Frequency-dependent attenuation of P and S waves in Yunnan region

We analyzed digital seismogram data of 5668 earthquakes that occurred in Yunnan region between July of 1999 and December of 2003. Among the 22 seismic stations, six ones, namely Baoshan, Yongsheng, Lijiang, Heqing, Yimen, and Luquan, were selected and their attenuations of P and S waves were measured by using the extended coda-normalization method. The six stations were classified into three regions according to their location, that is, Baoshan area, Yongsheng-Lijiang-Heqing (YLH) area and Luquan-Yimen (LY) area. The values of QS-1 and QP-1 are expressed as QS-1=0.00867 f-0.86, QP-1=0.01155 f-0.93, QS-1=0.01824 f-0.92, QP-1=0.02288 f-0.92, and QS-1=0.01647 f-0.91, QP-1=0.02826 f-0.97 in Baoshan, YLH, LY areas respectively. The attenuation of YLH, LY are closer to each other, however, Baoshan area is apparently lower. Comparing attenuation in the three areas with other areas of the world using the same method, it is suggested that the attenuations of P and S waves in YLH and LY areas are close to Kanto of Japan, but much higher than southeast of South Korea. The QS-1 and QP-1 in Baoshan area are slightly higher than southeastern South Korea. Furthermore, the results indicate that our QS-1 in Yunnan area is close to others by analyzing the coda attenuation.

Elastic modulus extraction based on generalized pre-stack PP–PS wave joint linear inversion

Joint PP–PS inversion offers better accuracy and resolution than conventional P-wave inversion. P-and S-wave elastic moduli determined through data inversions are key parameters for reservoir evaluation and fluid characterization. In this paper, starting with the exact Zoeppritz equation that relates P-and S-wave moduli, a coefficient that describes the reflections of P-and converted waves is established. This method effectively avoids error introduced by approximations or indirect calculations, thus improving the accuracy of the inversion results. Considering that the inversion problem is ill-posed and that the forward operator is nonlinear, prior constraints on the model parameters and modified low-frequency constraints are also introduced to the objective function to make the problem more tractable. This modified objective function is solved over many iterations to continuously optimize the background values of the velocity ratio, which increases the stability of the inversion process. Tests of various models show that the method effectively improves the accuracy and stability of extracting P and S-wave moduli from underdetermined data. This method can be applied to provide inferences for reservoir exploration and fluid extraction.

Inversion of teleseismic waves at Shidao Seismographic Station

Teleseismic datasets at the Shidao Seismographic Station, located in the northwestern South China Sea, are used to determine the earth anisotropy and the vertical distribution pattern of the shear wave velocity by inversion approaches. The rotated correction function is applied to analyzing high quality SeS records from five earthquakes at distance of 25°-35° to obtain shear wave splitting parameters of the lithosphere. The result from the deepest earthquake among the five events indicates that the polarization of the fast shear wave is N94°E, which means the direction of extensional stress or the moving of the upper mantle mass in Xisha Islands is nearly west to east and confirms that the crust in this region is a transitional one and the driving force beneath the crust is from the moving mass consistent with the Eurasian plate. The anisotropy effective thickness is estimated about 100 km based on the time delay of 1.3 s between the fast and slow shear waves. The receiver function is applied to analyzing high quality P wave records from nine earthquakes at distance of 20°- 60° to obtain the vertical distribution pattern of shear wave velocity beneath the station. The result indicates that the crust could be divided into three layers： the uppermost crust （5 km above） is a velocity gradient zone consisting of several small layers, where the shear wave velocity increases from 1.5 to 3.5 km/s gradually; the 5 - 16 km depth interval also consiss of several small layers of which the mean velocity is about 3.8 km/s; and the lower crust （ 16.0 - 26. 5 km） is an obvious low velocity layer with a velocity of about 3.6 km/s. The buried depth of the Moho discontinuity is 26.5 kin, the mean velocity of the layers beneath the Moho is about 4.7 km/s and there is an obvious low velocity layer just beneath the Moho. Moreover, analysis of the arrival time of converted waves and the swinging variation of velocity around the initial model suggests that smaller layers in the model maybe are not reliable but the low velocity layer between 16 and 26.5 km maybe is the real one that implies the plasticity of the lower crust.

Electrocardiographic findings in leads V_(3R) to V_(5R) in patients with a slurred or notched S wave in lead V_1

Background Right bundle branch block （RBBB） may present as slurred or notched S wave in lead V1. However, slurred or notched S wave may also represent slow conduction in the myocardium. Methods We retrospectively analyzed the QRS patterns in leads VgR to V5R in 7 patients with a slurred or notched S wave in lead V1. Results In the leads V3R to VSR, 6 patients showed incomplete or complete RBBB and 1 patient slurred or notched S wave. Conclusions In the majority of ECGs in a small patient series with slurred or notched S wave in lead V1, QRS morphology indicating incomplete or complete RBBB was present in leads V3R to VSR. A finding of fragmented QRS in these leads may indicate slow conduction in the myocardium.

Smooth constraint inversion technique in genetic algorithms and its application to surface wave study in the Tibetan Plateau

Smooth constraint is important in linear inversion, but it is difficult to apply directly to model parameters in genetic algorithms. If the model parameters are smoothed in iteration, the diversity of models will be greatly suppressed and all the models in population will tend to equal in a few iterations, so the optimal solution meeting requirement can not be obtained. In this paper, an indirect smooth constraint technique is introduced to genetic inversion. In this method, the new models produced in iteration are smoothed, then used as theoretical models in calculation of misfit function, but in process of iteration only the original models are used in order to keep the diversity of models. The technique is effective in inversion of surface wave and receiver function. Using this technique, we invert the phase velocity of Raleigh wave in the Tibetan Plateau, revealing the horizontal variation of S wave velocity structure near the center of the Tibetan Plateau. The results show that the S wave velocity in the north is relatively lower than that in the south. For most paths there is a lower velocity zone with 12-25 km thick at the depth of 15-40 km. The lower velocity zone in upper mantle is located below the depth of 100 km, and the thickness is usually 40-80 km, but for a few paths reach to 100 km thick. Among the area of Ando, Maqi and Ushu stations, there is an obvious lower velocity zone with the lowest velocity of 4.2-4.3 km/s at the depth of 90-230 km. Based on the S wave velocity structures of different paths and former data, we infer that the subduction of the Indian Plate is delimited nearby the Yarlung Zangbo suture zone.

Q value of anelastic S-wave attenuation in Yunnan region

The study of seismic attenuation property is a major subject in seismology. Seismic waves recorded by seismic stations （seismographs） contain source effect, seismic wave propagation effect, site response of seismic stations and instrumental response. The path effect of seismic wave propagation, site response of seismic stations and instrumental response must be taken out in the study of source property with seismic data. The path effect of seismic wave propagation （seismic attenuation） involves an important influential factor, the anelastic attenuation of medium, which is measured with quality factor Q, apart from geometric attenuation with the distance. As a basic physical parameter of the Earth medium, Q value is essential for quantitative study of earthquakes and source property （e.g. determination of source parameters）, which is widely used in earthquake source physics and engineering seismology.

S-wave velocity structure inferred from receiver function inversion in Tengchong volcanic area

Tengchong volcanic area is located near the impinging and underthrust margin of India and Eurasia plates. The volcanic activity is closely related to the tectonic environment. The deep structure characteristics are inferred from the receiver function inversion with the teleseismic records in the paper. The results show that the low velocity zone is influenced by the NE-trending Dayingjiang fault. The S-wave low velocity structure occurs obviously in the southern part of the fault, but unobviously in its northern part. There are low velocity zones in the shallow po-sition, which coincides with the seismicity. It also demonstrates that the low velocity zone is directly related to the thermal activity in the volcanic area. Therefore, we consider that the volcano may be alive again.

Study on Shear Wave Velocity Structure and Velocity Ratio Beneath Ordos Block and Its Eastern and Southern Margins

Using pure S wave fitting method, we studied the shear wave velocity structures under the Ordos block and its eastern and southern marginal areas. The results show that the velocity structure beneath Yulin station in the interior of Ordos block is relatively stable, where no apparent change between high and low velocity layers exists and the shear wave velocity increases steadily with the depth. There is a 12km thick layer at the depth of 25km under this station, with an S wave velocity (V S=3.90km/s) lower than that at the same depth in its eastern and southern areas (V S≥4.00km/s). The crust under the eastern margin of Ordos block is thicker than that of the Yulin station, and the velocity structures alternate between the high and low velocity layers, with more low velocity layers. It has the same characteristic as having a 10km-thick low velocity layer (V S=3.80km/s) in the lower crust but buried at a depth of about 35km. Moreover, we studied the V P/V S ratio under each station in combination with the result of P wave velocity inversion. The results show that, the average velocity ratio of the Yulin station at the interior of Ordos block is only 1.68, with a very low ratio (about 1.60) in the upper crust and a stable ratio of about 1.73 in the mid and lower crust, which indicates the media under this station is homogenous and stable, being in a state of rigidity. But at the stations in the eastern and southern margins of the Ordos block, several layers of high velocity ratio (about 1.80) have been found, in which the average velocity ratio under Kelan and Lishi stations at the eastern margin is systemically higher than that of the general elastical body waves (1.732). This reflects that the crust under the marginal areas is more active relatively, and other materials may exist in these layers. Finally, we discussed the relationship among earthquakes, velocity structures beneath stations and faults.

Interpretation of S-wave data from Tai′anXinzhou DSS profile and its relationship with Xingtai earthquakes

On the basis of S wave information from Tai′an Xinzhou DSS profile and with reference to the results from P wave interpretation, the 2 D structures , including S wave velocity V s, ratio γ between V p and V s; and Poisson′s ratio σ , are calculated; the structural configuration of the profile is presented and the relevant inferences are drawn from the above results. Upwarping mantle districts (V s≈4.30 km/s)and sloping mantle districts (V s≈4.50 km/s) of the profile with velocity difference about 4% at the top of upper mantle are divided according to the differences of V s , γ and σ in different media and structures, also with reference to the information of their neighbouring regions; the existence of Niujiaqiao Dongwang high angle ultra crustal fault zone is reaffirmed; the properties of low and high velocity blocks(zones) including the crust mantle transitionalzone and the boudary indicators of North China rift valley are discussed. A comprehensive study is conducted on the relation of the interpretation results with earthquakes. It is concluded that the mantle upwarps, thermal material upwells through the high angle fault, the primary hypocenter was located at the crust mantle juncture 30.0～33.0 km deep, and additional stress excited the M S=6.8 and M S=7.2 earthquakes at specific locations around 9.0 km below Niujiaqiao Dongwang, the earthquakes took place around the high angle ultra crustal fault and centered in the brittle media and rock strata with low γ and low σ values.

Inversion of waveforms from Xiangtang borehole seismic array for soil dynamic property

In order to understand the site soil response of the Xiangtang borehole seismic array under real strong ground motion, reveal the site response, verify the technique of borehole exploration, and improve the precision of in-situ test and laboratory test, this paper presents a new approach, which is composed of two methods. One is the layered site seismic response method, whose layer transform matrix is always real. The other is a global-local optimization technique, which uses genetic algorithm （GA）-simplex method. An inversion of multi-component waveforms of P, SV and SH wave is carried out simultaneously. By inverting the records of three moderate and small earthquakes obtained from the Xiangtang borehole array （2^#） site, the soil dynamic characteristic parameters, including P velocity, damping ratio and frequency-dependent coefficient b, which has not been given in previous literatures, are calculated. The results show that the soil S wave velocity of the Xiangtang 2^# borehole is generally greater than that obtained from the 1994 in-situ test, and is close to the velocity of the 3^# borehole, which is more than 200 m away from the 2^# borehole. Meanwhile, perceptible soil nonlinear behavior under peak ground motion of about 60×10^-2 m/s^2 is detected by the inversion analysis. The presented method can be used for studying the soil response of other borehole array sites.

Seismic waves in cracked two-phase anisotropic medium with fourfold rotation symmetry

In this paper, based on the propagation theorics of seismic waves in anisotropic medium and in cracked two-phase medium, the constitutive relations and dynamic equations of the propagation of seismic waves in Cracked twophase anisotropic medium with fourfold rotation symmetry have been derived, and the preliminary theoretical analysis have been made for plane wave as an example.

基于程函面波成像的华北克拉通岩石圈三维高分辨率S波速度结构研究

利用“中国地震科学台阵探测”项目Ⅱ期和Ⅲ期的流动地震台站以及中国区域地震台网中的部分固定台站的观测资料,采用程函面波成像方法获得了华北克拉通及周边区域10—120 s周期的瑞雷面波相速度分布和高分辨率的三维S波速度结构,并基于该速度模型估算了岩石圈厚度分布。结果显示,华北克拉通内部岩石圈厚度除了存在“西厚东薄”的一级分布特征外,还存在一些更小尺度的差异,包括鄂尔多斯地块内部岩石圈“南厚北薄”、鄂尔多斯地块周缘断陷带岩石圈显著的不均匀减薄以及燕山构造带与其南侧华北平原之间的显著差异等。山西断陷带北部与南部地区上地幔浅部(<100 km)存在不同程度的低速异常,它们被中部的高速异常区所分隔。在150 km以下深度从太行山南缘向北至山西断陷北缘存在一条NNE向展布的显著低速异常带,表明上地幔浅部南北部的低速异常在深部相连。结合已有的其它成像结果,我们推测这些低速异常起源于更深处(>200 km),并与由太平洋俯冲板块的滞留脱水导致的上地幔热物质上涌和小尺度地幔对流等密切相关。燕山构造带与华北平原的岩石圈结构存在明显差异,前者遭受的岩石圈破坏改造程度明显弱于后者,张家口—渤海地震带位于这两种不同壳幔结构的过渡带,地震活动较强,我们认为深部结构和热作用的显著差异,以及青藏高原远场挤压效应的共同作用是导致该区地震活动较强的主要原因。

Convective Schrodinger Equation:Insights on the Potential Energy’s Role to Wave Particle Decay

In this paper, we coupled the Quantum Mechanics conventional Schr&#246;dinger’s equation, for the particles, with the Maxwell’s wave equation, in order to study the potential’s role on the conversion of the electromagnetic field energy to mass and vice versa. We show that the dissipation (“conductivity”) factor and the particle implicit proper frequency are both related to the potential energy. We have also derived a new expression for the Schr&#246;dinger’s Equation considering the potential energy into this equation not as an ad hoc term, but also as an operator (Hermitian), which has the scalar potential energy as a natural eigenvalue of this operator.

Crustal S-wave velocity structure across the northeastern South China Sea continental margin: implications for lithology and mantle exhumation

The northeastern margin of the South China Sea (SCS), developed from continental rifting and breakup, is usually thought of as a non-volcanic margin. However, post-spreading volcanism is massive and lower crustal high-velocity anomalies are widespread, which complicate the nature of the margin here. To better understand crustal seismic velocities, lithology, and geophysical properties, we present an S-wave velocity (VS) model and a VP/VS model for the northeastern margin by using an existing P-wave velocity (VP) model as the starting model for 2-D kinematic S-wave forward ray tracing. The Mesozoic sedimentary sequence has lower VP/VS ratios than the Cenozoic sequence;in between is a main interface of P-S conversion. Two isolated high-velocity zones (HVZ) are found in the lower crust of the continental slope, showing S-wave velocities of 4.0–4.2 km/s and VP/VS ratios of 1.73–1.78. These values indicate a mafic composition, most likely of amphibolite facies. Also, a VP/VS versus VP plot indicates a magnesium-rich gabbro facies from post-spreading mantle melting at temperatures higher than normal. A third high-velocity zone (VP : 7.0–7.8 km/s;VP/VS: 1.85–1.96), 70-km wide and 4-km thick in the continent-ocean transition zone, is most likely to be a consequence of serpentinization of upwelled upper mantle. Seismic velocity structures and also gravity anomalies indicate that mantle upwelling/ serpentinization could be the most severe in the northeasternmost continent-ocean boundary of the SCS. Empirical relationships between seismic velocity and degree of serpentinization suggest that serpentinite content decreases with depth, from 43% in the lower crust to 37% into the mantle.

Forward modeling of ocean-bottom cable data and wave-mode separation in fluid–solid elastic media with irregular seabed

In marine seismic exploration, ocean-bottom cable techniques accurately record the multicomponent seismic wavefield; however, the seismic wave propagation in fluid–solid media cannot be simulated by a single wave equation. In addition, when the seabed interface is irregular, traditional finite-difference schemes cannot simulate the seismic wave propagation across the irregular seabed interface. Therefore, an acoustic–elastic forward modeling and vector-based P-and S-wave separation method is proposed. In this method, we divide the fluid–solid elastic media with irregular interface into orthogonal grids and map the irregular interface in the Cartesian coordinates system into a horizontal interface in the curvilinear coordinates system of the computational domain using coordinates transformation. The acoustic and elastic wave equations in the curvilinear coordinates system are applied to the fluid and solid medium, respectively. At the irregular interface, the two equations are combined into an acoustic–elastic equation in the curvilinear coordinates system. We next introduce a full staggered-grid scheme to improve the stability of the numerical simulation. Thus, separate P-and S-wave equations in the curvilinear coordinates system are derived to realize the P-and S-wave separation method.

Scattering of Tollmien-Schlichting waves as they pass over forward-/backward-facing steps

Forward-/backward-facing steps in boundary-layer flows are often seen in engineering applications, and they have potential impacts on laminar-turbulent transition through scattering of the oncoming instability modes(e.g., Tollmien-Schlichting(T-S) waves). This issue is studied in the present paper by applying a local scattering framework, which is a rather generic mathematical framework on describing the mode scattering process. In this framework, a high-Reynolds-number triple-deck formalism is employed, and a transmission coefficient, defined as the ratio of the asymptotic amplitude of the instability mode downstream of the step to that upstream, is introduced. Through the systematical study, it has been found that both the forward-and backward-facing steps have a destabilizing effect on the oncoming T-S waves in subsonic boundary layers, this effect increases with the height of the step and/or the frequency of the T-S wave, and a backward-facing step(BFS) always has a greater impact than a forward-facing step(FFS). These facts agree with most of the previous investigations.However, one numerical study(WORNER, A., RIST, U., and WAGNER, S. Humps/steps influence on stability characteristics of two-dimensional laminar boundary layer. AIAA Journal, 41, 192–197(2003)), which was based on an ad-hoc configuration, showed an opposite impact of an FFS. Through the investigation on the specific configuration, it is revealed that the wrong conclusion was drawn by misinterpreting the numerical results.

Tunnel effect of fractal fault and transient S-wave velocity rupture (TSVR) of in-plane shear fault

Transient S wave velocity rupture (TSVR) means the velocity of fault rupture propagation is between S wave velocity α and P wave velocity β . Its existing in the rupture of in plane ( i.e . strike slip) fault has been proved, but in 2 dimensional classical model, there are two difficulties in transient S wave velocity rupture, i.e ., initialization difficulty and divergence difficulty in interpreting the realization of TSVR. The initialization difficulty means, when v ↑ v R (Rayleigh wave velocity), the dynamic stress strength factor K 2(t) →+0, and changes from positive into negative in the interval ( v R, β ). How v transit the forbidden of ( v R, β )? The divergence difficulty means K 2(t) →+∞ when v ↓ β . Here we introduce the concept of fractal and tunnel effect that exist everywhere in fault. The structure of all the faults is fractal with multiple cracks. The velocity of fault rupture is differentiate of the length of the fault respect to time, so the rupture velocity is also fractal. The tunnel effect means the dynamic rupture crosses over the interval of the cracks, and the coalescence of the intervals is slower than the propagation of disturbance. Suppose the area of earthquake nucleation is critical or sub critical propagation everywhere, the arriving of disturbance triggers or accelerates the propagation of cracks tip at once, and the observation system cannot distinguish the front of disturbance and the tip of fracture. Then the speed of disturbance may be identified as fracture velocity, and the phenomenon of TSVR appears, which is an apparent velocity. The real reason of apparent velocity is that the mathematics model of shear rupture is simplified of complex process originally. The dual character of rupture velocity means that the apparent velocity of fault and the real velocity of micro crack extending, which are different in physics, but are unified in rupture criterion. Introducing the above mentioned concept to the calculation of K 2 (t) , the difficulty of initialization can be overcome, and the integral equation of triggering the initialization of TSVR is given quantitatively. By solving this integral equation, the lower limit of TSVR is 1.105 3 β , not β , and the divergence difficulty is overcome. TSVR is unstable solution, and may degenerate to sub Rayleigh wave velocity rupture immediately where the non critical condition can be measured. The results of this paper show that the initialization and continuum depends on the condition of earthquake nucleation in seismogenic area.