深部俯冲板片三维构造重建及其几何学、运动学研究——以汤加—克马德克地区俯冲板片为例

西南太平洋板块与澳大利亚板块之间的汤加—克马德克俯冲带,是研究地球动力学最重要的区域之一.本文研究根据MIT-P08地震数据,结合板块构造边界、地震活动分布、海岸地形数据等,基于GOCAD软件平台建立三维地震层析成像,对西南太平洋板块的汤加—克马德克俯冲板片进行三维解释.地震层析成像显示汤加—斐济地区地幔至少存在三个"高速"异常体.早期汤加—克马德克俯冲板片穿过地幔转换带,并进入下地幔,最大深度达1600 km.三维构造模型揭示了汤加—克马德克板片在深度600~800 km处存在断折形变,该俯冲板片去褶皱恢复后,测量其俯冲的最大位移达2600 km.汤加—克马德克板片开始快速俯冲的时间至少在30 Ma之前,平均移动速率约为68~104 mm /a.俯冲板片三维构造重建和恢复,可以有效揭示俯冲板片几何学、运动学,为研究深源地震成因、地球深部变化过程和动力学机制提供约束.     

大华北浅源地震与日本海西部及我国东北深震的关系

分析了大华北浅源地震与日本海西部及我国东北深震的关系,认为本世纪来日本海西部—我国东北深震经历了5个相对活跃期,大华北各地震区相应经历这5个活跃期的影响期。根据大华北M≧6级浅源地震与深震活动的相关性,建立了太平洋板块楔形俯冲带端部重大深震事件导致大华北浅源M≧6级地震发生的板块俯冲模型,应变波传播速度约94km/年,地表视速度约100km/年。重大深震事件突出、模型稳定性强,预测实验表明模型公式可做大华北地震监测参考。用本模型可以解释浅源地震迁移、各地震区地震与深震活动相关等现象。     

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

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

     

俯冲带断层粘滑运动机制数值模拟研究——以日本俯冲带为例

在新生代,西太平洋岛弧和边缘海地区已成为强烈的构造活动区。至晚第三纪以来,在西太平洋岛弧带上还发生了强烈的火山作用。在西太平洋岛弧的外侧发育着地球表面上最深的地沟带,即深海沟带。西太平洋板块向西俯冲,在日本海形成俯冲带,称为日本海俯冲带。日本俯冲带构成了西太平洋边缘俯冲带的重要组成部分。该俯冲带处于环太平洋火山地震带上,地质构造复杂,地震活动频繁,地震发生数量占全球20％。     

长白山天池火山地震活动与西北太平洋俯冲带内中深源地震关系的研究

根据长白山火山活动的地球动力学背景,综合分析西北太平洋俯冲带俯冲作用所引起的中深源地震与长白山火山地震活动及温泉流体化学气体释放变化趋势的关系。认为在北纬35°以北,西北太平洋板块俯冲作用引起的地幔对流的扰动作用下,长白山火山活动从1999年至目前大体经历了3个阶段:1999年6月至2002年5月为扰动起始阶段,2002年6月至2004年为扰动增强阶段,2005年至目前为扰动衰减阶段。在分析了3个阶段内长白山火山地震活动及温泉气体释放特征的基础上,讨论了长白山天池火山目前的活动状态。从日本海至中国东北深源地震区地震活动期和平静期存在的韵律性活动过程出发,对长白山天池火山未来喷发的可能性进行了分析     

中国东北的深源地震波形匹配检测及定位

中国东北珲春周边地区位于环太平洋地震带上,也是中国唯一存在的深源地震带.较大地震发生后常会有若干震级较小的余震发生,但在相同的地震震级情况下深源地震的余震一般比浅源地震的余震数量要少1~3个数量级,且在全球不同深震区的深震余震数量也存在显著差异.针对国际地震中心（ISC）2010年7月至2014年12月目录中给出的中国东北附近27次震源深度超过300 km的深源地震,我们首先利用区域固定地震台及NECsaids流动地震台阵的连续波形数据,选取已知地震事件作为模板,采用Match & Locate及Matched Filter方法进行波形互相关叠加分析来检测微小深震事件;然后对1966至2017年ISC目录中的东北地区深源地震进行双差重定位以提高震源位置的准确性,进一步分析深震活动与俯冲板片的关系.研究结果显示除ISC目录中给出的深震事件,我们未能检测出作为模板的深源地震的前震或余震活动,证明东北深源地震余震活动较少并不是由于台站分布有限而造成的漏检结果;重定位后震源延伸的角度与西太平洋板片在410~660 km地幔转换带内的俯冲角度较为一致,并且大部分深震震源位置位于俯冲板片中的亚稳态橄榄岩楔形区内部.结合双差定位结果、b值分析及前人研究成果,我们认为东北深源地震应不属于与俯冲非直接相关的"孤立地震",而是与西太平洋板块俯冲直接相关.     

