用面波方法研究上扬子克拉通壳幔速度结构

本文研究采用单台法和双台法提取了穿越上扬子的基阶面波相速度和群速度频散;通过对提取的面波群速度和相速度频散进行联合反演,得到的1-D SV速度模型显示上扬子块体下地壳S波速度与典型克拉通区域相当,其上地幔顶部80~170 km深处存在高速的岩石圈盖层,较AK135模型要快2%~3%,其岩石圈厚度约为180 km.在上扬子地区,径向各向异性集中分布在300 km以浅的岩石圈与软流圈部分,其中岩石圈部分SH波比SV波波速要快2%~4%,软流圈部分SH波比SV波波速要快3%~5%;Rayleigh波相速度方位各向异性分析结果显示,上扬子块体周期为25~45 s（大致相当于30~70 km深度范围内）的Rayleigh波相速度存在1.8%~2.7%不等的方位各向异性,其快波方向介于147°~174°.我们认为上扬子块体径向各向异性集中分布在岩石圈、软流圈部分,且各向异性随深度变化, 其岩石圈部分各向异性为大陆克拉通化的遗迹,软流圈部分各向异性与现今板块运动相关.     

地幔波震级的规范化标定与观测

本文基于瑞利面波的激发和传播理论,导出了零级地震的理论面波谱表达式;并以零级地震为参考地震,以零级地震所相应的面波谱和地震矩作为地幔波震级规范化标定的参照依据,进而提出了测定地幔波震级的新方案。利用CDSN-LPZ信道中的面波资料测定了1987—1990年间全球54次浅源大地震的地幔波震级M_m,并对Okal和Talandier(1989)用GEO-SCOPE-VLPZ及PAS-VLPZ信道中的基阶瑞利波所测1970—1987年间88次地震的结果作了校正。综合这两部分观测结果,得到M_m=1.00 1gM_0—12.13,由此测定的地幔波震级与中小地震的面波震级M_s趋势相一致,为解决面波震级饱和问题提供了一种可行的途径。     

中国东北地区地壳上地幔三维S波速度结构

收集了中国东北地区159个固定地震台2011年1月至2012年6月和27个流动地震台2011年1月至2011年6月间的垂向连续记录,根据噪声成像方法得到研究区（105°E-135°E, 39°N-52°N）较短周期（8~30 s）的瑞雷波群速度和相速度频散资料,再结合该区已有的天然地震长周期瑞雷波（36~145 s）的群速度频散资料,我们反演得到了中国东北地区200 km以浅深度范围内的三维壳幔S波速度结构,并得到了该区的岩石圈厚度分布图.结果表明：研究区中、下地壳S波速度结构的横向分布,在重力梯度带两侧有很大的不同,以东地区显示为大范围的高速,以西地区则呈现为大面积的低速;松辽盆地下方岩石圈地幔表现为显著的高速,岩石圈地幔底界面深度可能在90~100 km,薄的岩石圈盖层暗示东北地区的岩石圈可能发生了减薄;郯庐大断裂下方呈现出大范围的比较显著的低速特征,断裂下方上地幔顶部可能有热物质活动.     

用数字台网的面波资料研究中国大陆的地壳上地幔结构

本文采用适配滤波频时分析技术对中国数字地震台网记录的长周期面波资料进行处理,获得了路经中国大陆的23条波路径的基阶瑞雷面波群速度频散曲线,频散曲线的周期从10.5秒至150秒。使用HARKRIDER面波反演程序反演了各条波路径的平均地壳上地幔结构。所得模型直至230km。深处仍有较好的分辨。根据所得结果将中国大陆划分为四个构造单元。各构造单元之间差异明显且在某些单元内部横向变化也很显著。     

Love-Rayleigh wave incompatibility and possible deep upper mantle anisotropy in the Iberian peninsula

Love and Rayleigh wave phase velocities are analyzed with the goal of retrieving information about the anisotropic structure of the Iberian lithosphere. The cross-correlation method is used to measure the interstation phase velocities between diverse stations of the ILIHA network at periods between 20 and 120 s. Despite the 2-D structure of the network, the Love wave data are too few to enable an analysis of phase velocity azimuthal variations. Azimuthal averages of Love and Rayleigh wave phase velocities are calculated and inverted both in terms of isotropic and anisotropic structures. Realistic isotropic models explain the Rayleigh wave and short-period Love wave phase velocities. Therefore no significant anisotropy needs to be introduced in the crust and down to 100 km depth in the upper mantle to explain our data. A discrepancy is observed only at long periods, where the data are less reliable. Love wave data at periods between 80 and 120 s remain 0.15 km/s faster than predicted by isotropic models explaining the long-period Rayleigh wave data. Possibilities of biases in the measurements due to interferences with higher modes are examined but seem unlikely. A transversely isotropic model with 8% of S-wave velocity anisotropy in the upper mantle at depths larger than 100 km can explain the whole set of data. In terms of a classical model of mantle anisotropy, this corresponds to 100% of the crystals perfectly oriented in the horizontal plane in a pyrolitic mantle. This is a rather extreme model, which predicts at time delay between 0 and 2 seconds for split SKS.     

