地球物理学报  2021, Vol. 64 Issue (5): 1585-1596   PDF    
利用背景噪声和接收函数研究华北克拉通地壳结构
姜磊1, 丁志峰1,2, 高天扬1, 黄翔1,3     
1. 中国地震局地球物理研究所, 北京 100081;
2. 中国地震局震源物理重点实验室, 北京 100081;
3. 北京大学地球与空间科学学院, 北京 100871
摘要:收集华北克拉通地区188个宽频带流动台站观测资料进行处理.通过背景噪声面波数据和接收函数双重资料约束联合反演,得到了研究区沉积层厚度、地壳厚度及地壳S波速度结构.结果显示:(1)沉积盖层厚度与地质构造相对应,盆地区与隆起区分界明显.(2)研究区地壳厚度变化范围约29~46 km,自西向东逐渐变薄.(3)中、上地壳华北盆地S波速度偏高,可能与新生代以来多次沉降所造成的相对高的岩石强度有关;(4)下地壳S波速度显示研究区主要存在三个低速区,分别是唐山—天津周边、张北及太行山造山带地区;华北盆地存在显著高速异常,推测可能是由于华北盆地经历下地壳拆沉后,大规模的伸展作用相伴随的幔源基性铁镁质岩浆底侵至下地壳结晶所造成的.(5)多个发生过强震的区域表现出沉积层下方存在较大范围的(约10 km)高速体,并且高速体又被其下低S波速度包裹,壳内岩石强度的差异为应力积累及地震发生提供条件.
关键词: 联合反演      背景噪声      接收函数      华北克拉通      构造演化     
Crustal structure beneath the North China Craton from joint inversion of ambient noise and receiver function
JIANG Lei1, DING ZhiFeng1,2, GAO TianYang1, HUANG Xiang1,3     
1. Institute of Geophysics, China Earthquake Administration, Beijing 100081, China;
2. Key Laboratory of Seismic Observation and Geophysical Imaging, China Earthquake Administration, Beijing 100081, China;
3. School of Earth and Space Sciences, Peking University, Beijing 100871, China
Abstract: We collected waveforms from 188 broadband seismic stations deployed in the North China Craton (NCC). Based on these data, we obtained the sediment thickness, crustal thickness, and S-wave velocity structure of the crust in the study area through the joint inversion of ambient noise surface wave data and receiver function data. The results demonstrate: (1) The sediment thickness corresponds to the geological structure. (2) The crustal thickness changes from 29 km to 46 km, thinning sharply when moving northwest to southeast. (3) In the middle and upper crust, the North China Basin (NCB) shows slight high-speed anomalies, which may be related to the relatively high rock strength caused by multiple subsidence since the Cenozoic. (4) In the lower crust, the low velocity anomalies are found around the Tangshan-Tianjin, Zhangbei and Taihang Mountain orogenic belt. A significantly high-velocity anomaly is observed beneath the NCB, and we speculate that it may be caused by the underplating and crystallization of the basic mafic magma from the mantle, accompanied by large-scale extension after the lower crust detachment. (5) Located in several regions with strong earthquakes, there is a certain range of (~10 km) high-speed bodies under the sedimentary layer, where the high-speed bodies are wrapped by the low S-wave velocity under them. The difference of rock strength in the crust provides conditions for stress accumulation and earthquake occurrence.
Keywords: Joint inversion    Ambient noise    Receiver function    North China Craton    Tectonic evolution    
0 引言

华北克拉通是世界上最古老的克拉通之一,大量研究表明,大约18.5亿年前,华北克拉通由东、西两个太古宙陆块碰撞拼合而成,随后,华北克拉通基底经历了一段相对稳定的沉积过程,覆盖了较厚的古生代沉积物(Wang and Mo, 1995Zheng et al., 2017).进入中生代后华北克拉通开始遭受强烈破坏,在侏罗-白垩纪期间受制于古太平洋俯冲的燕山运动对华北克拉通的破坏起了重要作用(Wang and Li, 2008朱日祥等, 2011, 2012, 2019Wang et al., 2015a, 2017b),在该过程中以北西-南东向挤压为主,岩浆活动、火山活动频繁,形成了各种断层、褶皱.在新生代,受控于东侧太平洋板块俯冲和西侧印度—欧亚板块碰撞作用,发生了广泛的构造伸展变形和岩浆作用(Menzies et al., 2007刘俊来等,2008Li et al., 2018),孕育了大量的断陷及裂谷盆地.现今,华北克拉通周围被不同时期形成的造山带环绕(凌媛等,2017),北侧为中亚造山带,南侧为秦岭—大别造山带和祁连造山带.华北克拉通的破坏具有强烈的复杂性,即:东部块体破坏严重,西部块体相对保持稳定.诸多学者构建了多种模型来探讨其破坏机制,主要有地幔柱上涌模型、岩石圈或地壳拆沉模型等(乔彦超等,2013Wang et al., 2015b程华冬等,2016程华冬和黄金水,2017王永明等,2017).岩浆底侵和拆沉作用是两个相反的过程,主导了不同构造环境下的壳-幔演化过程(樊祺诚等, 2001, 2005),岩浆底侵作用主要是岩浆上涌导致岩石圈及地壳发生软化并在应力作用下造成破坏,而拆沉作用是由于重力失稳引起重力垮塌,造成地壳及岩石圈地幔断裂下沉破坏.然而何种机制起主导作用目前仍存在诸多争议.华北克拉通破坏的强烈非均匀性,表现在各块体在地貌特征,岩石物理、化学性质及深部结构存在强烈的非均匀性.整体来看,华北克拉通地壳厚度由西向东逐渐减薄(Wei et al., 2011危自根等,2015),并由地堑-地垒相间分布,该区大致可以被北东或北北东向构造带分为东、中、西部块体,块体边界大概与南北重力梯度带重合.华北克拉通东部是克拉通破坏最显著的区域,地质构造复杂,该区在新生代的裂陷作用下,发育多条张性断裂,并发生大规模沉降(Zheng et al., 2017).华北克拉通东部的华北盆地覆盖较厚的沉积层,可细分有冀中坳陷、沧县隆起、黄骅坳陷、埕宁隆起、济阳坳陷等次级构造单元(段永红等,2016),这些次级构造单元相间分布.与东部相比,中部和西部保持相对稳定,地壳及岩石圈较厚,沉积层较薄.地壳结构的差异反映了华北克拉通改造机制的差异,东部地壳整体破坏主要由中新生代以来东侧太平洋板块俯冲引起,而中、西部地壳改造可能是其非均匀性质构造带在长期演化过程中被连续的热-构造事件多次改造活化的叠加所造成的(陈凌等,2010危自根等,2015).华北克拉通也是我国地震频发地区,有史以来记录到23次7级以上的地震,这种强烈的地震构造活动也是其区别于其他稳定克拉通的主要特征.地震科学家们对地震的孕育条件及其发生的机制进行了探索,到目前为止仍是尚未解决的难题.前人研究表明在华北克拉通,脆性的上地壳是主要的孕震层也是一个多震层,上、下地壳分界面附近是地壳内部的脆韧转换带(段永红等,2016),由此可见地壳上地幔的非均匀性是导致地震产生的主要原因,所以对华北克拉通地壳上地幔精细结构的探测是我们揭示其孕震环境,探讨地震发生机理的重要手段.

