长白山火山下方地幔转换带中滞留的俯冲太平洋板块存在空缺吗?

针对近年来长白山火山下方地幔转换带中是否存在低波速异常指示的太平洋板块"空缺"而引起的不同科学认识的热烈辩论,本文主要回顾了我国东北地区地幔转换带的体波成像结果。使用相对走时残差的远震体波成像结果显示,长白山火山以西地幔转换带中存在低波速异常指示的太平洋板块"空缺";而使用绝对走时残差的区域成像和全球成像结果,尽管展示出长白山火山以西比以东略低的波速异常,但长白山火山以东至我国东北重力梯度带区域下方的地幔转换带均展示出明显的连续的高波速异常。在接收函数分析时,如果以全球平均值660km而非我国东北地区平均值670km作为基准,来分析660km间断面是抬升还是下沉;以全球平均值250km而非我国东北地区平均值260km作为基准,来分析地幔转换带是增厚还是减薄的话,则可以得到长白山火山以东至我国东北重力梯度带区域660km间断面下沉与地幔转换带增厚的认识。这种与绝对走时残差成像结果展示的地幔转换带为连续的高波速异常结果相一致的结果,说明太平洋板块俯冲前缘已由日本海沟抵达我国东北松辽盆地与大兴安岭交界处。结合高温高压实验、数值模拟与岩石地球化学研究结果,本文并不支持长白山火山以西的地幔转换带存在低波速异常指示的板块"空缺"和地幔转换带"减薄"的认识。长白山火山的深部起源与太平洋板块深俯冲至我国东北松辽盆地与大兴安岭交界处形成的"大地幔楔"结构动力学相关。     

郯庐断裂带及其邻区上地幔顶部Pn波速度与各向异性层析成像

本研究拾取了中国数字测震台网固定台站记录的2008-2016年2级以上地震事件中的27233条高质量Pn到时资料,反演得到了郯庐断裂带及其邻区上地幔顶部Pn波速度和各向异性结构模型.结果显示,研究区上地幔顶部Pn波速度结构存在强烈的横向不均匀性,速度异常形态与区域地质构造较为吻合.太行山造山带、鲁西隆起、大别造山带、苏鲁褶皱带、胶辽隆起和华北盆地南端等隆起区表现为低波速异常,而黄海北、南部盆地、渤海湾和华北盆地北部等凹陷区均为高波速异常.壳内强震主要发生在Pn低波速异常和高低波速异常的横向过渡地带,说明强震的发生与上地幔结构的横向变化之间存在有一定关联.郯庐断裂带两侧Pn波速度以郯城地震为界其东北侧和西南侧分别分布有与断裂带近平行的低波速异常条带,而西北侧和东南侧分别分布有高波速异常条带,各向异性快波方向近乎沿断裂带走向,可能由于上地幔热物质沿郯庐断裂带上涌形成低速异常后断裂带发生左旋平移运动所致.华北盆地内上地幔顶部Pn波速度结构和各向异性的明显变化,反映华北克拉通破坏过程中经历了地幔热物质上涌、莫霍面隆升以及岩石圈拆沉等复杂构造变形.     

基于多辛结构谱元法的保结构地震波场模拟

近年来构造高精度、高效且具有长时程跟踪能力的保结构算法已逐渐成为地震波模拟算法发展的重要方向之一. 本文基于谱元法(SEM)进行空间域离散结合新推导的三阶辛算法(NTSTO)进行时间域离散,构造了一种具有时-空保结构特性的新算法. 本文给出的多组数值试验对比结果表明,本算法无论在内存消耗、稳定性及计算耗时,还是长时程跟踪能力方面都有上佳的表现; 另外,本文给出的起伏地表多层介质模型的数值算例验证了该算法处理复杂几何形状和复杂介质时的有效性. 该多辛结构谱元法的发展将为长时程地震波传播的计算及模拟提供更为广泛而有效的选择.     

前郭MS5.8和松原MS5.7地震震源区地壳速度结构与孕震环境

