广东大亚湾周围地区地壳和上地幔结构的初步研究

在广东大亚湾周围地区,采用地震转换波法取得了该区地壳和上地幔结构的初步研究结果。在大亚湾地区的地壳中只能分出一个中间界面(界面C)和莫霍界面(M)。本区为两层地壳,地壳厚度为28—30公里。地壳和上地幔界面埋深变化不大,深断裂不发育,本区深部构造特征与地震活动性较高的华北地区相比有较大差别。     

Seismic scattering at the top of the mantle Transition Zone

We document strong seismic scattering from around the top of the mantle Transition Zone in all available high resolution explosion seismic profiles from Siberia and North America. This seismic reflectivity from around the 410 km discontinuity indicates the presence of pronounced heterogeneity in the depth interval between 320 and 450 km in the Earth’s mantle. We model the seismic observations by heterogeneity in the form of random seismic scatterers with typical scale lengths of kilometre size (10-40 km by 2-10 km) in a 100-140 km thick depth interval. The observed heterogeneity may be explained by changes in the depths to the α-β-γ spinel transformations caused by an unexpectedly high iron content at the top of the mantle Transition Zone. The phase transformation of pyroxenes into the garnet mineral majorite probably also contributes to the reflectivity, mainly below a depth of 400 km, whereas we find it unlikely that the presence of water or partial melt is the main cause of the observed strong seismic reflectivity. Subducted oceanic slabs that equilibrated at the top of the Transition Zone may also contribute to the observed reflectivity. If this is the main cause of the reflectivity, a substantial amount of young oceanic lithosphere has been subducted under Siberia and North America during their geologic evolution. Subducted slabs may have initiated metamorphic reactions in the original mantle rocks.     

Preliminary models of upper mantle P and S wave velocity structure in the western South America region

Upper mantle P and S wave velocities in the western South America region are obtained at depths of foci from an analysis of travel time data of deep earthquakes. The inferred velocity models for the Chile–Peru–Ecuador region reveal an increase of P velocity from 8.04 km/s at 40 km to 8.28 km/s at 250 km depth, while the S velocity remains almost constant at 4.62 km/s from 40 to 210 km depth. A velocity discontinuity (probably corresponding to the L discontinuity in the continental upper mantle) at 220–250 km depth for P and 200–220 km depth for S waves, with a 3–4% velocity increase, is inferred from the velocity–depth data. Below this discontinuity, P velocity increases from 8.54 km/s at 250 km to 8.62 km/s at 320 km depth and S velocity increases from 4.81 km/s at 210 km to 4.99 km/s at 290 km depth. Travel time data from deep earthquakes at depths greater than 500 km in the Bolivia–Peru region, reveal P velocities of about 9.65 km/s from 500 to 570 km depth. P velocity–depth data further reveal a velocity discontinuity, either as a sharp boundary at 570 km depth with 8–10% velocity increase or as a broad transition zone with velocity rapidly increasing from 560 to 610 km depth. P velocity increases to 10.75 km/s at 650 km depth. A comparison with the latest global average depth estimates of the 660 km discontinuity reveals that this discontinuity is at a relatively shallow depth in the study region. Further, a velocity discontinuity at about 400 km depth with a 10% velocity increase seems to be consistent with travel time observations from deep earthquakes in this region.     

Upper mantle heterogeneity in the northern part of Eurasia

