横向黏度变化的全地幔对流应力场初步研究

将地幔地震波速度异常转换为地幔横向黏度变化(达到3个数量级),在球坐标系下计算了瑞雷数为106、上边界为刚性、下边界为应力自由等温边界条件下的岩石层底部的地幔对流极型和环型应力场.结果表明,地幔对流极型应力场与地表大尺度构造具有良好的对应关系:俯冲带和碰撞带的应力呈现挤压状态,而洋中脊处的应力则呈现拉张状态.地幔对流环...     

基于震源机制解的鄂霍次克微板块东部俯冲带地区构造应力场反演

基于国际地震中心(ISC)提供的1970年1月～2016年12月期间的地震震源机制解,对鄂霍次克微板块东部俯冲带地区进行了应力张量反演,得到了日本海沟、千岛海沟和勘察加海沟3个俯冲带区域的构造应力场特征。研究结果显示:①海沟地区浅部区域(h100km)的水平主压应力轴与西北太平洋板块的俯冲方位一致,与海沟走向近似垂直,其洋壳一侧以拉张型应力状态为主,而陆壳一侧则以挤压型应力为主,且在弧后区域均存在拉张的应力状态;80～200km深度范围区域表现出双地震带"Ⅰ"型构造应力场特征。②日本海沟带由于俯冲角相对较小(相比于千岛海沟和勘察加海沟),水平方向沿NWW向延伸更远,大洋板块与上覆板块之间耦合更加强烈,逆冲型地震发生数量最多。③对于深部区域(h300km),千岛地区应力场表现出非均匀性特征,可能是由地幔阻力导致的;而勘察加地区应力场表现出拉张型,可能是因为俯冲板片的拉伸拖曳作用更强。     

克拉通岩石圈减薄与破坏机制的动力学数值模拟

大陆克拉通岩石圈的减薄是一种常见现象,但是其破坏并不常见.尽管前人针对这两种过程开展了不少研究,但是对其动力学机制尚缺乏统一认识.文章以大陆克拉通岩石圈,尤其是华北克拉通作为研究对象,基于前人的地质地球物理学观测以及动力学数值模拟研究,对克拉通岩石圈减薄与破坏的过程和机制进行了系统的分析和总结.其中,华北克拉通减薄和破坏的动力学机制限定为与大洋板块俯冲相关的两种主要模式,分别是:(1)"自下而上"过程,即俯冲板片扰动富水的地幔转换带,使其中流体上涌,软流圈物质受到扰动而发生对流,对上覆岩石圈底部进行"烘烤",并可能伴随流体/熔体-岩石反应,改变岩石圈地幔的矿物组成和流变学性质,导致岩石圈底部物质混合进入软流圈而发生减薄乃至破坏.(2)"自上而下"过程,强调大洋板片在克拉通岩石圈底部平俯冲并发生脱水,水化蚀变上覆岩石圈地幔,降低岩石圈地幔岩石流变特性;随后平俯冲转变为陡俯冲,伴随软流圈侧向上涌,被弱化的岩石圈地幔由于热侵蚀作用而进入软流圈,从而发生减薄乃至破坏.这两大类过程都是比较复杂的、多过程耦合的模式,都把大洋板块的俯冲及其流体活动作为上覆克拉通岩石圈减薄/破坏的主要驱动机制;不同之处在于俯冲大洋板片的位置,一种是在地幔转换带的平卧与滞留,另一种是沿着克拉通岩石圈底部的平俯冲.对于华北克拉通而言,东部岩石圈在中生代时期受到古太平洋板片的俯冲,板块俯冲带对华北克拉通之下的地幔和转换带物质的扰动及其伴随的流体活动,很大程度上控制着华北克拉通岩石圈的改造.为了更好的探讨和对比这两种模式,其核心关键点在于大洋俯冲带中水的活动.因此,文章重点对含水矿物相变及俯冲带和地幔转换带中水的迁移过程进行了探讨,分析了大洋俯冲带中可能的流体活动模型、各自存在的问题及其可行性.此外,为了全面认识克拉通岩石圈的减薄与破坏机制,还对其他可能的机制进行了总结,包含以下四组模式:拉张减薄模型、挤压增厚拆沉模型、大尺度地幔对流侵蚀模型以及地幔柱侵蚀模型.相对于前面介绍的两种大洋俯冲相关模型,这四种模型主要是由相对简单的单一过程和机制所致(分别是拉张、挤压、对流和地幔柱),可能是华北克拉通减薄与破坏的二级驱动机制.     