2011年5月中俄边界深震震源机制所展示的太平洋板块俯冲动力作用

太平洋板块向欧亚板块俯冲(常称为日本俯冲带)经日本海沟插入中国东北大陆下约600 km深处,到达吉林珲春一带,造成东北地区频繁发生深震活动.不少研究认为,太平洋板块俯冲是中国东北及邻区岩石圈减薄和伸展作用的重要动力来源.     

基于重力数据特征的东北地区深部动力机制解析

中国东北地区处于古亚洲构造域、蒙古-鄂霍茨克构造域和环太平洋构造域叠加作用最为显著的地区,是地学研究的热点区域.为了探析欧亚大陆下西太平洋板片的俯冲形态以及揭示该区域深部地球动力学机制,利用卫星重力数据通过预处理共轭梯度快速密度反演算法获得了包含东北地区、华北部分地区及日本海海域在内的研究区域上地幔三维密度结构,结合天然地震三维层析成像结果共同揭示太平洋板片的俯冲形态和深部动力机制.俯冲的太平洋板片在日本海沟处呈高密度异常,向西横向持续扩张,深度方向上有逐渐增加趋势.不连续的高密度体俯冲至地幔转换带（410~660 km）后继续水平西向俯冲,部分滞留板片可能进入下地幔;在大兴安岭断裂带下面转换带中同样发现水平分布的高密度体,推断是大兴安岭断裂带下方地幔岩石圈拆沉的结果,横向不均匀分布的俯冲板片边缘已抵至大兴安岭造山带附近,这对于研究东北地区深部动力学机制具有重要的意义.

     

苏门答腊俯冲带地区壳内密度结构研究

俯冲带上覆板片密度特征是认识俯冲及其引发深部岩浆过程的一个窗口.本文以苏门答腊俯冲带空间重力异常数据为基础,在2.5D密度结构剖面约束下,通过3D密度反演,获得了研究区3D密度结构分布.反演结果表明,俯冲板片角度和下倾极限深度均从研究区西北向东南方向逐渐增加;研究区上覆板片下地壳存在低密度异常,主要集中在东南部,分布范围也随深度而增加.分析认为这有可能是俯冲引发的洋壳脱水对上覆板片下地壳侵染所造成.此外,研究区东南部存在的地幔板片撕裂可能是造成该区下地壳低密度异常最为显著的另一个原因.3D密度切片揭示出下地壳密度异常具有沿NE方向延伸迹象,推断应与印澳板块在苏门答腊地区的斜向俯冲作用有关.本文还对研究区的地震分布特征进行了讨论.大部分的浅源地震集中在下倾极限附近,为脆性破裂或摩擦滑移所引发.震源深度大于200 km的地震基本分布在研究区中部和东南部,震源深度从西北向东南方向逐渐加大,这从另一角度为研究区东南部地幔存在板片撕裂的观点提供了佐证.

     

西北太平洋俯冲带东北地区壳幔结构研究进展

总结了近年来西北太平洋俯冲前端东北地区地壳的精细结构和上地幔间断面410和660km的研究成果.认为：该区莫霍面深约为30～39km.410km的间断面有30km的深度异常,这种异常是否在过渡带正常或含水相变展布厚度的内,还是于该区俯冲带的后退有关,还需做进一步的研究论证.在中国东北地区,日本海俯冲带向西倾斜,在660km的间断面上近似水平停滞汇集被大多数层析成像、接收函数、波形拟和以及深震研究所认同.然而,俯冲带具体的形式、大小范围,向下渗透到多深以及在局部地区表现的横向不连续性的看法并不一致.由此带来的动力学模式什么力来支撑着停滞的板片,在板片最终下沉到下地幔以前,有多少俯冲的板片能停滞在间断面之上,这些问题都还需要更合理的解释.     

Seismic evidence for a metastable olivine wedge in the subducting Pacific slab under Japan Sea

We apply a forward-modeling approach to high-quality arrival time data from 23 deep earthquakes greater than 400 km depth to investigate the detailed structure of the subducting Pacific slab beneath the Japan Sea. Our results show that a finger-like anomaly exists within the subducting Pacific slab below 400 km depth, which has a P-wave velocity 5% lower than the surrounding slab velocity (or 3% lower than that of the normal mantle), suggesting the existence of a metastable olivine wedge (MOW) in the slab. The MOW top and bottom depths are 400 and 560 km, respectively. The MOW is estimated to be about 50 km wide at 400 km depth and close to the slab upper boundary. At 560 km depth the MOW is located at about 25 km below the slab upper boundary. Most of the deep earthquakes are located in the MOW. Our results favor transformational faulting as the mechanism for deep earthquakes.     

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

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

