Upper mantle seismic anisotropy beneath a convergent boundary: SKS waveform modeling in central Tibet

This study used SKS waveforms from the International Deep Profiling of Tibet and the Himalayas (INDEPTH) III dataset and a new 2D method for modeling seismic waves in anisotropic media to construct an image of anisotropic structures beneath central Tibet. A preferred model revealed three-segment anisotropic structures in the upper mantle beneath the study region. Waveform modeling demonstrated that the anisotropy was mainly generated by the lithosphere but not the asthenosphere, and that an anisotropic model with a flatter axis of symmetry provides a more consistent interpretation of the observations than models having steeply dipping symmetry axes. A relatively low velocity zone may underlie or intermingle with the anisotropic structures in the northern portion of the region. Synthetic tests also indicate that variations in the elastic constants and depth extent of the anisotropy assumed by the calculations do not affect the general conclusions, although trade-offs exist among certain model parameters. The modeling results suggest that the complex seismic structures in central Tibet were associated with underthrusting of the Indian lithosphere beneath the Asian lithosphere; the inferred flat symmetry axis of the anisotropy was likely generated during this collision process. If this were not the case, the inherited anisotropy would exhibit a steeply dipping axis of symmetry, parallel to the direction of underthrusting.     

青藏高原东北缘壳幔各向异性研究:基于SKS和Pms震相分析

利用甘肃和青海两省固定宽频带地震台记录的远震波形资料,挑选高质量SKS震相,联合使用最小切向能量方法和旋转互相关方法获得230对高信噪比分裂参数;同时对接收函数中Pms震相随方位角的变化进行拟合,得到了24个台站的地壳各向异性分裂参数.整个区域SKS分裂快波方向均值为123°,Pms分裂快波方向均值为132°,且大部分区域SKS、Pms快波方向与地表构造走向相一致,说明青藏高原东北缘以岩石圈垂直连贯变形为主,地壳上地幔相互耦合.SKS、Pms分裂时差均值分别为1.0s和0.6s,显示地壳各向异性对于SKS分裂时差有较大贡献.昆仑断裂附近Pms、SKS分裂快波方向与昆仑断裂走向基本一致,说明昆仑断裂可能是岩石圈尺度深大断裂;而阿尔金断裂东缘二者快波方向显著差异意味着阿尔金断裂在东缘可能仅为地壳尺度的断裂.     

Backazimuthal Variations of Splitting Parameters of Teleseismic SKS Phases Observed at the Broadband Stations in Germany