大地电磁数据揭示中国东北地区牡丹江缝合带俯冲极性及活化特征

佳木斯地块和松嫩地块位于中亚造山带最东端,两者的汇聚受控于二叠纪至侏罗纪期间牡丹江洋的俯冲及闭合作用.牡丹江缝合带后期又受到古太平洋板块和太平洋板块俯冲作用的直接影响,因此是研究古缝合带俯冲极性和后期改造作用的理想地区.本文对跨越牡丹江缝合带的160km长大地电磁剖面数据进行三维反演,建立了缝合带及邻区的电阻率模型.结果揭示了牡丹江洋的俯冲极性和俯冲痕迹,并提供了(古)太平洋板块俯冲作用导致牡丹江缝合带活化的地电学证据.缝合带表现为复杂的高导结构,位于松嫩地块-佳木斯地块碰撞带下方的地壳尺度西倾高导体,代表牡丹江洋闭合的古俯冲带,表征了牡丹江洋西向俯冲的极性.同时,由地表地质定义的牡丹江断裂不能完全代表牡丹江缝合带的深部结构,对于该缝合带的定义应扩展至横向延伸约70km的整个高导区域.俯冲前缘弧下的固体导电矿物由深部被抬升至上地壳,其下与地幔连通的“烟囱”状高导结构代表了幔源物质的上侵通道,表明牡丹江缝合带受古太平洋板块和太平洋板块的俯冲作用影响发生活化,导致通道内原本冷却、结晶的物质发生重熔.因此,(古)太平洋板块的俯冲作用破坏了中亚造山带东段古碰撞带的岩石圈结构,其下方大尺度地壳高...     

西太平洋板片俯冲与后撤引起华北东部地幔置换并导致陆内盆-山耦合

华北中生代构造-岩浆活动频繁,深部岩石圈地幔性质发生变化,即克拉通发生活化作用.活化作用大致可分为三个阶段:(1)晚古生代至早侏罗世(至~170Ma),(2)中侏罗世至早白垩世早期(160~140Ma),(3)早白垩世至新生代(~140Ma以来).其中后两个阶段与古太平洋板片俯冲及后撤导致华北东部深部的岩石圈地幔置换并引起陆内浅部的盆山耦合过程是本文讨论的重点.在第一阶段,古亚洲洋俯冲和关闭引起华北北缘经历弧后拉张、碰撞挤压及碰撞后伸展等构造-岩浆活动,而且造成陆块边缘完整性的机械破坏和地幔性质的化学改造,成为后续软流圈物质上涌的通道和岩浆活动的优先发生区;受华南陆块俯冲的影响,华北南缘也发生了类似的过程.在第二阶段,蒙古鄂霍次克洋闭合及古太平洋板片俯冲剪切,引起华北北缘的两次近S-N向的挤压作用(燕山运动的A、B幕),近E-W向分布的陆缘盆地被晚中生代岩体和NE-SW断裂肢解为零星分布的盆岭省,岩浆作用由东北角向西迁移进入地块内部,同时郯庐断裂的性质由左行走滑转换为正断层,华北由早期的近S-N向的压扭性背景进入NW-SE向的弧后拉张阶段.第三阶段是华北克拉通破坏和岩石圈地幔增生的关键时期,深部难熔的克拉通型地幔被饱满的大洋型地幔置换,实现岩石圈大幅度减薄后的小幅增生增厚过程;浅部的表现是岩浆作用持续向东南迁移,陆内岩石圈薄弱带优先发生伸展变形,包括在早白垩世(140~110Ma)中部带侏罗纪逆冲断层反转为正断层、郯庐断裂的持续拉张引起中地壳拆离和大渤海湾盆地的沉降;晚白垩世至今(110Ma~),中部山带发生断陷作用形成汾渭盆地和沁水盆地,大渤海盆地内部断陷形成盆-山相间的地貌特征,苏鲁造山带则发育莱阳盆地等.华北克拉通规模小并发育陆内薄弱带,是克拉通容易破坏的内因.具这种特性的克拉通容易受周边多个俯冲构造域和上涌软流圈物质的共同影响.晚中生代(~160Ma)以来,华北克拉通破坏主要表现为周边块体的俯冲导致软流圈物质上涌、岩石圈减薄和浅部地壳滑脱,岩石圈薄弱带处(如中部山带)出现褶皱和逆冲,实现伸展背景下的局部挤压;俯冲板块后撤(~140Ma)则使上涌的软流圈回落形成岩石圈并实现地幔小幅增生置换(~125Ma)与伸展背景下浅部地壳断陷和成盆过程.因此,西太平洋板片俯冲和后撤是引起华北东部深部岩石圈地幔置换并导致陆内浅部盆-山耦合的外在动力来源,表明华北克拉通破坏是地块内部与地块边缘、深部过程与浅部盆-山耦合响应的综合地质记录,我们认为这也是燕山运动的本质.     

