Tomography for Q of the Eastern Section of the Tianshan Area from High-frequency Attenuation of S-wave

Based on the waveform data of 5,076 local earthquakes recorded at 25 stations in Xinjiang during the period from 2009 to 2014 and the observation reports provided by the Xinjiang Digital Seismic Network,a data set of 19,140 attenuation factors t*is obtained by fitting the high-frequency attenuation of S-wave spectra with a genetic algorithm. The spatial distribution of Q_S is determined by inverting the t*data with seismic tomography. The results show that the average Q0 in eastern Tianshan is 520,and there is a significant correlation between the Q_S value distribution or attenuation characteristics it disclosed and the surface structure of the study area. The Q_S value is lower in the intersection area of the mountain basin which is located on the north and south sides of the Tianshan Mountains,and the high Q_S distribution is more concentrated inside the Tianshan orogenic belt. The M≥6. 0 earthquakes have been basically located in the Low-Q_S region since 1900. 24 high heat flow points in eastern Tianshan are located at the north and south of Tianshan Mountains where low Q_S exists,indicating a negative correlation. In addition,there is a positive correlation between the velocity structure and the attenuation structure in the study area,which reflects the consistency of the 2-D attenuation structure with the velocity structure and the two-dimensional density structure.     

Anisotropic zoning in the upper crust of the Tianshan Tectonic Belt

The Tianshan Tectonic Belt is an intracontinental orogenic belt formed by continental convergence that has undergone long-term tectonic evolution. The reactivation that began during the Cenozoic Period has led to complex structural changes. The goals of this study are to review the seismic observational data obtained during 2009–2019 in the Xinjiang regional seismic network and analyze the anisotropy of the upper crust in the Tianshan area. Therefore, a shear-wave splitting system was adopted to collect and analyze shear-wave splitting parameters of 33 stations in the study area. The anisotropy of the upper crust of the Tianshan is spatially diverse, and the dominant polarization directions of fast shear-wave reflect the spatial variations of regional tectonic stress. In addition, the time delays of slow waves are proportional to the intensities of anisotropy in the upper crustal medium. The dominant polarization direction of the fast waves in the western segment of the North Tianshan Mountain,northwestern corner of the Tarim Basin, and northeastern edge of Pamir is consistent with the tectonic stress fields in the area. In the northern part of the Keping Block, the dominant polarization directions of the fast waves are consistent with the fault trends;however, they are at a high angle to the dominant directions of the regional tectonic stress field indicating that the anisotropy is affected by the faults in the area. The anisotropy of the eastern segment of the South Tianshan Mountains and the surrounding area of Urumqi are affected by the local stress field and fault structure. The polarization directions at some of the stations are subparallel to the directions of the regional principal stress. However, for other stations, the polarization directions are aligned with the neighboring faults. The polarization directions of the fast waves in most of the study area are consistent with the local tectonic stress fields. Thus, stress compression phenomena such as the Tarim Basin being thrusted and subducted between the Tianshan crust and the upper mantle due to the far field effects of the convergence between the Indian and Siberian plates are evident.Furthermore, the zoning of the time delays is distinct, and the time delays share an increasing trend from east to west in the North Tianshan and South Tianshan Mountain ranges. These results are consistent with the north-south convergence deformations across the Tianshan Mountains, where the deformation rate increased from east to west. The average values of time delays in northeastern Pamir are significantly higher than that in the other areas due to the occurrence of the most intensive tectonic movements suggesting that the anisotropy of the zone is significantly stronger than that of the other zones in the Tianshan Tectonic Belt. We successfully deciphered the seismic anisotropy in the upper crust and provided a comprehensive and systematic understanding of the dynamic mechanisms of the Tianshan Tectonic Belt.     

远震尾波地震矩的测定与全球剪切应力分布(英文)

Seismic moments are determined from coda waves of moderate and large earthquakes recorded by the broad-band SK instrument at large epicentral distance. These moments from coda waves are consistent with those given by Dziewonski and Woodhouse(1983), which come from GDSN digitized recorded data, the average of the differences being 0.24 order of magnitude. When the coda-obtained M0 values are compared to Ms, the relation obtained between log(Mo) and Ms is in good agreement with that given by Hanks and Kanamori (1979).The shear stress distribution on a global scale is studied by using the data of seismic moments and body wave magnitudes mb. It is found that the Western Circum-pacific region and the Hindu-kush and Assam areas, which are located in the thrust regions of the Indian plate, have very high shear stress together with relatively high seismicity. The rest of the seismic areas of the world have intermediate or low shear stress.     

