欧亚大陆东缘存在一个巨型深断裂系统——基于卫星重力的新发现

通过对卫星重力数据的精细处理,发现欧亚大陆东缘存在一个规则分布的重力异常组合,这个异常组合与中国大陆上的已知断裂严格对应。SinoProbe探测计划实施的反射地震探测出这些断裂下方存在Moho破裂,幔源物质上涌,剩余重力异常,莫霍面破裂,高密度的幔源物质上涌和深断裂之间存在密切联系。依照这个联系规律,对已有断裂的延伸补充,显示出欧亚大陆东缘存在一个巨型深断裂系统,主断裂带南起中国广州向北延伸直至鄂霍次克海,绵延逾3000 km。沿主断裂带东侧近似等间距的平行分布着9条北东向断裂,北东向断裂向东入海至大陆边缘。采用数值方法模拟了断裂系统的形成过程,在菲律宾板块和太平洋板块近NW方向构造力的作用下,欧亚大陆东缘产生NNE向走滑断层,其东侧生成彼此平行,间距大体相等的派生断层。断层形成过程中,Moho破裂,大量地幔物质上涌,形成中国东部中生代以来的巨量火山岩浆活动,控制了内生金属矿床的分布。这个断裂系统近代还在活动,导致了郯城8.5级大地震和长白山的火山喷发。     

天山造山带及邻区孕震应力场及断裂带活动性特征的数值实验分析

天山造山带构造环境复杂,活动断裂带和强震分布广泛,且主要分布于阿尔泰山、天山、西昆仑—帕米尔弧形构造带上,尤以天山地区最为集中.迄今为止,天山造山带地区的主要断裂带的活动特征与孕震应力场特征之间的动力学机理尚未有清晰的认识.本文以GPS等实际观测数据为约束,建立有限元数值模型,计算了研究区域地壳形变、应力/应变积累速率、弹性应变能密度以及库仑应力变化率等关键因素.模拟计算结果显示地表速度场与研究区域实际GPS观测值基本一致,且主要断裂带上弹性应变能密度分布与实际地震活动性也基本吻合,验证了数值模型和结果的可靠性.结合最新的观测和数值模拟结果分析发现,研究区的断层和地震活动性主要受控于近南北向的主压应力,与主要观测特征相一致.同时,帕米尔高原北部边界带—塔什库尔干断裂（TKF）、天山造山带南边界的东侧—迈丹断裂（MDF）、兴地断裂（XDF）库仑应力增大明显,在未来强震发生的可能性较高,应予密切关注.     

基于钻孔应变观测资料分析远震造成的动态库仑应力变化

已有证据充分证明了大地震远震触发作用的存在.为认识其机理,需要深入探究大地震产生的远场动态应力对断层面应力及断层地震活动性的影响.四分量钻孔应变仪作为超宽频地壳应力监测的重要手段,可以直接记录地震波传来时造成的动态水平应变和应力,这些资料对动态库仑应力变化及其影响的研究有着重要意义,而迄今尚未得到足够重视和应用.本文通过对高台、通化台的四分量钻孔应变仪记录到的2018年太平洋地区4次大地震事件数据进行分析,识别出了P波、S波和面波等不同震相,运用应力花瓣方法分析了它们的特点,计算了它们各自引起的动态应力变化,展示了利用钻孔应变实际观测资料定量研究不同震级大小、不同震中距离的远震地震波在台站位置处的地壳造成动态应力变化的可行性.研究发现,台站位置处动态应力变化的主应力轴方向大体与震中方位角一致,并讨论了对台站附近断层面上动态库仑应力的影响.定量计算结果表明,这4次地震在高台、通化台附近断层面产生的动态库仑应力变化峰值在数百帕量级,低于一般认为的动态触发阈值,也的确没有观测到被动态触发地震活动.但是本文的思路和方法,对今后更大远震可能产生的触发效应提供了研究基础.     

2020年第八届青藏高原东部构造与地球物理研讨会(8-WTGTP)顺利召开

第八届青藏高原东部构造与地球物理研讨会(8-WTGTP)暨中国地球物理学会中国大陆动力学专业委员会和中国地球物理学会固体地球物理委员会2020年学术年会,于2020年9月17日至11月24日顺利召开.由于疫情影响,今年的会议以网络直播的形式,分别在9月17日、9月24日、10月13日、10月28日、11月16日和11月24日召开6场学术报告会.会议由中国地球物理学会中国大陆动力学专业委员会、中国地球物理学会固体地球物理委员会、中国地震局地震预测重点实验室(中国地震局地震预测研究所)、中国科学院地球与行星物理重点实验室(中国科学院地质与地球物理研究所)、中国科学院计算地球动力学重点实验室(中国科学院大学)、中国地震学会构造地貌专业委员会、中山大学地球科学与工程学院和甘肃省地震学会联合主办(排名不分先后).     

华北地区水位下降是否会减缓气温上升——浅部地温影响的数值模拟分析