西北太平洋俯冲带日本本州至中国东北段应力场反演

本研究基于Global CMT提供的1196个1976年11月—2017年1月MW4.6地震矩心矩张量解,对西北太平洋俯冲带日本本州至中国东北段的应力场进行反演计算,得到了从浅表到深部俯冲带应力状态的完整分布.结果显示:俯冲带浅表陆壳一侧应力场呈现水平挤压、垂向拉伸状态,洋壳一侧的应力状态则相反,即近水平拉张、近垂向压缩.沿着俯冲板片向下,应力主轴逐渐向俯冲板片轮廓靠拢,其中位于双地震层(120km深度附近)之上的部分,主张应力轴沿俯冲板片轮廓展布而又比其更为陡倾;双地震层内的应力模式同典型I型双层地震带内的应力模式一致,即上层沿俯冲板片轮廓压缩、下层沿俯冲板片轮廓拉伸;双地震层之下,应力模式逐步转变为主压应力轴平行于俯冲板片轮廓.通观所研究的整个俯冲系统,水平面内主压和主张应力轴基本保持了与西北太平洋板片俯冲方向上的一致性,同经典俯冲板片的应力导管模型所预言的俯冲带应力模式相符;而主张应力轴在俯冲板片表面之下的中源地震深度范围内转向海沟走向,或许同研究区域横跨日本海沟与千岛海沟结合带,改变的浅部海沟形态致使完整俯冲板片下部产生横向变形有关.     

蒙古-贝加尔地区上地幔小尺度对流及 地球动力学意义

上地幔小尺度对流是控制区域地球动力学过程的主要机制之一,蒙古-贝加尔地区的一些区域动力学过程被认为与上地幔小尺度对流相关.本文目的在于利用重力资料研究蒙古-贝加尔地区的上地幔小尺度对流,并探讨其与构造动力学的关系.基于区域均衡重力异常与上地幔小尺度对流的相关方程,本文利用区域均衡重力异常资料反演了蒙古-贝加尔地区上地幔小尺度对流流场及作用于岩石层底部的应力场.结果显示,蒙古-贝加尔地区地幔流场及对流应力场呈现非常复杂的图像,流场及应力场分布与地表构造具有很好的相关性.西伯利亚地台和蒙古褶皱带下地幔流场和对流应力场均较弱,这与这些地区现今较弱的构造活动性是一致的.贝加尔裂谷区下存在地幔上升流,对流应力场呈拉张状态,但应力场的幅值较小(约8 MPa),表明地幔对流不是贝加尔裂谷开裂的主要控制因素.Hangay高原、阿尔泰和戈壁-阿尔泰下存在地幔上升流,对流应力场为拉张状态,这一方面可能构成Hangay高原隆升的深部动力机制,另一方面,也为Amurian板块西边界划分提供了动力背景.     