地幔波震级的规范化标定与观测

本文基于瑞利面波的激发和传播理论,导出了零级地震的理论面波谱表达式;并以零级地震为参考地震,以零级地震所相应的面波谱和地震矩作为地幔波震级规范化标定的参照依据,进而提出了测定地幔波震级的新方案。利用CDSN-LPZ信道中的面波资料测定了1987-1990年间全球54次浅源大地震的地幔波震级M_m,并对Okal和Talandier（1989）用GEO-SCOPE-VLPZ及PAS-VLPZ信道中的基阶瑞利波所测1970-1987年间88次地震的结果作了校正。综合这两部分观测结果,得到M_m=1.00 1gM₀-12.13,由此测定的地幔波震级与中小地震的面波震级M_s趋势相一致,为解决面波震级饱和问题提供了一种可行的途径。     

中国陆区S波分型与上地幔地震各向异性

本文利用中国大陆区共8个台站所观测到地震事件的S波数据，采用相关分析的手段以 观测数据重建的技术，进行了S波分裂参数－快慢波延迟时间与快波偏振方向的计算，并分析了S波分裂参数与中国陆区上地幔物质运动的相互关系，研究了结果表明，兰州台下的S波分裂时差最小，北京台下的最大，中国东部地区上地幔运动方向以西北向为主，西部地区以北偏西为优势。     

First-principles calculations of elasticity of minerals at high temperature and pressure

The elasticity of minerals at high temperature and pressure(PT) is critical for constraining the composition and temperature of the Earth's interior and understand better the deep water cycle and the dynamic Earth. First-principles calculations without introducing any adjustable parameters, whose results can be comparable to experimental data, play a more and more important role in investigating the elasticity of minerals at high PT mainly because of(1) the quick increasing of computational powers and(2) advances in method. For example, the new method reduces the computation loads to one-tenth of the traditional method with the comparable precise as the traditional method. This is extraordinarily helpful because first-principles calculations of the elasticity of minerals at high PT are extremely time-consuming. So far the elasticity of most of lower mantle minerals has been investigated in detail. We have good idea on the effect of temperature, pressure, and iron concentration on elasticity of main minerals of the lower mantle and the unusual softening in bulk modulus by the spin crossover of iron in ferropericlase. With these elastic data the lower mantle has been constrained to have 10–15 wt% ferropericlase, which is sufficient to generate some visible effects of spin crossover in seismic tomography. For example, the spin crossover causes that the temperature sensitivity of P wave at the depth of ~1700 km is only a fraction of that at the depth of ~2300 km. The disruptions of global P wave structure and of P wave image below hotspots such as Hawaii and Iceland at similar depth are in consistence with the spin crossover effect of iron in ferropericlase. The spin crossover, which causes anomalous thermodynamic properties of ferropericlase, has also been found to play a control role for the two features of the large low shear velocity provinces(LLSVPs): the sharp edge and high elevation up to 1000 km above core-mantle boundary. All these results clearly suggest the spin crossover of iron in the lower mantle. The theoretical investigations for the elasticity of minerals at the upper mantle and water effect on elasticity of minerals at the mantle transition zone and subducting slab have also been conducted extensively. These researches are critical for understanding better the composition of the upper mantle and water distribution and transport in the Earth's mantle. Most of these were static calculations, which did not include the vibrational(temperature) effect on elasticity, although temperature effect on elasticity is basic because of high temperature at the Earth's interior and huge temperature difference between the ambient mantle and the subducting slab. Including temperature effect on elasticity of minerals should be important future work. New method developed is helpful for these directions. The elasticity of iron and iron-alloy with various light elements has also been calculated extensively. However, more work is necessary in order to meet the demand for constraining the types and amount of light elements at the Earth's core.     

Crustal structure of the northeastern Tibetan plateau,the Ordos block and the Sichuan basin from ambient noise tomography