Application of the Levenberg-Marquardt Method in Improving Technology of the Seismic Influence Field

Attenuation relation of seismic intensity is fitted by using 152 complete isoseismic lines from 65 earthquakes which were greater than M_S4.0 from 1940 to 2015 in Inner Mongolia and neighboring regions. Meanwhile, based on the difference of land form and geological structure characteristics,the study area is the divided into eastern, western and central. The intensity attenuation relationships are established separately by using the LM (Levenberg-Marquardt) method and elliptic attenuation model. Comparisons are made by using the earthquake affecting field model of North China and computed results. The analytical study shows that the relation obtained represents the features of earthquake damage distribution in the areas well, and has positive meaning to guide loss assessments immediately after a destroyed earthquake.     

A wave propagation model with the Biot and the fractional viscoelastic mechanisms

Energy loss in porous media containing fluids is typically caused by a variety of dynamic mechanisms.In the Biot theory,energy loss only includes the frictional dissipation between the solid phase and the fluid phase,resulting in underestimation of the dispersion and attenuation of the waves in the low frequency range.To develop a dynamic model that can predict the high dispersion and strong attenuation of waves at the seismic band,we introduce viscoelasticity into the Biot model and use fractional derivatives to describe the viscoelastic mechanism,and finally propose a new wave propagation model.Unlike the Biot model,the proposed model includes the intrinsic dissipation of the solid frame.We investigate the effects of the fractional order parameters on the dispersion and attenuation of the P-and S-waves using several numerical experiments.Furthermore,we use several groups of experimental data from different fluid-saturated rocks to testify the validity of the new model.The results demonstrate that the new model provides more accurate predictions of high dispersion and strong attenuation of different waves in the low frequency range.     

Rayleigh wave phase velocities of South China block and its adjacent areas

Using records of continuous seismic waveforms from 609 broadband seismic stations in the South China Block and its adjacent areas in 2010–2012, empirical Green's functions of surface waves were obtained from cross-correlation functions of ambient noise data between these stations. High quality phase velocity dispersion curves of Rayleigh waves were obtained using time-frequency analysis. These interstation dispersion curves were then inverted to build Rayleigh wave phase velocity maps at periods of 6–50 s. The results of phase velocity maps indicate that phase velocities at 6–10 s periods are correlated with the geological features in the upper crust. Major basins and small-scale grabens and basins display slow velocity anomalies; while most of the orogenic belts and the fold belts display high velocity anomalies. With the gravity gradient zone along Taihang Mountain to Wuling Mountain as the boundary for the phase velocity maps at period of 20–30 s, the western area mainly displays low velocity anomalies, while the eastern side shows high velocity anomalies. Phase velocities in the eastern South China Block south to the Qinling-Dabie orogenic belt is higher than that in the eastern North China Block to the north, which is possibly due to the differences of tectonic mechanisms between the North China Craton and the South China Block. The phase velocities at periods of40–50 s are possibly related to the lateral variations of the velocity structure in the lower crust and upper mantle: The low-velocity anomalies in the eastern part of the Tibetan Plateau are caused by the thick crust; while the Sichuan Basin and the southern part of the Ordos Basin display distinct high-velocity anomalies, reflecting the stable features of the lithosphere in these blocks. The lateral variation pattern of phase velocities in the southern part of the South China Block is not consistent with the surface trace of the block boundary in the eastern Yunnan Province and its vicinities. The phase velocities in the Sichuan Basin are overall slow at short periods and gradually increase with period from the central part to the edge of the basin, indicating the features of shallower basement in the center and overall stable lithospheric mantle of the basin. The middle and upper crust of the southern Ordos Basin in the North China Block is heterogeneous, while in lower crust and the uppermost mantle the phase velocities mainly exhibit high anomalies. High-velocity anomalies are widespread at the middle of the Qinling-Dabie orogenic belt, as well as the areas in southeastern Guangxi with Caledonian granite explosion, but its detailed mechanism is still unclear.     

Constraining the characters of the upper mantle discontinuities beneath the NE margin of the Tibetan Plateau with a dense broadband seismic array

