南北地震带北段的远震P波层析成像研究

本文利用“中国地震科学台阵”探测项目在南北地震带北段布设的678个流动地震台站在2013年10月至2015年4月期间记录到的远震波形数据,经过波形互相关拾取到473个远震事件共130309条P波走时残差数据,通过远震层析成像研究获得了该区（30°N-44°N,96°E-110°E）下方0.5°×0.5°的P波速度扰动图像.结果显示,研究区下方P波速度结构显示强烈的不均一性和显著的分区、分块特征.岩石圈速度结构具有显著的东西差异：祁连、西秦岭和松潘甘孜地块组成的青藏东北缘地区显示明显的低速异常,而属于克拉通性质的鄂尔多斯地块和四川盆地则显示高速异常,表明东部克拉通块体对青藏高原物质的东向挤出起到了强烈的阻挡作用.阿拉善地块显示出弱高速和局部弱低速的异常并存的特征.阿拉善地块西部显示低速异常,而东部与鄂尔多斯相邻的地区显示高速异常,可能表明该地区的岩石圈的变形主要受到青藏高原东北缘的挤压作用.在鄂尔多斯和四川盆地之间的秦岭下方100~250 km深度上表现为明显的低速异常,表明该处可能存在软流圈物质的运移通道.鄂尔多斯北部的河套裂陷盆地下方在100~500 km深度内低速异常表现明显,说明该区有深部热物质上涌且至少来源于地幔过渡带.青藏东北缘上地幔显示低速异常且地幔过渡带中出现明显的高速异常,这种结构模式暗示了在青藏高原东北缘可能发生了岩石圈拆沉作用,而高速异常体可能是拆沉的岩石圈地幔.

     

Upper mantle structure of the Apulian plate from Rayleigh waves