We apply ambient noise tomography to significant seismic data resources in a region including the northeastern Tibetan plateau,the Ordos block and the Sichuan basin.The seismic data come from about 160 stations of the provincial broadband digital seismograph networks of China.Ambient noise cross-correlations are performed on the data recorded between 2007 and 2009 and high quality inter-station Rayleigh phase velocity dispersion curves are obtained between periods of 6 s to 35 s.Resulting Rayleigh wave phase velocity maps possess a lateral resolution between 100 km and 200 km.The phase velocities at short periods (20 s) are lower in the Sichuan basin,the northwest segment of the Ordos block and the Weihe graben,and outline sedimentary deposits.At intermediate and long periods (25 s),strong high velocity anomalies are observed within the Ordos block and the Sichuan basin and low phase velocities are imaged in the northeastern Tibetan plateau,reflecting the variation of crustal thickness from the Tibetan plateau to the neighboring regions in the east.Crustal and uppermost mantle shear wave velocities vary strongly between the Tibetan plateau,the Sichuan basin and the Ordos block.The Ordos block and the Sichuan basin are dominated by high shear wave velocities in the crust and uppermost mantle.There is a triangle-shaped low velocity zone located in the northeastern Tibetan plateau,whose width narrows towards the eastern margin of the plateau.No low velocity zone is apparent beneath the Qinling orogen,suggesting that mass may not be able to flow eastward through the boundary between the Ordos block and the Sichuan basin in the crust and uppermost mantle.     

中长周期数字化面波记录与中国东南地区地壳结构

本文使用适配滤波频时分析技术首次对中国数字地震台网的中长周期面波记录进行处理,获得穿过东南地区的82条勒夫波频散数据.使用随机反演理论,获得了东南沿海地区4°×4°网格的纯路径频散数据.这些频散的周期为1.95-68.27s,弥补了长周期面波所不能分辨的浅层结构.在网格反演的基础上,使用Harkrider的面波及演程序得出了中国东南地区的地壳和上地幔结构,浅部可分辨到1km,深部可达80km.在分辨率保证的前提下得出东南地区深至80km的三维剪切波速度结构.     

全球地震台网记录苏门答腊-安达曼地震的地幔波

介绍印度尼西亚苏门答腊-安达曼Ms8．7地震在全球地震台网（GSN）及所含“新一代中国数字地震台网”（NCDSN）记录的地幔面波。主要讨论本次大地震地幔瑞利波Ri,j（环球瑞利面波）无频散的3种可能原因,同时提出地幔妹泽波Mi,j的问题。     

中国西部三维速度结构及其各向异性

本文用覆盖中国的358条勒夫面波路径资料,研究了10.45-113.80s范围内中国西部的三维SH波速度结构.结果表明,各构造单元的SH波速度结构均有明显的差别.作为稳定块体的塔里木盆地,壳内重力分异程度较高,上、中、下地壳厚度差别小,壳内无明显的低速层,地壳平均速度比较小;上地幔低速层埋深大且层中速度大;区内横向变化小.构造活动区如天山、青藏高原,其突出的特征是下地壳厚度大且速度大,上地幔盖层速度值相当高.这与西伯利亚、印支板块的挤压有密切的关系.青藏高原东部及其北、东边缘地区壳内存在低速层,上地幔低速层埋深浅,一些地区存在壳幔过渡层.面波各向异性研究表明,青藏高原、天山及印支板块北缘下存在明显的各向异性,以构造边缘地区及上地幔低速层附近最为突出.印度板块、西伯利亚板块与中国大陆间的碰撞引起强大的水平压力和一定的下插作用,是造成青藏高原隆起、地壳增厚、天山隆起的最根本的因素,同时也促成壳幔中辉石、橄榄石的定向排列和物质运移,因而出现明显的各向异性现象.     

中国东部海域地壳-上地幔瑞利波速度结构研究

为了进一步了解中国东部沿海及相邻海域的地壳－上地幔结构特征,对该区域的构造演化历史、地震活动及深部构造等方面研究提供一些基础资料,利用３１个数字地震台记录的高质量瑞利波资料,采用一种新的混合路径频散的网格反演方法（Ｏｃｃａｍ方法）,对中国东部海域瑞利波群速度横向不均匀分布进行了初步研究．根据反演得到的１０-１５０ｓ共３６个中心周期的群速度分布特征,以及几个典型地点的剪切波速度结构的深度变化,对研究区域内各构造单元的划分以及它们在速度结构和上地幔低速层埋深等方面的特征进行了讨论。     

秦岭及周边地区瑞雷波方位各向异性