A two dimensional velocity model of the upper mantle has been compiled from a long-range seismic profile crossing the West Siberian young plate and the old Siberian platform. It revealed considerable horizontal and vertical heterogeneity of the mantle. A sharp seismic boundary at a depth of 400 km outlines the high-velocity gradient transition zone, its base lying at a depth of 650 km. Several layers with different velocities, velocity gradients and wave attenuation are distinguished in the upper mantle. They likewise differ in their inner structure. For instance, the uppermost 50–70 km of the mantle are divided into blocks with velocities from 7.9–8.1 to 8.4–8.6 km s^?1.Comparison of the travel-time curves for the Siberian long-range profile with those compiled from seismological data for Europe distinguished large-scale upper mantle inhomogeneities of the Eurasian continent and allowed for the correlation of tectonic features and geophysical fields. The velocity heterogeneity of the uppermost 50–100 km of the mantle correlates with the platform age and heat flow, i.e., the young plates of Western Europe and Western Siberia have slightly lower velocities and higher heat flows than the ancient East European and Siberian platforms. At greater depths (150–250 km) the upper mantle velocities increase from the ocean to the inner parts of the continent. The structure of the transition zone differs significantly beneath Western Europe and the other parts of Eurasia. The sharp boundary at a depth of 400 km, traced throughout the whole continent as the boundary reflecting intensive waves, transforms beneath Western Europe into a gradient zone. This transition zone feature correlates with positions of the North Atlantic-west Europe geoid and heat-flow anomalies.     

Deep structure of the East European platform according to seismic data

The paper presents a review and analysis of new seismic data related to the structure of the mantle beneath the East European platform. Analysis of observations of long-range profiles revealed pronounced differences in the structure of the lower lithosphere beneath the Russian plate and the North Caspian coastal depression. The highest P-velocities found at depths around 100 km are in the range 8.4–8.5 km s^?1. Deep structure of the Baltic shield is different from the structures of both these regions. No evidence of azimuthal anisotropy in the upper mantle was found. A distribution of P-velocity in the upper mantle and in the transition zone consistent with accurate travel-time data was determined. The model involves several zones of small and large positive velocity gradients in the upper mantle, rapid increases of velocity near 400 and 640 km depths and an almost constant positive velocity gradient between the 400 and 640 km discontinuities. The depth of the 640 km discontinuity was determined from observations of waves converted from P to SV in the mantle.     

华北地区上地幔及过渡带电性结构研究

采用远参考道和Robust技术,处理了华北地区14个地磁台站资料,得到了相干度超过0.8的地磁测深响应函数.并将其转换为大地电磁测深的响应函数,获取了105~107 s周期范围内的视电阻率和相位.应用ρ+理论对数据进行了一致性检验和反演,结果表明417km,850km深度附近可能存在电性间断面.同时采用基于一维最光滑模型的Occam反演方法得到了300~1000km范围的地幔电性结构,并与前人在其他地区的研究结果进行了对比.发现华北地区地幔过渡带的电导率在大兴安岭—太行山重力梯度带东西两侧表现不同,重力梯度带附近及西侧台站下方过渡带深度的电导率和北美的Tucson地区相当,而华北地区东部的电导率在地幔过渡带范围高出西侧约2~5倍,这很可能和太平洋板块的俯冲有关.     

华北克拉通东北部及邻区地壳和地幔转换带厚度研究

本文利用宽频流动台阵记录的远震波形资料和接收函数波动方程叠后偏移方法,获得了华北克拉通东北部边界及其邻近地区的地壳和地幔转换带的间断面结构图像.结果显示研究区域的地壳厚度存在显著的横向变化:以南北重力梯度带为界,西北部的兴蒙造山带地壳较厚(~40 km),东南部的燕山带、松辽盆地和辽东台隆地壳明显较薄(30~35 km).这有可能反映,研究区南北重力梯度带两侧地壳在中-新生代区域构造伸展过程中经历了不同程度的改造和减薄.地幔转换带成像结果显示,研究区410 km和660 km间断面结构存在横向差异.经度121°E-122°E之间,上地幔底部出现双重间断面,深度分别为660 km和690 km.经度122.5°E以东(北黄海地区),410 km间断面有5~20 km幅度的下沉,660 km间断面有5~15 km幅度的抬升;该地区地幔转换带厚度相对全球平均偏薄10~20 km,指示着该地区较热的上地幔底部温度环境.我们认为太平洋俯冲板块可能停滞在研究区119°E-122°E经度范围的地幔转换带中,但未延伸至118°E以西;而俯冲板块在124°E以东可能局部穿透了上地幔底部而进入下地幔,同时引起小尺度的地幔对流,导致北黄海地区下地幔物质的上涌.     

