中国地震  2019, Vol. 35 Issue (2): 286-294
南襄盆地及邻区地壳厚度与泊松比研究
余尚江, 成万里, 陈波, 卢亚, 陈贤, 黄恩贤     
河南省地震局, 郑州 450000
摘要:利用河南、湖北、陕西地震台网共27个固定台站记录到的2014~2016年远震波形的接收函数,通过叠加获得了南襄盆地及邻区26个台站下方的地壳厚度和泊松比,结合该区地形地貌特征、断裂分布和地震活动,得出以下结论:①南襄盆地及邻区的地壳厚度与地形变化关联紧密,呈现出由NE向SW方向增厚的趋势,地壳厚度为30~41km;②泊松比高值集中在房县-兴山附近,为0.31~0.34,这可能与区域内大量出露寒武-奥陶系的辉长岩、苏长岩等铁镁质基性岩相关;低值为0.18~0.24,分布在商南-西峡-淅川一带,与区域内存在大量的长英质岩有关,它们与南襄盆地内的新生代沉积地层形成整合接触关系。由此可推断,在淅川、西峡一带存在因介质的较大差异而导致地震射线在传播过程中出现反射的可能。
关键词南襄盆地    接收函数    地壳厚度    泊松比    
Study on Crustal Thickness and Poisson's Ratio in Nanxiang Basin and Its Adjacent Area
Yu Shangjiang, Cheng Wanli, Chen Bo, Lu Ya, Chen Xian, Huang Enxian     
Henan Earthquake Agency, Zhengzhou 450000, China
Abstract: In this paper, the receiver functions of the teleseismic waveforms from 2014 to 2016 recorded by 27 fixed seismic stations of Henan Digital Seismic Network, Hubei Digital Seismic Network, and Shanxi Digital Seismic Network were used to obtain the crust thickness and Poisson's ratio of the stations below the 26 seismic stations in the Nanxiang Basin and its adjacent area, which combined with the fracture distribution, seismicity, and topographic features of the area, has reached the following conclusions:(1)The crust thickness of the Nanxiang Basin and its adjacent areas is closely related to terrain change. It shows a trend of thickening from the northeast to the southwest, and the crust thickness is 30-41km.(2)The high value area of Poisson's ratio is located in the vicinity of Fangxian and Xingshan, the value of which is from 0.31 to 0.34. It may be associated with a large number of mafic rocks in the Cambrian-Ordovician, such as gabbro and sulphite. The range of low value is from 0.18 to 0.24, and the area distributes in Shangnan and Xixia, and is related to a large number of felsic rocks. They formed an integrated contact with Cenozoic sedimentary strata in the Nanxiang Basin. From this, it can be inferred that there is a possibility of reflection of seismic rays during the propagation process due to large differences in the medium between the Xichuan and Xixia areas.
Key words: Nanxiang Basin     Receiver functions     Crust thickness     Poisson's ratio    
0 引言

秦岭-大别构造带横贯中国大陆中部,分割着华北与华南2大地震带,与活跃的龙门山断裂带和郯庐断裂带紧密相连,地震活动具有频繁、震级小、震源较浅的特点,位于该构造带中部的南襄盆地,便具有这种特点。由于特殊的地质构造背景,南襄盆地在震相发育上有独特之处:当盆地内的淅川、西峡等地发生地震时,在相隔30km的南阳台(NY)的记录中,出现直达P波、S波的反射波,震相特征非常明显,而从震源往东,在南襄盆地内部的内乡、镇平、南阳、唐河等地震相正常,均未出现反射波震相(王林昌等,2011)。需要说明的是,这2个区域地震的震中距均为30~100km。在邻近台站,如卢氏台(LS)、平顶山台(PDS)等,记录到的南襄盆地震震相也存在这种现象。由震相基本理论可知,是否出现反射波应与发震区的地壳介质相关。若能探寻上述地区的地壳结构及物质组成信息,即可对解释南襄盆地发育特殊震相提供重要帮助。

