俯冲板块断离的深度

针对Davies等对俯冲板块断离判据的不足,提出了新的俯冲板块断离判据.该判据从大陆板块减速俯冲模型出发,计算了不同初速度和俯冲角度下大陆俯冲深度的变化情况.在此基础上,详细考查了影响俯冲板块热结构和强度结构的因素,计算了非热平衡状态下俯冲板块内部大陆板块与海洋板块过渡区"积分强度"的上下限,认为当俯冲入地幔的较轻陆壳在板块过渡区引发的拉张力位于这个"积分强度"上下限之间时,便可能发生俯冲板块断离.计算结果表明,一般情况下,俯冲板块断离的深度在85-135km之间,断离过程发生在大陆俯冲了1.3-4.0Ma之后.本文计算的俯冲板块断离深度可供超高压变质岩形成和折返的研究参考.     

上地幔俯冲板块的动力学过程:数值模拟

大洋板块俯冲到地幔转换带,进而可形成不同的形态：板块可以停滞在660 km不连续面,抑或穿过地幔转换带进入下地幔.这些不同的俯冲模式可进一步影响到海沟的运动.为更好地理解上地幔中俯冲板片的变形行为以及俯冲过程与海沟运动之间的关系,本文通过建立一系列高精度二维热-力学自由俯冲的数值模型,揭示了俯冲板块在上地幔中的变形方式及其与地幔转换带之间的相互作用过程.模拟结果显示,在俯冲板块与地幔转换带的相互作用过程中,其动力学过程可以分为以海沟后撤主导、海沟前进主导以及稳定型海沟等三种主要动力学类型.对于年龄较老,厚度较大的俯冲板块容易形成海沟后撤型俯冲,俯冲板块停滞在660 km不连续面.相反,年龄较小,塑性强度较小的板块容易形成海沟前进型俯冲,俯冲板块穿越660 km不连续面.

     

全球板块运动及其对地球表面积变化的影响

利用几个典型全球板块相对运动模型的欧拉矢量,计算沿这些模型板块边界段的两两板块之间的相对运动速度,据此进一步计算一定时间尺度内这些板块边界段的面积变化. 并依据已有对相应各边界段性质-包括张性、压性及走滑性质的研究,设定组成各边界段的两个板块如何分别承担其面积变化的情况,得到各模型中每个板块的面积变化和全球表面积变化总量. 结果显示,在近100万年时间尺度内,地球南半球表面积增大、北半球表面积减小,地球总表面积变化为2.9万～3.6万km~2. 如果假设地球表面积的增量完全由体积膨胀引起,这将对应地球半径160～250 m的增加. 以NUVEL-1模型为例,全球14个板块中,非洲、 南美与南极三个板块的面积持续增大,其余11个板块的面积在逐步减小;不过,非洲、南美与南极这三个板块面积增大部分之和,能够完全抵消掉其余11板块面积之减小,而且使全球总体表面积增加. 本文进一步讨论了一些其他因素,例如俯冲带前进与后退等影响地球表面积变化的情况.     

斜俯冲板块边界变形分配的力学分析

本文研究斜俯冲板块边界由俯冲的逆掩断层和平行于火山弧的走滑断层所切割出的岩石层窄条的宏观力学状态．分析计算表明,当断层强度小而构造压应力大时,海沟处板块俯冲形成的逆掩断层地震滑动方向与斜收敛板块运动方向将不一致,震源机制反映的主压应力和断层错动方向变得近乎垂直于海沟;板块间的强耦合和大陆板块边缘软弱带的存在为窄条滑动提供了有利条件．给定俯冲倾角和震源机制滑动角的观测值,可以计算逆掩和走滑断层上的有效摩擦系数．计算结果表明,两种断层上的有效摩擦系数一般都不大于０．２,反映或存在着高孔隙水压,或存在断层泥之类弱物质,或两者兼而有之．     

地球物理资料所揭示的南海东北部中生代俯冲增生带

