利用接收函数方法确定青藏高原东北缘近地表S波速度

公式推导和数值实验都表明远震接收函数直达波幅度对近地表S波速度比较敏感,可以用于确定近地表（如沉积层）的S波速度值.本文利用该方法分析了青藏高原东北缘区域数字地震台网36个台站的远震资料,得到了研究区近地表S波速度分布特征.研究结果表明：在深度5±2.5 km范围内,祁连山—河西走廊区域近地表S波速度均值较小;位于地形梯度带的甘东南区域近地表S波速度呈现较为复杂的分布.通过与该区域沉积层厚度结果对比,发现近地表S波速度与沉积层厚度呈现较好的负相关.本研究是对目前接收函数方法的有效补充,在探测研究近地表速度结构方面,相比人工地震方法其成本几乎可忽略且环保.在目前地震波观测数据海量增长的背景下,该方法对于确定不同区域近地表结构具有较大潜力.     

用接收函数方法反演青藏高原东北缘地壳结构

利用ASCENT计划于2007年布设在青藏高原东北缘的18个宽频带流动台站约一年的观测资料获得了2547个接收函数.使用H-?域的搜索算法,得到了14个台站下方的地壳厚度.对于数据质量较差的3个台站,通过Moho的Ps转换波的到时估算出地壳厚度值.计算结果表明,研究区的地壳结构复杂,Moho深度变化范围为40~60 km.海原断裂附近Moho模糊,而且较两侧明显地加深.秦祁地块由西向东Moho逐渐变浅,105°E以东,Moho平均深度约为45 km,以西则在50 km以上.结合面波研究结果推测,在105°E附近可能存在一个秦岭与祁连分界线.以Crust 2.0作为初始模型,把计算得到的地壳厚度值作为约束,用线性反演方法得到了15台站下方的S波速度结构.通过与研究区人工源地震测深结果比较,发现二者具有较好的一致性,表明反演结果可靠.综合分析认为:在青藏高原区,青藏高原与西秦岭过渡带和西秦岭北段均存在中下地壳低速层,推测它们可能与附近断裂和深部物质运移有关.从青藏高原内部到边缘,中下地壳速度逐渐升高.再结合研究区地壳的速度(尤其是低速层)的分布特点以及与地壳物质成分相关的泊松比的变化规律等综合分析,认为地壳流在青藏高原东北部的边缘地带可能不存在,地壳可能是通过在挤压方向上的缩短而加厚.     

青藏高原东北缘地壳S波速度结构与泊松比及其意义

利用甘肃地震台网16个台站记录的远震资料,采用最大熵谱反褶积方法,得到了各个台站的接收函数. 采用接收函数扫描法和线性反演方法对研究区的壳幔结构进行了研究,这两种接收函数方法得出的结果具有很好的一致性. 青藏高原东北缘地壳厚度变化剧烈,祁连块体为50～55 km、柴达木块体和河西走廊为45 km左右（合作台除外）,由北向南,Moho界面呈中央下凹的准对称状. 研究区地壳VP/VS介于1.66～1.85（σ=0.215～0.294,均值0.254）,其均值接近或略低于全球平均值;S波速度结构可见壳幔过渡带具有明显的突跳,结合其他地球物理学证据,推断该区可能不存在岩浆底侵作用和地壳部分熔融现象. 该区地壳VP/VS值与地壳厚度呈反相关关系,推断该区地壳的主要组成成分以中酸性岩石为主,其45～55 km厚的地壳可能主要是通过上地壳的叠置形成的.     

用宽频带远震接收函数研究青藏高原的地壳结构