华北地区上地幔温度及岩石圈厚度分布研究

上地幔温度是控制地幔流变性和动力学过程的关键参数之一.本文利用高分辨率S波地震层析成像速度结果,基于岩石温度与地震波速度的关系,研究了华北地区上地幔50~300 km深度范围内的温度分布和"热"岩石圈厚度.为了验证结果的可靠性,本文用计算的上地幔60 km深度处的温度作为底面约束,得到了相应的地表热流.计算地表热流与观测地表热流之间符合程度较好,相对误差大部分都在地表热流观测误差范围之内.通过对上地幔的温度分布进行分析,我们研究发现:(1)在上地幔浅部,温度与地表构造之间有很好的对应关系.在小于170 km的深度上,温度呈现东高西低的分布态势.温度较高的区域集中在东部的河淮盆地、渤海湾盆地、华北平原和中部陆块的交界处、西部鄂尔多斯高原北缘的银川―河套地堑以及阴山地区,同时,这些地区的岩石圈厚度也相应较薄,大约为80~100 km;(2)西部的鄂尔多斯高原是整个华北地区岩石圈地幔温度最低的地区,比东部地区低200~400 ℃,岩石圈厚度相应最厚,平均岩石圈厚度达到140~150 km,最厚处超过160 km.(3)在170 km以下的软流圈地幔部分,温度分布发生反转,西部温度高于东部,表明东、西部陆块在地质历史时期经历了不同的深部地幔动力学过程.     

The nature of the boundary between the upper and the lower mantle and its influence on convection

In the PREM seismic model, the boundary between the upper and the lower mantle is accepted at a depth of 670 km, where seismic velocities and density increase. However, until recently there was an obvious inconsistency in this model. The density increases abruptly, and the velocities, in addition to the jumps, have also the subsequent zones of increased gradient. The discontinuity between the upper and the lower mantle is related to the transition of olivine from the ringwoodite phase into the mixture of perovskite and magnesiowustite. However, in the pyrolyte model, the transition zone of the upper mantle consists not wholly of olivine, but partly of olivine (60%) and partly of garnet (40%). The latest data of the garnet measurement at high pressures show that it also experiences phase transition, being converted into magnesium perovskite with the impurity of calcium perovskite. In contrast to the sharp transition in olivine (within a depth interval of only 5 km), the transition in garnet is spread over the interval of depths of 660–710 km. In the widely used PREM and AK135 models, this additional transition corresponds to the zone of the increased gradient in seismic velocities, while in the density distribution it is included in the sharp transition of ringwoodite. Thus, the mineralogy data indicate the need for correction of the PREM and AK135 seismic models: the density jump at a depth of 660 km should be reduced by approximately a factor of two, and a subjacent layer with the increased density gradient should be added at the depth interval of 660–710 km. The phase transition in olivine hampers the mantle flows, although in garnet it accelerates them. Therefore, with an allowance for the smaller jump in density, the decelerating effect of the subducting plates, caused by the phase transition in olivine, decreases, and, furthermore, the effect of their acceleration, caused by the phase transition in garnet, is added. The decrease in the density jump by almost a factor of two will lead to essential changes in the results of the majority of recent works addressing the assessment of the deceleration of convection at the upper/lower mantle discontinuity on the basis of the PREM model.     

Tomographic structure of East Asia:Ⅱ.Stagnant slab above 660 km discontinuity and its geodynamic implications

P-wave arrival times of both regional and teleseismic earthquakes were inverted to obtain mantle structures of East Asia.No fast(slab) velocity anomalies was not find beneath the 660-km discontinuity through tomography besides a stagnant slab within the transition zone.Slow P-wave velocity anomalies are present at depths of 100-250 km below the active volcanic arc and East Asia.The western end of the flat stagnant slab is about 1 500 km west to active trench and may also be correlated with prominent surface topographic break in eastern China.We suggested that active mantle convection might be operating within this horizontally expanded "mantle wedge" above both the active subducting slabs and the stagnant flat slabs beneath much of the North China plain.Both the widespread Cenozoic volcanism and associated extensional basins in East Asia could be the manifestation of this vigorous upper mantle convection.Cold or thermal anomalies associated with the stagnant slabs above the 660-km discontinuity have not only caused a broad depression of the boundary due to its negative Clapeyron slope but also effectively shielded the asthenosphere and continental lithosphere above from any possible influence of mantle plumes in the lower mantle.     