Based on observations from a dense broadband seismic array located along the northeastern (NE) margin of the Tibetan Plateau in southeastern Gansu Province,we use receiver functions (RFs) to pick the arrival times of P-to-S converted waves and bin the traces in different grids according to the piercing points of the 410 and 660 km discontinuities in the upper mantle.The depths of the two discontinuities are estimated by the ray tracing method with the IASP91 velocity model and a 3-D tomography model.The results indicate the following:(1) The arrival times of the P410s and P660s converted phases are delayed by approximately 1 s than those predicted by the IASP91 model.The mantle transition zone (MTZ) is thicker than that in the global model.(2) The synchronous lags in the P410s and P660s arrival times are consistent with low-velocity anomalies in the upper mantle,which are believed to result mainly from the eastward migration of materials beneath the NE margin of the Tibetan Plateau.(3) Combined with previous tomography results,the depression of the'660'discontinuity and the thickened MTZ are somewhat consistent with the big mantle wedge (BMW) model.However,due to data limitations,more studies are required to explore the BMW in the future.     

测震数采仪记录的钻孔应变

A seismometer data acquisition unit has been used in the Changping seismic station to record the output of a strainmeter.The output of a strainmeter was sampled at a rate of 100/sec by seismometer acquisition from the original rate of 1 per minute.Plenty of high frequency sampled data was recorded.The minute value curve calculated from the seismometer acquisition are consistent with that of the original data sampled by the strain acquisition system.More complete waveforms were recorded with a higher sampling rate,and seismic phase parameters calculated by using higher sampling rate strain seismic waves are also in consistency with the results of its predecessors.Spectra of the strain seismic waves are compared with that of seismic waves recorded by a seismometer in the Shisanling seismic station,and their trends are almost the same.Besides,some lower frequency components still exist in strain seismic waves.     

Macroseismic intensity attenuation in Iran

Macroseismic intensity data plays an important role in the process of seismic hazard analysis as well in developing of reliable earthquake loss models. This paper presents a physical-based model to predict macroseismic intensity attenuation based on 560 intensity data obtained in Iran in the time period 1975–2013. The geometric spreading and energy absorption of seismic waves have been considered in the proposed model. The proposed easy to implement relation describes the intensity simply as a function of moment magnitude, source to site distance and focal depth. The prediction capability of the proposed model is assessed by means of residuals analysis. Prediction results have been compared with those of other intensity prediction models for Italy, Turkey, Iran and central Asia. The results indicate the higher attenuation rate for the study area in distances less than 70 km.     

Upper crust structure of eastern A’nyemaqên suture zone: Results of Barkam-Luqu-Gulang deep seismic sounding profile

Barkam-Luqu-Gulang deep seismic sounding profile runs from north of Sichuan Province to south of Gansu Prov- ince. It is located at the northeastern edge of Tibetan Plateau and crosses eastern A’nyemaqên suture zone. The upper crust structures around eastern A’nyemaqên suture zone and its adjacent area are reconstructed based on the arrival times of refracted Pg and Sg waves by using finite difference method, ray tracing inversion, time-term method and travel-time curve analysis. The results show that the depth variation of basement along profile is very strong as indicated by Pg and Sg waves. The basement rose in Zoigê basin and depressed in eastern A’nyemaqên suture zone, and it gradually rose again northward and then depressed. The results also indicate that eastern A’nyemaqên suture zone behaves as inhomogeneous low velocity structures in the upper crust and is inclined to- ward the south. Hoh Sai Hu-Maqên fault, Wudu-Diebu fault and Zhouqu-Liangdang fault are characterized by low velocity distributions with various scales. The distinct variation in basement depth occurred near Hoh Sai Hu-Maqên fault and Zhouqu-Liangdang fault, which are main tectonic boundaries of A’nyemaqên suture zone. Wudu-Diebu fault, located at the depth variation zone of the basement, possibly has the same deep tectonic back- ground with Zhouqu-Liangdang fault. The strongly depressed basement characterized by low velocity distribution and lateral inhomogeneity in A’nyemaqên suture zone implies crushed zone features under pinching action.     

Study of crustal structure with S—wave data from Maqen—Jingbian profile

2-D crustal velocity structure and vp/vs are obtained by processing and interpretation of S-wave data from Maqen-Jingbian deep seismic sounding(DSS)profile.The result shows that there exist obvious differences in 2-D S-wave velocity structure and vp/vs ratio structure along the profile.The S-wave velocities are low and vp/vs ration is high for the westem section of the profile and Haiyuan region,while they are normal for the middle and eastern sections.The changes in lithologic characters of two major anomalous zones are discussed according to lateral variation of S-wave velocity structure and vp/vs ratio structure.It is concluded that the development and occurrence of the Haiyuan strong earthquake is not only related to tectonic activities,but also to lithologic characters of the region.     

Research on the 3—D Seismic Structures in Qinghai—Xizang Plateau