在地震学研究中,可以利用地震P波和S波来推断地震波速度比vP/vS,进而能够测定出泊松比σ和地壳厚度,而泊松比和地壳厚度的给出就可以对研究区的物质组成及壳幔结构进行很好的约束(许卫卫等,2005)。泊松比是了解地球内部物质构成的一个非常重要的参数,其变化范围介于-1.0~0.5之间。泊松比与岩石组分之间有着密切的关系,绝大多数岩石的泊松比为0~0.5,低泊松比可能预示着地壳物质组成以长英质成分为主,高泊松比则表明富含铁镁质成分(嵇少丞等,2009)。若物质中存在流体或熔体则会提高岩石的泊松比,并且温度、压力等因素对泊松比的影响很小(Zandt et al,1995Christensen,1996Ji et al,2002)。另外,地壳厚度的给出可以大致反映该区域的壳幔结构。一般而言,同一地区地壳厚度基本恒定。若某个小区域地壳厚度增加,则可认为是地幔物质侵入了地壳,造成了地壳物质的增生;反之,地壳厚度变薄,则意味着该区域地壳发生了后期改造。

目前,使用远震波形的接收函数计算泊松比和地壳厚度的研究结果较多(吴庆举等,1998Salah et al,2008Yang et al,2011Tugume et al,2012Jo et al,2013黄海波等,2014高延光等,2014齐刚等,2015宫猛等,2015Wei et al,2016),这些结果主要是利用泊松比去反演地壳物质组成,同时利用地壳厚度去指示研究区地壳是否发生增生或拆沉作用。另外,也有部分结果(Salah et al,2011郑勇等,2013沈玉松等,2013钱辉等,2015龚辰等,2016邓嘉美等,2016)使用泊松比和地壳厚度的横向变化去追寻活动断裂和大地震孕育环境等,一般来说,泊松比或地壳厚度急剧变化区域很有可能成为孕育大地震和活动断裂的位置。

因此,本文旨在利用H-k叠加搜索的方法计算南襄盆地及邻区的泊松比和地壳厚度,进而试图对南襄盆地地壳结构和物质组成等信息进行约束,以期为南襄盆地发育特殊震相寻找可靠证据。

1 方法原理

基于宽频带地震计记录到的远震P波波形(包含壳-幔速度间断面的Ps转换波等信息),为有效地提取台站下方壳、幔结构的响应信息,消除震源时间函数和传播路径的影响,在对数据去除仪器响应、倾斜、均值及带通滤波之后,需用远震P波的垂直分量对其径向分量和切向分量作反褶积后得到时间序列,即P波接收函数。

在P波接收函数反映的远震记录中,经过一个速度间断面到达地表的初至波有多个研究所需要的后续震相,它们分别为Ps、PpPs、PpSs+PsPs震相。这几个续至波与初至P波之间的到时差分别为(Zhu et al,2000)

$ t_{\mathrm{P s}}=H *\left[\sqrt{\left(v_{\mathrm{P}} / v_{\mathrm{S}}\right)^{2}-p^{2} v_{\mathrm{P}}^{2}}-\sqrt{1-p^{2} v_{\mathrm{P}}^{2}}\right] / v_{\mathrm{P}} $ (1)
$ t_{\mathrm{PpPs}}=H *\left[\sqrt{\left(v_{\mathrm{P}} / v_{\mathrm{S}}\right)^{2}-p^{2} v_{\mathrm{P}}^{2}}+\sqrt{1-p^{2} v_{\mathrm{P}}^{2}}\right] / v_{\mathrm{P_S}} $ (2)
$ t_{\mathrm{PpSs}+\mathrm{PsPs}}=2 H * \sqrt{\left(v_{\mathrm{P}} / v_{\mathrm{S}}\right)^{2}-p^{2} v_{\mathrm{P}}^{2}} $ (3)

其中,H为地表到间断面之间的厚度;vpvS分别为该介质层的平均P波、S波速度;p为射线参数(对于远震,P波在接收区近垂直入射和地球半径已知的前提下,根据pr · sinλ/v可求射线参数p)。

给定地壳平均速度vP,在厚度H-波速比(k=vP/vS)平面上,上述2组曲线是2条斜率不同、相交于一点的曲线。该交点给出了这一介质层的厚度和波速比。由此可以在H-k平面上搜索2个震相叠加后的能量最大值,并得到间断面埋深厚度及波速比。由波速比与泊松比的关系即可求得测定台站下方区域内地壳的泊松比

$ \sigma=\frac{1}{2}\left\{1-\left[\left(v_{\mathrm{P}} / v_{\mathrm{S}}\right)^{2}-1\right]^{-1}\right\} $ (4)
2 数据处理计算