华北地区由于长期持续的地下水过量开采,导致了大面积地下水位大幅下降,引发地面塌陷、地下水质污染等一系列地质环境问题,这些现象早已为人们所熟知和关注.然而地下水位下降还会造成百米量级浅部地温及其梯度的变化,因此即使来自地球深部的大地热流密度没有变化,年度平均的从表浅部位通过地表实际传导进入大气的热流密度会减小,这是中外文献中尚未见讨论过的问题.我们通过数值模拟发现假定大地热流密度不变的条件下,华北数万平方公里地下水位下降会造成百米尺度内的地温降低,从而传入大气的热流密度降低40%以上,且会持续数百年以上的时间.这种长时间大范围的传导入大气的热流密度变化对环境会造成什么影响是一个十分值得关注的问题.这一预测在一定程度上得到了气象站地温观测数据的支持,但由于目前气象观测站只有3.2m深度范围内的地温资料,累计不超过5、60年,中间还有10余年的间断,而且表浅深度地温受地表多种因素的影响也较大,这些资料难以对我们关心的地下水位下降引起流入大气的热流密度变化这一问题提供直接确凿的数据来进行分析,因此今后有必要开展对地下数十乃至数百米地温进行持续精确的监测工作.     

气候模式模拟中海量数据的并行抽取与可视化分析诊断

在地球动力学和气候模拟等领域, 数值模拟产生的数据规模达到Tb至Pb量级。 实现这些海量数据的实时可视化和实时诊断分析面临很大的困难, 时空多尺度数据抽取可以解决这一瓶颈。 高精度数据的可视化结果展示需要高分辨率的显示设备, 并行的大屏幕显示技术是解决这一问题的有效手段。 地球科学数值模拟、 并行数据抽取和高分辨率显示都需要搭建高性能计算机集群。 本文在基于Lagrange插值的多维度、 多尺度、 多分辨率并行数据抽取算法的同时, 利用并行计算节点及LCD显示器, 基于Rocks cluster系统搭建起一个176核, 4×10×1024×1280分辨率的高性能计算模拟、 数据抽取和并行显示输出的集成平台, 并将该平台成功应用到气候模式模拟产生的海量数据的并行抽取和并行显示。     

固体潮应变各向异性研究

钻孔应变仪的观测过程中发现潮汐应变响应存在方向各向异性, 推测是受到台站附近断层影响导致的。 我们对此展开数值模拟研究, 讨论软弱断层对潮汐变形传播的影响。 我们的研究结果证实了软弱断层的存在确实能够产生能够观测到的各向异性现象。 同时我们还发现其各向异性变化对断层弹性参数较为敏感。 观测到的部分地震震前反映固体潮应变各向异性的潮汐因子"玫瑰图"变化有可能是由于台站附近断层在震前的性质变化引起的。     

汶川大地震地震波传播的谱元法数值模拟研究

汶川大地震震中位于震害活跃的龙门山断裂带上,其强度超过1976年唐山大地震．地震激发的地震波可以用来研究震源破裂过程及地球内部的层次结构．从弹性波理论基础出发,采用AK135理论地球模型,考虑地表地形、地球介质衰减及地球椭率等特性,利用谱元法结合高性能并行计算,分别对人工点源和复合源所激发的汶川大地震地震波的全球传播过程进行数值模拟重现．数值模拟结果显示了地震波在地球表面的传播形态．通过比较发现,复合源数值模拟结果较点源数值模拟结果更能体现汶川地震震源破裂过程的时空特性．另外,将两种震源的数值模拟结果分别与实际观测台站记录波形资料进行拟合对比,进一步定性地认识到多个点源组成的复合源的数值模拟结果较单个点源的数值模拟结果与实际观测资料的拟合程度更好,证明了汶川大地震破裂过程为多阶段破裂组成的一个复合破裂过程．     

虚拟川滇——基于千万网格并行有限元计算的区域强震演化过程数值模型设计和构建

利用合理的数值模型,对区域性的强震活动性分析,是近些年来地震学研究领域的一个热点问题.但无论从时间还是空间角度来讲,区域性强震演化的理论与数值模拟模型都是多尺度问题,对该问题的充分研究需要借助于大规模并行数值模拟技术.侧重于基于千万有限元网格技术,采用大规模并行有限元计算手段,对该问题的理论建模与并行数值算法实现,以及如何应用到对川滇地区区域强震演化过程的模拟进行了具体论述.本研究的关键目标是为中长期地震危险性估计提供一个大规模数值预报模拟平台系统.     

���������������¶�λ�������

????????λ???????????(111??120??E??35??42??N)1993??2004??6 771?ε???1.0??M??6.6???????λ???????λ???????????????????????????????????????????????;???o?????·???????????????????????????б???????????????????????????????????20 km????60 km???????25 km??34 km?????40 km;83%???????????λ??0??15 km??Χ?????????????λ???????????????(7??21 km)????????????????????????????????????????????????;???92%???????????λ??1??24 km??????24 km??????????????????????硣