长期以来华北克拉通的构造演化历程一直是热点问题,随着数字化地震观测系统的发展,利用各种地球物理手段针对华北地区地壳结构进行了大量的探测和研究,取得了显著的成果,包括接收函数、面波成像、人工地震测深、大地电磁测深(Yao et al., 2008房立华等, 2009, 2013陈凌等,2010Fang et al., 2010Zheng et al., 2011武岩等, 2011, 2014, 2018Tang et al., 2013Fu et al., 2015段永红等,2016凌媛等,2017王椿镛等,2017Ouro-Djobo等,2018Xu et al., 2018).近年来,由于华北地区科学台阵布设,为利用密集台阵开展壳幔结构成像提供了大量数据.然而,由于华北克拉通东部块体覆盖有较厚的沉积层,对该区地壳结构研究造成一定的困难,因而由不同的方法得到的地壳厚度,壳幔速度结构仍存在诸多差异.

近十年来,背景噪声成像方法得到广泛的应用.该方法是利用不同台站长时间的连续记录进行互相关得到经验格林函数,进而通过反演获得台站群(相)速度.很多学者利用背景噪声进行了诸多有意义的研究.例如,房立华等(2009, 2013)得到了华北克拉通瑞利面波(勒夫波)群(相)速度分布图.Feng等(2017)利用噪声干涉确定了华北地区的410 km和660 km不连续面.Tang等(2013)利用背景噪声和面波数据反演了华北地区壳幔结构,分析了华北盆地和鄂尔多斯块体的热、化学性质,并对其壳幔演化过程进行了探讨.然而,由于面波对S波速度结构敏感而对速度间断界面不敏感,故单独利用面波数据反演壳、幔速度结构有显著的非唯一性问题.因此很多学者开始尝试联合接收函数与面波频散数据反演S波速度结构(Shen et al., 2013Zheng et al., 2016郑晨等, 2016, 2018Wang et al., 2017a)并取得了很好的效果,这主要是利用接收函数对台站下方的速度界面的敏感性与面波数据对绝对速度的敏感性的互补.相比于地震面波,背景噪声对短周期瑞利面波有很好的反映,即能较好的反映中上地壳的S波速度(陈俊磊等,2019).因此我们利用噪声数据和接收函数联合反演能够更好地反映地壳浅部的速度结构特征.

本文所用数据源自于华北克拉通地区(111°E—120°E,37°N—42°N)布设的188个宽频带流动台站(图 1).我们通过对提取的背景噪声瑞利面波群速度和接收函数进行联合反演得到华北克拉通地壳S波速度结构,沉积层厚度及莫霍面深度,并对其构造演化过程进行分析.

图 1 研究区区域构造及台阵分布图 蓝色三角形代表流动台站位置,蓝色实线表示主要的断裂分布,红色实线为剖面K和L的位置,灰色实线为南北重力梯度带.主要次级构造单元有:JY:济阳坳陷,CN:埕宁隆起,HH:黄骅坳陷,CX:沧县隆起,JZ:冀中坳陷. Fig. 1 Map showing tectonic setting of the study area and station distribution The blue triangles are seismic stations, the blue solid lines denote major fracture distribution, the red solid line is the location of the cross sections K and L, and the gray line is north-south gravity gradient zone. The main secondary structural units are: JY: Jiyang depression, CN: Chengning uplift, HH: Huanghua depression, CX: Cangxian uplift, JZ: Jinzhong depression.
1 数据和方法 1.1 面波频散

我们收集了2007年1月到2007年12月华北科学台阵的188个宽频带流动台站(台站分布见图 1)的连续地震波形数据.地震计为CMG-3ESP,频带范围:50 Hz~60 s.根据Bensen等(2007)介绍的方法及步骤对垂直分量(经重采样:1 Hz、去仪器响应、带通滤波5~50 s、时域归一化)处理得到经验格林函数,筛选得到信噪比>7,台站间距>3倍波长的数据.利用姚华建等(2004)编制的多重滤波图像分析软件提取5~30 s瑞利面波群速度频散.得到频散曲线后,使用Ditmar和Yanovskaya(1987)Yanovskaya和Ditmar(1990)提出的面波层析成像方法反演得到的瑞利面波群速度分布,反演网格大小为0.25°×0.25°.从中提取各台站对应的瑞利面波群速度数据作为联合反演的面波频散数据.关于华北地区该数据已有诸多背景噪声瑞利面波成像的研究成果,如Fang等(2010)与本文使用数据一致且结果相似,因此本文不再罗列成像结果.

1.2 接收函数

采用来自于武岩等(2018)利用时间域反褶积计算得到的接收函数.对于不同震中距的接收函数,其Ps转换波与P波到时差不同,因此首先我们对其进行动校正到67°参考震中距并叠加得到平均接收函数(Yuan et al., 1997郑晨等,2018),然后截取前15 s的波形作为联合反演的接收函数数据.

1.3 联合反演

本文利用改进的NA算法(Sambridge,1999郑晨等, 2016, 2018)联合反演背景噪声面波和接收函数数据构建地壳S波三维速度结构,并得到沉积层厚度,莫霍面深度.该方法通过在模型空间进行随机搜索,计算200个模型对应的失配值,失配值由公式(1)计算,

(1)

式中,G表示接收函数及背景噪声面波数据,ω1ω2ω3ω4分别的代表计算失配值时各部分的权重,m是模型参数,L是Laplacian算子.反演过程中需要设置各部分权重,其中ω3ω4越大得到的结果越稳定但相应的数据失配值会增大,经过尝试我们分别将其设置为为1.0和0.3.对于ω1ω2,根据各台站波形拟合程度做出调整,即对于拟合较差的数据适当的调大其权重.

每次迭代选取失配值最低的5个模型,并再次均匀随机游走产生新的100个模型进行新一次迭代,整个过程共进行500次迭代.

反演过程中我们将台站下方划分为沉积层、地壳、上地幔三个部分,根据武岩等(2014, 据2018)对渤海湾地区沉积层厚度以及地壳厚度的研究,我么将反演沉积层厚度参数搜索范围设置为0~7 km,并且将莫霍面搜索范围设置为25~50 km.同时,固定各层vP/vS值,沉积层设置为2.0,地壳及地幔分别设定为1.75和1.795.在反演中利用样条曲线表示S波速度结构,将沉积层搜索厚度设置为0.1 km,地壳及上地幔顶部设置为2 km,其余上地幔部分设置为5 km(郑晨等,2016).各层密度根据Nafe-Drake速度-密度公式(Nafe and Drake, 1963)计算.