此次共收集了27个数字地震台站(图 1)的观测数据,其中,河南台网12个,湖北台网14个,陕西台网1个,覆盖面积达3.2万km2。我们从2014~2016年台站记录到的大量地震事件中挑选出震中距为30°~90°、具有清晰P波初动和高信噪比的远震波形。根据上述要求,共选出了145个地震事件,它们主要分布在环太平洋地震带上(图 1);然后根据地震发生时间对原始波形进行截取,截取过程中利用PREM模型计算理论到时以协助直达P波的识别,截取的时间区间为直达P波到时的前20s和后100s;使用SAC软件对截取的数据去除仪器响应、去倾斜、去均值,后使用频率为0.05~2.00Hz的带通滤波器进行滤波;在检查射线参数和各分量方位角后,将波形旋转至R、T、Z方位;利用时间反褶积法提取各远震事件的接收函数波形(Owens et al,1984Langston,1979);对提取结果做人工检查,剔除信噪比低、波形异常的接收函数,最终获得2577条有效记录,除河南许昌台(XC)有效接收函数记录数为35条外,其他台站均大于80条。

图 1 测震台站分布(a)及接收函数计算所用到的远震事件震中分布(b)

利用H-k叠加搜索法对这些台站下方的地壳厚度和波速比进行搜索,地壳平均P波速度设为6.3km/s(杨冬等,2008);叠加过程中转换S波震相Ps及其多次波PpPs、PpPs+PpSs的加权系数分别为0.7、0.2、0.1,地壳厚度控制在20~60km(王椿镛等,2017),波速比为1.5~2.0(Ji et al,2008)(图 2)。

图 2 部分台站接收函数叠加(a)及相应台站H/k叠加值(b) 图(a)、(b)为台站NY;图(c)、(d)为台站PDS;图(e)、(f)为台站XSH

除竹溪台未得到有效的接收函数外,此次通过H-k叠加搜索软件下的子程序iter_decon计算共获得了26个台站下方的地壳厚度H、波速比k及其相应的误差δH、δk(表 1)。由表 1可见,地壳厚度H高值区域主要集中在兴山-竹山-商南区域,厚度为38~41km;低值区域分布在平顶山-驻马店-南阳地区,厚度为30~32km。南襄盆地内部的南阳、唐河、丹江、襄樊分别为31.4、34.1、35.4、33.5km,厚度差为±2km,整体呈现东北薄、西南厚的特征,其误差δH均小于1.5km,反映数据基本真实可靠。

表 1 研究区台站下方的地壳厚度H、波速比k及泊松比σ

从波速比k来看,高值区域位于房县、兴山一带,其值为1.9~2.0;低值区位于商南附近,其值为1.60,差异较大,误差δk均小于0.07,数据基本可靠。

3 分析讨论 3.1 地壳厚度分布

由于中国大陆演化复杂,有关大陆地壳结构、厚度的研究一直是地学界的焦点。冯锐(1985)利用Parker-Oldenburg位场反演方法对中国大陆重力、沉积及地震测深等多种数据进行三维重力反演,得到中国地壳厚度分布结果,其中,位于太行-雪峰与秦岭-大别交界位置的南襄盆及邻区的地壳厚度为32~40km;随着中国大陆人工地震工作的持续开展,众多人工测深结果陆续问世,滕吉文等(2002)基于大量的人工地震数据计算出了东亚大陆及周边海域的地壳厚度,南襄盆地及邻区的厚度为32~42km,分布规律上展现出与地形起伏一致的东低西高的特点。与此同时,随着中国大陆数字地震台网的观测运行,使用天然地震记录反演地壳结构成为可能。Li等(2014)He等(2014)Wei(2016)使用接收函数方法对分布在中国大陆的宽频带地震波数据进行叠加搜索,得到了中国大陆地壳厚度分布,从其结果可知南襄盆地及邻区的地壳厚度为30~42km,呈现东北薄、西南厚的特征。王椿镛等(2017)利用接收函数方法计算出了大华北地区的地壳厚度与泊松比,其中,从黄淮地区、南襄盆地到鄂西地区地壳厚度由30km逐渐增加到42km。