为研究青藏高原的岩石层构造及其动力学过程,根据记录到的来自台站东北方向的大量宽频带远震林波波形资料,应用时间域的最大熵谱反褶积算法,得到了11个（全部）PASSCAL（地壳与岩石层的地震台阵研究计划）台站的接收函数．利用时间域广义线性反演的Jumping（跳动）算法,引入模型光滑度约束,并将合成地震图的Kennett完全算法及微分地震图的Randall快速算法用于接收函数的正演计算,由台站接收函数获得了各台站下方的一维S波速度分布．反演结果表明,青藏高原的Moho界面在班公—怒江缝合带附近存在明显的深度错断;在日喀则、拉萨、桑雄和安多等地的地壳内部,可以连续观测到三个显著的速度界面h1,h2和h3,其中h1和h2可以连续追踪到温泉、二道沟和不冻泉等地,而h3在班公—怒江缝合带以北消失;在日喀则、拉萨、桑雄、安多、二道沟和不冻泉等地有壳内低速层．关于班公—怒江缝合带附近Moho界面的错断现象,一个可能的解释是班公—怒江缝合带是印度地壳向欧亚下地壳挤入的前沿．     

用宽频带远震接收函数研究青藏高原的地壳结构

为研究青藏高原的岩石层构造及其动力学过程,根据记录到的来自台站东北方向的大量宽频带远震林波波形资料,应用时间域的最大熵谱反褶积算法,得到了11个（全部）PASSCAL（地壳与岩石层的地震台阵研究计划）台站的接收函数．利用时间域广义线性反演的Jumping（跳动）算法,引入模型光滑度约束,并将合成地震图的Kennett完全算法及微分地震图的Randall快速算法用于接收函数的正演计算,由台站接收函数获得了各台站下方的一维S波速度分布．反演结果表明,青藏高原的Moho界面在班公-怒江缝合带附近存在明显的深度错断;在日喀则、拉萨、桑雄和安多等地的地壳内部,可以连续观测到三个显著的速度界面h1,h2和h3,其中h1和h2可以连续追踪到温泉、二道沟和不冻泉等地,而h3在班公-怒江缝合带以北消失;在日喀则、拉萨、桑雄、安多、二道沟和不冻泉等地有壳内低速层．关于班公-怒江缝合带附近Moho界面的错断现象,一个可能的解释是班公-怒江缝合带是印度地壳向欧亚下地壳挤入的前沿．     

帕米尔东北缘地壳结构的P波接收函数研究

利用位于新疆帕米尔东北缘地带12个固定数字地震台和天山动力学Ⅱ期10个流动宽频带数字地震台记录的高质量远震波形数据,应用接收函数H-κ叠加方法研究了帕米尔东北缘的地壳厚度-泊松比特征和部分台站下方的壳内界面深度.研究发现:(1) 帕米尔东北缘的Moho面起伏变化剧烈,其总体分布呈现东薄西厚、南厚北薄的特征,由塔里木盆地向天山延伸,地壳厚度约从45 km加深到55 km,从塔里木盆地向西昆仑山延伸,地壳厚度约从45 km加深到69 km;(2)沿着天山动力学Ⅱ期剖面,位于塔里木盆地北缘台站的壳内间断面的深度约为13~16 km,向北进入天山南麓加深到20 km左右,继续向北进入南天山山区壳内间断面不明显,可能暗示塔里木盆地基底向北俯冲,俯冲距离可能到达南天山的山前;(3)研究区地壳泊松比变化复杂(约从0.20到0.31),显示地壳物质组成的复杂性和显著的不均匀构造;(4)整个研究区的地壳厚度和泊松比之间没有明显的相关性,但天山动力学Ⅱ期剖面的结果表明,从塔里木盆地北缘到西南天山,地壳厚度和泊松比之间存在反相关关系,意味着天山地壳的增厚可能主要是通过以长英质岩石为主要组成成分的上地壳叠置而成;(5) 研究区全部地震台地壳厚度与海拔高程的线性回归方程表明地壳厚度与海拔的相关性相对较弱(相关系数为0.66),天山动力学Ⅱ期10个台站的地壳厚度与海拔具有很好的相关性(相关系数为0.85),可能表明沿该剖面地壳整体上处于相对均衡的状态.     

利用接收函数反演青藏高原西部地壳S波速度结构