利用接收函数方法研究四川地区地壳结构

采用接收函数反演和共转换点(CCP)偏移叠加成像方法, 利用四川数字地震台网宽频带的52个区域固定地震台站和布设的两条52个宽频带流动地震观测台站的远震地震波形数据资料, 对四川地区地壳结构进行研究。 结果表明, 四川地区的Moho面深度在青藏高原和四川盆地差异明显, 在川西高原地区地壳厚度为52~68 km, 在川滇地块地壳厚度为50~60 km, 在中地壳内存在不连续的低速层分布; 而在四川盆地地壳厚度为38~45 km, 地壳内没有低速层存在。 Moho面深度从川西高原的60多公里至四川盆地的约40 km, 在二者的交界处龙门山断裂带下面, 存在厚度约30 km左右宽的下降过渡带, 说明其下的Moho面可能受断层影响, 结构比较复杂; 在高原地区的上地壳界面和下地壳上界面比四川盆地的相应界面深; 高原地区在中地壳的上部有不连续的低速层分布, 在松潘—甘孜地块的上地壳下部存在向南东运动的脆性推覆体, 在羌塘—理塘地块的上地壳下部存在向南东和南运动的脆性物质流动。     

基于波形拟合的中国东海地区410 km间断面附近速度结构研究

410 km间断面是地幔转换带的顶界面,对其速度结构和起伏形态开展地震学探测有助于认识地球内部物质组成和相关的地球动力学过程.本文选取了由中国数字地震台网记录到的位于琉球俯冲区的一个中源地震P波宽频带波形资料,利用三重震相波形拟合研究了中国东海地区410 km间断面附近的精细速度结构.结果表明：中国东海地区下方410 km间断面整体表现为一尖锐的速度界面且有8~15 km的小幅抬升;该间断面之上存在52~62 km厚的低速层,其P波速度降低0.5%~1.6%;440 km深度以下存在1.0%~3.0%的P波高速异常.结合前人在该地区的层析成像结果,我们推测该高速异常体与西太平洋俯冲板片在中国东海地区地幔转换带内的滞留有关;板片内水相E分解使得转换带内水含量增加,这引发了410 km间断面的抬升;410 km间断面之上的低速层应与含水矿物脱水导致的地幔橄榄岩部分熔融有关.     

Crust and upper mantle structures beneath Northeast China from receiver function studies

P-wave and S-wave receiver function analyses have been performed along a profile consisted of 27 broadband seismic stations to image the crustal and upper mantle discontinuities across Northeast China. The results show that the average Moho depth varies from about 37 km beneath the Daxing’anling orogenic belt in the west to about 33 km beneath the Songliao Basin, and to about 35 km beneath the Changbai mountain region in the east. Our results reveal that the Moho is generally flat beneath the Daxing’anling region and a remarkable Moho offset (about 4 km) exists beneath the basin-mountain boundary, the Daxing’anling-Taihang Gravity Line. Beneath the Tanlu faults zone, which seperates the Songliao Basin and Changbai region, the Moho is uplift and the crustal thickness changes rapidly. We interpret this feature as that the Tanlu faults might deeply penetrate into the upper mantle, and facilitate the mantle upwelling along the faults during the Cenozoic era. The average depth of the lithosphere-asthenosphere boundary (LAB) is ~80 km along the profile which is thinner than an average thickness of a continental lithosphere. The LAB shows an arc-like shape in the basin, with the shallowest part approximately beneath the center of the basin. The uplift LAB beneath the basin might be related to the extensive lithospheric stretching in the Mesozoic. In the mantle transition zone, a structurally complicated 660 km discontinuity with a maximum 35 km depression beneath the Changbai region is observed. The 35 km depression is roughly coincident with the location of the stagnant western pacific slab on top of the 660 km discontinuity revealed by the recent P wave tomography.     