根据本文计算结果,通过插值得到南襄盆地及邻区的地壳厚度分布(图 3(a))。结果显示,南襄盆地及邻区的地壳厚度与地形变化间呈正相关,总体呈现出由NE向SW方向增厚的趋势,这与研究区处于中国大陆二、三级阶梯有密切的联系。根据厚度分布,大致可分为3个梯度区,其中,第1梯度区(Ⅰ区)为蓝色区域,主要对应黄淮平原和南襄盆地部分地区,其地壳相对较薄,基本为30~33km,横向变化不大,显示地壳最薄处为平顶山台所在的下方地壳,厚度为30.0km,这与王椿镛等(2017)的计算结果“河淮地区为华北平原地壳厚度最薄的地区,仅为30km”一致;第2梯度区(Ⅱ区)为米黄色区域,主要对应华北平原向东秦岭-大别构造带过渡区域,地壳逐渐增厚,厚度为33~36km,镶嵌在东秦岭山脉中间的南襄盆地厚度基本属于这一范围;第3梯度区(Ⅲ区)为红黄区域,主要对应鄂西北、陕南等东秦岭山脉地区,地壳相对最厚,厚度为36~41km。其中,兴山台所在区域的地壳为全区最厚,达到41.0km,接近第二级阶梯地壳厚度42~44km的范围。

图 3 南襄盆地及邻区地壳厚度分布(a)及泊松比分布(b)

① 淅川-西峡震群:2008~2018年发生在淅川、西峡等地ML≥2.0地震,共16次

3.2 泊松比分布

一般而言,泊松比作为岩石的弹性力学参数,是探索地球内部物质成分和物理状态的直接依据。Ji(2008)基于各种岩石密度和P波、S波速度随围压变化的波速数据,计算了岩石圈中主要种类岩性的泊松比,他认为:①在长英质岩石,如花岗岩、花岗闪长岩以及闪长岩中,石英含量的增加会使岩石的泊松比减小,而长石含量的增加会使岩石的泊松比增加;②铁镁质岩如辉长岩、辉绿岩以及基性片麻岩具有较高的泊松比,主要原因就是这些岩石不含石英,而含大量基性斜长石;③新鲜的橄榄岩较之辉长岩、辉绿岩和基性片麻岩具有较低的泊松比,但橄榄岩的泊松比实际上随其蛇纹石化程度的增加而明显地增加,原因是蛇纹石具有特别高的泊松比。

根据波速比计算结果进行2次计算,得到了研究区的泊松比,通过插值得到泊松比的空间分布(图 3(b))。如图 3(b)所示,泊松比高值集中在房县-兴山附近,其中,房县台所在地区的泊松比为0.34,为研究区最高,同时在房县周边的兴山、十堰等地也显示出高泊松比,指示该地区地壳物质广泛存在铁镁质岩体,这与该区大量出露的寒武-奥陶系的辉长岩、苏长岩等铁镁质基性岩恰好吻合。低值区主要位于商南-西峡-淅川一带,最低值区位于商南地区,其值为0.18,与其邻近的内乡、丹江、卢氏等地泊松比也为0.18~0.24,表明该区的商南-西峡-淅川一带地壳中包含了较多的石英成分,即说明以花岗岩、花岗闪长岩等为代表的长英质岩体在该区大量存在。这一推断得到了野外地质证据的验证(王晓霞等,2015):在商丹断裂(F3)以南的刘岭群中出现绢云母石英片岩、石英岩等,以及更靠南侧的陡岭群中出露新元古代混合片麻岩、混合花岗质片麻岩等,基本处于泊松比低值区内;处于朱夏断裂(F2)和商丹断裂(F3)之间的秦岭杂岩群在商南-西峡出露的花岗岩以古生代Ⅰ型为主,岩石类型为片麻状花岗闪长岩、黑云母花岗闪长岩和石英闪长岩等,也在计算得到的泊松比低值区内;位于栾川-南召断裂(F1)和朱夏断裂(F2)之间的二郎坪群规模较大,包含了古生代黑云母花岗闪长岩,部分位于泊松比低值区内。

值得注意的是,上述低值区揭示着长英质岩体自商南经西峡一直向东至少延伸到淅川境内,这一点可从南阳凹陷的钻探记录中得到验证:在盆地西侧的镇平浅1井中钻取到陡岭群岩芯,在镇平浅2井中钻到刘岭群的岩石芯样,它们均富含长英质岩石(刘绍龙等,1998)。另外,基于各构造岩相带均受到商丹、朱夏等NW向断裂围限的事实,以及对盆地两侧构造岩相铅同位素的研究表明,盆地东南侧的桐柏地区是秦岭造山带的构造东延部分,桐柏断裂为商丹断裂在盆地东侧的延续(Zhang et al,1997张宏飞等,1999)。由此可推断,西侧的二郎坪群、秦岭群、陡岭群和刘岭群等构造岩相可能沿断裂向东延伸到盆地内部,在淅川、西峡一带与盆地的沉积地层形成某种接触关系,因介质的较大差异,进而导致地震射线在传播过程中出现反射。

