中国台湾地区-菲律宾群岛深部纵横波结构成像及其构造意义

使用地震层析成像方法反演了中国台湾地区-菲律宾群岛的1197126条P波震相和1217821条S波震相,首次同时得到了该地区从地表到100km的纵横波速度结构。成像结果揭示了沿着马尼拉海沟向东俯冲的欧亚板块从中国台湾西南部到吕宋岛南端构造形态上的变化;中国台湾中央山脉由于受到造山运动的影响,地壳厚度可达55~60km,而山脉两侧的地壳厚度多为20~35km;此外,成像结果还提供了菲律宾海板块在中国台湾东北部及琉球海沟下方北向俯冲的地球物理学证据;研究区域南部的菲律宾群岛由于同时受到两侧欧亚板块和菲律宾海板块的双向俯冲,使得岛内的岩石圈变形严重,岛弧岩浆和地震活动十分发育;群岛东侧的菲律宾海板块在东吕宋海槽的活动性较弱,很可能是本哈姆海台的碰撞和俯冲所造成的,而菲律宾海板块在菲律宾海沟的俯冲活动则十分强烈,随着深度的增加逐渐向南发展。研究表明,板块俯冲造成了中国台湾地区-菲律宾群岛的地壳及上地幔具有较强的不均一性,这不仅孕育了大量的地震和火山活动,同时也对该地区的地质构造产生了深远的影响。     

Vp/Vs ratios in the Yellowstone National Park region, Wyoming

In this paper we study the variation of V_p/V_s and Poisson's ratio (δ) in the Yellowstone National Park region, using earthquakes which were well recorded by a local seismic network. We find that the average V_p/V_s value within the geothermally active Yellowstone caldera is about 7% lower than in the area outside the caldera. Within the caldera itself there may be a further 2–7% reduction of V_p/V_s in the hydrothermally active Norris Geyser Basin, the Upper and Lower Geyser Basins, and the Yellowstone Lake and Mud Volcano regions. After considering various possible causes for V_p/V_s changes, such as geologic and structural differences, thermal effects, partial melting, and hydrothermal activity, we conclude that the most plausible explanation for the observed V_p/V_s reduction is the presence of hot-water at temperatures and pore-pressures near the water steam transition in the caldera geothermal reservoirs.     

Use of mode-converted waves in marine seismic data to investigate the lithology of the sub-bottom sediments in Isfjorden,Svalbard

The compression wavefield is efficiently converted to shear-wave energy at post-critical angles in areas of high impedance contrast at the sea floor. We have analysed mode-converted shear waves in a data set acquired with a hybrid marine/land geometry in Isfjorden, Svalbard. Through a kinematic 2D ray-tracing modellingV _p/V_s ratios for part of the uppermost 5km of the crust are obtained. Low values (V _p /V _s =1.65) are tentatively associated with the section of Devonian sandstones which appears to attain a minimum thickness of 1.5km below 3 km depth about 10km west of Kapp Thorden.     

Seismic anisotropy of upper mantle in eastern China

Based on the polarization analysis of teleseismic SKS waveform data recorded at 65 seismic stations which respectively involved in the permanent and temporary broadband seismograph networks deployed in eastern China, the SKS fast-wave direction and the delay time between the fast and slow shear waves at each station were determined by use of SC method and the stacking analysis method, and then the image of upper mantle anisotropy in eastern China was acquired. In the study region, from south to north, the fast-wave polarization directions are basically EW in South China, gradually clockwise rotate to NWW-SEE in North China, then to NW-SE in Northeast China. The delay time falls into the interval [0.41 s, 1.52 s]. Anisotropic characteristics in eastern China indicate that the upper mantle anisotropy is possibly caused by both the collision between the Indian and Eurasian Plates and the subduction from the Pacific and Philippine Sea Plates to the Eurasian Plate. The collision between two plates made the crust of western China thickening and uplifting and the material eastwards extruding, and then caused the upper mantle flow eastwards and southeastwards. The subduction of Pacific Plate and Philippine Sea Plate has resulted in the lithosphere and the asthenosphere deformation in eastern China, and made the alignment of upper mantle peridotite lattice parallel to the deformation direction. The fast-wave polarization direction is consistent with the direction of lithosphere extension and the GPS velocity direction, implying that the crust-upper mantle deformation is possibly a vertically coherent deformation. Supported by Special Project for the Fundamental R & D of Institute of Geophysics, China Earthquake Administration (Grant No. DQJB06B06), Special Program of the Ministry of Science and Technology of China (Grant No. 2006FY110100), China Digital Earthquake Observation Network Project “North China Seismic Array”, and National Natural Science Foundation of China (Grant Nos. 40334041 and 40774037)     