华南陆缘高热流区的壳幔温度结构与动力学背景

华南陆缘是我国重要的矿产、地热资源区.晚中生代以来,在太平洋板块西向俯冲,地幔热对流活动共同作用下,该区出现多期岩浆-热事件和大规模爆发式成矿作用.在前人研究基础上,本文利用地表热流观测资料、地震剪切波资料、重力位球谐系数,计算了壳-幔温度结构,分析了动力学背景.计算结果表明：华南陆缘东南沿海地带,地壳10 km以浅温度达200℃以上,居里点温度475℃,莫霍面平均温度550℃.地壳浅层较热,花岗岩中放射性元素衰变放热是地壳浅层地下水热活动的重要热源,但地壳总体温度不高,为"冷壳热幔"型热结构.地幔中,90 km深度,温度950~1250℃;120 km深度,温度1050~1400℃;150 km深度,温度1200~1450℃;220 km深度,温度1500~1700℃."热"岩石圈底界深度在110~150 km之间,西深东浅.岩石圈内,地幔应力场为挤压-伸展相间格局;岩石圈之下,地幔应力场为一个以南昌为中心、长轴NE-SW向的椭圆.分析认为,晚中生代以来,太平洋板块的西向俯冲,导致华南陆缘在区域性SE向地幔对流背景上叠加局域性不稳定热扰动,在175~85Ma期间,上地幔物质向上流动,形成不同的岩浆活动高峰期.同时,岩石圈地幔受俯冲洋壳流体的影响,含水量高,黏度小,在地幔流切向应力场作用下,岩石圈底界由西向东"波浪"状减薄.现今岩石圈之下仍具备地幔小尺度热对流温度条件,但除地表浅层外,地壳整体温度不高,岩石圈构造稳定.     

帕米尔—兴都库什地区构造应力场反演及拆离板片应力形因子特征研究

从Global CMT目录搜集了1976年1月至2016年6月之间的震源深度大于70km的255个震源机制解,用阻尼应力反演方法,分70～160km和170～310km两个深度,计算了帕米尔—兴都库什地区的构造应力场;同时以10km为间隔计算了兴都库什地区深度介于70～310km之间的应力形因子.得到以下初步结论:兴都库什板片向下俯冲和帕米尔地区断裂带的横向拉张,可能是导致应力场不同的原因.兴都库什俯冲带与帕米尔俯冲带碰撞,导致交汇地区(37°N—37.5°N)的应力场参数突变.兴都库什俯冲板片受到深部温度、压力等因素,出现薄弱面进而形成拆离板片.其脱离了主俯冲板片的束缚后,重力的上下拉张作用导致空区附近张轴倾伏角接近90°,拆离板片俯冲至上地幔不连续面,导致板片部分熔融进而应力形因子随着深度变小.而拆离板片受到地幔挤压其内部发生破碎,其压应力轴由西部的NS到东部NW-SE方向偏转及纵向张应力轴倾伏角变小.     

大陆岩石圈、地幔底部异常体与地幔对流相互作用的数值模拟

在地球表层存在着占地表面积约30%的具有低固有密度、高黏度的大陆岩石圈.由于其特殊的物理化学性质,大陆岩石圈通常不直接参与下方的地幔对流,但其与地幔对流格局有着重要的相互影响.大量研究显示,在中太平洋和非洲的下地幔底部,存在着两块占核幔边界（CMB）面积约20%的高密度热化学异常体（由于其剪切波速度较低,常称作低剪切波速度省（LSVPs））.LSVPs的演化既受地幔对流的影响,同时也影响地幔物质运动的格局和动力学过程.本文系统研究了存在大陆岩石圈,下地幔LSVPs的地幔对流模型.模拟结果显示：（1）当大陆体积较小时,其边缘常伴随着俯冲,大陆区域地幔常处于下涌状态,其上地幔温度较低,大陆岩石圈在水平方向处于压应力状态.随着大陆体积的增大,大陆边缘的俯冲逐渐减弱,大陆区域地幔由下涌转为上涌,其上地幔温度较高,大陆岩石圈水平方向处于拉应力状态.（2） 岩石圈与软流圈边界（LAB）在大陆下方较深,温度较低;在海洋区域较浅,温度较高.随着大陆体积的增大,陆洋之间LAB深度、温度的差异逐渐减小.（3）大陆区域地幔底部LSVPs物质的丰度与大陆的体积呈正相关.当大陆体积较小时,大陆下方的LSVPs丰度比海洋区域少.随着大陆体积的增大,大陆下方LSVPs的丰度逐渐增大.（4）海洋地区地表热流高,且随时间波动大,大陆地区地表热流低,随时间波动较小;LSVPs区域的核幔边界热流低.     