4 结论

(1) 南襄盆地及邻区的地壳厚度与地形变化关联紧密,呈现出由NE向SW方向增厚的趋势,地壳厚度为30~41km。

(2) 泊松比高值区集中在房县-兴山附近,其值为0.31~0.34,可能与区域内大量出露寒武-奥陶系的辉长岩、苏长岩等铁镁质基性岩相关;低值为0.18~0.24,分布在商南-西峡-淅川一带,与区域内存在大量的长英质岩有关,它们与南襄盆地内的新生代沉积地层形成整合接触关系。由此可推断,在淅川、西峡一带存在因介质的较大差异而导致地震射线在传播过程中出现反射的可能。

考虑到地震反射的物理机制,我们认为应该存在一个规模较大的岩层带,并且岩层带与盆地沉积形成某种接触关系。然而,仅通过本文的泊松比数据无法进一步验证,需要搜集此区域的地球物理勘探、地质等证据,以进一步证实上述推断。

致谢: 衷心感谢审稿人对本文的热心指导;感谢中国科学院南海海洋研究所黄海波博士在波形数据处理方面的耐心指导;感谢广东省地震局沈玉松博士对本文GMT绘图提供帮助。
参考文献
邓嘉美、陈佳、高琼等, 2016, 云南宾川地区的地壳厚度和泊松比分布特征研究, 地震, 36(3): 99-108. DOI:10.3969/j.issn.1000-3274.2016.03.010
冯锐, 1985, 中国地壳厚度及上地幔密度分布(三维重力反演结果), 地震学报, 7(2): 143-157.
高延光、李永华, 2014, 中国东北-华北地区地壳厚度与泊松比及其地质意义, 地球物理学报, 57(3): 847-857.
龚辰、李秋生、叶卓等, 2016, 远震P波接收函数揭示的张家口(怀来)-中蒙边境(巴音温多尔)剖面地壳厚度与泊松比, 地球物理学报, 59(3): 897-911.
宫猛、李信富、张素欣等, 2015, 利用接收函数研究河北及邻区地壳厚度与泊松比分布特征, 地震, 35(2): 34-42. DOI:10.3969/j.issn.1000-3274.2015.02.004
黄海波、郭兴伟、夏少红等, 2014, 华南沿海地区地壳厚度与泊松比研究, 地球物理学报, 57(12): 3896-3906. DOI:10.6038/cjg20141204
嵇少丞、王茜、杨文采, 2009, 华北克拉通泊松比与地壳厚度的关系及其大地构造意义, 地质学报, 83(3): 324-330.
刘绍龙、刘宁、薛玉玮等, 1998, 南襄盆地隐伏区前中生界展布规律及其大地构造意义, 河南石油, 12(4): 1-6.
齐刚、陈棋福, 2015, 太行山与燕山交汇部位的地壳厚度与泊松比分布特征, 地球物理学报, 58(9): 3239-3250.
钱辉、李海兵、Mechie J等, 2015, 接收函数H-k叠加方法在龙门山台阵的应用, 地球物理学进展, 30(1): 91-98.
沈玉松、康英、徐果明, 2013, 广东及其邻域的地壳厚度和泊松比分布, 中国地震, 29(2): 210-218. DOI:10.3969/j.issn.1001-4683.2013.02.003
滕吉文、曾融生、闫雅芬等, 2002, 东亚大陆及周边海域Moho界面深度分布和基本构造格局, 中国科学:D辑, 32(2): 89-100.
王椿镛、吴庆举、段永红等, 2017, 华北地壳上地幔结构及其大地震深部构造成因, 中国科学:地球科学, 47(6): 684-719.
王林昌、王勉、王晓谦等, 2011, 南阳盆地特殊震相研究, 高原地震, 23(3): 30-33. DOI:10.3969/j.issn.1005-586X.2011.03.006
王晓霞、王涛、张成立, 2015, 秦岭造山带花岗质岩浆作用与造山带演化, 中国科学:地球科学, 45(8): 1109-1125.
吴庆举、曾融生, 1998, 用宽频带远震接收函数研究青藏高原的地壳结构, 地球物理学报, 41(5): 669-679. DOI:10.3321/j.issn:0001-5733.1998.05.010
许卫卫、郑天愉, 2005, 渤海湾盆地北西盆山边界地区泊松比分布, 地球物理学报, 48(5): 1077-1084. DOI:10.3321/j.issn:0001-5733.2005.05.014
杨冬、王典、刘殿秘等, 2008, 南襄盆地及邻域壳幔结构地球物理特征与油气条件分析, 世界地质, 27(1): 48-52. DOI:10.3969/j.issn.1004-5589.2008.01.009
张宏飞、张利、高山等, 1999, 桐柏地区变质杂岩和侵入岩类Pb同位素组成特征及其地质意义, 地球科学:中国地质大学学报, 24(3): 269-273.
郑勇、葛粲、谢祖军等, 2013, 芦山与汶川地震震区地壳上地幔结构及深部孕震环境, 中国科学:地球科学, 43(6): 1027-1037.
Christensen N I, 1996, Poisson's ratio and crustal seismology, J Geophys Res, 101(B2): 3139-3156. DOI:10.1029/95JB03446
He R Z, Shang X F, Yu C Q, et al, 2014, A unified map of Moho depth and VP/VS ratio of continental China by receiver function analysis, Geophys J Int, 199(3): 1910-1918. DOI:10.1093/gji/ggu365
Ji S C, 2008, Deformation mechanisms, rheology, and seismic properities of rocks, 539, Beijing: Geological Publishing House.
Ji S C, Wang Q, Xia B, 2002, Handbook of seismic properties of minerals, rocks and ores, 630, Montreal: Polytechnic International Press.
Jo E, Hong T K, 2013, VP/VS ratios in the upper crust of the southern Korean Peninsula and their correlations with seismic and geophysical properties, J Asian Earth Sci, 66: 204-214. DOI:10.1016/j.jseaes.2013.01.008
Langston C A, 1979, Structure under Mount Rainier, Washington, inferred from teleseismic body waves, J Geophys Res, 84(B9): 4749-4762.
Li Y H, Gao M T, Wu Q J, 2014, Crustal thickness map of the Chinese mainland from teleseismic receiver functions, Tectonophysics, 611: 51-60. DOI:10.1016/j.tecto.2013.11.019
Owens T J, Zandt G, Taylor S R, 1984, Seismic evidence for an ancient rift beneath the Cumberland Plateau, Tennessee:a detailed analysis of broadband teleseismic P waveforms, J Geophys Res, 89(B9): 7783-7795.
Salah M K, Seno T, 2008, Imaging of VP, VS, and Poisson's ratio anomalies beneath Kyushu, southwest Japan:Implications for volcanism and forearc mantle wedge serpentinization, J Asian Earth Sci, 31(4/5/6): 404-428.
Salah M K, Şahin S, Aydin U, 2011, Seismic velocity and Poisson's ratio tomography of the crust beneath East Anatolia, J Asian Earth Sci, 40(3): 746-761. DOI:10.1016/j.jseaes.2010.10.021
Tugume F, Nyblade A, Julià J, 2012, Moho depths and Poisson's ratios of Precambrian crust in East Africa:Evidence for similarities in Archean and Proterozoic crustal structure, Earth Planet Sci Lett, 355/356: 73-81. DOI:10.1016/j.epsl.2012.08.041
Wei Z G, Chen L, Li Z W, et al, 2016, Regional variation in Moho depth and Poisson's ratio beneath eastern China and its tectonic implications, J Asian Earth Sci, 115: 308-320. DOI:10.1016/j.jseaes.2015.10.010
Yang W C, Qu C, Yu C Q, 2011, Crustal Poisson's ratio anomalies in the eastern part of North China and their origins, Geosci Front, 2(3): 313-321.
Zandt G, Ammon C J, 1995, Continental crust composition constrained by measurements of crustal Poisson's ratio, Nature, 374(6518): 152-154. DOI:10.1038/374152a0
Zhang H F, Gao S, Zhang B R, et al, 1997, Pb isotopes of granitoids suggest Devonian accretion of Yangtze(South China)craton to North China craton, Geology, 25(11): 1015-1018. DOI:10.1130/0091-7613(1997)025<1015:PIOGSD>2.3.CO;2
Zhu L P, Kanamori H, 2000, Moho depth variation in southern California from teleseismic receiver functions, J Geophys Res, 105(B2): 2969-2980. DOI:10.1029/1999JB900322