通过对各台站逐个进行反演发现,位于基岩上的台站反演的接收函数和面波频散数据与实际数据有很好的拟合(图 2a),但是部分台站,特别是华北盆地地区反演效果相对较差,这主要是由于接收函数数据和面波频散数据分辨率的不同,接收函数对沉积层十分敏感,而面波频散数据却无法有效的反映沉积层(Shen et al., 2013),即面波频散高频成分的速度是地壳浅层速度的平均效应,其无法反映沉积层极低的剪切波速度.对此类台站我们适当的截掉小部分高频成分使我们的结果更为准确(图 2b2c).

图 2 台站A103(a)、A605(b)、K046(c)联合反演结果 (a1, b1, c1)中红色实线失配值最低的S波速度结构,黑色背景为200组次优解;(a2, b2, c2)中黑色实线为真实接收函数,红色实线为反演拟合的接收函数;(a3, b3, c3) 中蓝色星形代表真实面波频散数据,红色实线为反演拟合的面波频散曲线. Fig. 2 Joint inversion results under station A103 (a)、A605 (b)、K046 (c) The red solid line in the left figure represents the best vS model with smallest misfit value while the black background for the best 200 models; In top right figure, receiver function data are shown in black solid curve and predicted waveform is in dashed line. In the bottom right, the blue stars represent real surface wave dispersion data and the solid lines denote predicted surface wave dispersion.
2 结果

通过插值,我们得到了研究区沉积层厚度,莫霍面深度及地壳S波速度结构.图 3a展示了研究区沉积层厚度分布,整个研究区沉积层厚度范围在0~7 km,以太行山为界,两侧反映出了不同的沉积层厚度分布.其中太行山西侧及燕山造山带区域沉积层厚度较小小于1 km,但是在山西地堑带沉积层厚度较大,约为2~3 km. 东侧华北盆地沉积层明显较厚大于3 km,局部地区可达到5~6 km.可以看出沉积盖层厚度与地质构造相对应,坳陷区沉积层厚度较大,隆起区沉积层厚度较小.

图 3 (a) 研究区沉积层厚度及(b)研究区地壳厚度 Fig. 3 Map of sedimentary thickness (a) and crustal thickness (b) in study region

地壳厚度成像结果与前人结果较为一致(嘉世旭等,2009葛粲等,2011Tang et al., 2013武岩等,2018),从图 3b可以看出地壳厚度变化范围约29~46 km,由西向东逐渐减薄.同样以太行山为界,华北克拉通东部华北盆地地壳厚度小于34 km,燕山造山带东侧地壳厚度小于36 km.而华北克拉通中、西部地壳厚度较大约为36~46 km,在山西地堑莫霍面局部隆起约36~39 km.研究区地壳最厚的地区位于清水河附近大于43 km,地壳最薄地区位于保定附近及渤海湾周边地区,莫霍面轻微隆起.

研究区0~5 km范围内的平均S波速度结构(图 4a)显示,华北盆地存在显著的低速异常,其中沧县隆起区S波速度相对略高于其两侧的冀中坳陷和黄骅坳陷.在南北重力梯度带以西地区,清水河—朔州地区呈现显著的高速异常,太行山和燕山造山带区域也显示高速特征,而山西地堑区域S波速度略低于造山带.整体来看,在0~5 km范围内,华北克拉通东部与中、西部地区速度分界明显,且速度结构基本与地表地质和构造带相对应,坳陷区速度相对较低,隆起区速度相对较高.

图 4 联合反演获得的研究区地壳平均S波速度分布 (a) 地壳顶部(0~5 km)平均S波速度分布;(b) 中、上地壳(5~20 km)平均S波速度分布;(c) 下地壳(20 km~莫霍面)平均S波速度分布. Fig. 4 Average of S-wave velocity distribution in the crustal from the joint inversion in study region (a) Average S-wave velocity distribution in uppermost crust (0~5 km); (b) Average S-wave velocity distribution in the middle and upper crust (5~20 km); (c) Average S-wave velocity distribution in the lower crust (20 km~Moho).

图 4b展示了研究区5~20 km范围的平均S波速度结构,该范围主要反映了研究区沉积层以下的中、上地壳的S波速度结构.在这一深度范围内的S波速度结构仍然基本与区内主要的构造带相对应,但是与结晶基底以上的S波速度结构存在显著的差异.研究区内主要的两条造山带(太行山造山带、燕山造山带)显示较低的S波速度.在研究区内可以非常清楚的识别出三个高速异常区域,分别是清水河—朔州—大同地区、天津—唐山地区以及华北盆地地区.

下地壳S波速度结构(图 4c)显示研究区下地壳存在多个明显的高速异常区域,分别是华北盆地区域、北京北部区域、清水河—朔州及大同区域.而在张北、太行山造山带及天津—唐山周边地区显示出低速异常,特别是天津—唐山周边地区低速异常显著.与中、上地壳速度结构相比,研究区下地壳在天津—唐山周边地区发生了明显的速度反转.

3 讨论

华北克拉通破坏不仅是发生在岩石圈地幔的破坏,而且也造成其下地壳甚至整个地壳的破坏.目前,对于华北克拉通破坏机制主要有拆沉作用和热-化学侵蚀作用.华北克拉通历史上发生过多次构造活动,尤其是中、新生代以来,发生了燕山运动和古太平洋板块的俯冲作用,使得华北克拉通地壳和岩石圈发生强烈的破坏,其表现在沉积层厚度、地壳厚度、S波速度结构等方方面面.南北重力梯度带是华北克拉通一条重要的分界线,其两侧经历了不同的构造演化过程.在本文的结果中,无论是沉积层厚度还是地壳厚度或者地壳S波速度结构都明显以该重力梯度带为界,其两侧表现出截然不同的特征.目前对该区已经有不少的研究成果,我们将对本文所得到的结果与前人的结果进行对比分析.

华北克拉通东部覆盖的巨厚沉积层给该区的研究带来一定的困难,诸多学者利用不同的方法对该区的沉积层厚度进行了研究(罗艳等,2008武岩等, 2014, 段永红等,2016).我们的结果与这些研究结果有着良好的对应,华北盆地沉积层较厚约为3~7 km,其中在冀中坳陷东侧可达到5~6 km.山西地堑带沉积层厚度大约2~3 km,太行山造山带与燕山造山带结晶基底几乎出露地表(沉积层厚度<1 km).研究区地壳厚度由西向东逐渐减薄,以南北重力梯度带为界,两侧厚度剧烈变化.诸多研究表明古太平洋大约于180 Ma年前开始向欧亚大陆之下俯冲(Wang et al., 2015a, 2017b),东部严重破坏减薄与古太平洋的俯冲及其引起的弧后拉张作用有重要关联(朱日祥等, 2011, 2012, 2019Wang et al., 2018).