板块运动与地震各向异性及应力场的相关性统计分析

利用收集到的各种来源共计7 959组的地震各向异性观测数据和21 750组应力场数据,结合板块绝对运动模型计算给出的各板块的运动规律,分别统计分析了板块运动与地震各向异性及应力场的相关性,并对板块运动对地震各向异性及应力场特征产生的影响进行了分析． 统计结果表明,阿拉伯、 加勒比、 胡安德富卡、 北美、 纳兹卡、 太平洋和南美板块上地震各向异性与板块运动均具有较好的相关性,而非洲、 南极洲、 澳大利亚、 欧亚、 印度和菲律宾板块上二者的相关性则相对较差. 讨论分析发现,板块运动拖动软流圈流动、 橄榄岩晶格优选方位、 化石各向异性和地幔流动或岩石圈流动等因素均在一定程度上控制并影响着地震各向异性与板块运动的一致性． 而板块基底拖曳力、 洋脊推力、 浮力作用和碰撞及俯冲作用等多种因素共同制约了板块运动与应力场的相关性,使得非洲、 可可斯、 欧亚、 胡安德富卡、 北美、 纳兹卡、 菲律宾和南美板块上二者的相关性较好,其它板块上其相关性则较差． 对于俯冲带地区,由于俯冲机制的复杂性和软流圈、 岩石圈地幔流动方向的不确定性,其板块运动与地震各向异性及应力场的相关性图像表现复杂,需要结合具体的俯冲带构造进行近一步研究．      

Three-dimensional mantle lithosphere deformation at collisional plate boundaries: A subduction scissor across the South Island of New Zealand

The continental plate collision across the South Island of New Zealand is highly oblique (dextral) and bounded by oppositely verging ocean plate subduction zones. As such, the region can be considered as a type of ‘subduction scissor’. Within this tectonic context, we use three-dimensional computational geodynamic models to consider how convergent mantle lithosphere can be modified by scissor and strike–slip effects. Bounding subduction at both ends of the continental collision causes flow of the descending mantle lithosphere in the direction along strike of the model plate boundary, with thinning in the centre and thickening towards the subduction zones that bifurcates the continental mantle lithosphere root. With dipping bounding subduction, the mantle lithosphere root takes on a more complex morphology that folds over from one subduction polarity to the other, but remains as a continuous feature as it folds under the collision zone. In the absence of bounding subduction, the plate convergence causes a linear (along strike) mantle lithosphere root to develop. A rapid strike–slip motion between the converging plates transfers material in the plate boundary-parallel direction and tends to blur out features that develop in this direction—such as descending viscous instabilities. The along-strike variations in the morphology of the mantle lithosphere root that develop in the models—viz., thickening of the root towards the subduction edges, thinning in the centre—are consistent with recent, albeit poorly constrained, geophysical interpretations of the large-scale lithospheric structure of the South Island. We speculate that this reflects the nature of the evolution of the South Island collision as a limited continental segment of the plate boundary that it is dominated and guided by adjacent well-developed/developing ocean plate subduction.     

Deep dynamical processes in the central-southern Qinghai-Tibet Plateau—Receiver functions and travel-time residuals analysis of north Hi-Climb