根据新处理的重、磁和广角地震图件的解释, 推测在南海北部从台西南盆地到深海盆北缘存在一条大致NE45°走向的中生代俯冲增生带. 主要依据包括: 台西南-中沙东布格重力异常总梯度峰值带在强度和规模上都与马尼拉海沟俯冲增生带引起的总梯度峰值带相近; 该峰值带与海底地形和新生代构造都斜交, 指示前新生代构造; 该峰值带被NW向断层左行错断成雁列状, 符合中生代区域应力场特征; 其西北方在陆架区有与之大致平行的高磁异常带, 指示中生代火山弧. 此外, 海底地震仪资料显示在俯冲增生带位置上出现北倾陡坎、海底地震仪和长电缆反射地震剖面都显示相应部位出现双莫霍面叠置, 也可作为佐证. 这段中生代俯冲增生带的发现正好填补了欧亚大陆东南缘晚中生代俯冲增生带在南海东北部的一段空白.     

基于随温度变化的热系数模拟板块俯冲动力学过程

热传导系数和热膨胀系数是影响板块俯冲动力学过程的两个重要参数. 由于地球介质的不均匀性,热系数也会随深度发生变化.然而,这种变化在地球动力学模拟研究中往往被忽略.本文针对随温度变化的热传导系数和热膨胀系数, 模拟板块俯冲的动力学过程,分析热系数、黏度对板块俯冲形态的影响及其对应的地幔对流特征.结果表明,依温度变化的热传导系数和热膨胀系数会影响地幔温度及黏度分布,进而改变板块的俯冲角度;黏度是控制板块俯冲动力学演化过程的重要因素;地幔对流受黏度结构的影响,呈现分层对流及局部多个对流环等多种不同形态的对流场特征.     

斜俯冲板块边界变形分配的力学分析

本文研究斜俯冲板块边界由俯冲的逆掩断层和平行于火山弧的走滑断层所切割出的岩石层窄条的宏观力学状态．分析计算表明，当断层强度小而构造压应力大时，海沟处板块俯冲形成的逆掩断层地震滑动方向与斜收敛板块运动方向将不一致，震源机制反映的主压应力和断层错动方向变得近乎垂直于海沟；板块间的强耦合和大陆板块边缘软弱带的存在为窄条滑动提供了有利条件．给定俯冲倾角和震源机制滑动角的观测值，可以计算逆掩和走滑断层上的有效摩擦系数．计算结果表明，两种断层上的有效摩擦系数一般都不大于０．２，反映或存在着高孔隙水压，或存在断层泥之类弱物质，或两者兼而有之．     

海山的倾斜俯冲对上覆板块变形的影响

伴随洋壳的俯冲,驼伏其上的海山会导致上覆板块的强烈变形.为解释该构造变形特征,本文运用物理模拟实验的方法,着重分析海山的斜向俯冲对上覆板块变形的影响,并将模拟结果与正向俯冲过程进行对比.实验结果显示:海山开始进入俯冲,前缘楔体的增生会被阻止,同时楔体被抬升并出现脱顶构造,未被海山破坏的楔体会出现后冲断层的激活,后冲断层轴平行于海山的俯冲方向.海山进一步俯冲,突起项部发育一系列张扭性质的微断裂和走滑性质的共轭断裂,尾随突起之后的楔体由于重力会产生正断层系统.比起正向俯冲,斜向俯冲过程中所产生的后逆冲体、海山两侧的叠瓦状逆冲推覆构造都出现不对称分布,断裂和微断裂束的走向不规则散开,后冲断层的轴向及海山俯冲过后在楔体上产生的凹槽的轨迹都不断斜向迁移,且凹槽两侧的地势不一致等.最后利用文中的物理模拟结果,很好的解释了马尼拉海沟中段俯冲构造的构造特征,同时对其他俯冲大陆边缘的构造解释具有指导意义.     

西太平洋板块俯冲运动与中国东北深震带

中国东北地区有一条震源深度达600公里的深震带,从震源剖面投影可以明显地看出,它是西太平洋板块俯冲的结果,俯冲的角度约为26°。 根据板块学说,建立了一个板块俯冲模型,主要考虑地幔物质和岩石层板块的热传导,计算了板块俯冲到一定深度时的温度分布。结果表明,西太平洋板块俯冲到中国东北地区深度达600公里时,其中心温度约为1200℃,仍比周围地幔物质的温度低得多,因而能产生弹性断裂,发生地震。 震源机制的结果表明,中国东北地区深震的主压应力轴方位为93°—113°,正好迎着板块俯冲的方向,仰角为27°—28°,与板块俯冲的角度大体一致。以上结果说明中国东北深震是西太平洋板块俯冲到中国东北大陆之下造成的。     