—SKS phases observed at broadband stations in Germany show significant shear-wave splitting. We have analyzed SKS and SKKS phases for shear-wave splitting from 13 stations of the German Regional Seismic Network (GRSN), from 3 three-component stations of the Gräfenberg array (GRF) and from one Austrian station (SQTA). The data reveal strong differences in the splitting parameters (fast direction φ and delay time δt from a single event at various stations as well as variations at the individual stations for events with different backazimuths. The backazimuthal variations of the splitting parameters at some stations can be explained by two-layer anisotropy models with horizontal symmetry axes. The best resolved two-layer model is the GRA1 model (upper layer φ = 40°, δt = 1.15s; lower layer φ = 115°, δt = 1.95s). The upper layer can be attributed to the lithosphere. Because of the magnitude of the delay time of the upper layer, the lower layer must lie within the asthenosphere. At other stations splitting parameters are consistent with an anisotropic one-layer model for the upper mantle. Stations near the Bohemian Massif show fast directions near EW. Throughout NE Germany the directions are oriented NW/SE. The reason for this direction is probably the nearby Tornquist-Teisseyre line. The observed fast axes are subparallel to this prominent Transeuropean suture zone. At stations in southern Germany near the Alps we observed ENE/WSW directions. Below some stations we also found indications of inclined anisotropic layers.     

海拉尔CTBTO地震台阵下方小尺度非均匀体研究

利用国际禁核组织在海拉尔布设的小孔径地震台阵的16个深远地震波资料,分离出了P波波场中不相干的尾波.频率域中对波场扰动的研究表明,在0.5～3.0 Hz的频率范围内,观测到的波场扰动可以利用P波在随机介质中的散射进行解释.不同远震事件得到的稳定结果表明,在该台阵下方34 km厚的地壳和110 km厚的岩石层中可能存在相关长度为2.0～7.4 km的散射体.     

Azimuthal anisotropy in lithosphere on the Chinese mainland from observations of SKS at CDSN

Ａｚｉｍｕｔｈａｌ　ａｎｉｓｏｔｒｏｐｙ　ｉｎ　ｌｉｔｈｏｓｐｈｅｒｅ　ｏｎ　ｔｈｅ　Ｃｈｉｎｅｓｅ　ｍａｉｎｌａｎｄ　ｆｒｏｍ　ｏｂｓｅｒｖａｔｉｏｎｓ　ｏｆ　ＳＫＳ　ａｔ　ＣＤＳＮ（郑斯华）（高原）Ａｚｉｍｕｔｈａｌａｎｉｓｏｔｒｏｐｙｉｎｌｉｔ...     

基于PpSms、P660S、SKS震相分裂研究上海地区地壳、地幔各向异性

基于上海地震台阵记录的远震利用台阵处理技术研究了PpSmS、P660S及SKS震相的分裂。根据三震相的分裂参数(快慢波延时时间及快波偏振方向)的差异对上海地震台阵地区地壳、上地幔及下地幔介质的各向异性情况进行了简要分析。     

中天山及邻区S波分裂研究及其动力学意义

本文利用天山及其邻区布设的37个宽频带地震台站记录到的远震波形数据,分别采用最小能量法和旋转相关法对SKS和SKKS波震相进行了偏振分析,计算出了台站下方介质的S波分裂参数：快波的偏振方向（φ）和慢波延迟时间（δt）.本文研究结果表明：中天山内部大多数台站的各向异性快波方向呈NEE-SWW向,与天山走向平行,慢波时间延迟为0.4～1.7 s,这是塔里木、哈萨克斯坦的南北双向俯冲及其导致的天山地区岩石圈地幔南北向缩短变形的直接反映.本文研究发现中天山北部部分台站下方地震各向异性快波方向与慢波延时随方位角呈现规律性的变化,可能暗示该区上地幔各向异性不能仅用单层水平各向异性这一简单模式来解释.75°E以西的天山地区台站下方S波快波方向和延时具有强烈的横向变化,可能与研究区下方存在的小规模对流有关.中天山不同地段地震台站下方各向异性明显不同,进一步证实了天山地区构造变形的复杂性.     

Study on the Anisotropy of the Underground Medium in Hohhot

In this paper, the data of earthquake events of magnitude M_S6.0 and above produced in Hohhot Seismic Station from 2008 to 2015 and the data of M_L ≥ 1.0 seismic events from 2015 to 2016 in Horinger Seismic Station and the surrounding mobile stations in southern Hohhot are selected. Using Splitlab and SAM software, the spin-correlation method, the least-energy method and the cross-correlation coefficient method are used to analyze the teleseismic and near-seismic phases (SKS, S). The results of this study are in good agreement with the results previously obtained by other researchers. The study of teleseismic SKS splitting reflects the characteristics of the anisotropy of the upper mantle beneath Hohhot, that is, the anisotropy of the upper mantle shows NW, which reflects "fossil" Anisotropy, mainly in the continental structure of stable units and preserves the history of mantle deformation information. The crustal anisotropy reflected by the near-earthquake S-wave splitting study is similar to that of the active fault zone, trending NE as a whole and is consistent with the tectonic stress field background of the northeastern margin of Ordos block.     

D″ anisotropy inverted from shear wave splitting intensity

The D″ layer, located at the bottom of the mantle, is an active thermochemical boundary layer. The upwelling of mantle plumes, as well as possible plate subduction in the D″ layer, could lead to large-scale material transformation and mineral deformation, which could result in significant seismic anisotropy. However, owing to limited observations and immense computational cost, the anisotropic structures and geodynamic mechanisms in the D″ layer remain poorly understood. In this study, we proposed a new inversion method for the seismic anisotropy in the D″ layer quantitatively with shear wave splitting intensities. We first proved the linearity of the splitting intensities under the ray-theory assumption. The synthetic tests showed that, with horizontal axes of symmetry and ray incidences lower than 30º in the D″ layer (typical SKS phase), the anisotropy is well resolved. We applied the method to the measured dataset in Africa and Western Europe, and obtained strong D″ anisotropy in the margins of the large low shear-wave velocity provinces and subducting slabs. The new method makes it possible to obtain D″ anisotropy, which provides essential constraints on the geodynamical processes at the base of the mantle.     

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)     

Seismic anisotropy beneath the Chinese Mainland: Constraints from shear wave splitting analyses

