Double-difference tomography of P- and S-wave velocity structure beneath the western part of Java, Indonesia

West Java in the western part of the Sunda Arc has a relatively high seismicity due to subduction activity and faults.In this study,double-difference tomography was used to obtain the 3D velocity tomograms of P and S waves beneath the western part of Java.To infer the geometry of the structure beneath the study area,precise earthquake hypo・center determination was first performed before tomographic imaging.For this,earthquake waveform data were extracted from the regional Meteorological,Climatological,Geophysical Agency(BMKG)network of Indonesia from South Sumatra to Central Java.The P and S arrival times for about 1,000 events in the period April 2009 to July 2016 were selected,the key features being events of magnitude>3,azimuthal gap<210°and number of phases>8.A nonlinear method using the oct-tree sampling algorithm from the NonLinLoc program was employed to determine the earthquake hypocenters.The hypocenter locations were then relocated using double-difference tomography(tomoDD).A significant reduction of travel-time(root mean square basis)and a better clustering of earthquakes were achieved which correlated well with the geological structure in West Java.Double-difference tomography was found to give a clear velocity structure,especially beneath the volcanic arc area,i.e.,under Mt Anak Krakatau,Mt Salak and the mountains complex in the southern part of West Java.Low velocity anomalies for the P and S waves as well as the vp/vs ratio below the volcanoes indicated possible partial melting of the upper mantle which ascended from the subducted slab beneath the volcanic arc.

The strain seismograms of P- and S-waves of a local event recorded by four-gauge borehole strainmeter

At a sampling rate of 100 samples per second,the YRY-4 four-gauge borehole strainmeters（FGBS） are capable of recording transient strains caused by seismic waves such as P and S waves or strain seismograms. At such a high sampling rate, data from the YRY-4 strainmeters demonstrate fairly satisfactory self-consistency. The strain tensor seismograms demonstrate the senses of motion of P waves, that is, the type of seismic wave travels in the direction of the maximum normal strain change. The observed strain patterns of S waves significantly differ from those of P waves and should contain information about the source mechanism. Spectrum analysis shows that the strain seismograms are consistent with conventional broadband seismograms from the same site.

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.

Modeling of Surface Waves in a Fluid Saturated Poro-Elastic Medium under Initial Stress Using Time-Space Domain Higher Order Finite Difference Method

In this present context, mathematical modeling of the propagation of surface waves in a fluid saturated poro-elastic medium under the influence of initial stress has been considered using time dependent higher order finite difference method (FDM). We have proved that the accuracy of this finite-difference scheme is 2M when we use 2nd order time domain finite-difference and 2M-th order space domain finite-difference. It also has been shown that the dispersion curves of Love waves are less dispersed for higher order FDM than of lower order FDM. The effect of initial stress, porosity and anisotropy of the layer in the propagation of Love waves has been studied here. The numerical results have been shown graphically. As a particular case, the phase velocity in a non porous elastic solid layer derived in this paper is in perfect agreement with that of Liu et al. (2009).

Propagation of P-and S-waves in initially stressed gravitating half space

The present paper contributes in studying the phase velocities of P- and S-waves in a half space subjected to a compressive initial stress and gravity field. The density and acceleration due to gravity vary quadratically along the depth. The dispersion equation is derived in a closed form. It is shown that the phase velocities depend not only on the initial stress, gravity, and direction of propagation but also on the inhomogeneity parameter associated with the density and acceleration due to gravity. Various particular cases are obtained, and the results match with the classical results. Numerical investigations on the phase velocities of P- and S-waves against the wave number are made for various sets of values of the material parameters, and the results are illustrated graphically. The graphical user interface model is developed to generalize the effect.

Assessments of Elastic Anisotropy of Banded Amphibolite as a Function of Cleavage Orientation Using S- and P-Wave Velocity