中国西部及邻区岩石圈S波速度结构面波层析成像

本文利用瑞利波群速度频散资料和层析成像方法,研究了中国西部及邻近区域(20°N—55°N,65°E—110°E)的岩石圈S波速度结构.结果表明这一地区存在三个以低速地壳/上地幔为特征的构造活动区域:西蒙古高原—贝加尔地区,青藏高原,印支地区.西蒙古高原岩石圈厚度约为80 km,上地幔低速层向下延伸至300 km深度,说明存在源自地幔深部的热流活动.缅甸弧后的上地幔低速层下至200 km深度,显然与印度板块向东俯冲引起俯冲板片上方的热/化学活动有关.青藏高原地壳厚达70 km,边缘地区厚度也在50 km以上并且具有很大的水平变化梯度,与高原平顶陡边的地形特征一致.中下地壳的平均S波速度明显低于正常大陆地壳,在中地壳20~40 km深度范围广泛存在速度逆转的低速层,这一低速层的展布范围与高原的范围相符.这些特征说明青藏高原中下地壳的变形是在印度板块的北向挤压下发生塑性增厚和侧向流动.地幔的速度结构呈现与地壳显著不同的特点.在高原主体和川滇西部地区上地幔顶部存在较大范围的低速,低速区范围随深度迅速减小;100 km以下滇西低速消失,150 km以下基本完全消失.青藏高原上地幔速度结构沿东西方向表现出显著的分段变化.在大约84°E以西的喀喇昆仑—帕米尔—兴都库什地区,印度板块的北向和亚洲板块的南向俯冲造成上地幔显著高速;84°E—94°E之间上地幔顶部速度较低,在大约150~220 km深度范围存在高速板片,有可能是俯冲的印度岩石圈,其前缘到达昆仑—巴颜喀拉之下;在喜马拉雅东构造结以北区域,存在显著的上地幔高速区,可能阻碍上地幔物质的东向运动.川滇西部岩石圈底界深度与扬子克拉通相似,约为180 km,但上地幔顶部速度较低.这些现象表明青藏高原岩石圈地幔的变形/运动方式可能与地壳有本质的区别.     

应用层析成像技术研究华北地壳速度结构

利用华北地区１６９个地震台站记录的１９８４￣１９９５年的１６９６个地震的Ｐ波到时资料,重建了华北地区地壳三维速度图像,不同深度的速度图像表明：地壳中地震波速度的横向不均匀性显著,上地壳的速度分布与地有已知构造特征相联系;太行山山前断裂带是切割到地壳底部的深断裂带,２５ｋｍ深度普遍存在速度间断面;地壳底部起伏明显;强震孕育发生与其周围介质的不均匀性和土地幔的隆起或凹陷密切相关。     

中国东部地区的壳-幔过渡带结构

莫霍面是地壳和上地幔的边界,但莫霍面并不是一个简单的"面",而是一个反映地壳和地幔物质交换、相互作用等动力学意义的"过渡带".本文综合深地震反射、宽角地震折射和高温高压岩石物理实验结果,确定壳-幔过渡带的地震P波速度变化范围为6.8~7.5 km·s^-1.在克拉通等构造活动稳定地区壳-幔过渡带内的速度梯度强且壳-幔过渡带厚度薄,而在造山带等构造活动区域壳-幔过渡带内的速度梯度弱且壳-幔过渡带厚.中国东部地区的壳-幔过渡带的平均厚度约为5~10 km,在四川盆地下方最薄（<5 km）,而在华北克拉通中部造山带下方的壳-幔过渡带最厚（~30 km）.综合地球化学结果,华北中部巨厚壳-幔过渡带主要是幔源岩浆的底侵作用和堆晶作用而形成.     

鲁西隆起重磁异常特征及其构造活动性分析