相对于宽阔的腹地,青藏高原西部南北向宽度仅约600 km,却记录了印度和欧亚板块汇聚的深部过程及其响应.本文用22台宽频带流动地震台站在西缘构建了一条南北向探测剖面(~80°E,TW-80试验).利用接收函数反演剖面下方S波速度结构,综合西部已有的宽频带探测结果,分析认为:印度板块向北俯冲可能已到达班公湖-怒江缝合带附近,俯冲过程中下地壳发生榴辉岩化;喀拉昆仑断裂带、班公湖-怒江缝合带、阿尔金断裂带均为切穿地壳的深断裂,莫霍面发生错断;喀拉昆仑断裂带和龙木错断裂带之间的中上地壳没有发现连续的S波低速体,说明可能缺乏解耦层,支持青藏高原西部地壳为整体缩短增厚模式.     

青藏高原东北缘地壳三维速度结构

本文用1980—2000年M≥1.5的2 032个天然地震事件的38 052个〖AKP-〗、〖AKS-〗、Pm、Sm、Pn和Sn震相到时及人工地震测深给出的Moho面形态资料,利用地震层析技术反演了32°~40°N, 100°~108°E区域内地壳地震波速度结构.从层析成像图象中可以得到,本区的地壳可分成4个层位.第1层（埋深约在0~3 km）为沉积层, 速度梯度约为0.2 s-1;第2层（埋深约在3~17 km）为上地壳, 其顶部速度梯度约为0.1 s-1, 下部速度横向变化较大且存在低速块体;第3层（埋深约在17~36 km）为中地壳, 速度梯度约为0.03 s-1;第4层（埋深约在36 km—Moho）为下地壳, 是一个契形层,总的趋势是西厚东薄,青藏高原较厚逐渐向鄂尔多斯地块和扬子准地台方向变薄,各处的地震波速度梯度不尽相同.     

青藏高原东北缘—鄂尔多斯地块地壳上地幔S波速度结构

1999～2000年从青海玛沁到陕西榆林,横跨青藏高原东北缘和鄂尔多斯布设了一条由47台宽频带数字地震仪组成的长约1000km的流动地震台阵观测剖面.利用记录到的远震体波波形资料和接收函数方法获得了剖面下0～100km深度的地壳和上地幔S波速度结构.结果表明,沿观测剖面地壳结构显示了明显的分块特征; 地壳厚度自东向西由40km增加到64km左右;在海原地震带下方和西秦岭断裂以西到日月山断裂之间的区域Moho间断面结构复杂;在1920年海原震区及其西侧,上地壳存在明显的低速层,在该地区的绝大部分地震分布在该低速层东边界偏向高速区一侧;祁连山东缘Moho面有约4km的深度间断,壳内向西逐渐减薄的低速层内有大量微震发生,沿祁连山的逆冲加走滑的构造运动在深度上已经穿透了Moho面;在玛沁断裂和日月山断裂之间,上地壳存在厚度很大的低速层,同时该区域下地壳也明显加厚.研究结果表明,青藏高原东北缘与鄂尔多斯地块之间的过渡带地壳变形强烈,地壳结构较为破碎,这与该地区地震频发相一致.     

用远震接收函数反演上海及其邻区地壳速度结构

通过提取江苏、浙江及上海数字地震台网台站远震记录的接收函数,反演了这些台站下方地壳的S波速度结构,同时计算出了地壳的厚度及泊松比。初步分析发现,该地区莫霍面起伏不大,江苏地区地壳厚度大约是30km,地壳中存在4～5km厚的低速夹层,位于地面以下15～22km之间。浙江地区地壳内部不存在低速区,并且莫霍面深度自北向南逐渐加深。     

利用双差成像方法反演青藏高原东北缘及其邻区地壳速度结构

利用区域尺度双差层析成像方法,使用2009年1月至2017年2月的近震资料,对青藏高原东北缘及其邻区内记录到的地震事件进行震源位置和三维速度结构的联合反演.重新定位后震源空间位置得到明显改善,浅层的地震波速与地形和沉积层厚度对应较好,研究区地震主要发生在河西走廊过渡带的低泊松比区域.本文将研究区分为五个区域并分别对其层析成像结果进行了讨论,结果显示研究区不同地块之间地壳结构变化明显,地壳物质整体呈酸性,青藏高原东北缘地壳增厚可能主要发生在中下地壳.

     