研究区内浅层(0~5 km)S波速度结构主要受到低速的沉积层影响,其速度分布特点与我们所得到沉积层厚度分布结果相吻合,并基本与地表地质构造带相对应,华北盆地低速特征明显,山西地堑速度相对于造山带偏低,华北盆地与山西地堑厚的沉积层是其呈现相对低速的主要原因.同时,在华北盆地中的次级构造带的隆起区与坳陷区也存在差异,速度分界较为明显,房立华等(2013)也得到了相似的结果,这也反映了不同次级构造单元结晶基底埋深的差异性.

研究区沉积层以下中、上地壳(5~20 km),造山带和盆地区的速度结构相对于结晶基底以上的速度结构显示出不同的特征,在华北盆地地区显示高速异常,而太行山造山带及燕山造山带显示低速异常,自新生代以来,整个华北盆地产生张裂和不均匀的沉降(曾融生等,1991),华北盆地中、上地壳相对高速可能与沉降所造成的相对高的岩石强度相关.清水河—朔州及大同地区显示高速异常并一直延伸到下地壳,该区高速异常在其他学者相关研究中也都有体现(房立华等, 2009, 2013Wang et al., 2012Tang et al., 2013),该区玄武岩广泛出露,推测可能与玄武质岩浆冷却凝固有关.

下地壳是壳幔物质交换的重要纽带,本文结果显示S波速度结构较为复杂,并且研究区内不同块体间速度对比相比于中、上地壳更为强烈.华北盆地区域显示强烈的高速异常,Tang等(2013)Fu等(2015)也得到了华北盆地下地壳显示高速异常的结果.然而一些研究成果显示华北盆地地区在下地壳为低速异常,他们将其解释为岩浆侵入造成下地壳大范围部分熔融,从而导致下地壳低速,但是泊松比研究显示华北东部地区内没有出现异常高的泊松比(ν>0.31), 暗示地壳内不存在大规模的部分熔融(Owens and Zandt, 1997),并且华北东部块体南部泊松比值比较低,这被认为是由拆沉作用造成的(许卫卫等,2005刘俊来等,2008武岩等,2018).岩石学证据也表明华北克拉通东部伸展与拆离作用强烈,存在大规模的变质核杂岩(Darby et al., 2004Liu et al., 2006).结合泊松比特征及岩石学证据,我们推测拆沉作用可能是华北克拉通东部地壳强烈破坏的主导因素,该区下地壳高速异常可能是由于华北盆地在古太平洋俯冲挤压后经历不同程度的下地壳拆沉作用,随后大规模的伸展作用相伴随的幔源基性铁镁质岩浆底侵至下地壳结晶所造成的. 太行山造山带及以西地区显示低速异常,武岩等(2018)研究显示太行山山前断裂的泊松比可达到0.30以上,Ji等(2009)研究显示张北—张家口—怀来、保定—繁峙及太行山一带为岩浆底侵区域,与我们得到的下地壳低速异常区形态基本吻合(图 4b),根据其华北克拉通伸展分区模式,推论在太行山隆起区下地壳存在岩浆底侵造成区域部分熔融.唐山地震区位于燕山造山带与华北盆地的交界部位,我们的结果显示该区下地壳存在显著的低速异常,暗示该区存在地幔物质上涌通道,对于该区下地壳低速异常在其他学者的相关研究中也都有体现,地幔热物质向地壳运移被认为是导致该区应力累积和发生地震的重要因素(曾融生等,1991刘启元等,2007段永红等,2016).

为了进一步了解研究区壳内高、低速体与构造之间的关系,我们沿K测线(图 5a)和L测线(图 5b) 绘制两条速度结构垂直剖面图.王未来等(2009)利用双台法得到的面波频散曲线与接收函数联合反演得到了K测线0~100 km的S波速度结构,与其相比我们利用噪声面波频散数据包含更短的周期,对测线下方地壳浅部有更好的分辨能力.结果显示测线下方速度变化较好的勾画了沉积层和莫霍面形态,由西向东沉积层厚度增大、莫霍面深度减小.在20 km以上范围内,S波速度较低,主要反映中、上地壳的速度特征.在下地壳内存在多个低速体(20~40 km),其中张北附近地区下地壳低速体最为明显,岩石学和地球化学观测显示在张北地区玄武岩捕掳体广泛出露(Zhang et al., 2018),三维导电性模型显示该区下地壳高导电性异常(Zhang et al., 2016),这些证据意味着张北附近地区存在强烈的壳幔热化学作用造成下地壳部分熔融.在唐山—唐海、三河—宝坻、商都以西及繁峙地区也存在一定范围的下地壳低速体,在前文中已指出这些区域可能存在一定范围的岩浆底侵作用或者存在幔源物质上涌通道,从而导致S波速度低速异常.

图 5 沿图 1中L剖面(a)和K剖面(b)S波速度模型垂直切片 绿色五角星圆圈从左至右依次表示:张北—尚义6.2级地震,三河—平谷8.0级地震和唐山7.8级地震. Fig. 5 Cross sections of vS model along L profile (a) and K profile (b) in Fig. 1 From left to right, the earthquake marked in green stars are: Zhangbei-Shangyi M6.2 earthquake, Sanhe-Pinggu M8.0 earthquake, and Tangshan M7.8 earthquake.

唐山地震区一直是地球物理学家们研究的重点区域,梅世蓉等(1999)根据唐山地震区层析成像结果提出震源区的中、上地壳存在较大规模的高速体,而四周速度较低的孕震模式.在这种模式下,中、上地壳由于岩石强度存在差异在受到水平应力场和垂向应力场的作用下更容易积累应力进而产生破裂.图 5a显示,在唐山地震区表现出在沉积层下存在较大范围的高速体,并且该高速体又被其下低S波速度包裹的特征,表明唐山地区符合这种孕震模式.除了唐山地震区外,张北地区及太行山与华北盆地块体交界地带下方速度结构也存在这种特征,历史上在这两个区域附近也曾分别发生过1998年张北6.2级地震和1679年三河—平谷8级地震(震源深度不详),该特征是否具有广泛的代表性仍有待进一步研究.