南阿拉斯加地壳及上地幔结构成像研究

通过反演562891个纵波和156321个横波走时数据,第一次同时获得了阿拉斯加地区地壳及上地幔的纵波与横波速度以及泊松比图像,为更好地认识阿拉斯加地区的深部地震结构、太平洋板块与亚库塔特板块的俯冲几何形态提供了科学依据.成像结果表明P波和S波速度图像与泊松比结构具有很好的一致性,强的高速度和低泊松比异常沿着阿拉斯加俯冲带延伸至200 km深度,该高速度和低泊松比异常体与俯冲带的地震空间分布吻合,因此,我们认为该高速体为俯冲的太平洋板块和亚库塔特板块.从地震空间分布发现,大部分大地震（M>6.5）发生在高速度与低速度异常交界处,可能反映了俯冲板块之间强耦合作用.在俯冲带的地幔楔显示出广泛的低速度和高泊松比异常,并且这些异常与岛弧火山的位置相对应,这与大洋板块俯冲所形成的岩浆入侵作用有关.研究结果表明在南阿拉斯加俯冲带,俯冲板块的俯冲角度从兰格尔块体下方的平坦变成在布里斯托尔湾下方的陡峭,这与亚库塔特板块俯冲在兰格尔块体下方和太平洋板块俯冲在布里斯托尔湾下方有关.在基奈半岛和科迪亚克岛连接处的上地幔位置存在强烈的低速与高泊松比异常体,使该处的大洋俯冲板块变薄.这一现象可能与亚库塔特板块和太平洋板块相互碰撞作用以及软流圈强烈的上升流入侵有关.     

Dynamics of thinning and destruction of the continental cratonic lithosphere: Numerical modeling

Thinning of the cratonic lithosphere is common in nature, but its destruction is not. In either case, the mechanisms for both thinning and destruction are still widely under debate. In this study, we have made a review on the processes and mechanisms of thinning and destruction of cratonic lithosphere according to previous studies of geological/geophysical observations and numerical simulations, with specific application to the North China Craton (NCC). Two main models are suggested for the thinning and destruction of the NCC, both of which are related to subduction of the oceanic lithosphere. One is the “bottom-up” model, in which the deeply subducting slab perturbs and induces upwelling from the hydrous mantle transition zone (MTZ). The upwelling produces mantle convection and erodes the bottom of the overriding lithosphere by the fluid-melt-peridotite reaction. Mineral compositions and rheological properties of the overriding lithospheric mantle are changed, allowing downward dripping of lithospheric components into the asthenosphere. Consequently, lithospheric thinning or even destruction occurs. The other is the “top-down” model, characterized by the flat subduction of oceanic slab beneath the overriding cratonic lithosphere. Dehydration reactions from the subducting slab would significantly hydrate the lithospheric mantle and decrease its rheological strength. Then the subduction angle may be changed from shallow to steep, inducing lateral upwelling of the asthenosphere. This upwelling would heat and weaken the overriding lithospheric mantle, which led to the weakened lithospheric mantle dripping into the asthenosphere. These two models have some similarities, in that both take the subducting oceanic slab and relevant fluid migration as the major driving mechanism for thinning or destruction of the overriding cratonic lithosphere. The key difference between the two models is the effective depth of the subducting oceanic slab. One is stagnation and flattening in the MTZ, whereas the other is flat subduction at the bottom of the cratonic lithosphere. In the NCC, the eastern lithosphere was likely affected by subduction of the Izanagi slab during the Mesozoic, which would have perturbed the asthenosphere and the MTZ, and induced fluid migration beneath the NCC lithosphere. The upwelling fluid may largely have controlled the reworking of the NCC lithosphere. In order to discuss and analyze these two models further, it is crucial to understand the role of fluids in the subduction zone and the MTZ. Here, we systematically discuss phase transformations of hydrous minerals and the transport processes of water in the subduction system. Furthermore, we analyze possible modes of fluid activity and the problems to explore the applied feasibility of each model. In order to achieve a comprehensive understanding of the mechanisms for thinning and destruction of cratonic lithosphere, we also consider four additional possible dynamic models: extension-induced lithospheric thinning, compression-induced lithospheric thickening and delamination, large-scale mantle convection and thermal erosion, and mantle plume erosion. Compared to the subduction-related models presented here, these four models are primarily controlled by the relatively simple and single process and mechanism (extension, compression, convection, and mantle plume, respectively), which could be the secondary driving mechanisms for the thinning and destruction of lithosphere.     