吉林-日本区深震特征及板块俯冲图像分析

根据USCS地震资料：分析了吉林-日本深震区地震的深度特征．结果表明：吉林-日本深震区地震位于太平洋板块与欧亚大陆板块的交界带上．其北支地震活动强于南支．同时也揭示了该区地震深度、地震频次以及地震能量之间的关系．指出其最容易发生强震的地震层位在570-580km左右。由深度资料推断板块间的作用方式，太平洋板块向欧亚大陆的挤压是一个由浅入深的过程，在不同的部位其俯冲角度与俯冲距离是不一样的，一般来说，倾角越大，俯冲越陡．水平延伸距离越短；反之，倾角越小．俯冲越平缓，水平延伸距离越长。     

三维板块几何形态对大陆深俯冲动力学的制约

大陆深俯冲及超高压变质作用是大陆动力学的重要研究内容,前人进行了系统的地质、地球物理观测以及数值模拟研究.然而,自然界中大陆板块的俯冲、碰撞及造山过程大部分具有明显的沿走向的差异性,这种典型的三维特征可能很大程度上依赖于会聚大陆板块的初始几何学和运动学特征.本文采用三维高分辨率的动力学数值模拟方法,建立了方形大陆板块和楔形大陆板块两种不同的俯冲-碰撞模型,并且俯冲大陆板块侧面与大洋俯冲带相邻.数值模拟结果揭示大洋板块可以持续地俯冲到地幔之中,而大陆板块俯冲到一定深度处,其前端的俯冲板块将发生断离,并进而造成残余的大陆板块俯冲角度的减小.方形大陆俯冲板块的断离深度约为150km,而楔形大陆俯冲板块的断离深度较大,约250~300km,这很大程度上取决于俯冲带中大洋板块的牵引力和大陆板块的负浮力之间的竞争关系.同时,无论方形还是楔形大陆板块俯冲模型中,板块断离后,侧向的大洋俯冲板块仍可以拖曳约60~70km宽的大陆边缘岩石圈持续向下俯冲,揭示了新西兰东部的洋-陆空间转换俯冲带的动力学机制.并且,数值模型与喜马拉雅造山带和秦岭—大别—苏鲁造山带进行了对比,进而对其高压-超高压岩石空间展布沿走向的差异性特征和机制提供了一定的启示.     

Character of sediments entering the Costa Rica subduction zone: Implications for partitioning of water along the plate interface

Abstract   Sediments deposited off the Nicoya Peninsula advect large volumes of water as they enter the Costa Rica subduction zone. Seismic reflection data, together with results from Ocean Drilling Program Leg 170, show that hemipelagic mud comprises the upper ∼135 m of the sediment column (ranging from 0 to 210 m). The lower ∼215 m of the sediment column (ranging from 0 to 470 m) is pelagic carbonate ooze. We analyzed samples from 60 shallow (<7 m) cores to characterize the spatial variability of sediment composition on the incoming Cocos Plate. The bulk hemipelagic sediment is 10 wt% opal and 60 wt% smectite on average, with no significant variations along strike; the pelagic chalk contains approximately 2 wt% opal and <1 wt% smectite. Initially, most of the water (96%) in the subducting sediment is stored in pore spaces, but the pore water is expelled during the early stages of subduction by compaction and tectonic consolidation. Approximately 3.6% of the sediment's total water volume enters the subduction zone as interlayer water in smectite; only 0.4% of the water is bound in opal. Once subducting strata reach depths greater than 6 km (more than 30 km inboard of the subduction front), porosity drops to less than 15%, and temperature rises to greater than 60°C. Under those conditions, discrete pulses of opal and smectite dehydration should create local compartments of fluid overpressure, which probably influence fluid flow patterns and reduce effective stress along the plate boundary fault.     

Helium and carbon isotopes in fluorites: implications for mantle carbon contribution in an ancient subduction zone