菲律宾群岛深部速度结构成像与双向俯冲板片构造特征

菲律宾群岛受到欧亚板块、菲律宾海板块和印度-澳大利亚板块的碰撞作用,地质环境复杂,构造因素多样.尽管近几年来已经有了少数关于该区域层析成像的研究,但这些研究的区域主要集中在马尼拉海沟、吕宋岛及中菲律宾地区,而关于群岛周围其他海沟和南菲律宾地区的讨论相对较少.到目前为止,还没有同时获得过关于菲律宾群岛深部纵、横波速度结构的研究,本次研究通过反演155779条P波震相和59642条S波震相,同时获得了菲律宾群岛从地表至150 km深度的纵、横波速度结构.地震层析成像结果表明该地区的壳幔速度结构具有较强的不均一性,地壳内部存在着广泛的低速异常,而表征俯冲板块的高速异常则沿着群岛周边的海沟展布.南海块体在马尼拉海沟中段的俯冲角度和俯冲活动性比南段小;菲律宾海板块在东吕宋海槽南段微弱的俯冲作用很有可能同本哈姆海台的碰撞有关.菲律宾群岛大部分M_W>6.0的强震沿着各个板块的边界发生,体现出菲律宾海板块同欧亚板块之间的强耦合作用,群岛西侧的南海块体在马尼拉海沟16°N-20°N之间呈现出的弱耦合状态可能跟北吕宋地区的拉张应力环境有关,南海块体在16°N以南的地区同上覆块体之间的耦合作用较强;群岛东侧的菲律宾海板块在14°N以北的地区没有强震发生,它与菲律宾群岛之间的耦合程度从北向南逐渐增强,在12°N以南的地区要强于12°N以北的地区;此外苏禄海盆和菲律宾构造带之间也存在着强耦合关系.     

DEEP VELOCITY STRUCTURE AND ITS TECTONIC IMPLICATIONS IN SOUTHERN PHILIPPINES