Based on the recording data from the analogue and broadband digital seismic stations in and around Qinghai-Xizang (Tibet)Platean,the three dimensiomal 3-D) seismic velocity stroctures in Qinghai-Xizang Plateau were obtained by using the regional body wave tomography and surface wave tomography.The results from these two tomography methods have similar characteristics for P-and S-wave velocity structures in crust and upper mantle.They show that there are remarkahle low velocity zones in the upper crust of L hasa block in the southern Qinghai-Xizang Plateau and the lower crust and upper mantle of Qiangtang block in the northern Qinghai-Xizang Plateau.These phenomena may be related to the different steps of collision process in southern and northern Qinghai-Xizang Plateau.     

Seismic isolation effect of lined circular tunnels with damping treatments

The Havriliak-Negami model for dynamic viscoelastic material behavior and Biot＇s theory of poroelasticity are employed to develop an exact solution for three-dimensional scattering effect of harmonic plane P-SV waves from a circular cavity lined with a multilayered fluid-filled shell of infinite length containing viscoelastic damping materials and embedded within a fluid-saturated permeable surrounding soil medium. The analytical results are illustrated with numerical examples where the effects of liner/coating structural arrangement, viscoelastic material properties, liner-soil interface bonding condition, seismic excitation frequency, and angle of incidence on the induced dynamic stress concentrations are evaluated and discussed to obtain representative values of the parameters that characterize the system. It is demonstrated that incorporating viscoelastic damping materials with a low shear modulus in the constrained layer configuration is an efficient means of enhancing the overall seismic isolation performance, especially for near-normally incident seismic shear waves where the amplitudes of induced dynamic stresses may be reduced by up to one-third of those without isolation in a relatively wide frequency range. Some additional cases are considered and good agreements with solutions available in the literature are obtained.     

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.     

Research on the Q-value,Site Response and Seismic Source Parameters in Jiangsu and Its Adjacent Areas

Based on 49 digital seismograms recorded by 73 seismic stations in the Jiangsu Telemetered Seismic Network,the paper uses Atkinson's method to calculate the inelastic attenuation coefficient of the Jiangsu area. We find that the frequency-dependent Q in the Jiangsu region is Q( f) = 272. 1·f~(0. 5575). We also use Moya's method to invert the 63 stations' site responses. The results show that the site responses of the 25 stations in Jiangsu are approximately 1 at a range between 1Hz and 20 Hz, which is consistent with their basements on rocks. The response curves of the site responses of the 14 underground stations are similar to each other. Their site responses show an amplification at low frequencies and minimization at high frequencies. The calculation of the Brune model on the waveform data of M_L≥2. 5 earthquakes from Jiangsu Digital Seismic Network between October 2010 and May 2015 in terms of seismic source parameters of 58 seismic waves shows that there are good correlations between seismic magnitude and other source parameters such as seismic moment, source radius and corner frequency, while the correlations between seismic magnitude and stress drop,and stress drop and source radius are not so good.     

Seismogram Analysis and Fitting of Three Earthquakes,West Papua,Recorded at Stations ENH and KMI,China

The S wave velocity structure of the earth below Eastern Southeast Asia has been investigated by analyzing the seismogram from surface wave to multiple depth waves in the time domain and three Cartesian components simultaneously. The wave passes across the front area of subduction zone between the Philippine plate and the Asian plate. The main data are waveform comparisons, instead of the arrival times. The synthetic seismogram is calculated using the GEMINI method. The synthetic seismogram constructed by PREMAN global earth model deviates greatly from the measured one. To solve this problem, corrections are needed for the β speed structure. Corrections cover the gradient change of βh, which turns from negative to positive in upper mantle layers as in the PREMAN, change of earth crust depth and change of zero order coefficients of β velocity function in all earth mantle layers. So, the fitting is obtained, as well as the arrival time or the waveform of Love and Rayleigh surface waves, the S wave and the repetitive depth waves ScS2 and ScS3. This result reveals that the Southeast Asia, being stretched due to tectonic release, has a mantle in some parts with negative anomaly of S wave velocity and vertical anisotropy in all earth mantle layers.     

双相介质中地震波衰减的物理机制