青藏高原东北缘远震P波走时层析成像研究

利用青藏高原东北缘地区固定地震台网2010年4月至2015年3月期间记录的远震事件,采用多道波形互相关方法（Multi-Channel Cross-Correlation）拾取了10697个有效P波相对走时残差数据,进而采用FMTT（Fast Marching Teleseismic Tomography）方法获得了青藏高原东北缘上地幔400 km深度范围内的P波速度结构,结果显示：秦祁地块下面存在深达70 km的高速异常,阻断了青藏高原块体中下地壳低速层向东北方向的延伸;40~140 km深度范围内,四川西南部存在一个低速区,该低速区穿过龙门山断裂带进入到四川盆地内部;祁连山造山带东部低速异常区从地壳一直延伸到上地幔400 km处,表明这里可能存在一个上地幔到地壳间的热流通道;松潘—甘孜地块分布大面积的低速异常区,而鄂尔多斯地块西南缘相对速度较高,这与青藏高原为软块体、介质密度低和鄂尔多斯块体为硬块体、介质密度高相吻合.

     

利用接收函数研究兴蒙造山带阿巴嘎地区的地壳结构

本文利用内蒙阿巴嘎地区布设的38个宽频带地震台站记录到的远震数据,采用P波接收函数共转换点叠加方法（CCP）揭示台站下方Moho面起伏形态,并利用H-κ方法进一步得到地壳厚度和壳内平均波速比值.结果显示,研究区地壳厚度为35~44 km,均值约为40 km,西南部的鄂尔多斯盆地边缘地壳较厚,东北部的阿巴嘎火山群地区地壳显著变薄.研究区地壳平均波速比值在1.70~1.87之间,均值为1.76,其中阿巴嘎火山地区波速比值明显偏高.CCP叠加结果显示研究区Moho界面较平缓,但在缝合带附近存在明显的变化.我们推测,新生代阿巴嘎火山地区薄的地壳和高波速比值可能是由于火山活动底侵作用引起上地幔铁镁质物质侵入下地壳所致.     

利用远震接收函数探测四川盆地及周边地区的地壳结构

青藏高原东南缘的龙门山断裂两侧具有陡峭的地形特征,在约50~100 km的水平距离内,地形高程从2000 m增加到4000 m,该区强烈的壳幔变形特征及地球动力学模式一直是研究的热点问题.本文从四川地区49个固定台站记录的远震资料提取了P波接收函数,获得了四川盆地及周边的地壳厚度和泊松比,并以此构建反演的初始模型.在线性反演的基础上,引入了分别拟合低频和高频接收函数的两步反演技术,用以反演台站下方的地壳S波速度结构.数字试验表明,该方法可以有效抑制接收函数反演的不唯一性,为了得到最优解,最后用Bootstrap重采样技术估计解的不确定性.结果表明,四川盆地的地壳厚度在40~46 km,松潘-甘孜块体北部的地壳厚度为46~52 km,而南部增厚到50~60 km.从四川盆地向西跨过龙门山断裂,地壳厚度增加了10~15 km.在四川盆地及周边地区,地壳泊松比在0.26~0.32之间,呈块体分布特征,高泊松比（0.28~0.32）主要沿龙门山断裂以及安宁河-小江断裂分布.地壳S波速度结构表明,来自青藏高原中部的中下地壳低速层可能受到了坚硬的四川盆地阻挡,改变原来的运动方向并沿龙门山断裂展布,由于低速层的囤积导致该区地形陡峭和下地壳增厚.

     

利用P波接收函数研究青藏高原东南缘地幔转换带结构

研究青藏高原东南缘的深部结构对于理解印度-欧亚板块的碰撞机理和青藏高原的形成演化具有重要的科学意义.本研究对布设在研究区域内566个固定和流动地震台站的波形资料进行了处理,获得77853条高质量P波接收函数,应用接收函数共转换点（CCP）叠加技术获得了研究区域下方精细的地幔转换带间断面起伏形态及转换带厚度变化图像.结果表明：研究区域南北方向上具有两个明显的转换带增厚异常区,南侧异常区位于滇中次级块体与印支块体下方,可能是新特提斯洋板片与上部印度板块间断离并部分滞留在转换带底部的结果;北侧川西地区异常增厚可能与上方岩石圈拆沉并降至转换带有关;腾冲火山起源可能是板块俯冲过程中发生断离造成软流圈物质部分熔融,湿热物质上涌所致.

     