3-D V_P and V_S images of southern Philippines at the 0~100km depths are generated by inverting a large number of travel-time data from the International Seismological Centre(1960-2017)through seismic tomography method. The results show lateral variation exists in the crust and upper mantle:High V_P and V_S anomalies emerge in mid-west Mindanao and Bohol Island, which might be caused by the combined action of huge magmatism and ophiolite accretion in the lower crust; low velocity anomalies of the upper mantle in the west of Mindanao are consistent with locations of volcanoes on the surface. It, thus, could be inferred that the low velocity anomaly is closely related to magmatic activity. The dense earthquake distribution along plate margin extending to 100km coincides with the strong activity of the Philippine Sea Plate which is located in the northeast and southeast of Mindanao. Relative weak activity of Sulawesi Sea Basin is presented simultaneously. The subduction of the Philippine Sea Plate is mostly concentrated in the crust and the top of the uppermost mantle. Our tomographic images show that lateral heterogeneities exist in the crust and uppermost mantle of the southern Philippines. Low V_P and V_S anomalies emerge in Philippine Trench and Cotabato Trench, in contrast, high V_P and V_S anomalies appear in shallow crust of land area where a large number of earthquakes and magmatic activities develop. This may reflect strong tectonic processes between the Philippine Sea Plate and Philippine Mobile Belt. Low V_P and V_S anomalies in the crust of eastern Mindanao coinciding with the location of volcanoes on the surface may show partial melting of crust material caused by dehydration of the subducting Philippine Sea Plate. Such a similar phenomenon can be also seen in the south of Negros Island and Cotabato Trench. Thus we infer that active tectonic behaviors are constrained within the crust of the Philippine Sea Plate, Sulu Sea Basin and Sulawesi Basin.Low V_P and V_S anomalies of the mantle in the mid-west of Mindanao island are associated with magmatic activity which may be caused by a collision between the east and west part of Mindanao at 5Ma. The fracture system in the west of Mindanao provides the possible passage ways of mantle hot material upwelling, coinciding with the model of geothermal distribution in this area. According to the geochemical analysis, ophiolite observed in Sanbaoyan and the western part of Mindanao could indicate material composition from crust to upper mantle on Eurasian continental margin which may show the evidence of rapid expansion environment of mid-ocean ridge. High V_P and V_S anomalies in the mantle of northeast and southeast of Mindanao coinciding with the distribution of massive earthquake along boundaries show a well agreement with the shape of the Philippine Sea Plate. Dense earthquake distribution in south Mindanao at 100km shows the Philippine Sea Plate has strong activity and stress accumulation in the upper mantle. On the contrary, the seismicity in southwest Mindanao and Cotabato Trench reduces rapidly at the depth from 50km to 100km, revealing weak subduciton and stress release of Sulawesi Basin in the mantle.     

Receiver Functions for Three-layer Media

We have extended the H-k stacking method of receiver functions applicable to a three-layer model, which is useful in studying detailed crustal structures. We have demonstrated its application with two sample sites in Taiwan, making use of travel times of converted phases from the direct P waves of teleseismic events as the P reflects and refracts at different discontinuities in the crust. This three-layer extension allows a closer examination on the crust, as well as the relationship between crustal layers and the associated V_p/V_s ratios. Data were processed using the multiple-taper correlation technique to obtain the radial receiver functions (RRFs). The relative time delays of the converted phases measured from the RRFs were used to estimate the depths of the crustal discontinuities. Results not only yield the depths to principal crustal layers, including the Moho, but also give the corresponding ranges of V_p/V_s ratio which are comparable with findings from other tomographic studies.     

Deep tectonic setting of the 2008 Wenchuan <Emphasis Type="Italic">M</Emphasis><Subscript>s</Subscript>8.0 earthquake in southwestern China

Teleseismic P-wave receiver functions at 20 broadband seismic stations in the Longmenshan fault zone (LMFZ) and its vicinity were extracted, and the crustal thickness and the P- and S-wave velocity ratio were calculated by use of the H-k stacking algorithm. With the results as constraints, the S-wave velocity structures beneath each station were determined by the inversion of receiver functions. The crustal structure of the Rear-range zone is similar to that of the Songpan-Garze Block, whereas the velocity structure of the Fore-range zone resembles that of Sichuan Basin, implying that the Central Principal Fault of LMFZ is the boundary between the eastern Tibetan Plateau and the Yangtze Block. Lower velocity zone exists in lower crust of the Songpan-Garze Block and the central-southern segment of the Rear-range zone, which facilitates the detachment of the material in upper and middle crust. Joint analysis of the receiver functions and the Bouguer gravity anomalies supports the thesis on the detachment-thrust mode of the LMFZ. A double-detachment pattern is suggested to the tectonic setting in the Songpan-Garze Block. The upper detachment occurs at the depth of 10-15 km, and represents a high-temperature ductile shear zone. There is a lower detachment at the depth of about 30 km, below which the lower crust flow exists in the eastern Tibetan Plateau. Interpretation of the Bouguer gravity anomalies indicates that the Sichuan Basin is of higher density in upper and middle crust in comparison with that of the Songpan-Garze Block. The LMFZ with higher density is the result from the thrusting of the Songpan-Garze Block over the Sichuan Basin. In the lower crust, higher P velocity and higher density in the Sichuan Basin are related to more rigid material, while lower S velocity and lower density in the Songpan-Garze Block are related to the softened and weakened material. The higher density block beneath the Sichuan Basin obstructs the eastward flow of lower crustal material from the Tibetan Plateau, which is driven by the compression of northward movement of Indian Plate. The eastward movement of upper and middle crustal material is also obstructed by the rigid Yangtze Block, resulting in the stress concentrated and accumulated along the LMFZ. When the stress releases sharply, the Wenchuan M _s8.0 earthquake occurs. Supported by the National Natural Science Foundation of China (Grant Nos. 40334041, 40774037) and Joint Foundation of Earthquake Science (Grant No. 1040062)     