为了解2013年11月23日吉林前郭M_S5.8和2018年5月28日吉林松原M_S5.7地震震源区的地壳速度结构及其孕震环境,本研究收集了东北地区CEA固定台站、NECESSArray流动台站和五大连池WAVESArrary流动台站的连续波形数据,使用背景噪声成像方法,获得了前郭地震和松原地震震源区地壳S波速度结构.结果显示,这两个地震主要发生在相对高波速异常区,而震源区下方存在明显低波速异常,且该低速层呈北北东走向,并有向北延展的趋势.结果表明,震源区相对高波速异常区易于积累能量孕育地震,而震源区下方的低波速异常可能代表了流体作用.这种流体作用会降低断层面有效正应力从而触发地震.这种流体作用可能和太平洋板块深俯冲至我国东北地区（包括前郭地震和松原地震震源区）下方的地幔转换内形成"大地幔楔"结构与动力学密切相关.在"大地幔楔"结构中,由于地幔转换带中滞留板块脱水作用和地幔角流作用,容易形成热湿物质上涌,进而导致松辽盆地中北部岩石圈物质拆沉、地幔流体入侵至中下地壳、然后作用于断裂带,从而导致了两个震源区中强度地震的发生.     

最优化辛格式广义离散奇异核褶积微分算子地震波场模拟

将波动方程变换至Hamilton体系,构造了一种新的保结构算法,即最优化辛格式广义褶积微分算子(OSGCD). 在时间离散上,首先引入了Lie算子设计二级二阶辛格式,基于最小误差原理得到了优化的辛格式. 在空间离散上,引入广义离散奇异核褶积微分算子计算空间微分,提出了一种有效方法优化GCD并得到了稳定的算子系数. 针对本文发展的新方法,给出了OSGCD稳定性条件. 在数值实验中,将OSGCD与多种方法比较,从精度和计算效率两方面分析了OSGCD的计算优势,计算结果也表明OSGCD长时程以及非均匀介质中地震波模拟亦具有较强能力.     

中国东北地区SKS分裂的上地幔各向异性结构与动力学

本研究收集了中国东北地区2008—2016年九年时间内207个固定地震台站和127个NECESSArray流动地震台站的波形资料,利用SKS波分裂的最小切向能量网格搜索方法获得了243个台站的有效分裂结果.研究结果显示,尽管研究区各向异性快波方向基本以NW-SE向为主,但无论是在快波方向上还是快慢波时间延迟上不同构造单元内部与不同构造单元之间均存在着较大差别.大兴安岭造山带北部的各向异性快波方向自北向南由NNE-SSW向转变为NNW-SSE向,在中部以NW-SE向为主,而南部自北而南由NE-SW向逐渐转变为近E-W向;松辽盆地的各向异性快波方向在北部自西向东主要表现为由NNW-SSE向逐渐转变为NW-SE向,在中部自西向东由NE-SW向转变为近E-W向,而在南部既有NE-SW向又有NW-SE向;佳木斯地块各向异性方向由西部的NW-SE转变为东部的NNW-SSE,同时快慢波时间延迟逐渐变大;长白山造山带北部自北向南由NW-SE向逐渐转变为近E-W向,中部各向异性快波方向为NNW-SSE向,且快慢波时间延迟较大,而南部以NW-SE向为主;燕山造山带的各向异性快波方向主要沿E-W向分布,基本平行于燕山造山带的走向.这些结果说明,尽管复杂的各向异性快波方向与局部岩石圈拆沉和热物质上涌有关,但更重要是与"大地幔楔"中物质水平流等动力过程密切相关,也有待将来结合更多地震资料如面波不同深度的特征各向异性进行分析.在阿巴嘎火山群、哈拉哈火山群、长白山火山、龙岗火山和镜泊湖火山区及五大连池火山区等特殊构造区的周边地区,各向异性快波方向围绕这些构造区随方位均发生明显变化,暗示了火山区下方热物质上涌可能影响了"大地幔楔"中的软流圈物质水平流方向.