青藏高原东南缘地壳速度结构及云南漾濞MS6.4地震深部孕震环境

青藏高原东南缘岩石圈变形强烈、地震活动频繁,对其深部结构进行研究有助于提高对其演化及强震发震机理等问题的认识.本研究利用青藏东南缘固定和流动观测台站记录的地震P波走时数据,采用体波走时层析成像方法得到研究区地壳、上地幔顶部三维P波速度结构.结果显示,峨眉山大火成岩省内带范围呈明显的高速异常,推测为二叠纪时期地幔柱活动残留在地壳内的基性和超基性幔源物质.川西北次级块体和大火成岩省内带东西两侧存在低速带,可能是壳内部分熔融存在和中下地壳流动的证据:高原物质向南运移时受大火成岩省壳内高速体阻挡后分为两支,东支沿安宁河断裂—则木河断裂—小江断裂延伸,西支沿红河断裂向上地壳运移并逐渐穿过红河断裂.云南漾濞M_S6.4地震序列全部位于西支低速通道之上,推测构造块体SE向运动和地壳流作用使得应力在上地壳震源区进一步集中,共同驱动断裂活动导致漾濞地震发生.     

The crust and upper mantle discontinuity structure beneath Alaska inferred from receiver functions

In this study, three receiver function stacking methods are used to study the detailed crust and upper mantle structure beneath south-central Alaska. We used teleseismic waveform data recorded by 36 stations in the Broadband Experiment Across the Alaska Range (BEAAR) and 4 permanent stations in Alaska. H − κ stacking method using P-to-S converted wave and its multiply reflected waves between the Earth's surface and the Moho discontinuity is adopted to estimate the crustal thickness (H) and average crustal V_P/V_S ratio (κ) in this region. The receiver function results for 24 stations show that the crustal thickness under Alaska ranges from 26.0 to 42.6 km with an average value of 33.8 km, and the V_P/V_S ratio varies from 1.66 to 1.94 with an average value of 1.81 which corresponds to an average Poisson's ratio of 0.277 with a range from 0.216 to 0.320. High Poisson's ratios under some stations are possibly caused by partial melting in the crust and the uppermost mantle. Common converted point (CCP) stacking results of receiver functions along three lines show clear Moho and slab images under this subduction zone. The depths of the slab from our CCP stacking images are consistent with those estimated from the Wadati–Benioff Zone (WBZ). In the area between two stations DH2 (147.8°W, 63.3°N) and DH3 (147.1°W, 63.0°N), a Moho depth offset of about 10 km is found by both the H − κ and CCP stacking techniques. Common depth point (CDP) stacking of receiver functions shows not only the 410-, 520- and 660-km discontinuities, but also significant variations (−30 to 15 km) in the transition zone thickness under the southwest and southeast parts of the study region. The transition zone becomes thinner by 20–30 km, indicating that the temperature there is 150–200 K higher than that of the normal mantle.     

利用远震接收函数研究黄海周边地区的地壳上地幔结构

利用宽频带流动台站(YSBSN)记录的远震波形数据和远震接收函数方法,反演了黄海东、西两侧地壳上地幔的S波速度结构.结果表明,莫霍面深度在30~38 km之间变化,位于中方一侧的JNN台下方地壳厚度最大,可以归因于华北板块和扬子板块的碰撞;韩方一侧的地壳厚度自北向南逐渐变厚,但仍然难以厘定朝鲜半岛南部潜在碰撞带的位置,这些问题的解决需要更大范围的流动台站观测.由于部分台站位于巨厚的沉积层和多孔的火山岩之上,与浅部构造的相关性使得接收函数表现出较大振幅的混响,从而影响了来自深部结构的转换震相.