Crust and uppermost mantle structure of the Ailaoshan-Red River fault from receiver function analysis

S-wave velocity structure beneath the Ailaoshan-Red River fault was obtained from receiver functions by using teleseismic body wave records of broadband digital seismic stations. The average crustal thickness, Vp/Vs ratio and Poisson’s ratio were also estimated. The results indicate that the interface of crust and mantle beneath the Ailaoshan-Red River fault is not a sharp velocity discontinuity but a characteristic transition zone. The velocity increases relatively fast at the depth of Moho and then increases slowly in the uppermost mantle. The average crustal thickness across the fault is 36―37 km on the southwest side and 40―42 km on the northeast side, indicating that the fault cuts the crust. The relatively high Poisson’s ratio (0.26―0.28) of the crust implies a high content of mafic materials in the lower crust. Moreover, the lower crust with low velocity could be an ideal position for decoupling between the crust and upper mantle.     

Lg coda Q0 value and its relation with the tectonics in Chinese mainland and adjacent regions

We have collected 432 vertical component records from 45 stations of new CENC （China Earthquake Network Center） in Chinese mainland and adjacent regions. These records were used to calculate Q0 （Q at 1Hz） and η values of Lg coda from each station by the stack spectral ratio （SSR） method. Then the tomography method was applied to obtaining lateral variation of Q0 and η values in Chinese mainland and adjacent regions. The result indicates that Q0 value varies between 150 and 600 in the studied areas. Yunnan, southwest Sichuan, and northwest Myanmar show the lowest Q0 value （Q0〈240） and the crust of these regions is characterized by complicated crack and strong hydrothermal activity. The highest Q0 value （Qo〉510） exists in the border of southern Mongolia, Alxa and Ordos block. The η value varies between 0.45 and 0.75 in Chinese mainland and its adjacent regions.     

Fault deformation anomaly and intermediate and short-term prediction of the Jingtai M s=5.9 earthquake

The time-space distribution characteristics of fault deformation anomaly in the near-source region and its outlying zone in the seismogenic process of the Jingtai M _s=5.9 earthquake occurred on June 6, 2000 in Gansu Province is studied preliminarily. The distribution scope of fault deformation anomaly before the earthquake is wide, the anomaly shape is complicated and the pattern anomalous zone of fault deformation (strain) information index is obvious. The shape and amplitude of fault deformation anomaly in different regions differ significantly, which is closely related with the tectonic location of anomaly. The fault deformation anomaly of α, β, and γ phases along the western segment of Haiyuan fault zone shows the process from the quasi-linearity to non-linearity of fault movement in the near-source region, matches the high-value anomalous area of fault deformation (strain) information index, and reflects the high strain accumulation in the seismogenic region. However, the anomaly of abrupt jump and cusp with a large amplitude occurred in the areas far from the earthquake, such as Liupanshan fault zone which is the tectonic convergent section does not reflect the strain accumulation of its location, maybe it is a sign that the regional tectonic stress field is strengthened in the seismogenic process. Based on the above-mentioned facts and combined with the preliminary summary of experiences and lessons in the intermediate and short-term prediction of the Jingtai M _s=5.9 earthquake, we study and explore the application of fault deformation anomaly to earthquake judgment. Foundation item: National Key Basic Research Development Program (G1998040703 and G1998040705), and State Scientific and Technological Project of the “Ninth Five-Year Plan” (96-913-09-01-02-03 and 96-913-09-02-02-03), China.     