除了使用前人提取的面波频散曲线外, 还从秦岭及周边地区布设的59个宽频流动地震台的数据和中国地震局及各省局台网的数据中筛选出的地震事件波形和台站间噪声互相关格林函数中提取瑞雷波群速度频散曲线. 利用二维面波层析成像反演获得了瑞雷波周期为10—50 s的各向同性群速度及方位各向异性分布. 结果显示: 周期为10 s的各向同性群速度和方位各向异性分布与各构造单元存在明显的对应关系; 周期为10—50 s的面波在四川盆地和鄂尔多斯盆地内的快波方向多为近NS向. 与前人研究结果不同的是, 本文得到的秦岭、 大巴山构造带周期为10—50 s的面波快波方向均与山脉走向近似平行, 且与SKS波分析得到的快波方向一致. 这表明秦岭和大巴山之下整个岩石圈的快波方向都与山脉走向平行, 预示着秦岭和大巴山整个地壳, 甚至岩石圈发生了类似的变形. 由于四川盆地和鄂尔多斯盆地面波快波方向与SKS波结果差别较大, 推测青藏高原隆升扩展对这两个盆地的地壳基本无影响, 但对其岩石圈上地幔却产生了重大影响.      

黑龙江地区背景噪声面波群速度层析成像

本文选取黑龙江省台网及吉林、内蒙古部分台站共44个固定地震台站连续2年的地震数据,利用两个台站间互相关方法获得瑞利面波经验格林函数,再用自适应时频分析方法提取群速度频散曲线,最后利用传统面波层析成像方法反演获得黑龙江地区周期为8～40 s范围内群速度分布图像.反演结果揭示黑龙江地区地壳及上地幔群速度结构存在明显横向不均...     

Subducted slabs beneath the eastern Indonesia-Tonga region: insights from tomography

Tomographic images of mantle structure beneath the region north and northeast of Australia show a number of anomalously fast regions. These are interpreted using a recent plate tectonic reconstruction in terms of current and former subduction systems. Several strong anomalies are related to current subduction. The inferred slab lengths and positions are consistent with Neogene subduction beneath the New Britain and Halmahera arcs, and at the Tonga and the New Hebrides trenches where there has been rapid rollback of subduction hinges since about 10 Ma. There are several deeper flat-lying anomalies which are not related to present subduction and we interpret them as former subduction zones overridden by Australia since 25 Ma. Beneath the Bird’s Head and Arafura Sea is an anomaly interpreted to be due to north-dipping subduction beneath the Philippines-Halmahera arc between 45 and 25 Ma. A very large anomaly extending from the Papuan peninsula to the New Hebrides, and from the Solomon Islands to the east Australian margin, is interpreted to be the remnant of south-dipping subduction beneath the Melanesian arc between 45 and 25 Ma. This interpretation implies that a flat-lying slab can survive for many tens of millions of years at the bottom of the upper mantle. In the lower mantle there is a huge anomaly beneath the Gulf of Carpentaria and east Papua New Guinea. This is located above the position where the tectonic model interprets a change in polarity of subduction from north-dipping to south-dipping between 45 and 25 Ma. We suggest this deep anomaly may be a slab subducted beneath eastern Australian during the Cretaceous, or subducted north of Australia during the Cenozoic before 45 Ma. The tomography also supports the tectonic interpretation which suggests little Neogene subduction beneath western New Guinea since no slab is imaged south of the New Guinea trench. However, one subduction zone in the tectonic model and many others, that associated with the Trobriand trough east of Papua New Guinea and the Miocene Maramuni arc, is not seen in the tomographic images and may require reconsideration of currently accepted tectonic interpretations.     

On asymptotics of low-frequency oscillations of the liquid core: 2. Frequency dependence of the inertial core-mantle coupling

New, unique information on the inertial and dissipative coupling of the liquid core and the mantle has been retrieved from modern high-precision (radiointerferometer and GPS) data on tidal variations in the rotation velocity and nutation of the Earth. Comparison of theoretical and observed data provided new estimates for the dynamic flattening of the outer liquid and the inner solid cores, mantle quality factor, viscosity of the liquid core, and electromagnetic coupling of the liquid core and the mantle [Molodensky, 2004, 2006]. As was shown in the first part of the paper [Molodensky, 2008] (further referred to as [I]), generation of eddy flows in Proudman-Taylor columns, whose orientation is controlled by the topography of the liquid core-mantle boundary, should be taken into account for correct estimation of the inertial coupling (see formulas (8) and (34) in [I]). The range of periods within which this effect plays a significant role is determined by the decay time of these flows. This time is estimated in the paper for the case where dissipation is related to viscous friction at the core-mantle boundary or with the electromagnetic coupling of the liquid core and the mantle. Because of significant uncertainties in modern data on the viscosity of the liquid core, the magnetic field intensity at the core-mantle boundary, and the electrical conductivity of the lower mantle, the dissipative coupling of the liquid core and the mantle cannot be calculated as yet. However, as shown in the paper, the decay time of eddy flows is connected with the attenuation time of subdiurnal free nutation and with the liquid core viscosity. This enables the estimation of the frequency dependence of the dissipative coupling in a fairly wide range. It is shown that the range of periods for which relations (8) and (34) in [I] are valid encompasses the best-studied length-of-day variations and, therefore, these relations are applicable to analysis of the majority of modern data.     