High-frequency seismic attenuation is conventionally attributed to anelastic absorption. In this paper, I present three studies on high-frequency seismic attenuation and propose that the physical mechanism results from the interference of elastic microscopic multiple scattering waves. First, I propose a new theory on wave propagation in a two-phase medium which is based on the concept that the basic unit for wave propagation is a nano- mass point. As a result of the elasticity variations of pore fluid and rock framework, micro multiple scattering waves would emerge at the wavelength of the seismic waves passing through the two-phase medium and their interference and overlap would generate high- frequency seismic attenuation. Second, I present a study of the frequency response of seismic transmitted waves by modeling thin-layers with thicknesses no larger than pore diameters. Results indicate that high-frequency seismic waves attenuate slightly in a near-surface water zone but decay significantly in a near-surface gas zone. Third, I analyze the seismic attenuation characteristics in near-surface water and gas zones using dual-well shots in the Songliao Basin, and demonstrate that the high-frequency seismic waves attenuate slightly in water zones but in gas zones the 160-1600 Hz propagating waves decay significantly. The seismic attenuation characteristics from field observations coincide with the modeling results. Conclusions drawn from these studies theoretically support seismic attenuation recovery.     

Investigation on regional attenuation of Vrancea (Romania) intermediate-depth earthquakes

This paper aims at investigating possible regional attenuation patterns in the case of Vrancea(Romania) intermediate-depth earthquakes.Almost 500 pairs of horizontal components recorded during 13 intermediate-depth Vrancea earthquakes are employed in order to evaluate the regional attenuation patterns.The recordings are grouped according to the azimuth with regard to the Vrancea seismic source and subsequently,Q models are computed for each azimuthal zone assuming similar geometrical spreading.Moreover,the local soil amplification which was disregarded in a previous analysis performed for Vrancea intermediate-depth earthquakes is now clearly evaluated.The results show minor differences between the four regions situated in front of the Carpathian Mountains and considerable differences in attenuation of seismic waves between the forearc and backarc regions(with regard to the Carpathian Mountains).Consequently,an average Q model of the type Q(f) = 115×f~(1.25) is obtained for the four forearc regions,while a separate Q model of the type Q(f) = 70×f~(0.90) is computed for the backarc region.These results highlight the need to evaluate the seismic hazard of Romania by using ground motion models which take into account the different attenuation between the forearc/backarc regions.     

Seismic ground motion relationship in southern China based on stochastic finite-fault model

The characteristics of seismic ground motions in southern China are difficult to determine statistically due to a lack of strong ground motion data. In this study, a stochastic finite-fault ground motion model was adopted to simulate the seismic ground motions at bedrock for southern China, based on parameters derived from small and medium earthquakes that have occurred in the region. From these, the response spectra was estimated. A set of ground motion attenuation relations hipswas then developed based on simulated peak ground motions and response spectral parameters through regression, which would be applicable for use in engineering practice. Through comparisons, it was demonstrated that the proposed ground motion relationships are generally consistent with those obtained from other reported ground motion attenuation models for southern China.     

Investigation of geothermal structure of the Sulawesi,using gravity and magnetic method

The Sulawesi Sea and Sulawesi Island are located in the western Pacific area where volcanic activity,plate subduction,and seismic activity are very active.The Sulawesi basin formed during the Middle Eocene-Late Eocene and nearly half of the Eocene oceanic crust has subducted below the North Sulawesi Trench.The Sulawesi Island was spliced and finalized in the Early Pliocene-Pleistocene during volcanic activity and is recently very active.This area is an optimal location to study volcanic geothermal conditions and subduction initiation mechanisms in the southern part of the western Pacific plate margin,which are important in geothermal and geodynamic research.In this study,we combined 133 heat flow data with gravity and magnetic data to calculate the Moho structure and Curie point depth of the Sulawesi Sea and periphery of the Sulawesi Island,and analyze the distribution characteristics of the geothermal gradient and thermal conductivity.The results show that the average depths of the Moho and Curie surfaces in this area are 18.4 and 14.3 km,respectively,which is consistent with the crustal velocity layer structure in the Sulawesi Basin previously determined by seismic refraction.The average geothermal gradient is 4.96°C(100 m)-1.The oceanic area shows a high geothermal gradient and low thermal conductivity,whereas the land area shows a low geothermal gradient and high thermal conductivity,both of which are consistent with statistical results of the geothermal gradient at the measured heat flow points.The highest geothermal gradient zone occurs in the transition zone from the Sulawesi Sea to Sulawesi Island,corresponding to the spreading ridge of the southward-moving Sulawesi Basin.Comprehensive gravity,magnetic,and geothermal studies have shown a high crustal geothermal gradient in the study area,which is conducive to the subduction initiation.The northern part of the Palu-koro fault on the western side of Sulawesi is likely the location where subduction initiation is occurring.During the process of moving northwest,the northern and eastern branches of Sulawesi Island have different speeds;the former is slow and the latter is fast.These branches also show different deep tectonic dynamic directions;the northern branch tilts north-up and the eastern branch tilts north-down.