Time domain solution for multiple scattering and the coda envelopes

This article summarizes work on multiple scattering based on models of media with randomly distributed scatterers. The scatterers are isotropic and statistically uniform. Measuring distance in terms of mean-free pathL _s and time in terms of the mean-free timesL _s/V, whereV is the velocity of scattered waves, we have more convenient dimensionless distance and time. It can be shown that after the dimensionless time equals 0.65 energy contributed from multiple scattering becomes predominant. Thus the later coda reflects the effect of multiple scattering rather than single scattering. Treating the seismic record, including starting and tail parts, as a whole, the diffusion theory predicts that at a dense distribution of scatterers and a small distance between source and receiver, codas reflect mainly intrinsicQ _i. Of course, this conclusion is coincident with the presumption of the diffusion theory,Q _s>Q _i. However, from a new integral equation of multiple scattering, which deals with the scattered waves and primary waves separately, the conclusion is similar but clearer. This article quotes the new expression for coda energy in two-dimensional space. It shows that if the receiver is close to the source, the coda decay reflects only intrinsicQ _i, then as the distance increases, effects of scatteringQ _s, are involved in the decay feature. The theoretical plots of coda decay show that it seems in most cases in the earthQ _i should not be smaller than one tenth ofQ _s.Project Sponsored by the Joint Earthquake Science Foundation of China.     

Focal depths for moderate-sized aftershocks of the Wenchuan <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>8.0 earthquake and their implications

Sliding-window cross-correlation method is firstly adopted to identify sPn phase, and to constrain focal depth from regional seismograms, by measuring the time separation between sPn and Pn phases. We present the focal depths of the 17 moderate-sized aftershocks (M _S⩾5.0) of the Wenchuan M _S8.0 earthquake, using the data recorded by the regional seismic broadband networks of Shaanxi, Qinghai, Gansu, Yunnan and Sichuan. Our results show focal depths of aftershocks range from 8 to 20 km, and tend to cluster at two average depths, separate at 32.5°N, i.e., 11 km to the south and 17 km to the north, indicating that these aftershocks are origin of upper-to-middle crust. Combined with other results, we suggest that the Longmenshan fault is not a through-going crustal fault and the Pingwu-Qingchuan fault may be not the northward extension of the Longmenshan thrust fault. Supported by the National Natural Science Foundation of China (Grant Nos. 40604009 and 40574040) and Special Project for the Fundamental R & D of Institute of Geophysics, China Earthquake Administration (Grant No.DQJB08B20)     

四川、甘肃地区VS30经验估计研究

目前我国建筑工程抗震设计规范中对于工程场地条件的判断依据主要是地表以下20m深度范围内土层的等效剪切波速,简称VS20。相比之下,国外应用较广的是地表以下30m深度范围内的等效剪切波速,简称VS30。这种差别导致国内科研工作者在应用国外的地震工程、工程抗震模型时经常遇到对场地条件描述不准确的困难。为了解决这个问题,本文根据147个四川、甘肃地区国家强震动台站20m左右深度的钻孔剪切波速数据,利用延拓方法、场地分类统计方法以及基于地形特征的VS30估计方法研究各台站VS30与VS20的经验关系,对比发现基于速度梯度延拓的结果最为可取。参考国际上通用的Geomatrix Classification场地分类标准,最终得到四川、甘肃地区各类场地的平均VS30,此结果可以为缺乏钻孔数据的工程场地的VS30估计提供参考。     

Application of Fisher method to discriminating earthquakes and explosions using criterion m b/ M S

We try to give a quantitative and global discrimination function by studying m _b/M _S data using Fisher method that is a kind of pattern recognition methods. The reliability of the function is also analyzed. The results show that this criterion works well and has a global feature, which can be used as first-level filtering criterions in event identification. The quantitative and linear discrimination function makes it possible to identify events automatically and achieve the goal to react the events quickly and effectively. Contribution No.05FE3018, Institute of Geophysics, China Earthquake Administrtion.     