Probability of radial anisotropy in the deep mantle

It is well established that the Earth's uppermost mantle is anisotropic, but observations of anisotropy in the deeper mantle have been more ambiguous. Radial anisotropy, the discrepancy between Love and Rayleigh waves, was included in the top 220 km of PREM, but there is no consensus whether anisotropy is present below that depth. Fundamental mode surface waves, for commonly used periods up to 200 s, are sensitive to structure in the first few hundred kilometers and therefore do not provide information on anisotropy below. Higher mode surface waves, however, have sensitivities that extend to and below the transition zone and should thus give insight about anisotropy at greater depths, but they are very difficult to measure. We previously developed a new technique to measure higher mode surface wave phase velocities with consistent uncertainties. These data are used here to construct probability density functions of a radially anisotropic Earth model down to approximately 1500 km. In the uppermost mantle, we obtain a high probability of faster horizontally polarized shear wave speed, likely to be related to plate motion. In the asthenosphere and transition zone, however, we find a high probability of faster vertically polarized shear wave speed. To a depth of 1500 km in the lower mantle, we see no significant shear wave anisotropy. This is consistent with results from laboratory measurements which show that lower mantle minerals are anisotropic but LPO is unlikely to develop in the pressure–temperature conditions present in the mid-mantle.     

天山造山带地区瑞利面波相速度与方位各向异性

天山造山带作为世界上陆内最大的造山带之一,现今地震活动频繁,造山运动强烈,是开展陆内造山和内陆地震活动研究的天然试验场.本文利用整个天山造山带地区国内及国际台网的108个地震台站连续三年的背景噪声资料,提取了8~50 s周期的瑞利面波相速度频散曲线,并构建了整个天山造山带地区的二维瑞利面波相速度与方位各向异性分布图像.结果表明：浅部结构与地表的地质构造单元具有较大的相关性.低波速异常主要分布于沉积层厚度较大的盆地地区,而高波速异常主要分布于构造活动比较活跃的山脉地区.东天山地区中下地壳存在比较弱的低波速异常,而塔里木盆地和准噶尔盆地汇聚边缘的上地幔区域则表现为明显的高波速异常,各向异性快波方向呈现近NS向的特征,暗示着塔里木盆地和准噶尔盆地的岩石圈已经俯冲至东天山的下方.中天山地区的中下地壳至上地幔区域均呈现为明显的低波速异常,且各向异性快波方向变化比较复杂,表明中天山地区的整个岩石圈结构已经弱化,热物质上涌可能对介质的方位各向异性有一定的影响.西天山及帕米尔高原的上地幔区域存在低波速异常,各向异性表现为NW-SE方向,可能与欧亚板块的大陆岩石圈南向俯冲有关.塔里木盆地内部存在相对弱的低波速异常,推测塔里木盆地可能已经受到上涌的地幔热物质的侵蚀和破坏.     

Crustal and upper mantle structures of the brazilian coast

The Rayleigh wave phase and group velocities in the period range of 24–39 sec, obtained from two earthquakes which occurred in northeastern brazil and which were recorded by the Brazilian seismological station RDJ (Rio de Janeiro), have been used to study crustal and upper mantle structures of the Brazilian coastal region. Three crustal and upper mantle models have been tried out to explain crustal and upper mantle structures of the region. The upper crust has not been resolved, due basically to the narrow period range of the phase and group velocities data. The phase velocity inversions have exhibited good resolutions for both lower crust and upper mantle, with shear wave velocities characteristic of these regions. The group velocity data inversions for these models have showed good results only for the lower crust. The shear wave velocities of the lower crust (3.86 and 3.89 km/sec), obtained with phase velocity inversions, are similar to that (=3.89 km/sec) found byHwang (1985) to the eastern South American region, while group velocity inversions have presented shear velocity (=3.75 km/sec) similar to that (=3.78 km/sec) found byLazcano (1972) to the Brazilian shield. It was not possible to define sharply the crust-mantle transition, but an analysis of the phase and group velocity inversions results has indicated that the total thickness of the crust should be between 30 and 39 km. The crustal and upper mantle model, obtained with phase velocity inversion, can be used as a preliminary model for the Brazilian coast.