Teleseismic receiver functions and travel-time residuals along the north Hi-Climb broadband seismic array in the central-southern Qinghai-Tibet Plateau show that the lithosphere structures in the central and western Qinghai-Tibet Plateau are different. In the central Qinghai-Tibet Plateau, the Indian Plate is northward subducted beneath the Qiangtang block and arrives at the greatest depth beneath the central-southern Qiangtang block. The delaminated Indian lithospheric slab remains beneath the central Lhasa block to a depth possibly greater than that of the upper interface of the mantle transform zone. In the western Qinghai-Tibet Plateau, the Indian lithospheric plate is gently northward subducted and may have arrived to the south of Tarim plate. Due to the resistance from the gently northward subduction of the Indian mantle lithosphere in the western Qinghai-Tibet Plateau, the upwelling mantle material be-neath the Qiangtang block moves mostly toward the east to bring about the lateral eastward flow of the deep mantle hot material in the central Qinghai-Tibet Plateau.     

华北克拉通及东邻西太平洋活动大陆边缘地区的P波速度结构：对岩石圈减薄动力学过程的探讨

本文通过地震层析成像研究获得了华北克拉通及其东邻地区（30°N-50°N,95°E -145°E）1°×1°的P波速度扰动图像.结果显示,在西太平洋俯冲带地区,上地幔中西倾的板片状高速异常体与其上方的低速异常区构成俯冲带与上覆地幔楔的典型速度结构式样.俯冲板片高速体在约300~400 km深度范围内被低速物质充填,暗示俯冲板片可能发生了断离.在华北克拉通地区的上地幔中发现三个东倾排列的高速异常带.在此基础上,本文构建了华北克拉通及其东邻西太平洋活动大陆边缘地区的上地幔速度结构模式图,并据此探讨克拉通岩石圈减薄与西太平洋活动大陆边缘的深部动力学联系.本文认为,太平洋板片的俯冲（断离）,触发热地幔物质上涌并在上覆地幔楔中形成对流,使克拉通岩石圈受到改造（底侵与弱化）.随着俯冲板片后撤,地幔楔中的对流场以及对岩石圈改造的影响范围均随之东移,最终导致华北克拉通岩石圈自下而上、从西向东分三个阶段依次拆沉减薄.这一模式能很好地解释现今克拉通岩石圈自西向东呈台阶状减薄的深部现象.     

俯冲带耦合作用对苏门答腊地区应变场影响的三维数值模拟

采用有限元方法模拟了俯冲带耦合作用对巽他弧及其邻区的影响.根据模拟结果,对比GPS、地震和地质学观测数据,定量分析了苏门答腊及其周边地区的应变强度和主应变方向的分布特征,据此探讨了该区构造特征、地震发生模式与耦合面积之间的关系.模型由具有黏弹性性质的岩石圈和软流圈上地幔组成,其中岩石圈包括了大陆岩石圈和大洋岩石圈以及俯冲至上地幔中的俯冲板片.研究结果如下:(1) 通过对不同俯冲带耦合面积模拟,发现苏门答腊前弧伴随耦合面积的增加应变强度逐渐增大,而增大的应变强度又影响了其周边地区的应变分布特征,因此整个苏门答腊前弧呈现出明显的分段性,这与该区地震破裂模式有较好的对应.(2)苏门答腊北部地区主应变方向与南部相比存在一定的差异,该差异是俯冲带的俯冲方向、俯冲速度、俯冲形态以及不同区域间耦合面积共同作用的结果.(3)虽然苏门答腊2004年地震主震区处于弱耦合状态,但从本文模拟的结果中可以看到,在俯冲作用下该区依然存在垂直向下的位移,这为地震激发海啸提供了有利的构造环境.     

华北中西部和青藏高原东北缘上地幔各向异性变形特征

基于华北中西部和青藏高原东北缘3个流动台阵共480个台站新得到的远震XKS(SKS、SKKS和PKS)波分裂结果,并结合研究区已得到的987个台站的分裂结果,获得了高分辨率的上地幔各向异性图像.分析表明,鄂尔多斯块体的时间延迟较小,反映了其稳定性和弱的各向异性变形特征,可能保留了古老克拉通根的"化石"各向异性,但其靠近...     

汤加—克马德克俯冲带现今非均匀应力场特征及其动力学意义