We use earthquakes recorded by the China National Seismic Network from 2015 to 2019 and measure shear wave splitting parameters of SKS to study the anisotropic characteristics beneath the mainland of China. In general, the fast directions change from nearly E-W in western China (northwest China and Qinghai-Tibetan Plateau) to nearly N–S in central China (Ordos and Sichuan-Yunnan), and then turn to approximately E-W in eastern China (North and South China). The delay times of slow wave in eastern China are about 1.0–1.7 ​s, larger than those in central and western China (about 0.6–1.0 ​s). In addition, the fast directions in eastern China are highly consistent with the plate motion direction and horizontal GPS velocities with respect to Eurasia, indicating that the observed anisotropy is mainly from the asthenosphere which is strongly coupled to the overlying lithosphere. However, the fast directions in western China are mostly in accord with the strike of the surface structures (such as faults), possibly due to the directional arrangement of crystal lattices caused by shear deformation under tectonic activities.     

阿尔金断裂西部邻区的上地幔各向异性研究

本文利用横波分裂方法对北京大学于田流动台阵记录的SKS震相进行分析,获到了阿尔金断裂西部及邻区的上地幔各向异性参数.分析结果显示,快波偏振方向在整个研究区基本呈近E-W向,与研究区内阿尔金断裂的走向几乎一致,分裂延迟时间在0.93~1.20s之间.综合研究区附近前人横波分裂研究结果,我们认为,在印度和欧亚大陆板块碰撞作用下,青藏高原北部上地幔软流圈物质向北流动,遇到塔里木盆地"克拉通"较厚岩石圈阻挡并发生了旋转,向东西两侧流动,导致在青藏高原和塔里木盆地边界地带软流圈上地幔橄榄岩中晶格沿近E-W向优势排列.这一模式显示阿尔金断裂可能是一个岩石圈尺度的大型走滑断裂:它既控制近地表的上地壳构造运动,同时也影响了上地幔软流圈物质的流动.另外,在向塔里木盆地内部延伸的台站也观测到显著的各向异性和近EW向的快波偏振方向.这些结果表明塔里木盆地"克拉通"岩石圈的中、下部分在南部边界被青藏高原北部上地幔软流圈流动"热侵蚀"而损失一部分,导致青藏高原软流圈向东西两侧的流动已经延伸到塔里木盆地内部.本文的研究结果揭示克拉通岩石圈"活化"不仅可以在垂直方向发生(如,岩石圈拆沉或软流圈上涌导致的热侵蚀),也可以在水平方向上发生,即软流圈的水平流动对克拉通岩石圈边界的热侵蚀作用.     

Passive Seismology and Deep Structure in Central Italy

—In the last decade temporary teleseismic transects have become a powerful tool for investigating the crustal and upper mantle structure. In order to gain a clearer picture of the lithosphere-asthenosphere structure in peninsular Italy, between 1994 and 1996, we have deployed three teleseismic transects in northern, central, and southern Apennines, in the framework of the project GeoModAp (European Community contract EV5V-CT94–0464). Some hundreds of teleseisms were recorded at each deployment which lasted between 3 and 4 months. Although many analyses are still in progress, the availability of this high quality data allowed us to refine tomographic images of the lithosphere-asthenosphere structure with an improved resolution in the northern and central Apennines, and to study the deformation of the upper mantle looking at seismic anisotropy through shear-wave splitting analysis. Also, a study of the depth and geometry of the Moho through the receiver function technique is in progress. Tomographic results from the northernmost 1994 and the central 1995 teleseismic experiments confirm that a high-velocity anomaly (HVA) does exist in the upper 200–250 km and is confined to the northern Apenninic arc. This HVA, already interpreted as a fragment of subducted lithosphere is better defined by the new temporary data, compared to previous works, based only on data from permanent stations. No clear high-velocity anomalies are detected in the upper 250 km below the central Apennines, suggesting either a slab window due to a detachment below southern peninsular Italy, or a thinner, perhaps continental slab of Adriatic lithosphere not detectable by standard tomography. We found clear evidence of seismic anisotropy in the uppermost mantle, related to the main tectonic processes which affected the studied regions, either NE–SW compressional deformation of the lithosphere beneath the mountain belt, or arc-parallel asthenospheric flow (both giving NW–SE fast polarization direction), and successive extensional deformation (～E–W trending) in the back-arc basin of northern Tyrrhenian and Tuscany. Preliminary results of receiver function studies in the northern Apennines show that the Moho depth is well defined in the Tyrrhenian and Adriatic regions while its geometry underneath the mountain belt is not yet well constrained, due to the observed high complexity.     