通过对鲁西隆起区重磁资料的分析和反演计算,研究了沂沭断裂带、齐河—广饶断裂带、聊城—兰考断裂带、丰沛断裂带以及地块内部断裂的重磁异常、莫霍面和居里面深度特征,并讨论了鲁西隆起的地质构造特征和构造活动性.结果显示:鲁西隆起基底广泛出露,沉积层主要分布在由断裂下降盘控制的凹陷内,区内断裂深度达20 km以上,其中蒙山断裂深入至上地幔,控制了蒙山金伯利岩型金刚石矿的产出;鲁西隆起区莫霍面深度为30—35 km,整体呈向西开口的箕形,地块中部地壳厚度较厚,除西侧地壳呈阶梯状增厚外隆起地块四周地壳逐渐减薄;居里面深度介于20—33 km之间,中部地区较深,为整体稳定的地块,断裂带分布位置对应于居里面梯度带;地震活动主要集中于断裂带与莫霍面梯度带交会区以及断裂带上的居里面突变区.      

Phase transformations in the earth's mantle and subducting slabs: Implications for their compositions, seismic velocity and density structures and dynamics

Abstract Phase transformations in model mantle compositions and those in subducting slabs have been reviewed to a depth of 800 km on the basis of recent high-pressure experimental data. Seismic velocity and density profiles in these compositions have also been calculated using these and other mineral physics data. The nature of the seismic velocity and density profiles calculated for a pyrolite composition was found to generally agree with those determined by seismic observations (e.g. PREM). The locations of the seismic discontinuities at 400 and 670 km correspond almost exactly to the depths where the transformations of the olivine component to denser phases take place. Moreover, the steep gradients in the seismic velocity/density profiles observed between these depths are qualitatively consistent with those expected from the successive transformations in the complementary pyroxene-garnet component in the pyrolite composition. Further, the calculated seismic velocity and density values agree well with those observed in the upper mantle and mantle transition region within the uncertainties attached to these calculations and observations. Pyrolite or peridotite compositions are thus most likely to represent the composition of the mantle above 670 km depth, although some degrees of chemical heterogeneity may exist in the transition region. The observed sharp discontinuous increases of seismic velocities and density at this depth may be attributed either to the phase transformation to a perovskite-bearing assemblage in pyrolite or to chemical composition changes. Density profiles in subducted slabs have been calculated along adequate geotherms assuming that the slabs are composed of the former oceanic crust underlain by a thicker harzburgitic layer. It is shown that the former oceanic crust is substantially less dense than the surrounding pyrolite mantle at depths below 670 km, while it is denser than pyrolite in the upper mantle and the transition region. The subducted former oceanic crust may be trapped in this region, forming a geochemically enriched layer at the upper mantle-lower mantle boundary. Thick and cool slabs may penetrate into the lower mantle, but the chemically derived buoyancy may result in strong deformation and formation of megalith structures around the 670 km seismic discontinuity. These structures are consistent with those detected by recent seismic tomography studies for subduction zones.     

中国东部随县─启东地带上地幔结构研究

随县─启东人工地震探测剖面长７００ｋｍ，穿越郯城─庐江深大断裂带．提取了该区莫霍界面以下上地幔中的反射界面信息，获得了上地幔的分层结构．结果表明，该区上地幔中有两个深度分别为７６ｋｍ和１００ｋｍ的界面．７６ｋｍ深度以下为一低速层，１００ｋｍ深度以下则为速度更低的低速层．该区岩石层厚度为７６ｋｍ，与华北地区相似，具有薄岩石层性质．     

中国东部随县─启东地带上地幔结构研究

随县─启东人工地震探测剖面长７００ｋｍ，穿越郯城─庐江深大断裂带．提取了该区莫霍界面以下上地幔中的反射界面信息，获得了上地幔的分层结构．结果表明，该区上地幔中有两个深度分别为７６ｋｍ和１００ｋｍ的界面．７６ｋｍ深度以下为一低速层，１００ｋｍ深度以下则为速度更低的低速层．该区岩石层厚度为７６ｋｍ，与华北地区相似，具有薄岩石层性质．