The concentrations of helium and carbon in fluorite associated with Cretaceous to Neogene (90–13 Ma) granitic magmatism in the Japanese arc have been measured. Concentrations of Li, U, Th and Gd were measured to correct for secondary generated ³He. The CO₂/³He of fluorites are almost uniform (1.5×10¹⁰–4×10¹⁰) and in fair agreement with the range of present island arc volcanic gases. The calculated mantle C contribution in the Mesozoic subduction zone appear to have been identical to the present one (7–19%) indicating that the C flux from the mantle in supra-subduction zone environments has remained fairly constant during the past 70 million years.     

A long-term rock uplift rate for eastern Crete and geodynamic implications for the Hellenic subduction zone

The island of Crete in the forearc of the Hellenic subduction zone has a rugged topography with local relief exceeding 2 km. Based on the elevation of marine shorelines, rates of rock uplift during the Late Holocene were previously estimated to range between 1 and 4 mm/a in different parts of the island. These rates may, however, not be representative for longer timescales, because subduction earthquakes with up to 9 m of vertical coseismic displacement have affected Crete in the Late Holocene. Here we use a well preserved sequence of marine terraces near Kato Zakros in eastern Crete to determine the rate of rock uplift over the last ∼600 ka. Field investigations and topographic profiles document a flight of more than 13 marine bedrock terraces that were carved into limestones of the Tripolitza unit. Preliminary age constraints for the terraces were obtained by ¹⁰Be exposure dating of rare quartz-bearing sandstone clasts, which are present on some terraces. The ¹⁰Be ages of these samples, which have been corrected for an inherited nuclide component, yielded exposure ages between ∼100 ka and zero. Combined with geomorphologic evidence the two oldest ¹⁰Be ages suggest that the terraces T₄ and T₅, with shoreline angles at an elevation of ∼68 and ∼76 m above sea level, respectively, formed during the marine isotope stage 5e about 120 ka ago. The correlation of the higher terraces (T₆ to T₁₃) with regional sea-level highstands indicates sustained rock uplift at a rate of ∼0.5 m/ka since at least ∼600 ka. As normal faulting has dominated the tectonics of Crete during the last several million years, upper crustal shortening can be ruled out as a cause for rock uplift. We argue that the sustained uplift of the island results from the continuous underplating of sediments, which are transferred from the subducting African plate to the base of the crust beneath Crete.     

The seismogenic zone of subduction thrust faults

Abstract Subduction thrust faults generate earthquakes over a limited depth range. They are aseismic in their seaward updip portions and landward downdip of a critical point. The seaward shallow aseismic zone, commonly beneath accreted sediments, may be a consequence of unconsolidated sediments, especially stable-sliding smectite clays. Such clays are dehydrated and the fault may become seismogenic where the temperature reaches 100--150°C, that is, at a 5--15 km depth. Two factors may determine the downdip seismogenic limit. For subduction of young hot oceanic lithosphere beneath large accretionary sedimentary prisms and beneath continental crust, the transition to aseismic stable sliding is temperature controlled. The maximum temperature for seismic behavior in crustal rocks is ～ 350°C, regardless of the presence of water. In addition, great earthquake ruptures initiated at less than this temperature may propagate with decreasing slip to where the temperature is ～ 450°C. For subduction beneath thin island arc crust and beneath continental crust in some areas, the forearc mantle is reached by the thrust shallower than the 350°C temperature. The forearc upper mantle probably is aseismic because of stable-sliding serpentinite hydrated by water from the underthrusting oceanic crust and sediments. For many subduction zones the downdip seismogenic width defined by these limits is much less than previously assumed. Within the narrowly defined seismic zone, most of the convergence may occur in earthquakes. Numerical thermal models have been employed to estimate temperatures on the subduction thrust planes of four continental subduction zones. For Cascadia and Southwest Japan where very young and hot plates are subducting, the downdip seismogenic limit on the subduction thrust is thermally controlled and is shallow. For Alaska and most of Chile, the forearc mantle is reached before the critical temperature, and mantle serpentinite provides the limit. In all four regions, the seismogenic zones so defined agree with estimates of the extent of great earthquake rupture, and with the downdip extent of the interseismic locked zone.     

Receiver functions in northeast China - implications for slab penetration into the lower mantle in northwest Pacific subduction zone