四川芦山7.0级地震前中小地震P轴方位角CV值的变化

通过对2003年1月1日—2013年4月1日芦山地震前震源区中小地震震源机制解的分析,发现不同阶段的震源机制解在一定程度上反映了强震孕育过程中构造应力场随时间的变化。震源区中小地震的P轴方位角C_V值在芦山M7.0地震发生前有一个上升-下降-上升的过程,只是相比于汶川8.0级地震前C_V值的下降-上升过程经历了更长的时间,这表明四川芦山M7.0地震的孕育经历了长时间的应力积累,与许多研究结果相一致。2007年1月1日—2014年4月1日C_V值空间分布的非均匀性特征在龙门山断裂带南段有显著的增强与减弱过程,对于发震地点可能有一定的指示意义。     

Tomographic imaging of seismic velocity structure in and around the Onikobe volcanic area, northeastern Japan: implications for fluid distribution

We estimated three-dimensional P- (V_p) and S-wave velocity (V_s) and V_p/V_s structures in and around the Onikobe volcanic area, northeastern Japan, by local travel time tomography. We used travel time data from source and receiver pairs located within and outside the study area, which plays an important role in obtaining the optimum ray coverage and in elucidating the deeper structure more accurately. Detailed information on deeper structures is essential for imaging the complete volcanic system from the magmatic source zone through areas of shallow hydrothermal circulation. More than 50 000 travel time data for the P-waves and 35 000 for the S-waves were used to image the velocity structure. Our results show the following dominant features: (1) two conduits in the upper crust with low V_p and low V_s indicative of H₂O-rich fluid pathways: one lying beneath Naruko volcano, the other beneath the focal area of the 1962 Northern Miyagi earthquake (M6.5); (2) an underlying broad region in the lower crust with low V_p, low V_s and high V_p/V_s, suggestive of a zone of partial melt, from which the fluids in (1) are derived; and (3) low V_p/V_s areas near the surface of the Sanzugawa and Onikobe calderas, suggesting a diffuse vapor-saturated cap.     

2022年门源MS6.9和2016年门源MS6.4地震序列比较分析

2022年1月8日青海省海北州门源县发生M_S6.9地震,震中距离2016年1月21日门源M_S6.4地震震中约33km,两次门源地震均发生在冷龙岭断裂附近,但在震源机制、主发震断层破裂过程及地震序列余震活动等方面显著不同。针对两次门源地震序列的比较分析,对研究冷龙岭断裂及其附近区域强震序列和余震衰减特征等具有重要研究意义。通过对比分析2022年门源M_S6.9地震和2016年门源M_S6.4地震余震的时空演化特征,发现二者在震源过程和断层破裂尺度上存在明显差异,前者发震断层破裂充分,震后能量释放充分,余震丰富且震级偏高;而后者发震断层未破裂至地表,余震震级水平偏低。综合分析两次门源地震序列表现出来的差异性,认为其可能与地震发震断层的破裂过程密切相关,且同时受到区域构造环境的影响。     

Secular variations ing 1 0 component of geomagnetic field and its origin

The secular variations (SV) ing ₁ ⁰ of geomagnetic field, caused by the interaction between the geomagnetic field and the radial fluid-flow of the outer-core, are discussed with both statistical and analytic methods. When the value ofg ₁ ⁰ s ₂ ⁰ is relatively high compared with other terms, SV ing ₁ ⁰ is characterized by exponential change. When the effect of non-dipole field is notable, SV ing ₁ ⁰ shows complex features. Especially when the value of g ₁ ⁰ is close to zero, SV ing ₁ ⁰ is completely determined by the random process of non-dipole fields. Project supported by the National Natural Science Foundation of China (Grant No. 49429405) and by the Chinese Academy of Sciences.