Linking the Alps and Apennines subduction systems: New constraints revealed by high-resolution teleseismic tomography

We report a new model of the upper mantle structure beneath Italy obtained by means of P-wave teleseismic tomography. Besides the recent and remarkable development of the Italian Seismic Network, a high model resolution has been achieved improving the inversion method upon the ACH method used in previous investigations and picking high quality arrival times with the Multi-Channel Cross-Correlation technique. The finer details of our Vp model yield new insights into the heterogeneous structure of the Adria continental lithosphere involved in the collision between the Africa and Europe plates. A wide low Vp anomaly located in the northern Adria mantle, facing the Alpine high Vp slab, supports the idea that the Adria lithosphere has been hydrated and thinned during the Alpine subduction. We argue that this mantle softening may have played a key role in favoring the subsequent delamination of the Adria lithosphere in the northern Apennines. We hypothesize that delamination of continental lithosphere previously thinned in a back-arc setting may be considered a key process to favor subduction polarity reversal and recycling of continental material into the mantle circulation. Conversely, in the central-southern Apennines, the velocity structure is consistent with the existence of a deeper oceanic slab that flattens at the base of the upper mantle, in agreement with the widely accepted geodynamic evolution of the central Mediterranean by slab retreat and back-arc spreading. The oceanic slab is discontinuously detached from the surface plate, suggesting a different structure of the Adria lithosphere, which resists subduction instead of favoring delamination.     

未知OBS水平方位的SKS分裂参数测量

宽频带海底地震仪(OBS)探测是研究海底深部结构最有效的地球物理方法,其中SKS分裂参数测量已被广泛应用于上地幔各向异性的研究.传统的SKS分裂参数测量需要知道确切的OBS水平分量方位,然而OBS的投放过程是自由式的,其两个水平分量在海底的方位是不确定的.本文通过SKS快慢波的运动学特性,发展了一种在未知OBS水平分量方位的情况下直接求取SKS分裂参数的方法,分别用合成地震图和实际观测到的远震资料对该方法进行了数值检验,表明其是一种有效的获得SKS分裂参数的方法,而且相比于传统的先确定OBS水平分量方位,再进行SKS分裂的方法,此方法能够同时获得SKS分裂参数和OBS水平分量方位,应用简便,且减少了累积误差,获得的SKS分裂参数可靠性更高.     

The use of teleseismic P-wave complexity for seismic event screening – results determined from GRF and GERESS array data

Seismic data recorded at the broad-band teleseismic GRF array and theshort-period regional GERESS array, which is a designated IMS primarystation, are analyzed to determine the effectiveness of teleseismic P-wave complexity for the purpose of seismic event screening within theframework of Comprehensive Nuclear-Test-Ban Treaty verification. For theGRF array, seismic waveform data from nearly 200 nuclear explosions havebeen recorded since its installation in the late 1970's, which were studiedalong with several thousand earthquakes from the last few years.Additionally, we investigated teleseismic P wave complexity for a similarnumber of earthquakes recorded at GERESS. However, owing to itsoperation starting in 1991, only a limited number of nuclear explosionseismograms are available for study.For nuclear explosions, complexity does not exceed levels of 0.3 except fora number of events from the Nevada Test Site recorded only at the GRFarray and located at a large distance where PcP may interfere with the initialP wavelet. Since all events with complexity at GRF larger than 0.3 areexclusively located on Pahute Mesa within the Nevada Test Site,near-source geology or topography must play a dominant role for theseincreased complexity values, while PcP may not contribute significantly tothe high-frequency energy measured by the complexity parameter.Although many earthquakes show complexities below this level, for morethan 25% of the earthquakes investigated the complexities determined arelarger than 0.7, thus showing distinctly larger values than nuclearexplosions. Therefore, this percentage may be screened as earthquakes fromall seismic events detected. As currently only about half of the eventsdetected by the global IMS network are screened out based on focal depthand the m _b :M _s criterion, teleseismic P-wavecomplexity may contribute significantly to the task of seismic eventscreening.     

蒙古中部地区SKS分裂研究的初步结果