Suboceanic earthquake location and seismic structure in the Kanto district,central Japan

We present a combined method, using sP depth-phase data and double-difference arrival times, to determine the precise hypocenter locations of earthquakes that occur under the Pacific Ocean outside of the area covered by the land-based seismic network. We assess the effectiveness of the combined method using a data set of P- and S-wave arrival times and sP depth phase from suboceanic earthquakes recorded by both land-based seismic stations and offshore seismic stations (OFS). The hypocenters of the offshore earthquakes relocated using the combined method are consistent with those determined using the standard location method and OFS data. The differences in the hypocenters relocated by the two methods are less than 4 km. We applied the method to the subduction region that underlies the Kanto district, central Japan, and located a large number of earthquakes that occurred beneath the Pacific Ocean. We then determined the detailed 3D seismic velocity structure by inverting a large number of arrival times of P- and S-waves and sP depth phase from the relocated earthquakes in the study region. High-velocity anomalies related to the cold subducting Pacific slab and low-velocity anomalies related to the hot mantle wedge are clearly imaged. Beneath active volcanoes, low-velocity zones are visible from the surface to a depth of 100 km, reflecting fluids released by dehydration of the subducting Pacific slab. Strong lateral heterogeneities are revealed on the upper boundary of the Pacific slab beneath the forearc region. The low-velocity areas under the offshore region are associated with low seismicity and weak interplate coupling. A low-velocity layer is imaged along the upper boundary of the Philippine Sea slab in the northern part of Kanto district, which may reflect dehydration of the slab. Our tomographic images indicate that the overlaying Philippine Sea plate has effects on the spatial distribution of active volcanoes related to the subducting Pacific slab in the study region.     

Kamchatka subduction zone,May 2013: the Mw 8.3 deep earthquake,preceding shallow swarm and numerous deep aftershocks

A sequence of 98 teleseismically recorded earthquakes occurred off the east coast of Kamchatka at depths between 10-90 km around latitude 52.5°N and longitude 160°E on May 16–23, 2013. The swarm occurred along the northern limit of the rupture area of the 1952 M_w 9.0 great Kamchatka earthquake, the fifth largest earthquake in the history of seismic observations. On May 24, 2013 the strongest deep earthquake ever recorded of M_w 8.3 occurred beneath the Sea of Okhotsk at a depth of 610 km in the Pacific slab of the Kamchatka subduction zone, becoming the northernmost deep earthquake in the region. The deep M_w 8.3 earthquake occurred down-dip of the shallow swarm in a transition zone between the southern deep and northern shallow segments of the Pacific slab. Several deep aftershocks followed, covering a large, laterally elongated part of the slab. We suppose that the two described earthquake sequences, the May 16–23 shallow earthquake swarm and the May 24–28 deep mainshock-aftershock series, represent a single tectonic event in the Pacific slab having distinct properties at different depth levels. A low-angle underthrusting of the shallow part of the slab recorded by the shallow earthquake swarm activated the deep part; this process induced the deep mainshock-aftershock series only three days after the swarm. The domain of the subducting slab activated by the May 2013 earthquake occurrence was extraordinarily large both down-dip and along-strike.     

东北亚边缘地区地幔过渡带内滞留太平洋板片上界面的三重震相研究

东北亚地区下方地幔过渡带与太平洋俯冲滞留板片的相互作用对于区域深部物质运移和区域构造演化具有重要影响.基于中国国家数字台网记录的发生于鄂霍次克海与千岛群岛地区的两个深源地震的宽频带波形资料, 我们发现了与660-km间断面上方一个特殊界面有关的清晰的新P波三重震相波形; 通过对观测波形进行波形拟合以及搜索分析, 本文获得了东北亚边缘地区下方660-km间断面附近的精细速度结构.结果表明: 东北亚边缘下方的660-km间断面之上存在一高速异常层, 其具有尖锐的速度异常上界面, 深度介于455~510 km之间, P波速度异常达2%~4.5%;与此同时, 该地区下方660-km间断面整体速度跃变量较小, 且存在0~15 km的下沉; 660-km间断面下沉与高速层上界面的形态具有很强的相关性, 均表现为西浅东深、南浅北深的特征.结合前人研究结果, 我们推测受日本—千岛海沟回撤速率差异的影响, 西北太平洋俯冲板片对上地幔底部间断面作用的差异是造成660-km间断面整体形态南北差异的主要原因; 而俯冲板片在地幔过渡带内的滞留并向西展平堆积使得板片上界面西侧较东侧更浅; 滞留板片上界面的速度异常特征显示, 俯冲沉积物可能到达地幔过渡带且经脱水相变形成了含黄玉的矿物.