华北克拉通尤其是东部地区地壳结构复杂,地堑-地垒构造交替分部,强烈的不均匀性是导致地震发生的直接原因.面波和接收函数是深入了解地壳结构的有效手段.在研究区范围内发生过多次7级以上大震并造成了巨大的人员财产损失,尤其是唐山地震.通过开展地壳精细结构的研究有利于我们探索其孕震机制.我们利用华北地区188个台站连续1年的波形数据互相关并用频时分析得到噪声瑞利面波数据并联合接收函数反演得到了华北克拉通沉积层厚度、莫霍面深度及地壳S波速度结构,并与前人的结果进行对比分析,讨论了华北克拉通破坏机制及其孕震环境.具体总结如下:

(1) 华北克拉通沉积层厚度坳陷带和隆起带分界明显,其中华北盆地存在巨厚的沉积层,并且在冀中坳陷东部达到最大(5~6 km),山西地堑也显示较厚的沉积层.地壳由东向西逐渐增厚(29~46 km),与中、西部地区相比,东部华北盆地地壳厚度<34 km,表示该区域可能经历了强烈的伸展减薄作用,严重破坏减薄与古太平洋的俯冲引起的弧后拉张作用有重要关联.

(2) S波速度结构显示,在地壳顶部速度特征与地形地貌、地壳厚度存在对应关系,华北盆地、山西地堑受到沉积层影响表现为低速异常.在中、上地壳范围内华北盆地的高速异常可能与该区新生代以来多次沉降有关.清水河—朔州及大同区域显示高速异常并延伸到下地壳,可能与玄武质岩浆冷却凝固有关.

(3) 在下地壳范围内,华北盆地地区显示显著高速异常.该区下地壳高速异常可能是由于华北盆地在古太平洋俯冲挤压后经历不同程度的下地壳拆沉作用,随后大规模的伸展作用相伴随的幔源基性铁镁质岩浆底侵至下地壳结晶所造成的.张北及太行山造山带表现出低速异常,暗示岩浆侵入造成下地壳部分熔融.唐山地震区及周边显示显著的低速异常,该区地震多发可能与地幔物质的运移有关.

(4) 剖面显示多个发生过强震的区域表现出沉积层下方存在较大范围的(约10 km)高速体,并且高速体又被其下低S波速度包裹的特征,壳内岩石强度的差异为应力积累及地震发生提供条件.该特征是否具有广泛的代表性仍有待进一步研究.