为了进行蒙古中部地区的深部构造观测与模型研究, 中国地震局地球物理研究所和蒙古科学院天文与地球物理研究中心合作, 在该地区架设了60台宽频带流动地震台. 这些台站分布在被额尔齐斯—中蒙古—额尔古纳断裂隔断的加里东、 海西两个构造域. 本文选取了2011年8月—2012年7月部分台站记录的震中距在85°—135°且M_W≥5.0的远震事件, 通过SplitLab软件来进行SKS分析. 使用最小能量法、 旋转相关法和最小特征值法进行处理, 取得了蒙古中部地区各向异性的初步结果. 结果显示: 快慢波延时在0.8—2.0 s之间. 西北部区域靠近杭爱山脉与肯特山脉的区域构造走向均为ENE--WSW方向, 而快波方向大多为N40°W, 与区域构造走向近似垂直.根据HS3-NUVEL-1A板块模型计算的该地区板块绝对运动方向与快波偏振方向呈顺时针约15°夹角. 然而, 有个别台站表现出了两种模式的快波分裂, 其中一种与当地ENE--WSW向的断裂方向近似一致.      

On the <Emphasis Type="Italic">P</Emphasis> and <Emphasis Type="Italic">S</Emphasis> Receiver Functions Used for Inverting the One-Dimensional Upper Mantle Shear-Wave Velocities

The teleseismic P receiver functions are customarily inverted to attain the seismic velocities beneath a seismic station. Surface wave dispersion data are often added to reduce the effect of the non-uniqueness. The combination of P receiver function and surface wave works well in resolving the structures in the crust and uppermost mantle, but is less effective in characterizing greater (lithosphere and asthenosphere) depths due to the interference from crustal multiples. A solution to this problem is jointly to model teleseismic S receiver functions with surface wave and P receiver functions. This study adopts a fast, one-dimensional (1-D) inversion scheme. To avoid the effect of multidimensional structures away from the seismic station, we eliminate multiples that reverberate between the surface and interfaces below a restriction depth (RD), as well as S-to-P conversions below an inversion depth (ID). P-to-S conversions off the interfaces above the half-space and S-to-P conversions above the ID and multiples above the RD are properly modelled. This approach favours ray paths travelling close to stations and is, therefore, more suitable for 1-D inversions. We perform numerical experiments with and without noise and highlight the advantages of a joint receiver function and surface wave analysis.     

青藏高原东北部岩石圈有效弹性厚度及其各向异性

青藏高原东北部是中国大陆构造环境特殊的主要构造域,毗邻青藏高原羌塘地块、塔里木盆地、四川盆地和华北克拉通,属于不同构造类型块体俯冲、碰撞及陆内汇聚的结合部,在中国大陆形成与演化的历史中扮演着重要角色.岩石圈有效弹性厚度(T_e)及其各向异性与岩石圈流变性、力学结构紧密相连,研究青藏高原东北部的岩石圈T_e及其各向异性将为我们认识大陆岩石圈的流变性及动力学过程提供重要信息.本文基于Fan小波相关性分析法,运用布格重力和地形资料获得了青藏高原东北部岩石圈T_e及其各向异性二维分布的详细信息.研究结果表明研究区域内T_e的分布范围在5～100 km之间;松潘-甘孜地块、祁连山造山带和龙门山地区的T_e较薄(5 kme<40 km)、各向异性较强;而周缘的断裂带、缝合带的T_e值都较低,其中龙门山断裂带T_e只有5～20 km,且南、北两段各向异性存在明显差异.内部的若尔盖盆地T_e值略显高值,说明其是仍保留有刚性的块体;北缘的柴达木盆地下伏为古生代的地壳,被认为是古老的克拉通碎片,T_e较大(50 kme<90 km),显示为轮廓分明的刚性块体.并且我们发现研究区域内T_e的各向异性轴垂直于大的块体边界.通过比较T_e各向异性与SKS波的快波偏振方向、Rayleigh面波方位角各向异性的相互关系我们推测阿拉善地块各向异性源自地幔橄榄岩晶格的优势取向,岩石圈变形趋于垂直连贯变形模式;柴达木盆地各向异性源于历史构造事件残留在岩石圈中的"化石"各向异性;松潘-甘孜地块各向异性源自物质的侧向流动.     

鄂尔多斯块体北缘呼包鄂地区岩石圈各向异性初步研究

选取呼包鄂区域11个宽频带地震台站记录的156个远震事件,使用剪切波分裂软件,对所选远震波形数据进行处理,求取并分析震中距85°—110°范围内同时被各地震台站记录的远震震相SKS的各向异性参数,即快波偏振方向及慢波到时差,结果表明:呼包鄂区域地下深层介质存在各向异性,快波偏振优势方向为ESE,慢波延迟时间为0.87—...