As most rocks are of an anisotropic nature, the measurement elastic modulus of anistropic rocks is of general interest. Nevertheless, uniaxial compression test is common method to measure the dynamic elastic constants of anisotropic rocks;the use of ultrasonic pulse test is attractive, because the test is non-destructive and easy to apply. This study aimed to demonstrate the influence of orientation of foliation planes of banded amphibolite rocks on the compressional (Vp), shear wave (Vs) velocities propagating and elastic modules using ultrasonic pulse test. The result showed that the planes of foliation have a major effect on the wave velocity, where the Vp and Vs were taken parallel to the foliation plane show higher values than those obtained in the other directions (β = 30。, 60。 and 90。). The preliminary conclusions are developed concerning that the elastic modulus is vary continuously as a function of cleavage orientation with respect to the direction of wave propagations, where Poisson’s ratio having the smallest relative change. The highest values of Young’s modulus and shear modulus are observed for foliation dip angles of 0? and the lowest values are for foliation dip angles of 90。. This indicates that the observed intrinsic anisotropy and the close relations of the directional dependent seismic anisotropy to the foliation planes are mainly a result of crystallographic preferred orientation of major minerals (e.g. horn- blende and elongated quarts grains).

Quantum Mechanics: Harmonic Wave-Packets, Localized by Resonant Response in Dispersion Dynamics

From a combination of Maxwell’s electromagnetism with Planck’s law and the de Broglie hypothesis, we arrive at quantized photonic wave groups whose constant phase velocity is equal to the speed of light c = ω/k and to their group velocity dω/dk. When we include special relativity expressed in simplest units, we find that, for particulate matter, the square of rest mass , i.e., angular frequency squared minus wave vector squared. This equation separates into a conservative part and a uniform responsive part. A wave function is derived in manifold rank 4, and from it are derived uncertainties and internal motion. The function solves four anomalies in quantum physics: the point particle with prescribed uncertainties;spooky action at a distance;time dependence that is consistent with the uncertainties;and resonant reduction of the wave packet by localization during measurement. A comparison between contradictory mathematical and physical theories leads to similar empirical conclusions because probability amplitudes express hidden variables. The comparison supplies orthodox postulates that are compared to physical principles that formalize the difference. The method is verified by dual harmonics found in quantized quasi-Bloch waves, where the quantum is physical;not axiomatic.

Waves induced by a submerged moving dipole in a two-layer fluidof finite depth

The waves induced by a moving dipole in a two-fluid system are analytically and experimentally investigated. The velocity potential of a dipole moving horizontally in the lower layer of a two-layer fluid with finite depth is derived by superposing Greens functions of sources (or sinks). The far-field waves are studied by using the method of stationary phase. The effects of two resulting modes, i.e. the surface- and internal-wave modes, on both the surface divergence field and the interfacial elevation are analyzed. A laboratory study on the internal waves generated by a moving sphere in a two-layer fluid is conducted in a towing tank under the same conditions as in the theoretical approach. The qualitative consistency between the present theory and the laboratory study is examined and confirmed.

Diffraction of plane P waves by a canyon of arbitrary shape in poroelastic half-space (Ⅰ): Formulation

This paper presents an indirect boundary integration equation method for diffraction of plane P waves by a two-dimensional canyon of arbitrary shape in poroelastic half-space. The Green＇s functions of compressional and shear wave sources in poroelastic half-space are derived based on Biot＇s theory. The scattered waves are constructed using the fictitious wave sources close to the boundary of the canyon, and magnitude of the fictitious wave sources are determined by the boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, the comparison between the degenerated solutions of single-phased half-space and the well-known solutions, and the numerical stability of the method.

Numerical Simulation of Wave Field near Submerged Bars by PLIC-VOF Model

Investigating the wave field near structures in coastal and offshore engineering is of increasing significance. In the present study, simulation is done of the wave profile and flow field for waves propagating over submerged bars, using PLIC-VOF （Pieeewise Linear Interface Construction） to trace the free surface of wave and finite difference method to solve vertical 2D Navier-Stokes （N-S） equations. A comparison of the numerical results for two kinds of submerged bars with the experimental ones shows that the PLIC-VOF model used in this study is effective and can compute the wave field precisely.

The S to P convert wave from the bottom of sediment basin in the near-field seismic records