Seismic studies of the subducting lithosphere and the upper mantle discontinuities in the northwest Pacific subduction zone beneath Japan and northeast China have suggested contrary subduction scenarios. There was little consensus on the issue whether the subducting slab penetrates the upper mantle discontinuities into the lower mantle or it is deflected atop of the 660-km discontinuity over several hundred kilometers. We calculate receiver functions from a recent seismic broadband station network located in northeast China and find topographic variations of the upper mantle discontinuities. A deeper-than-normal 660-km discontinuity is observed over an area of 400 km and it coincides with the stagnant slab imaged by seismic tomography. The 660-km discontinuity is locally depressed by more than 35 km and the transition zone is thickened by more than 20 km in the east of the region where it encounters the slab. These observations provide evidence of the slab accumulating in the mantle transition zone and locally penetrating into the lower mantle.     

西北太平洋俯冲地区410-km间断面上覆低速层探测

自20世纪90年代首次探测到410-km间断面上覆低速层以来,全球多个俯冲带和大陆克拉通地区都陆续发现了该低速层结构.对其特性及形成机理的探讨是深部地幔结构、物性和动力学研究的热点问题.本文聚焦于西北太平洋俯冲地区410-km间断面上覆低速层的探测及特性研究上.通过对发生于日本北海道地区两个中等深度地震区域波形资料的分析,利用三重震相波形拟合方法获得了我国东北及日本海西北部下方410-km间断面附近的P波速度结构.速度模型明确显示,410-km间断面上方存在厚~47±14 km,异常值~2%的低速层,横向展布近700 km.结合区域地震层析成像、矿物高温高压物理实验及动力学模拟结果,我们否定了"从下至上"的上涌热物质导致410-km间断面上覆低速层的模型;认为较老且快速俯冲的太平洋板块在地幔过渡带顶部脱水导致硅酸盐矿物的部分熔融,由于熔体密度较大能够稳定存在于410-km间断面之上,从而产生了观测到的横向展布较广的410-km间断面上覆低速层结构.     

马尼拉俯冲带北段增生楔前缘构造变形和精细结构

马尼拉俯冲带是南海的东部边界,记录了南海形成演化的关键信息,同时也是地震和海啸多发区域.本文利用过马尼拉俯冲带北段的高分辨率多道地震剖面,分析了研究区内海盆和海沟的沉积特征,精细刻画了区内增生楔前缘的构造变形、结构以及岩浆活动特征.研究区内增生楔下陆坡部分由盲冲断层、构造楔和叠瓦逆冲断层构成,逆冲断层归并于一条位于下中新统的滑脱面上,滑脱面向海方向的展布明显受到增生楔之下埋藏海山和基底隆起的影响;上陆坡的反射特征则因变形强烈和岩浆作用而难以识别;岩浆活动开始于晚中新世末期并持续至第四纪.马尼拉俯冲带北段增生楔的形成时间早于16.5 Ma,并通过前展式逆冲向南海方向扩展;马尼拉俯冲带的初始形成时间可能在晚渐新世,而此时南海海盆扩张仍在持续.南海东北缘19°N-21°N区域为南海北部陆坡向海盆的延伸,高度减薄的陆壳的俯冲造成马尼拉海沟北段几何形态明显地向东凹进.

     

阿留申—阿拉斯加俯冲带及周边地区地幔过渡带结构研究

利用国家测震台网固定台站和"中国地震科学台阵探测"项目在南北地震带北段布设的宽频带流动台阵记录到的极远震事件,通过SS前驱波震相研究,获得了阿留申—阿拉斯加俯冲带东段及邻区下方410 km和660 km间断面的埋深和起伏形态特征.为增强对SS前驱波震相的识别,我们采用了时差校正和共反射点叠加分析.叠加结果显示,毗邻阿留申俯冲带的白令海、阿拉斯加半岛、以及阿拉斯加中南部和东部地区下方,410 km和660 km间断面的埋深基本呈正相关关系,因而具有正常的过渡带厚度.这表明在阿留申—阿拉斯加俯冲带东段,北太平洋板块还没有俯冲到地幔过渡带深度范围内.其次,在阿拉斯加西部地区下方,660 km间断面出现明显下沉,而上覆的410 km间断面埋深接近于全球平均值,从而导致过渡带明显加厚.据此,我们推测在阿拉斯加西部地区下方地幔过渡带底部可能存在库拉残留板块.