汤加—克马德克俯冲带是太平洋板块向澳大利亚板块俯冲碰撞的动力作用区,是全球俯冲带动力学研究的热点区域.本研究基于EHB地震目录,对汤加—克马德克俯冲带(18.5°S—28.5°S)区域进行平面拟合,得到该范围内俯冲带走向约为196°,倾角约为48°;利用该俯冲带研究区域内Global CMT目录,对不同位置、不同深度进行区域应力张量反演,得到汤加—克马德克俯冲带研究区内精细的应力图像.结果显示:(1)俯冲带浅部(60~300km)应力结构非均匀特征明显,主应力轴倾伏角变化多样,并且最大主压应力轴方位在24°S左右发生明显偏转,我们推测这可能与洋底构造路易斯维尔海链俯冲有关;(2)中部(300~500km)最大主压、主张应力轴由北向南逐渐发生偏转,这可能与由北向南流动的地幔流对俯冲板片产生推挤作用有关,并且这种推挤作用向南逐渐减弱;(3)深部(500~700km)最大主压应力轴沿俯冲方向分布;(4)本文的结果还发现了主俯冲带深部西侧"偏移"板片与主俯冲带应力结构不同,表明"偏移"板片与主俯冲带是分离的.     

冲绳板块应力场数值模拟及其动力学特征div>

冲绳板块位于菲律宾海板块向欧亚板块俯冲形成的西太平洋边缘活动带上,构造应力场图像及其动力学机制表现得相当复杂．采用伪三维有限元方法,以WSM2008 观测应力场数据的应力取向和应力型两方面指标作为主要约束,对冲绳板块构造应力场进行了数值模拟．通过对计算结果的分析,对模型涉及的各种作用力作出了估计．在此基础上,对冲绳板块岩石层的状态,以及该地区的板块动力学特征进行了探讨,并得到了以下一些初步认识:① 软流层静压推力控制着该地区构造应力场的基本形态;② 冲绳海槽的演化过程,例如该地区的岩石层减薄与其下地幔流的上升等,也在很大程度上影响了该地区的板内应力场空间分布特征;③ 琉球俯冲带边界力的作用是分段的,不同区段作用力对板内应力场的影响有所不同．      

Geodynamic model of the evolution of the Pacific Ocean

The present Pacific Ocean differs significantly in its structure and evolution from the expanding Atlantic Ocean. The Pacific is asymmetric. Its mid-ocean ridge is located not along its median line but is closer to South America and adjoins North America. The Pacific is surrounded by a ring of subduction zones but has marginal seas only at its Eurasian margins. After the breakup of Pangea, the Atlantic began to open and the Pacific began to close. This paper examines the evolution of the Pacific Ocean and, in particular, the formation mechanisms of its present structures. Numerical modeling of the long-term drift of a large continent is performed, with the initial position of the continent corresponding to the state after the breakup of the supercontinent. At first the continent, driven by the nearest descending mantle flow, begins to approach a subduction zone. Since the mantle flows beneath a large continent have different directions, its velocity is a few times lower than that of the mantle flows near the subduction zone. As a result, a zone of extension arises at the active continental margin and a fragment is broken off from the continent; this fragment rapidly moves away and stops above the descending mantle flow as in a trap. A marginal sea forms at the active continental margin. The continent continues its slow movement toward the subduction zone. The oceanic lithosphere, which earlier sank vertically, begins to descend obliquely. This evolutionary stage corresponds to the present position of Eurasia. The modeling shows how the interaction of the continent with the mantle causes the subduction zone to roll back toward the ocean. Subsequently, the continent nevertheless catches up with the subduction zone, and they move together for a while. The marginal sea then closes and high compressive stresses arise at the active continental margin. This state corresponds to the present position of South America. During the subsequent drift, the continent together with the subduction zone reaches the mid-ocean ridge and partially overrides it. This state corresponds to North America, which was the first to break off from Pangea and passed through the stages of both Eurasia and South America. The large and slowly moving Eurasia, which formed only at the time of Pangea, is still in the first evolutionary stage of the Pacific Ocean closure.