From the near-field records of aftershock of Datong earthquakes in October 1989 and March 1991,an extra phase between P and S arrivals is found. High-precision epicenter location shows that some of these records are obtained with the incidental angle less than the critical angle. This excludes the possibility that the extra wave phase is a refractive wave from ground surface. Particle motion analysis shows that the characteristic of the extra wave is similar to that of wave, so it is possible that the extra phase is an S to P convert wave from the bottom of sediment basin. Suppose a low velocity layer covers a high velocity basement. Successful simulation by the synthetic seismogram confirms that the extra phase is an S-P convert wave from the interface of the basin bottom. Modifying the depth of interface at each ray path to match the waveform, we obtain an interface distribution in space. In this way, a brief bottom image could be shown, and the Datong basin has a 'V' shaped bottom.

Pure S-waves in land P-wave source VSP data

常规土地垂直地震介绍(VSP ) 探索通常为接收装置使用 P 波浪来源和三部件的地震检波器，强调 P -- 并且变换 S 波浪。以前的研究证明从在表面的可控制的地震来源的两炸毁地上凿穿射击和垂直颤动将生产相对强壮的纯 P 波浪和更弱的纯 S 波浪。大泊松的比率差别的接口在强壮的播送变换 S 波浪的形成上有积极影响。由从来源的纯 S 波浪的比较分析并且在去的 S 波浪下面变换了，当时，我们相信纯 S 波浪的主要频率通常比纯 P 波浪低主要频率在变变换下面， S 波浪接近 P 波浪的。我们学习了零偏移量并且从陆地 P 波浪的偏移量 VSP 数据采购原料。结果证明纯 S 波浪通常在波浪紧张与差别在这些数据存在。S 波浪速度能从 P 波浪来源被获得零偏移量的 VSP 数据。最后，我们讨论 VSP P 的联合申请的明亮的未来 -- 并且 S 波浪和在 P 波浪的 S 波浪的完整的使用采购 VSP 数据。

Estimation of Shallow S-Wave Velocity Structure of Two Practical Sites from Microtremors Array Observation in Tangshan Area

Microtremors array observation for estimating S-wave velocity structure from phase velocities of Rayleigh and Love wave on two practical sites in Tangshan area by a China-US joint group are researched.The phase velocities of Rayleigh wave are estimated from vertical component records and those of Love wave are estimated from three-component records of microtremors array using modified spatial auto-correlation method.Haskell matrix method is used in calculating Rayleigh and Love wave phase velocities,and the shallow S-wave velocity structure of two practical sites are estimated by means of a hybrid approach of Genetic Algorithm and Simplex.The results are compared with the PS logging data of the two sites,showing it is feasible to estimate the shallow S-wave velocity structure of practical site from the observation of microtremor array.

A new method of testing pile using dynamic P-S-curve made by amplitude of wave train

A new method of detecting the vertical bearing capacity for single-pile with high strain is discussed in this paper. A heavy hammer or a small type of rocket is used to strike the pile top and the detectors are used to record vibra- tion graphs. An expression of higher degree of strain (deformation force) is introduced. It is testified theoretically that the displacement, velocity and acceleration cannot be obtained by simple integral acceleration and differential velocity when long displacement and high strain exist, namely when the pile phase generates a whole slip relative to the soil body. That is to say that there are non-linear relations between them. It is educed accordingly that the force P and displacement S are calculated from the amplitude of wave train and (dynamic) P-S curve is drew so as to determine the yield points. Further, a method of determining the vertical bearing capacity for single-pile is dis- cussed. A static load test is utilized to check the result of dynamic test and determine the correlative constants of dynamic-static P(Q)-S curve.

Breaking of a Langmuir wave in cold electron-positron-ion plasmas

The space-time evolution of a given density perturbation in cold homogeneous electron- positron-ion plasma is investigated with an assumption of infinitely massive ions by employing a numerical calculation method. The phase-mixing time and wave-breaking time can be effectively distinguished with this method. It is found that an increase of the ratio of equilibrium ion density to equilibrium electron density can attenuate plasma oscillations, leading to a delay in wave breaking. The dependence of the phase-mixing and wave-breaking times on the amplitude of the initial perturbation is also discussed.