Phase velocities of Rayleigh waves for the Adriatic Sea area are obtained in the period range 25–190 sec along the path (l'Aquila-Trieste) AQU-TRI and 20–167 sec along the path (Trieste-Bari) TRI-BAI.The phase velocities are systematically higher than the known values for the surrounding regions. The data inversion indicates the presence of a lithosphere typical of stable continental areas with clear high-velocity lid (V _s 4.6 km/sec) overlying a well developed low velocity zone (V _s 4.2 km/sec).P. F. Geodinamica C.N.R., Roma Pubbl. N. 189.     

应用远震有限频率层析成像反演首都圈上地幔速度结构

本文选用首都圈数字地震台网2003年9月~2005年12月记录的300多个远震事件的波形资料,采用分频带多道互相关方法得到三个不同频段的P波相对走时数据共18499个,计算了每个频段的走时灵敏度核,应用有限频率层析成像反演得到首都圈地区的上地幔三维P波速度结构模型.利用检测板估计了反演结果的分辨率,并与射线层析成像方法的结果进行了比较,说明了反演结果的可靠性.研究结果表明,各构造单元具有明显不同的速度结构特征,其差异可到150 km深:燕山隆起区表现高速;太行山隆起区整体以低速为主并存在小范围高速块体;华北盆地、渤海湾下浅层上地幔中存在大范围的强低速异常,其顶面在50~70 km,可视为软流圈顶面的埋深,这一结果说明华北盆地、渤海湾下岩石圈明显减薄;张家口—蓬莱断裂带是上地幔浅部速度结构的变异带,也是岩石圈减薄的边界带,区内大部分强震都发生在该构造带上,由此看来该带上强震的发生不仅与地壳结构的不均匀性有关,还可能有较深的构造背景.     

Three-dimensional structure of the crust and upper mantle beneath the south China area from long period surface waves

In the present study the long period surface wave records of 238 wave-paths from 79 earthquakes within China and its adjacent regions received by 30 seismic network stations are measured by using the improved match-filtering frequency-time analysis technique and the grid dispersion inversion method to obtain the rayleigh pure-path dispersion values for 147 slant grids of 4° × 4° in this area, then a three-dimensional shear wave velocity model of the crust and upper mantle beneath south China area to a depth of 170 km is inversed. It is found that there are obvious differences among the main structural units, and there are also certain differences among the subordinate elements even in the individual unit. The crustal thickness of this area is ranging from 30 to 43 km, and is getting thicker gradually from the east to the west. The average shear velocity of crust is ranging form 3.48 to 3.68 km/s with the lowest in the northeast part and highest in the west part. No obvious crustal low velocity layer of large scale is detected. There exist upper mantle low velocity zones in the most of south China area with the starting depth ranging from 75 to 106 km. The lowest shear velocity within the low velocity zones is about 4.28–4.38 km/s. Despite of the existing of upper mantle low velocity zones beneath the most of south China area, the interfaces between the important layers are quite clear, the variation of the bedding surfaces is very gentle, and the lateral changes measured in a larger scale of the underground structure are rather small. It may indicate that the crustal and upper mantle structure of the main part of south China area belongs to the relatively stable structure of the continental blocks except for the fringe areas such as the fold-faulted region in the west part and the fault system along the southeastern coast which may belong to the tectonically active area. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,15, 159–167, 1993. This subject is supported by the National Science Foundation of China.     

基于P波三重震相的华南地区上地幔速度结构研究

华南块体是研究太平洋板块俯冲和岩石圈减薄机制等问题的最佳场所之一.本文基于中国地震观测台网和大型流动台阵记录到的震中距10°~30°之间的两个中深源地震P波记录,利用三重震相波形拟合技术,获得了中扬子克拉通和华夏地块上地幔高精度P波速度结构.研究结果表明:(1)中扬子克拉通过渡带底部存在高速异常,系太平洋俯冲板块的滞留体.俯冲的板块并没有进入下地幔,660-km间断面下沉约11 km,与后尖晶石相变的克拉伯龙斜率为负有关.而华夏地块过渡带底部并无明显高速异常,接近全球平均模型;(2)整个华南块体,410-km间断面上方普遍存在低速层,主要与上地幔部分熔融有关,与IASP91相比P波速度减小了1.38%~2.29%;(3)在研究区域内,中扬子克拉通和华夏地块都存在岩石圈减薄(80 km),推测可能与太平洋板块俯冲和快速回撤导致的岩石圈拆沉有关.且华夏地块减薄程度较明显,下伏软流圈速度较低,说明其上地幔强度较弱、温度较高.另外,中扬子克拉通过渡带中存在一个较宽的速度梯度带,可能与520-km间断面有关,其具体成因有待进一步研究.     

Upper mantle anisotropy beneath central and southwest Japan: An insight into subduction-induced mantle flow

Analysis of seismic anisotropy in the crust and mantle wedge above subduction zones gives much information about the dynamic processes inside the Earth. For this reason, we measure shear wave polarization anisotropy in the crust and upper mantle beneath central and southwestern Japan from local shallow, intermediate, and deep earthquakes occurring in the subducting Pacific slab. We analyze S phases from 198 earthquakes recorded at 42 Japanese F-net broadband seismic stations. This data set yields a total of 980 splitting parameter pairs for central and southwestern Japan. Dominant fast polarization directions of shear waves obtained at most stations in the Kanto–Izu–Tokai areas are oriented WNW–ESE, which are sub-parallel to the subduction direction of the Pacific plate. However, minor fast polarization directions are oriented in NNE–SSW directions being parallel to the strike of the Japan Trench, especially in the north of Izu Peninsula and the northern Tokai district. Generally, fast directions obtained at stations located in Kii Peninsula and the Chubu district are oriented ENE–WSW, almost parallel to the Nankai Trough, although some fast directions have NW–SE trends. The fast directions obtained at stations in northern central Honshu are oriented N–S. Delay times vary considerably and range from 0.1 to 1.25 s depending on the source depth and the degree of anisotropy along the ray path. These lateral variations in splitting character suggest that the nature of anisotropy is quite different between the studied areas. Beneath Kanto–Tokai, the observed WNW–ESE fast directions are probably caused by the olivine A-fabric induced by the corner flow. However, the slab morphology in this region is relatively complicated as the Philippine Sea slab is overriding the Pacific slab. This complex tectonic setting may induce lateral heterogeneity in the flow and stress state of the mantle wedge, and may have produced NNE–SSW orientations of fast directions. The ENE–WSW fast directions in Kii Peninsula and the Chubu district are more coherent and may be partly induced by the subduction of the Philippine Sea plate. The N–S fast directions in northern central Honshu might be produced by the trench-parallel stretching of the wedge due to the curved slab at the arc–arc junction.     

Structure of the upper mantle in Asia from phase and group velocities of Rayleigh waves

Dispersion curves of phase velocities of Rayleigh waves are determined by the method of frequency-time analysis in a range of periods of 10–200 s from data of 43 interstation traces in Central Asia. Because the joint use of phase and group velocities significantly decreases the uncertainty in the determination of S wave velocity structures, the same traces were used for calculating group velocities from tomographic reconstructions obtained in [Yanovskaya and Kozhevnikov, 2003, 2006] and determining average velocity structures along these traces. The velocity structures were calculated by the Monte Carlo and linear inversion methods, which gave consistent results. Using velocity values obtained at fixed depths by the 2-D tomography method, lateral variations in velocities at these depths were estimated, which allowed us to construct smoothed vertical velocity structures at some points in the region. The resulting structures were used as initial approximations for constructing local velocity structures solely from previously obtained local dispersion curves of group velocities in the area (32°–56°N, 80°–120°E). Based on these structures, we mapped the lateral distribution of velocity variations at upper mantle depths of 75–400 km and along three vertical profiles. The inferred velocity variations are in good agreement with data on the tectonics of the region.     

Upper mantle velocities on the Northern Cocos Plate

Refraction lines on the Northern Cocos Plate between the Orozco and Clipperton Fracture Zones have been used to determine upper mantle velocities over the plate. The velocities range from 7.50 to 8.43 km/sec. Azimuthal variations are found near the rise crest with low velocities parallel to the rise crest. The low velocities increase with age, lessening the observed azimuthal dependence away from the rise crest. A low-velocity zone is found in the mantle and may extend over a considerable portion of the plate near the rise crest. A 7.1-km/sec basal crustal layer is also observed and makes up a substantial portion of the crustal thickness.     

基于电性结构模型的青藏高原东缘上地幔热结构研究

研究青藏高原东缘地区的深部物质结构对于理解青藏高原的隆升及扩张机制具有重要的科学意义.本文将青藏高原东缘实测大地电磁测深剖面反演所得的岩石圈电性结构模型与高温高压岩石物理实验测得的上地幔矿物和熔融体导电性定量关系相结合,通过Hashin-Shtrikman（HS）边界条件建立上地幔电导率与温度、熔融百分比等参数的定量关系,在此基础上计算得到了青藏高原东缘上地幔热结构及熔融百分比分布模型.研究结果表明在青藏高原东缘地区通过大地电磁测深方法所探测到的上地幔低阻体可以解释为由高温作用所产生的局部熔融区域.其中,松潘—甘孜地块上地幔高导体对应的温度介于1300~1500℃之间,熔融百分比可高达10%,支持前人将松潘—甘孜地块内部的低阻体解释为局部熔融的观点.龙门山断裂带以东、四川盆地西缘的上地幔高导体温度介于1200~1400℃之间,熔融百分比介于1%~5%左右,表明扬子克拉通的西缘可能正在经历一定程度的活化作用.龙门山断裂带下方的上地幔高阻体温度介于1100℃附近,基本没有发生局部熔融,具有较冷的刚性块体特征,与该区域频发的地震活动相吻合.四川盆地东部的扬子上地幔温度介于800~900℃之间,没有发生局部熔融,符合古老稳定的克拉通块体的基本特征.

     

A 3D model of the central asia mantle from dispersion of Rayleigh-wave phase velocities

Frequency-time analysis was applied to records of Rayleigh surface waves due to teleseismic earthquakes (M ≥ 6.0) to obtain dispersion curves of fundamental mode phase velocities for paths between 51 pairs of digital seismic stations in the Asian IRIS networks; the range of periods was 10–200 s. For each of these pairs we derived 1D-shear-velocity sections that fit the corresponding dispersion curves and image the integral earth structure down to depths of ∼650–700 km beneath each of the interstation seismic paths, which traverse the area of study in different directions. These structures were used by Yanovskaya (2001) to derive a 3D-shear-velocity model by 2D tomography, imaging the distribution of the larger horizontal inhomogeneities in the Central Asia mantle for depths of 50 to 600 km, and to construct 2D-velocity sections for five lines passing through the major tectonic features of Central Asia.     

Apparent phase velocities and fundamental-mode phase velocities of Rayleigh waves

Rayleigh wave dispersion data usually appear in the form of multimodal spectra for a layered model structure. The number of dispersion modal curves depends on the number of layers in the model. The measured dispersion velocities from the frequency–wavenumber (F–K) space, however, may not represent the true phase velocities of the fundamental-mode dispersion curve, but apparent phase velocities. The present study discusses how multimode curves are generated in the multichannel analysis of surface waves (MASW) method and the cause of the apparent velocity produced by the F–K method. Results from a field trial demonstrate that apparent phase velocities fail to reveal thin layers with low velocities. A better agreement of the inverted model with the geotechnical record is obtained by using the data points extracted from the fundamental-mode curve of the MASW spectral image.     

Deep structure of the bohemian massif from phase velocities of Rayleigh and Love waves

Summary Phase velocities of Rayleigh and Love waves have been measured between the broadband seismic stations KHC (Kaperské Hory, South Bohemia) and KSP (Ksi, Lower Silesia), a profile that nearly coincides with the Interactional DSS Profile VII. The data for both wave types were separately inverted into models of shear-wave velocity versus depth. Novotný's model KHKS 82[1] for the DSS Profile VII was used as a start model. While the crustal section of Novotný's model is compatible with both of our data sets, our Rayleigh-wave data require smaller shear-wave velocities, on the average by 0.24 km/s, in the top 180 km of the mantle. The average difference between Novotný's model and our Love-wave model in that depth range is only 0.06 km/s. If our identification of the observed Love waves as the fundamental mode is correct, this result indicates the presence of polarization anisotropy in the uppermost mantle.

---

u a mu u a ¶rt; unu uuu mauu (an , a u) u S (, ua uu). u - S nuuum mmmm ¶rt;a¶rt; nu VII . a ¶rt; u mun m¶rt; umua ¶rt;u auumu mu nn mu. am mam¶rt;u a unaa ¶rt; S 82 m[1] ¶rt; nu VII . a am ¶rt;u m aua uu aau au ¶rt;a, ¶rt;a au ¶rt;a ¶rt; mm uu m nn — ¶rt; a 0,24 / 180 uma amuu. ¶rt;a aua ¶rt; ¶rt; m u a ¶rt; ¶rt; a m ¶rt;uanau n¶rt;mam 0,06 /. u aa u¶rt;muuau a¶rt;a a a ¶rt;ama ¶rt;a m nau, m mm mam u¶rt;ummm numuu nuau aumnuu amuu.

     

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.     

Electrical structure of the upper mantle beneath Central Europe: Results of the CEMES project

In the years 2001–2003, we accomplished the experimental phase of the project CEMES by collecting long-period magnetotelluric data at positions of eleven permanent geomagnetic observatories situated within few hundreds kilometers along the south-west margin of the East European Craton. Five teams were engaged in estimating independently the magnetotelluric responses by using different data processing procedures. The conductance distributions at the depths of the upper mantle have been derived individually beneath each observatory. By averaging the individual cross-sections, we have designed the final model of the geoelectrical structure of the upper mantle beneath the CEMES region. The results indicate systematic trends in the deep electrical structure of the two European tectonic plates and give evidence that the electrical structure of the upper mantle differs between the East European Craton and the Phanerozoic plate of west Europe, with a separating transition zone that generally coincides with the Trans-European Suture Zone.     

The velocity structure of the upper mantle in the transition zone from the East European Platform to Western Europe from seismic noise data

The surface wave tomography from ambient seismic noise recorded at stations in Western Europe (WE) and on the East European Platform (EEP) revealed the structure of the crust and upper mantle in the transitional zone from the Precambrian platform to the younger geological units in Western Europe. The Tornquist-Teisseyre Line separating these structures is clearly traced as a transition zone from the high velocities beneath EEP to the low velocities beneath WE in the crust and upper mantle, which extends to a depth of 150?C170 km. Below 200 km the relationship between the velocities beneath EEP and WE becomes the opposite. A similar relationship between the velocities in the upper mantle down to a depth of 300 km is observed on the southern boundary, where EEP borders on the northern segment of the Alpine-Himalayan seismic belt.     

Velocity structure of uppermost mantle beneath China continent from Pn tomography

39473 Pn travel times are inverted to tomographically image both lateral variation and anisotropy of uppermost mantle velocities beneath China continent. The result indicates that the overall average Pn velocity of uppermost mantle in the studied region is 8.0 km/s and the regional velocity fluctuation varies from ?0.30 km/s to +0.35 km/s. Pn velocities higher than 8.2 km/s are found in the regions surrounding Qingzang Plateau, such as Junggar Basin, Tarim Basin, Qaidam Basin and Sichun Basin. Pn velocities slightly lower than the average are found in western Sichuan and Yunnan, Shanxi Graben and Bohai Bay region. A Pn velocity as low as 7.8 km/s may exist in the region striding the boundary between Guangxi and Guangdong provinces. In general, Pn velocity in tectonically stable region like cratonic platform tends to be high, while that in tectonically active region tends to be low. The regions in compressive setting usually show higher Pn velocity, while extensional basins or grabens generally display lower one. Anisotropy of Pn velocity is seen in some regions. In the southeastern region of Qingzang Plateau the directions of fastest Pn velocity show a rotation pattern, which may be related to southeastward escape of the plateau material due to the collision and compression of Indian Plate to Asia along Himalaya arc. Notable anisotropy also exists around Bohai Bay region, likely indicating crustal extending and possible magma activity therein.     

利用剪切波分裂研究内蒙古阿巴嘎地区上地幔各向异性

内蒙古阿巴嘎地区壳幔经历强烈变形,岩石圈变形机制尚不明确.利用布设在研究区的32个流动地震台站所记录到的远震剪切波数据,测量得到120对各向异性参数和113个无效分裂结果.结果表明,研究区快慢波延迟时间变化范围为0.4~1.4 s,平均0.77±0.21 s;各向异性快波方向变化范围为N101°E-N45°W.其中一组快波偏振方向为N82.0°E±12.3°,与区域内断裂走向平行,反映地幔矿物晶格定向排列;另一组快波方向集中位于华北克拉通内部,平均为N146.8°E±9.5°,平行于早白垩纪岩石圈伸展变形方向,推测由残留在岩石圈中的化石各向异性所引起.在研究区北部部分台站,只观测到无效分裂而没有观测到有效分裂结果,可能存在局部热地幔物质上涌.

     

Crust and upper mantle structure of the Bohemian Massif from the dispersion of seismic surface waves

Summary The phase velocity dispersion of Rayleigh waves for the Moxa-Vienna (MOX-VIE) and Moxa-Kaperské Hory (MOX-KHC) profiles, and of both Rayleigh and Love waves for the Kaperské Hory-Ksi (KHC-KSP) profile have been measured and inverted into models of shearwave velocity vs. depth. The three paths cross, respectively, the central part of the Bohemian Massif, its western margin, and the Bohemian Pluton and Cretaceous. For the MOX-VIE profile mean and lower crustal shear wave velocities of 3.7 and 3.9 km/s, respectively, a mean Moho depth of 34 km, and no existence of a low-velocity layer in the lower crust were found. The model obtained for the MOX-KHC profile is characterized by a slightly lower velocity in the lower crust (3.8 km/s), by a slightly lower Moho depth (32 km), and by the appearance of a weak low-velocity channel between 55 and 140 km. The crustal section of the final model for the KHC-KSP profile agrees well with the KHKS82 model derived by Novotný from results of DSS along international profile VII. Our final Rayleigh-wave model has significantly lower shear-wave velocities down to 215 km in the mantle. A systematic difference of 0.18 km/s between the average velocities of Rayleigh and Love waves has been revealed for the depth range from 30 to 215 km. Since almost no contamination of the fundamental Love mode with higher modes has been observed, and since the investigated structure hardly contains an unresolved system of thin, alternately low- and high-velocity layers, the cause of the difference is evidently polarization anisotropy of the upper mantle beneath the Bohemian Massif. It is recommended that the discussed investigations should be supplemented with data from the fan of KSP-GRF (Gräfenberg Array, Germany) paths and from the KHC-BRG (Berggiesshübel, Germany) profile.