致谢  武岩博士提供了接收函数数据.郑晨及王兴臣博士在数据处理期间给予了指导.审稿专家及编辑给出了修改意见.中国地震局地球物理研究所“中国地震科学探测台阵数据中心”为本研究提供了地震波形数据.在此一并致谢.
References
Bensen G D, Ritzwoller M H, Barmin M P, et al. 2007. Processing seismic ambient noise data to obtain reliable broad-band surface wave dispersion measurements. Geophysical Journal International, 169(3): 1239-1260. DOI:10.1111/j.1365-246X.2007.03374.x
Chen J L, Zheng D C, Long F. 2019. Joint inversion of receiver function and ambient noise: research and application. Progress in Geophysics (in Chinese), 34(3): 862-869. DOI:10.6038/pg2019CC0218
Chen L, Wei Z G, Cheng C. 2010. Significant structural variations in the Central and Western North China craton and its implications for the craton destruction. Earth Science Frontiers (in Chinese), 17(1): 212-228.
Cheng H D, Huang J S, Fu R S. 2016. Numerical simulation on convective thinning of the cratonic lithosphere. Chinese Journal of Geophysics (in Chinese), 59(4): 1293-1308. DOI:10.6038/cjg20160412
Cheng H D, Huang J S. 2017. Numerical simulation on the detachement of eclogitic lower crust. Science China Earth Sciences (in Chinese), 47(1): 82-94. DOI:10.1360/N072016-00019
Darby B J, Davis G A, Zhang X H, et al. 2004. The newly discovered Waziyu metamorphic core complex, Yiwulv shan, western Liaoning Province, Northwest China. Earth Science Frontiers, 11(3): 145-155.
Ditmar P G, Yanovskaya T B. 1987. Generalization of Backus-Gilbert method for estimation of lateral variations of surface wave velocities. Izvestiya, Academy of Sciences, USSR. Physics of the Solid Earth, 23(6): 470-477.
Duan Y H, Wang F Y, Zhang X K, et al. 2016. Three-dimensional crustal velocity structure model of the middle-eastern North China Craton (HBCrust1.0). Science China Earth Sciences, 59(7): 1477-1488. DOI:10.1007/s11430-016-5301-0
Fan Q C, Sui J L, Liu R X, et al. 2001. Eclogite facies garnet-pyroxenolite xenolith in Hannuoba area: New evidence of magma underplating. Acta Petrologica Sinica (in Chinese), 17(1): 1-6.
Fan Q C, Zhang H F, Sui J L, et al. 2005. Magmatic underplating and crust-mantle boundary in Hannuoba-evidence from xenolith petrology and geochemistry. Science in China: Earth Science (in Chinese), 35(1): 1-14.
Fang L H, Wu J P, Lu Z Y. 2009. Rayleigh wave group velocity tomography from ambient seismic noise in North China. Chinese Journal of Geophysics (in Chinese), 52(3): 663-671. DOI:10.1002/cjg2.1388
Fang L H, Wu J P, Ding Z F, et al. 2010. High resolution Rayleigh wave group velocity tomography in North China from ambient seismic noise. Geophysical Journal International, 181(2): 1171-1182. DOI:10.1111/j.1365-246X.2010.04571.x
Fang L H, Wu J P, Wang W L, et al. 2013. Love wave tomography from ambient seismic noise in North-China. Chinese Journal of Geophysics (in Chinese), 56(7): 2268-2279. DOI:10.6038/cjg20130714
Feng J K, Yao H J, Poli P, et al. 2017. Depth variations of 410 km and 660 km discontinuities in eastern North China Craton revealed by ambient noise interferometry. Geophysical Research Letters, 44(16): 8328-8335. DOI:10.1002/2017GL074263
Fu Y Y, Gao Y, Li A B, et al. 2015. Lithospheric shear wave velocity and radial anisotropy beneath the northern part of North China fromsurface wave dispersion analysis. Geochemistry, Geophysics, Geosystems, 16(8): 2619-2636. DOI:10.1002/2015GC005825
Ge C, Zheng Y, Xiong X. 2011. Study of crustal thickness and Poisson ratio of the North China craton. Chinese Journal of Geophysics (in Chinese), 54(10): 2538-2548. DOI:10.3969/j.issn.0001-5733.2011.10.011
Ji S C, Wang Q, Salisbury M H. 2009. Composition and tectonic evolution of the Chinese continental crust constrained by Poisson's ratio. Tectonophysics, 463(1-4): 15-30. DOI:10.1016/j.tecto.2008.09.007
Jia S X, Zhang C K, Zhao J R, et al. 2009. Crustal structure of the rift-depression basin and Yanshan uplift in the northeast part of North China. Chinese Journal of Geophysics (in Chinese), 52(1): 99-110.
Li S Z, Zhao S J, Liu X, et al. 2018. Closure of the Proto-Tethys Ocean and Early Paleozoic amalgamation of microcontinental blocks in East Asia. Earth-Science Reviews, 186: 37-75. DOI:10.1016/j.earscirev.2017.01.011
Ling Y, Chen L, Wei Z G, et al. 2017. Crustal S-velocity structure and radial anisotropy beneath the southern part of central and western North China Craton and the adjacent Qilian Orogenic Belt from ambient noise tomography. Science China Earth Sciences, 60(10): 1752-1768. DOI:10.1007/s11430-017-9092-8
Liu J L, Guan H M, Ji M, et al. 2006. Late Mesozoic metamorphic core complexes: new constraints on lithosphere thinning in North China. Progress in Natural Science, 16(6): 633-638. DOI:10.1080/10020070612330045
Liu J L, Davis G A, Ji M, et al. 2008. Crustal detachment and destruction of the North China craton: constraints from Late Mesozoic extensional structures. Earth Science Frontiers (in Chinese), 15(3): 72-81. DOI:10.1016/S1872-5791(08)60063-9
Liu Q Y, Wang J, Chen J H, et al. 2007. Seismogenic tectonic environment of 1976 great Tangshan earthquake: results given by dense seismic array observation. Earth Science Frontiers (in Chinese), 14(6): 205-213. DOI:10.1016/S1872-5791(08)60012-3
Luo Y, Chong Y J, Ni S D, et al. 2008. Moho depth and sedimentary thickness in Capital region. Chinese Journal of Geophysics (in Chinese), 51(4): 1135-1145.
Mei S R, Xue Y, Yin J Y. 1999. Relationship between the earthquake sequences of Tangshan and Xingtai and the three-dimensional velocity structure: Discussing on predicting strong earthquakes of swarm type. Acta Seismologica Sinica (in Chinese), 21(2): 159-165. DOI:10.3321/j.issn:0253-3782.1999.02.006
Menzies M, Xu Y G, Zhang H F, et al. 2007. Integration of geology, geophysics and geochemistry: a key to understanding the North China Craton. Lithos, 96(1-2): 1-21. DOI:10.1016/j.lithos.2006.09.008
Nafe J E, Drake C L. 1963. Physical properties of marine sediments. //Hill M N ed. The Sea. New York: Interscience, 3: 794-813.
Ouro-Djobo S B, Wei W B, Ye G F, et al. 2018. Experiments of magnetotelluric observation network on North China and lithospheric conductivity structure from fast imaging method. Chinese Journal of Geophysics (in Chinese), 61(6): 2508-2524. DOI:10.6038/cjg2018L0695
Owens T J, Zandt G. 1997. Implications of crustal property variations for models of Tibetan plateau evolution. Nature, 387(6628): 37-43. DOI:10.1038/387037a0
Qiao Y C, Guo Z Q, Shi Y L. 2013. Thermal convection thinning of the North China Craton: Numerical simulation. Science China Earth Sciences, 56(5): 773-782. DOI:10.1007/s11430-013-4588-3
Sambridge M. 1999. Geophysical inversion with a neighbourhood algorithm-I. Searching a parameter space. Geophysical Journal International, 138(2): 479-494. DOI:10.1046/j.1365-246X.1999.00876.x
Shen W S, Ritzwoller M H, Schulte-Pelkum V. 2013. A 3-D model of the crust and uppermost mantle beneath the Central and Western US by joint inversion of receiver functions and surface wave dispersion. Journal of Geophysical Research: Solid Earth, 118(1): 262-276. DOI:10.1029/2012JB009602
Tang Y C, Chen Y J, Zhou S Y, et al. 2013. Lithosphere structure and thickness beneath the North China Craton from joint inversion of ambient noise and surface wave tomography. Journal of Geophysical Research: Solid Earth, 118(5): 2333-2346. DOI:10.1002/jgrb.50191
Wang C Y, Wu Q J, Duan Y H, et al. 2017. Crustal and upper mantle structure and deep tectonic genesis of large earthquakes in North China. Science China Earth Sciences, 60(5): 821-857. DOI:10.1007/s11430-016-9009-1
Wang H Z, Mo X X. 1995. An outline of the tectonic evolution of China. Episodes, 18(1-2): 6-16. DOI:10.18814/epiiugs/1995/v18i1.2/003
Wang W L, Wu J P, Fang L H. 2009. Crust and upper mantle S-wave velocity structure beneath Tanghai-Shangdu seismic array profile. Chinese Journal of Geophysics (in Chinese), 52(1): 81-89.
Wang W L, Wu J P, Fang L H. 2012. High resolution Rayleigh wave phase velocity tomography in northern North China. Geophysical Journal International, 189(1): 647-658. DOI:10.1111/j.1365-246X.2012.05381.x
Wang X C, Li Y H, Ding Z F, et al. 2017a. Three-dimensional lithospheric S wave velocity model of the NE Tibetan Plateau and western North China Craton. Journal of Geophysical Research: Solid Earth, 122(8): 6703-6720. DOI:10.1002/2017JB014203
Wang Y, Li H M. 2008. Initial formation and mesozoic tectonic exhumation of an intracontinental tectonic belt of the northern part of the Taihang Mountain belt, Eastern Asia. Journal of Geology, 116(2): 155-172. DOI:10.1086/529153
Wang Y, Zhou L Y, Zhao L J, et al. 2015a. Tectonic transformations in the north of eastern China during 170~150 Ma: Causal linkage to the rapid formation of the paleo-Pacific plate. //Anderson T H, Didenko A N, Johnson C L, et al eds. Late Jurassic Margin of Laurasia-A Record of Faulting Accommodating Plate Rotation. London: Geological Society of America, 513(18): 561-587.
Wang Y, Zhou L Y, Luo Z H. 2017b. Kinematics and timing of continental block deformation from margins to interiors. Terra Nova, 29(4): 253-263. DOI:10.1111/ter.12272
Wang Y, Sun L X, Zhou L Y, et al. 2018. Discussion on the relationship between the Yanshanian Movement and cratonic destruction in North China. Science China Earth Sciences, 61(5): 499-514. DOI:10.1007/s11430-017-9177-2
Wang Y M, Huang J S, Zhong S J. 2015b. Episodic and multistaged gravitational instability of cratonic lithosphere and its implications for reactivation of the North China Craton. Geochemistry, Geophysics, Geosystems, 16(3): 815-833. DOI:10.1002/2014GC005681
Wang Y M, Li S Z, Li X Y, et al. 2017. Destruction mechanisms of the north China Craton: a review from numerical simulations. Marine Geology and Quaternary Geology (in Chinese), 37(4): 137-150. DOI:10.16562/j.cnki.0256-1492.2017.04.009
Wei Z G, Chen L, Xu W W. 2011. Crustal thickness and Vp/Vs ratio of the central and western North China Craton and its tectonic implications. Geophysical Journal International, 186(2): 385-389. DOI:10.1111/j.1365-246X.2011.05089.x
Wei Z G, Chu R S, Chen L. 2015. Regional differences in crustal structure of the North China Craton from receiver functions. Science China Earth Sciences, 58(12): 2200-2210. DOI:10.1007/s11430-015-5162-y
Wu Y, Ding Z F, Zhu L P. 2011. Crustal structure of the North China Craton from teleseismic receiver function by the Common Conversion Point stacking method. Chinese Journal of Geophysics (in Chinese), 54(10): 2528-2537. DOI:10.3969/j.issn.0001-5733.2011.10.010
Wu Y, Ding Z F, Zhu L P. 2014. Sedimentary basin structure of the Bohai Bay from teleseismic receiver functions. Acta Seismologica Sinica (in Chinese), 36(5): 837-849. DOI:10.3969/j.isn.0253-3782.2014.05.008
Wu Y, Ding Z F, Wang X C, et al. 2018. Crustal structure and geodynamics of the North China Craton derived from a receiver function analysis of seismic wave data. Chinese Journal of Geophysics (in Chinese), 61(7): 2705-2718. DOI:10.6038/cjg2018L0244
Xu W W, Zheng T Y. 2005. Distribution of Poisson's ratios in the northwestern basin mountain boundary of the Bohai Bay Basin. Chinese Journal of Geophysics (in Chinese), 48(5): 1077-1084.
Xu X B, Zhao L, Wang K, et al. 2018. Indication from finite-frequency tomography beneath the North China Craton: The heterogeneity of craton destruction. Science China Earth Sciences, 61(9): 1238-1260. DOI:10.1007/s11430-017-9201-y
Yanovskaya T B, Ditmar P G. 1990. Smoothness criteria in surface wave tomography. Geophysical Journal International, 102(1): 63-72. DOI:10.1111/j.1365-246X.1990.tb00530.x
Yao H J, Xu G M, Xiao X, et al. 2004. A quick tracing method based on image analysis technique for the determination of dual stations phase velocities dispersion curve of surface wave. Seismological and Geomagnetic Observation and Research (in Chinese), 25(1): 1-8.
Yao H J, Beghein C, Van Der Hilst R D. 2008. Surface wave array tomography in SE Tibet from ambient seismic noise and two-station analysis-Ⅱ. Crustal and upper-mantle structure. Geophysical Journal International, 173(1): 205-219. DOI:10.1111/j.1365-246X.2007.03696.x
Yuan X H, Ni J, Kind R, et al. 1997. Lithospheric and upper mantle structure of southern Tibet from a seismological passive source experiment. Journal of Geophysical Research: Solid Earth, 102(B12): 27491-27500. DOI:10.1029/97JB02379
Zeng R S, Zhu L P, He Z Q, et al. 1991. A seismic source model of the large earthquakes in north China extensional basin and discussions on the genetic processes of the extensional basin and earthquakes. Acta Geophysica Sinica (in Chinese), 34(3): 288-301.
Zhang C, Yao H J, Liu Q Y, et al. 2018. Linear array ambient noise adjoint tomography reveals intense crust-mantle interactions in North China craton. Journal of Geophysical Research: Solid Earth, 123(1): 368-383. DOI:10.1002/2017JB015019
Zhang H Q, Huang Q H, Zhao G Z, et al. 2016. Three-dimensional conductivity model of crust and uppermost mantle at the northern Trans North China Orogen: Evidence for a mantle source of Datong volcanoes. Earth and Planetary Science Letters, 453: 182-192. DOI:10.1016/j.epsl.2016.08.025
Zhen J X, Song X D, Zhu L P. 2013. Crustal and uppermost mantle S velocity structure under Hi-CLIMB seismic array in central Tibetan Plateau from joint inversion of surface wave dispersion and receiver function data. Tectonophysics, 584: 209-220. DOI:10.1016/j.tecto.2012.08.024
Zheng C, Ding Z F, Song X D. 2016. Joint inversion of surface wave dispersion and receiver functions for crustal and uppermost mantle structure in Southeast Tibetan Plateau. Chinese Journal of Geophysics (in Chinese), 59(9): 3223-3236. DOI:10.6038/cjg20160908
Zheng C, Ding Z F, Song X D. 2018. Joint inversion of surface wave dispersion and receiver functions for crustal and uppermost mantle structure beneath the northern north-south seismic zone. Chinese Journal of Geophysics (in Chinese), 61(4): 1211-1224. DOI:10.6038/cjg2018L0443
Zheng D, Li H Y, Shen Y, et al. 2016. Crustal and upper mantle structure beneath the northeastern Tibetan Plateau from joint analysis of receiver functions and Rayleigh wave dispersions. Geophysical Journal International, 204(1): 583-590. DOI:10.1093/gji/ggv469
Zheng T Y, Duan Y H, Xu W W, et al. 2017. A seismic model for crustal structure in North China Craton. Earth and Planetary Physics, 1(1): 26-34. DOI:10.26464/epp2017004
Zheng Y, Shen W S, Zhou L Q, et al. 2011. Crust and uppermost mantle beneath the North China Craton, northeastern China, and the Sea of Japan from ambient noise tomography. Journal of Geophysical Research: Solid Earth, 116(B12): B12312. DOI:10.1029/2011JB008637
Zhu R X, Chen L, Wu F Y, et al. 2011. Timing, scale and mechanism of the destruction of the North China Craton. Science China Earth Sciences, 54(6): 789-797. DOI:10.1007/s11430-011-4203-4
Zhu R X, Xu Y G, Zhu G, et al. 2012. Destruction of the North China Craton. Science China Earth Sciences, 55(10): 1565-1587. DOI:10.1007/s11430-012-4516-y
Zhu R X, Xu Y G. 2019. The subduction of the west Pacific plate and the destruction of the North China Cration. Science China Earth Sciences, 62(9): 1340-1350. DOI:10.1007/s11430-018-9356-y
陈俊磊, 郑定昌, 龙飞. 2019. 基于接收函数与背景噪声联合反演的研究与应用. 地球物理学进展, 34(3): 862-869. DOI:10.6038/pg2019CC0218
陈凌, 危自根, 程骋. 2010. 从华北克拉通中、西部结构的区域差异性探讨克拉通破坏. 地学前缘, 17(1): 212-228.
程华冬, 黄金水, 傅容珊. 2016. 克拉通岩石圈对流减薄的数值模拟. 地球物理学报, 59(4): 1293-1308. DOI:10.6038/cjg20160412
程华冬, 黄金水. 2017. 榴辉岩下地壳拆沉的数值模拟. 中国科学: 地球科学, 47(1): 82-94. DOI:10.1360/N072016-00019
段永红, 王夫运, 张先康, 等. 2016. 华北克拉通中东部地壳三维速度结构模型(HBCrust1.0). 中国科学: 地球科学, 46(6): 845-856. DOI:10.1360/N072015-00447
樊祺诚, 隋建立, 刘若新, 等. 2001. 汉诺坝榴辉岩相石榴辉石岩-岩浆底侵作用新证据. 岩石学报, 17(1): 1-6.
樊祺诚, 张宏福, 隋建立, 等. 2005. 岩浆底侵作用与汉诺坝现今壳-幔边界组成: 捕虏体岩石学与地球化学证据. 中国科学D辑: 地球科学, 35(1): 1-14.
房立华, 吴建平, 吕作勇. 2009. 华北地区基于噪声的瑞利面波群速度层析成像. 地球物理学报, 52(3): 663-671.
房立华, 吴建平, 王未来, 等. 2013. 华北地区勒夫波噪声层析成像研究. 地球物理学报, 56(7): 2268-2279. DOI:10.6038/cjg20130714
葛粲, 郑勇, 熊熊. 2011. 华北地区地壳厚度与泊松比研究. 地球物理学报, 54(10): 2538-2548. DOI:10.3969/j.issn.0001-5733.2011.10.011
嘉世旭, 张成科, 赵金仁, 等. 2009. 华北东北部裂陷盆地与燕山隆起地壳结构. 地球物理学报, 52(1): 99-110. DOI:10.3969/j.issn.1672-7940.2009.01.020
凌媛, 陈凌, 危自根, 等. 2017. 利用噪声层析成像研究华北克拉通中-西部南段及邻区祁连造山带地壳S波各向异性速度结构. 中国科学: 地球科学, 47(10): 1202-1219. DOI:10.1360/N072017-00070
刘俊来, Davis G A, 纪沫, 等. 2008. 地壳的拆离作用与华北克拉通破坏: 晚中生代伸展构造约束. 地学前缘, 15(3): 72-81. DOI:10.3321/j.issn:1005-2321.2008.03.005
刘启元, 王峻, 陈九辉, 等. 2007. 1976年唐山大地震的孕震环境: 密集地震台阵观测得到的结果. 地学前缘, 14(6): 205-213. DOI:10.3321/j.issn:1005-2321.2007.06.026
鲁来玉, 何正勤, 丁志峰, 等. 2009. 华北科学探测台阵背景噪声特征分析. 地球物理学报, 52(10): 2566-2572. DOI:10.3969/j.issn.0001-5733.2009.10.015
罗艳, 崇亚军, 倪四道, 等. 2008. 首都圈地区莫霍面起伏及沉积层厚度. 地球物理学报, 51(4): 1135-1145. DOI:10.3321/j.issn:0001-5733.2008.04.022
梅世蓉, 薛艳, 尹京苑. 1999. 唐山、邢台地震序列特征与三维速度结构的关系——兼论强震群型地震的预测问题. 地震学报, 21(2): 159-165. DOI:10.3321/j.issn:0253-3782.1999.02.006
Ouro-Djobo S B, 魏文博, 叶高峰, 等. 2018. 华北大地电磁测深阵列观测实验与岩石圈导电性快速成像模型. 地球物理学报, 61(6): 2508-2524. DOI:10.6038/cjg2018L0695
乔彦超, 郭子祺, 石耀霖. 2013. 数值模拟华北克拉通岩石圈热对流侵蚀减薄机制. 中国科学: 地球科学, 43(4): 642-652.
王椿镛, 吴庆举, 段永红, 等. 2017. 华北地壳上地幔结构及其大地震深部构造成因. 中国科学: 地球科学, 47(6): 684-719. DOI:10.1360/N072016-00168
王未来, 吴建平, 房立华. 2009. 唐海-商都地震台阵剖面下方的地壳上地幔S波速度结构研究. 地球物理学报, 52(1): 81-89.
王永明, 李三忠, 李玺瑶, 等. 2017. 华北克拉通破坏动力学机制: 数值模拟的启示. 海洋地质与第四纪地质, 37(4): 137-150. DOI:10.16562/j.cnki.0256-1492.2017.04.009
危自根, 储日升, 陈凌. 2015. 华北克拉通地壳结构区域差异的接收函数研究. 中国科学地球科学, 45(10): 1504-1514.
武岩, 丁志峰, 朱露培. 2011. 利用共转换点叠加方法研究华北地区地壳结构. 地球物理学报, 54(10): 2528-2537. DOI:10.3969/j.issn.0001-5733.2011.10.010
武岩, 丁志峰, 朱露培. 2014. 利用接收函数研究渤海湾盆地沉积层结构. 地震学报, 36(5): 837-849. DOI:10.3969/j.isn.0253-3782.2014.05.008
武岩, 丁志峰, 王兴臣, 等. 2018. 华北克拉通地壳结构及动力学机制分析. 地球物理学报, 61(7): 2705-2718. DOI:10.6038/cjg2018L0244
许卫卫, 郑天愉. 2005. 渤海湾盆地北西盆山边界地区泊松比分布. 地球物理学报, 48(5): 1077-1084. DOI:10.3321/j.issn:0001-5733.2005.05.014
姚华建, 徐果明, 肖翔, 等. 2004. 基于图像分析的双台面波相速度频散曲线快速提取方法. 地震地磁观测与研究, 25(1): 1-8. DOI:10.3969/j.issn.1003-3246.2004.01.001
曾融生, 朱露培, 何正勤, 等. 1991. 华北盆地强震的震源模型兼论强震和盆地的成因. 地球物理学报, 34(3): 288-301. DOI:10.3321/j.issn:0001-5733.1991.03.004
郑晨, 丁志峰, 宋晓东. 2016. 利用面波频散与接收函数联合反演青藏高原东南缘地壳上地幔速度结构. 地球物理学报, 59(9): 3223-3236. DOI:10.6038/cjg20160908
郑晨, 丁志峰, 宋晓东. 2018. 面波频散与接收函数联合反演南北地震带北段壳幔速度结构. 地球物理学报, 61(4): 1211-1224. DOI:10.6038/cjg2018L0443
朱日祥, 陈凌, 吴福元, 等. 2011. 华北克拉通破坏的时间、范围与机制. 中国科学: 地球科学, 41(5): 583-592.
朱日祥, 徐义刚, 朱光, 等. 2012. 华北克拉通破坏. 中国科学: 地球科学, 42(8): 1135-1159.
朱日祥, 徐义刚. 2019. 西太平洋板块俯冲与华北克拉通破坏. 中国科学: 地球科学, 49(9): 1346-1356. DOI:10.1360/N072018-00282