v_p/v_s ESTIMATION IN TIBETAN CRUST FROM INVERSION OF SURFACE WAVE DISPERSIONS

We estimated crustal v p/ v s ratio of Tibetan Plateau by combined inversion of Love and Rayleigh wave dispersion data. It is developed by us that the joint inversion methods using both Love and Rayleigh wave dispersion data. Thickness and S\|wave velocity of each sub layer are taken from Love wave dispersion data, then P\|wave velocity structure was deduced using Rayleigh wave dispersion data. Densities of sub layers were estimated by the empirical relationships between seismic velocity and rock density. Having S\| and P\|wave velocities, v p/ v s ratio is calculated for each sub layer. Six sub layers in crust of Tibetan Plateau has been identified, which are 0～8km, 8～30km, 30～40km, 40～62km, 62～68km and 68～75km respectively. The S\|wave velocity structure of the Plateau is 3 13, 3 32, 3 15, 3 92, 3 45 and 3 87 km/s for each sub layer; and P\|wave velocities are 6 00, 6 10, 5 72, 6 35, 6 78 and 6 64km/s respectively v p/ v s ratios in sub layers are 1 92, 1 84, 1 82, 1 62, 1 96 and 1 72; and corresponding Poisson ratios are 0 31, 0 29, 0 28, 0 19, 0 32 and 0 24. Our result on Poisson ratios of Tibetan crust was supported by seismic waveform modelling by Rodgers and Schwartz (1998).

Experimental Investigation of Lower Hybrid Wave Coupling in HT-7

Experiments on lower hybrid wave （LHW） coupling were investigated in the HT- 7 tokamak. Good coupling of LHW plasma has been demonstrated at different conditions in the HT-7 tokamak. Relevant results have proved that LHW-plasma coupling is affected by the phase difference between adjacent waveguides. Furthermore, the edge density around the grill and relevant coupling can be adjusted by changing the plasma line average density or the gap value between the LH grill and the last closed flux surfaces （LCFS）. It is found that the coupling of LHWs becomes poor when the edge density around the LH grill is large enough in the HT-7 tokamak, and that coupling remains good with a proper edge density. With increasing LHW power, it is also found that the reflection coefficients （RCs） increase due to non-linear effects under conditions of low edge recycling, but can decrease under high edge recycling. The edge density depends mainly on the competition between the ponderomotive force （PMF） and the edge recycling intensity in the HT-7 tokamak.

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.

Propagation of Love waves in non-homogeneous substratum over initially stressed heterogeneous half-space

The paper studies the propagation of Love waves in a non-homogeneous substratum over an initially stressed heterogeneous half-space. The dispersion equation of phase velocity is derived. The velocities of Love waves are calculated numerically as a function of kH and presented in a number of graphs, where k is the wave number, and H is the thickness of the layer. The case of Gibson＇s half-space is also considered. It is observed that the speed of Love waves is finite in the vicinity of the surface of the half-space and vanishes as the depth increases for a particular wave number. It is also observed that an increase in compressive initial stresses causes decreases of Love waves velocity for the same frequency, and the tensile initial stress of small magnitude in the half-space causes increase of the velocity.

Analysis of wave propagation in a functionally graded nanobeam resting on visco-Pasternak＇s foundation

Wave propagation analysis for a functionally graded nanobeam with rectangular cross-section resting on visco-Pasternak＇s foundation is studied in this paper. Timoshenko＇s beam model and nonlocal elasticity theory are employed for formulation of the problem. The equations of motion are derived using Hamilton＇s principals by calculating kinetic energy, strain energy and work due to viscoelastic foundation. The effects of various parameters such as wavenumber, non-homogeneous index, nonlocal parameter and three parameters of foundation are performed on the phase velocity of the nanobeam. The obtained results indicate that some parameters such as non-homogeneous index, nonlocal parameter and wavenumber have significant effect on the response of the system.