Crustal attenuation in the southern Tibetan Plateau by reverse two-station Lg-wave Q value tomography
-
摘要: 首次基于2017—2019年西藏自治区区域台网27个宽频带固定台站记录的757次地震的波形资料，利用反双台法开展了青藏高原南部地区1 Hz的Lg波Q值层析成像研究。研究中采用3.5—2.4 km/s的速度窗截取了1 981条Lg波，计算得到13 543条路径上的Q值，测试了1°×1°和0.5°×0.5°网格下的棋盘格恢复情况，得到了0.5°×0.5°分辨率的Lg波Q0值层析成像。反演结果显示：青藏高原南部地壳整体的Lg波呈高衰减、低Q值，与P波速度负异常、地热分布及东部的两条裂谷系对应良好，因此推断青藏高原南部地壳存在广泛的熔融物质；两条可能存在的流体-熔融物质通道中，主通道位于亚东—谷露裂谷与桑日—错那裂谷之间，副通道沿雅鲁藏布江缝合带分流而出。此外，还对亚东—谷露裂谷两侧熔融物质的分布差异予以分析，结果表明，印度板块与欧亚板块碰撞前端存在不同的动力学演化模式，亚东—谷露裂谷以西符合缩短增厚理论，以东符合“水泵”模式。Abstract: Based on the seismic waveform data of 757 earthquakes recorded by 27 broadband stations of Tibet Autonomous Region seismic network from 2017 to 2019, it is for the first time that the reverse two-station method has been applied to 1 Hz Lg-wave Q value tomography study in the southern Tibetan Plateau. In this research, total 1 981 Lg-waves were intercepted from the velocity window of 3.5−2.4 km/s and the Q values of 13 543 paths were calculated. After testing the checkerboard recovery of 1°×1° and 0.5°×0.5° grids respectively, we got the tomography of the Lg-wave Q0 value of the southern Tibetan Plateau with 0.5°×0.5° resolution. The inversion result shows that there exhibits high attenuation and low Q values of Lg-wave in the southern Tibetan Plateau crust, which is highly consistent with the negative anomaly of P-wave velocity, the geothermal distribution and two rift valleys in the east of the Tibet. Therefore it is inferred that there may be widespread molten material in the crust of the southern Tibetan Plateau and two fluid-melting channels. The main channel is located between the Yadong-Gulu rift and the Sangri-Cona rift, and the secondary channel flows out along the Yarlung Zangbo suture zone. By analyzing the differences of molten material distribution on both sides of the Yadong-Gulu rift, it is considered that there are different dynamical evolution models in the front-end of collision between Indian Plate and Eurasian Plate, the dynamical evolution to the west of Yadong-Gulu rift is in accordance with the theory of shortening and thickening, and that to the east of Yadong-Gulu rift is consistent with the “pump” mode.
-
Keywords:
- Tibetan Plateau /
- Lg wave attenuation （QLg） /
- reverse two-station method /
- tomography
-
引言
华北克拉通形成于太古代至古元古代，是世界上最古老的克拉通之一（图1）。华北克拉通在中生代以前较为稳定（朱日祥等，2011）。根据地幔捕虏体的证据，中生代以来华北克拉通东部岩石圈重新活化，发生了明显的减薄和破坏，失去了克拉通性质（Zhu et al，2012b）；其西部的鄂尔多斯地块内部长期以来缺乏岩浆活动及中强地震（马瑾，2009），基本保持了克拉通稳定性，表明华北克拉通东部和西部经历了差异性的动力学过程（Bao et al，2013）。新生代以来，受到印度—欧亚板块碰撞的远场效应（Tapponnier，Molnar，1977）及西太平洋俯冲（Yin，2010）的影响，鄂尔多斯地块周缘形成了一系列断陷盆地，现今仍处于拉张状态（盛书中等，2015），且强震活动频繁（韩竹军等，2003）。
![]()
图 1 研究区构造背景及地震活动性Figure 1. Tectonic background and seismicity in the study area鄂尔多斯地块东缘的山西断陷带长约1 000 km （Clinkscales et al，2020），是世界上最大的新生代大陆裂谷之一（Xu et al，2020）。大陆裂谷的出现可能预示着岩石圈尺度的改造和破坏（Thybo，Nielsen，2009）。但与世界上其它的大陆裂谷相比，山西断陷带的火山活动水平相对较低，主要集中在北部地区（Xu，Ma，1992；Xu et al，2020），其中第四纪仍在活动的大同火山的喷发产物主要为碱性玄武岩和拉斑玄武岩（Xu et al，2005）。地球化学研究认为大同火山是了解华北克拉通岩石圈结构剧烈变化过程的关键（Xu et al，2005）。然而目前对大同火山地区的岩浆上涌模式及动力学过程尚存在较大争议，关于其形成机制的主要模型有大地幔楔模型（Huang，Zhao，2006）、地幔柱模型（Lei，2012）以及青藏高原东北缘软流圈流动模型（Li et al，2018）等。
鄂尔多斯地块及周边地区已开展过一系列的层析成像研究（Tian et al，2009；Lei，2012；Jiang et al，2013；Lü，2019）。尽管大多数成像结果都揭示出鄂尔多斯下方存在高速异常，地块东北部呈现出明显的低速特征，但由于这些研究大多采用分布稀疏的固定地震台站和少量的流动地震台站的观测资料，其结果显示的异常体范围和延伸深度存在较大的差异，影响了我们对该地区动力学过程的理解。例如，一些研究认为鄂尔多斯下方的高速异常下延深度为200—400 km，保留了较厚的岩石圈根（Huang，Zhao，2006；Tian et al，2009；Zhao et al，2009；Lei，2012；Jiang et al，2013；毛慧慧等，2016；Li et al，2018），而少数面波成像研究则认为鄂尔多斯地块下方高速异常下延深度仅为120—150 km （An et al，2009；Tang et al，2013；Yao et al，2020），岩石圈可能已经被局部改造。这些深部结构成像的差异，均会对该地区岩石圈活化位置及成因机制的认识产生影响。
中国地震局地球物理研究所自2016年以来在鄂尔多斯及周边地区开展了密集流动地震台阵观测，平均台间距约为35 km。本文拟联合使用密集流动台站与固定地震台站的观测资料，采用远震P波层析成像方法获得鄂尔多斯及周边地区的深部结构，并对其构造意义进行分析讨论，以探索该地区深部构造信息及动力学过程。
1. 数据和方法
本文收集了鄂尔多斯地块及周边地区943个台站记录到的404个远震事件（30°—90°，M≥5.5）（图2），共计15万3 963条观测数据。使用的台站包括中国地震台阵探测第三期的472个流动台站（2016年1月至2019年1月）和中国地震台网的471个固定台站（2016年1月至2019年1月）（郑秀芬等，2009）。中国地震台阵探测第三期流动台站覆盖了鄂尔多斯地块及中部地块，台间距约为35 km，极大改善了此前鄂尔多斯地块中东部及华北克拉通中部稀疏的台站分布。
![]()
图 2 台站（a）及地震震中（b）分布图图（a）中实线为构造边界，虚线为南北重力梯度带，倒三角形表示固定台站，三角形表示中国地震科学台阵三期流动台阵；图（b）中圆点为地震，星形为研究中心Figure 2. Distribution of stations （a） and epicenters （b）In Fig. （a），solid lines delineate the tectonic units，the dashed line is the north-south gravity gradient zone，the inverted triangles indicate permanent stations，and the triangles indicate the temporary stations of phase III of the ChinArray Project；In Fig. （b），the dots denote the earthquakes，and the star denotes the center of studied region我们使用远震体波层析成像方法来获取研究区内的速度结构。对于选取的远震波形，首先去均值、去线性趋势和带通滤波（0.02—1.0 Hz），然后利用波形互相关方法计算初至P波走时（Rawlinson，Kennett，2004）。为避免发震时刻误差、震源位置及震源区速度异常的影响，采用相对走时残差进行反演。以往的相对走时残差计算方法中，首先在IASP91 （Kennett，Engdahl，1991）模型下利用射线追踪计算得到理论走时，其与实际观测走时的差值即为走时残差，每个台站的走时残差减去所有台站的平均走时残差即可得到远震相对走时残差。
由于IASP91模型是一个平均层状模型，并不包含地表的沉积层速度结构，因此会对计算出的走时残差产生不利影响。从CRUST1.0模型（Laske et al，2012）提取出沉积层厚度及速度，用射线在沉积层内的走时与在IASP91模型下沉积层厚度内的走时之差来校正走时残差，使得走时残差中包含了射线在沉积层内的走时更为准确。此外，利用传统的相对走时残差计算方法，在台站分布具有较高不均匀性时对成像结果可能有不利影响。为减轻这种影响，本文使用了一种新的算法（Wang et al，2021）：将研究区划分为0.5°×0.5°的网格（图2a），先计算每个网格内台站的平均走时残差，然后计算所有网格平均走时残差的平均值，最后用每个台站的走时残差减去所有网格平均走时残差的平均值，以此作为相对走时残差，使得每个台站具有等权重的平均相对走时残差。图3a显示了原始走时残差与相对走时残差的对比，可见原始走时残差的分布明显不均，负值多于正值。原始走时残差的2σ范围（95.5%）为−2.97—3.19 s，其中σ代表正态分布中的标准差。相比之下，相对走时残差则呈现出了良好的正态分布，其2σ范围为−1.44—1.42 s。
![]()
图 3 原始走时残差与相对走时残差的对比（a）及光滑和阻尼因子测试（b）Figure 3. Comparison of the raw travel time residuals with relative ones （a） and trade-off curves for the norm of solution and RMS （root-mean-square） travel time residual （b）获得相对走时残差后，采用Zhang和Thurber （2003，2006）发展的双差层析成像程序反演P波三维速度结构。首先通过射线追踪获取IASP91模型下每一个远震在研究区底面或四个侧面的入射点。与通常的双差层析成像采用的地震观测走时数据不同，我们将所有到达同一个地震台站的射线归为一组，以此得到每一个台站的射线组和相应入射点坐标，然后将原有的地表地震台站设定为震源位置，入射点设为观测台站，形成新的观测数据。
模型参数化设置横向网格节点为0.5°×0.5°，为避免模型边缘误差的影响，模型在经纬度方向（98.5°E—124°E，26.5°N—46°N）分别向外延伸。深度方向网格节点分别为−100，0，50，120，190，270，350，440，530，620，710，800，900，1 000 km，其中−100，800，900，1 000 km仅用于插值计算，并不呈现在最终的成像结果中。初始速度模型考虑了华北地区利用接收函数获得的莫霍深度横向变化，地壳内部采用两层模型，选取IASP91模型作为上地幔深度的初始速度模型。
在反演过程中，利用LSQR （least squares QR-decomposition）方法来求解大型观测方程（Paige，Saunders，1982），以获得研究区最终的速度结构。在此过程中需要引入阻尼因子和光滑因子进行约束。经过反复测试，最终选定阻尼因子为60，光滑因子为200 （图3b）。
2. 结果
2.1 分辨率检验
我们利用检测板方法对成像结果进行分辨率测试，以检验成像结果的分辨率。在每相邻网格节点分别赋予3%和−3%的速度扰动，并添加标准差为0.1 s的随机误差，计算合成走时。采用与实际观测数据相同的计算方法得到理论走时，然后反演速度结构。通过检测板反演结果可以直观地评估研究区内的分辨率信息。
图4为0.5°×0.5°网格检测板结果，可见所得速度图像在研究区的绝大部分地区都较好地恢复了输入的速度异常，表明本文的成像结果可以分辨出横向尺度约0.5°的速度异常，浅部的分辨率受到了台站分布的影响。中国地震科学台阵三期流动台站极大地提高了鄂尔多斯地块及中部地块的台站密度，有效地改善了浅部的分辨效果（图4a）。
![]()
图 4 不同深度h下0.5°×0.5°网格检测板结果Figure 4. Results of checkerboard test for 0.5°×0.5° grid spacing with different depths h（a） h＝50 km；（b） h＝120 km；（c） h＝190 km；（d） h＝270 km；（e） h＝350 km；（f） h＝440 km；（g） h＝530 km；（h） h＝620 km；（i） h＝710 km2.2 三维P波速度结构
图5展示了不同深度处P波速度横切片，图6给出了纵剖面的P波速度结果。
![]()
图 5 研究区不同深度h处的横切片P波速度结果Figure 5. P-wave velocity images in map view at the different depths h in the studied area（a） h＝50 km；（b） h＝120 km；（c） h＝190 km；（d） h＝270 km；（e） h＝350 km； （f） h＝440 km；（g） h＝530 km；（h） h＝620 km；（i） h＝710 km![]()
图 6 纵剖面位置（a）及其P波速度结果（b−e）图（b）−（e）中，虚线代表410 km及660 km间断面Figure 6. Location of the four vertical cross sections （a） and P-wave velocity results along the sections （b−e）The dashed lines represent the 410 km and 660 km discontinuities in Figs. （b）−（e）在地壳及上地幔深度处，观测到的速度异常与已知的地质学地块边界有较好的一致性（图5a，b）。在鄂尔多斯地块下方存在显著的高速异常，与鄂尔多斯地块的克拉通属性相吻合。高速异常在西部和东部的深度约为180 km，在中部地区可达300 km （图6c），表明在鄂尔多斯地块中部仍存在巨厚的克拉通核；而在东北部，高速异常仅延伸至150 km （图6e）。在50 km深度处高速体的北缘与鄂尔多斯地块边界非常吻合（图5a）。在50—440 km深度，鄂尔多斯东北部表现为低速异常，该低速异常可延伸至中部地块中部及大同火山下方（图5）。
华北克拉通中部地块的南部上地幔低速异常并不明显，与前人结果形成鲜明对比（Lei，2012；Jiang et al，2013；Guo，Chen，2017）。北部尤其是大同火山下方的上地幔中存在强烈的低速异常，其扩展方向大致为NNW−SSE，这些低速异常的分布位置与该处的断陷及火山群相对应（图5a）。在50 km深度处，该低速体向北可延伸至鄂尔多斯地块东北部地区，向南可延伸至山西断陷带中部，向东则到达了华北克拉通东部地块。大同火山下方的低速异常可一直延伸至地幔转换带深度（图6b）。在50—270 km深度，山西断陷带中部表现为高速异常（图5），且该高速体在西侧与鄂尔多斯地块相连。渭河盆地附近在50 km深度表现为高速异常（图5a），而运城盆地附近在50 km深度存在显著的局部低速异常，这可能意味着其断陷作用已超过莫霍面深度。鄂尔多斯东南缘在50 km深度表现出高低速相间的复杂结构，可能预示了该地区复杂的岩石圈演化。
华北克拉通东部地块在上地幔表现为大范围的低速异常，这与中生代以来东部地块遭受强烈破坏与改造，以及中、新生代火成岩广泛侵入与喷发相符合，且这些低速异常与大同火山下方的低速异常相连（图5）。华北克拉通东部的河淮盆地在50—120 km深度大部分地区表现为弱高速（图5a，b），与该地区较低的岩石圈地幔密度以及地表热流相一致（Xia et al，2020）。其内部地震相对较弱且分布均匀，表明其现今的岩石圈结构可能相对稳定。河淮盆地在190—350 km深度存在局部的高速异常（图5c，d，e），推测可能为岩石圈拆离体。紧邻河淮盆地北部的渤海湾盆地在50—270 km深度呈现为强烈的低速异常（图5a，b）。河淮盆地与渤海湾盆地上地幔速度异常的明显对比，与河淮盆地晚中生代岩浆活动较弱且地温场较低、渤海湾盆地晚中生代岩浆活动强烈且地温场较高的特征相符合（刘池洋等，2018）。在440—530 km深度，东部地块东部出现明显的高速异常，可能为太平洋滞留板片（图5f，g）。
2.3 恢复性试验测试
为了进一步检验本文成像结果的可靠性，我们进行恢复性分辨率测试。采用的步骤与检测板测试完全一致，仅输入模型有所差异（图7a）。图7b和7d为输入模型，青藏高原东北缘下方设置了−3%的低速异常，鄂尔多斯地块下方设置了3%的高速异常，东部地块下方为−3%的低速异常，低速异常下方的高速异常代表了太平洋滞留板片。图7c和7e为反演后得到的输出模型，可以看到尽管速度异常值及形态存在一定差异，但输出模型仍很好地恢复了输入模型中的高低速异常，表明了成像结果的可靠性。
![]()
图 7 恢复性分辨率试验结果（a） 异常输入分布，图中HV1和HV2代表3%高速异常，LV1和LV2代表−3%低速异常，分别用蓝色和橘色方框表示，且LV2和HV2在116°E以东重合；（b，d） 输入模型；（c，e） 输出模型Figure 7. Restoring resolution test results（a） The distribution of input anomalies，where HV1 and HV2 denote 3% high anomalies，and LV1 and LV2 denote −3% low anomalies，noting that LV1 and LV2 overlap to the east of 116°E；（b，d） The input models；（c，e） The output models3. 讨论
3.1 与已有成像结果的对比
本文成像结果表现出与前人结果大体上的相似性，但在许多细节上存在差异。例如，与大多数人的结果相似，本文成像结果显示在鄂尔多斯地块下方存在显著的高速异常，鄂尔多斯东北部及东部地区主要以低速异常为主。但不同于Huang和Zhao （2006）的鄂尔多斯地块下方高速异常向西倾斜至300 km的结果，本文结果显示高速异常具有向南倾斜的特征，其中东北部厚约150 km，西部和东部厚约180 km，中部局部地区可达300 km。这一结果与Lei （2012）获得的约400 km厚的高速异常也存在差异。毛慧慧等（2016）的远震P波层析成像结果显示高速异常厚度约为150—200 km。而面波层析成像结果显示的高速异常厚度相对体波层析成像结果普遍偏小，例如：An等（2009）的结果显示鄂尔多斯下方岩石圈厚度约为130 km；Li等（2018）的结果显示鄂尔多斯下方高速异常厚度约为200 km，并向东逐渐减薄至约150 km；Yao等（2020）的结果显示鄂尔多斯岩石圈厚度为150—170 km，在鄂尔多斯东北缘的厚度为80—130 km。
一些远震成像结果在地幔转换带中观测到了高速异常，并将其解释为太平洋滞留板片，但其分布形态存在争议。本文成像结果显示的滞留板片厚度相对前人结果（Huang，Zhao，2006；Tian et al，2009；Lei，2012）较薄。前人结果显示大同火山下方普遍呈低速特征，本文结果与Huang和Zhao （2006）结果相似，均观察到向西倾斜的低速异常，但本文结果中的低速异常西缘下延深度相对较浅，在410 km间断面附近（图6b）。
这些结果的差异性可能是由于使用的成像方法以及数据的不同所致。前人在华北克拉通地区主要使用固定台以及一些小尺度流动台站，成像分辨效果受到了台站分布密度的影响。本文成像研究利用了分布于鄂尔多斯地块及中部地块的中国地震台阵探测第三期流动台站，能够有效地显示出该地区的精细结构。
3.2 鄂尔多斯地块岩石圈特征
成像结果揭示鄂尔多斯地块下方的岩石圈厚度存在较大的横向差异。鄂尔多斯地块下方的高速异常代表了相对稳定的克拉通岩石圈（图5a-e和图6c-e），高速异常在中部局部地区可达300 km。本文结果观察到的巨厚岩石圈根与鄂尔多斯地块内部弱各向异性的特征（Zhao et al，2013）相符合。
鄂尔多斯下方的高速异常体在东北部地区减薄至大约150 km （图6e），其下方存在显著的低速异常，该低速异常向东与大同火山地区及华北盆地下方的低速异常相连。这些大范围的低速异常可能代表了太平洋滞留板片上方热的大地幔楔，以及板片前缘的局部对流。
Dong等（2014）利用大地电磁方法在鄂尔多斯地块北部及河套断陷的上地幔中发现一大尺度高导体，本文观察到的低速异常位置与其相近。对鄂尔多斯地块东北缘的晚中新世玄武岩成分分析表明，其来源于上地幔，可能是软流圈与岩石圈地幔相互作用的产物（邹和平等，2010）。鄂尔多斯地块东北缘地震活动水平较强，在杭锦旗、和林格尔和大同火山区存在新生代玄武岩出露（图1）。这些证据均表明鄂尔多斯地块东北缘上地幔中存在热物质上涌，可能预示着该地区的克拉通改造。故本文推测，受鄂尔多斯地块东北缘地幔热物质上涌的影响，地块东北部及北部的岩石圈受到了加热和改造，导致了岩石圈的减薄，但岩石圈仍保留相对稳定的特征。
鄂尔多斯地块西部的高速异常厚度约为180 km （图6c），岩石圈厚度比中部明显减薄，在190—270 km深度上表现出局部弱低速异常，在350 km深度表现出明显的低速异常（图5c−e）。青藏高原东北缘上地幔表现出低速异常（图5a−e）。剪切波分裂结果显示，青藏高原东北缘至鄂尔多斯地块西南缘的各向异性较强，变化相对连续，与鄂尔多斯地块中部较小的各向异性形成鲜明对比（Yu，Chen，2016；Chang et al，2017；Liu et al，2020，2021）。接收函数研究揭示，青藏高原东北缘的岩石圈厚度约为90 km （Shen et al，2015）。自上新世开始，青藏高原的急剧抬升扩展至其东北缘（Tapponnier et al，2001）。GPS结果显示地壳运动速度从青藏高原东北缘至鄂尔多斯地块有明显的降低（Niu et al，2005）。上述证据表明，具有较厚的岩石圈根和较大的力学强度的鄂尔多斯地块在很大尺度上阻止了青藏高原向东北部横向扩展。鄂尔多斯地块西缘的地震活动性较高（图1），主压应力方向为NE−SW向（盛书中等，2015）。地热研究表明，青藏高原上地幔温度相对鄂尔多斯地块较高（Sun et al，2013）。鄂尔多斯地块西部与其中部地区相比具有相对较薄的岩石圈厚度（图6c），并且其下方存在明显的低速异常（图5d，e），我们推测青藏高原上地幔的低速热物质受横向挤出作用的影响，向外扩张至鄂尔多斯西缘地区，导致了该地区的岩石圈被部分改造和减薄。
鄂尔多斯地块相对太古宙克拉通有较高的热流值（62 mW/m2）。热流主要来源于地球内部（Pollack et al，1993），可以作为判断克拉通是否保持稳定性的重要依据（He，2015）。研究表明，太古宙克拉通通常具有较低热流值，平均为40 mW/m2 （Ballard，Pollack，1987），明显低于全球大陆热流均值65 mW/m2 （Pollack et al，1993）。鄂尔多斯地块的热流值明显高于正常的太古宙克拉通，这可能意味着已经存在部分失稳及破坏。在本文的成像结果中，鄂尔多斯地块的岩石圈厚度沿横向变化较大，在中部可达300 km，在西部和东部减薄至约180 km，在东北部及北部甚至减薄至约150 km。鄂尔多斯地块内部岩石圈厚度较大的横向差异性、较高的热流值以及地表存在的差异性隆升等表明鄂尔多斯地块的稳定性可能正在逐步减弱。
3.3 鄂尔多斯地块周边地区的上地幔结构
山西断陷带是世界上最大的新生代大陆裂谷之一，也是中国东部著名的地震构造带，位于东部已经被改造的岩石圈与西部相对稳定的鄂尔多斯地块的转换带处（Bao et al，2013；Xu et al，2020）。山西断陷带中段表现出高速特征，南部和北部为低速异常。南部的局部低速异常深度较浅（图5a−d），表明可能存在局限于对地壳尺度和上地幔顶部附近的改造，该低速异常与Yin等（2017）在该地区观测到的低电阻异常的分布形态具有一致性。北部的低速异常范围较大，深度较深（图5a−c）。地球化学结果显示山西断陷带北部与东部地块的新生代玄武岩组分存在相似性，这些组分主要来源于软流圈（Xu et al，2004），表明山西断陷带北部与东部地块的岩石圈可能经历了相似的动力学过程。
大同火山位于山西断陷带北部。本文成像结果显示大同火山地区上地幔顶部呈现为强烈的低速异常，这与宋晓燕等（2020）的结果存在一致性（图5）。大同火山及周围区域在50—440 km深度存在大范围的低速异常，270 km以下深度的低速异常主要出现在华北克拉通中部地块及西部地块的北部，而190 km以上深度的低速异常则主要出现在华北克拉通东部地块以及中部地块。
有研究认为华北克拉通中部地块下方的低速异常可能来源于深部地幔柱（Zhao et al，2009；Lei，2012）。地幔柱起源于上下地幔边界或核幔边界，在本文成像结果中未观测到来自深部的地幔上涌以及类似地幔柱形态的低速异常。热点是地幔柱的地表表现（Morgan，1971；Wilson，1973），但该地区也缺乏热点轨迹证据。此外，华北克拉通的火成岩中也没有表现出与新生代地幔柱活动相关的确切证据（Menzies et al，2007）。大陆洋岛玄武岩的成因可能与地幔柱存在联系（Morgan，1981）。但在华北克拉通地区，Yang和Zhou （2013）的地球化学研究结果认为大陆洋岛玄武岩的产生不需要地幔柱模型来支撑，辉石的熔融物即可以解释洋岛玄武岩的组分特征。因此，本文成像结果中的低速异常可能难以用地幔柱模型来解释。
在华北地区最东部的440—530 km深度处，地幔转换带中存在一显著的高速异常（图5f，g）。前人成像结果显示该位置附近及以东存在太平洋滞留板片（Huang，Zhao，2006；Tian et al，2009）。接收函数结果中，华北克拉通在约115°E以东的地区地幔转换带较厚（Chen，Ai，2009）。结合本文440—530 km成像结果（图5f，g），我们推测在转换带中的高速体可能是太平洋滞留板片，其前缘可达116°E附近，基本位于南北重力梯度带以东。滞留板片的厚度约为150 km （图6c），对比前人的成像结果（Huang，Zhao，2006；Wei et al，2012）明显较薄，并且位置主要分布在地幔转换带中上部。
在地幔转换带底部，太平洋滞留板片下方存在显著的低速异常（图5h，i），这种滞留板片下方存在低速异常的现象在其它研究中也被揭示（Tang et al，2014；Liu et al，2017）。尽管An等（2009）认为板片下方的低速异常可以逃逸并导致软流圈地幔中出现热物质，但在本文的纵剖面结果中未发现该低速热异常有明显的上涌通道（图6b，c）。
在本文的成像结果中，太平洋滞留板片前方和上方的上地幔中存在大范围的低速异常（图6b，c）。板片前方，特别是鄂尔多斯地块北部上地幔的大范围低速异常成因目前尚不明确（图5a−f）。位于太平洋板片前方和上方的低速异常可能代表了上地幔熔融物。晚中生代以来，由于古太平板片的俯冲及后撤，华北克拉通东部地块经历了大规模岩浆活动以及伸展变形（Zhu et al，2012a）。东部地块上地幔温度相对较高（Sun et al，2013）。我们推测，东部地块上地幔的低速异常可能与晚中生代以来东部伸展环境下深部地幔上涌及减压熔融存在联系。
同时，这些大范围的低速异常也不能排除是受到了太平洋滞留板片潜在脱水的影响。尽管有一些研究认为，华北克拉通东部岩石圈目前的水含量较低（Xia et al，2019），然而这种低含水量可能是受到了岩浆作用对岩石圈加热的影响（Yang et al，2008）。对于软流圈及以下深度的上地幔水含量目前尚无定论。女山地区新生代玄武岩中的橄榄岩具有相对高的水含量，且该地区的氢氧同位素研究表明交代流体来自于海洋板块（Xia et al，2004；郝艳涛等，2005）。这意味着该地区的岩浆活动可能受到了太平洋板片脱水的影响。赤峰地区的地球化学特征也与太平洋滞留板片存在联系（Wang et al，2015）。一些研究认为，现今观测到的太平洋滞留板片前缘俯冲年龄约为30 Ma，滞留在转换带中的时间不超过10—20 Ma，而更古老的板片已沉入下地幔（Liu et al，2017；Tao et al，2018）。俯冲引起的地幔流数值模拟结果表明太平洋板块的俯冲扰动可以促使地幔转换带释放流体（Wang et al，2016）。流体的加入会导致地幔黏度降低（Hirth，Kohlstedt，2003），可能加速上地幔中产生对流。对流可以向岩石圈底部传递热量（He，2015）。这种温度升高及对流的现象同样存在于Guo等（2016）的热化学研究中。滞留板片脱水、板片前缘的对流及温度升高可能导致了华北克拉通东部地块、中部地块及鄂尔多斯地块北部上地幔大范围低速异常。He （2015）的数值模拟研究认为，地幔对流引起岩石圈减薄需要数十个百万年。尽管现今的太平洋板片俯冲引发的对流可能向华北克拉通的岩石圈底部传递了热量，但由于实际俯冲时间可能较短，因此可能不足以使华北克拉通岩石圈发生大范围明显的减薄破坏。
华北克拉通东部和大同火山区周边的岩石圈在晚中生代至新生代发生了减薄（Xu et al，2005；邹和平等，2010；Ai et al，2020）。地球化学研究显示，大同火山区的玄武岩以及汉诺坝中新世玄武岩具有相似的软流圈来源（Xu et al，2005）。Xu等（2020）认为大同盆地较高的沉积速率可能与8—10 Ma以来该地区下方软流圈上涌及被动裂谷作用相关。我们推测，俯冲板片前缘的局部对流向岩石圈提供了热量。在大同火山地区新生代伸展作用背景下（Zhang et al，1998；Xu et al，2005；邹和平等，2010；Yao et al，2020），岩石圈或软流圈的熔融物沿着软弱带上涌并形成了包括大同火山在内的火山群。这些部分熔融物在大地电磁（Zhang et al，2016）及地球化学（Xu，2007）结果中也有显示。
4. 结论
本文采用研究区内943个台站记录到的404个远震事件的15万3 963个走时数据，使用远震P波层析成像方法获取了华北克拉通地区高分辨率的上地幔速度结构。
本文结果显示：鄂尔多斯地块下方存在平均厚度约为200 km的显著高速异常，在局部地区高速异常可达300 km深，表明鄂尔多斯地块总体上保持了克拉通特性；高速异常在东北部减薄至约150 km，推测热的地幔上涌对鄂尔多斯地块东北部的岩石圈进行了加热和改造，导致了克拉通岩石圈的局部活化；鄂尔多斯地块西部高速异常延展深度约180 km，其下方存在低速异常，推测其上地幔可能受到了青藏高原东北缘深部热物质扩展的影响，岩石圈地幔已受到了部分改造。
华北克拉通东部地块、中部地块及鄂尔多斯地块东北部的上地幔表现出大范围的低速异常，这可能与太平洋板块俯冲脱水以及俯冲板片前缘的局部对流有关。对流向岩石圈底部传递热量，在该地区的伸展背景下，岩石圈或软流圈的熔融物质沿着软弱带上涌并形成了包括大同火山在内的火山群。
中国地震局地球物理研究所地震科学国际数据中心为本研究提供了地震波形数据，中国地震台网中心提供了地震目录，作者在此表示衷心的感谢。
-
![]()
图 1 青藏高原南部地区缝合带和地体的分布（蓝色矩形为本文研究区域）
Figure 1. Suture zones and terranes distribution in the southern Tibetan Plateau（The blue box is the research area of this study）
![]()
图 2 反双台法几何路径示意图
（a） 理想条件下；（b） 实际条件下
Figure 2. Schematic diagram for the geometry of the reverse two-station method
（a） An ideal recording geometry；（b） A more practical geometry
![]()
图 3 本文反双台法计算所用台站及地震事件的分布
Figure 3. Distribution of seismic stations and earthquakes calculated by reverse two-station method in this paper
![]()
图 4 Pn波与Lg波窗口的拾取实例
Figure 4. An example showing the Pn-waves and the window of Lg-waves
![]()
图 7 1°×1° （a）和0.5°×0.5° （b）网格下的棋盘格测试结果
Figure 7. The checkerboard test results with gridding of 1°×1° （a） and 0.5°×0.5° （b）
![]()
图 8 青藏高原南部地区Lg波Q0值层析成像
Figure 8. Lg-wave Q0 value tomography of the southern Tibetan Plateau
![]()
图 9 青藏高原南部熔融通道（红色实线区域）的推断示意图
Ⅰ ，Ⅱ ，Ⅲ和Ⅳ代表四个由雅鲁藏布江缝合带和亚东—谷露裂谷划分的地块
Figure 9. The inference of melting channels in the southern Tibetan Plateau
The area delineated by red solid lines denotes melting channels，and the symbols of Ⅰ ，Ⅱ ，Ⅲ and Ⅳ represent four blocks separated by the Yarlung Zangbo suture and Yadong-Gulu rift
-
艾印双,郑天愉. 1997. 青藏高原地震活动及其构造背景[J]. 地球物理学进展,12(2):30–40. Ai Y S,Zheng T Y. 1997. Seismic activity in Tibetan Plateau and its tectonic implication[J]. Progress in Geophysics,12(2):30–40 (in Chinese).
白嘉启,梅琳,杨美伶. 2006. 青藏高原地热资源与地壳热结构[J]. 地质力学学报,12(3):354–362. doi: 10.3969/j.issn.1006-6616.2006.03.010 Bai J Q,Mei L,Yang M L. 2006. Geothermal resource and crustal thermal structure of the Qinghai-Tibet Plateau[J]. Journal of Geomechanics,12(3):354–362 (in Chinese).
才巴央增,赵俊猛. 2018. 藏南裂谷系的研究综述[J]. 地震研究,41(1):14–21. doi: 10.3969/j.issn.1000-0666.2018.01.002 Caibayangzeng,Zhao J M. 2018. A summary of researches on southern Tibet rift system[J]. Journal of Seismological Research,41(1):14–21 (in Chinese).
冯昭贤,赵文津. 1997. INDEPTH与INDEPTH-MT项目简介[J]. 现代地质,11(3):363–365. Feng Z X,Zhao W J. 1997. Simple introduction to project INDEPTH and subproject NDEPTH-MT[J]. Geoscience,11(3):363–365 (in Chinese).
何静,吴庆举,李永华,雷建设. 2017. 天然地震Lg波衰减研究进展及其在中国大陆地区的应用[J]. 地球物理学进展,32(2):466–475. doi: 10.6038/pg20170204 He J,Wu Q J,Li Y H,Lei J S. 2017. Developments of earthquake Lg-wave attenuation study and its application in the continental China[J]. Progress in Geophysics,32(2):466–475 (in Chinese).
李仕虎,黄宝春,朱日祥. 2012. 青藏高原东南缘构造旋转的古地磁学证据[J]. 地球物理学报,55(1):76–94. doi: 10.6038/j.issn.0001-5733.2012.01.008 Li S H,Huang B C,Zhu R X. 2012. Paleomagnetic constraints on the tectonic rotation of the southeastern margin of the Tibetan Plateau[J]. Chinese Journal of Geophysics,55(1):76–94 (in Chinese). doi: 10.1002/cjg2.1702
李永华,吴庆举. 2007. 中国地学热点研究区几个地学问题的探讨[J]. 国际地震动态,(9):11–19. doi: 10.3969/j.issn.0253-4975.2007.09.003 Li Y H,Wu Q J. 2007. Study on some “Hot Spot” areas of geo-sciences in China[J]. Recent Developments in World Seismology,(9):11–19 (in Chinese).
李娱兰. 2019. 青藏高原地震活动性及冈底斯成矿带东段上地幔顶部Pn波速度结构[D]. 北京: 中国地震局地球物理研究所: 1–8. Li Y L. 2019. Seismicity of the Qinghai-Tibet Plateau and Velocity Structure of Pn Wave on the Top of the Upper Mantle in the Eastern Gangdise Metallogenic Belt[D]. Beijing: Institute of Geophysics, China Earthquake Administration: 1–8 （in Chinese）.
刘建华,刘福田,阎晓蔚,胥颐,郝天珧. 2004. 华北地区Lg尾波衰减研究:Lg尾波Q0地震成像[J]. 地球物理学报,47(6):1044–1052. doi: 10.3321/j.issn:0001-5733.2004.06.017 Liu J H,Liu F T,Yan X W,Xu Y,Hao T Y. 2004. A study of Lg coda attenuation beneath North China:Seismic imaging Lg coda Q0[J]. Chinese Journal of Geophysics,47(6):1044–1052 (in Chinese).
吴中海,叶培盛,王成敏,张克旗,赵华,郑勇刚,尹金辉,李虎侯. 2015. 藏南安岗地堑的史前大地震遗迹、年龄及其地质意义[J]. 地球科学:中国地质大学学报,40(10):1621–1642. Wu Z H,Ye P S,Wang C M,Zhang K Q,Zhao H,Zheng Y G,Yin J H,Li H H. 2015. The relics,ages and significance of prehistoric large earthquakes in the Angang graben in south Tibet[J]. Earth Science:Journal of China University of Geosciences,40(10):1621–1642 (in Chinese). doi: 10.3799/dqkx.2015.147
喻成,乔学军,王伟,史永明. 2014. 亚东—谷露裂谷带与块体运动的特征[J]. 大气测量与地球动力学,34(2):36–40. Yu C,Qiao X J,Wang W,Shi Y M. 2014. Characteristics of crust deformation in Yadong-Gulu rift and its surrounding areas with GPS data[J]. Journal of Geodesy and Geodynamics,34(2):36–40 (in Chinese).
张戈铭,李细兵,郑晨,宋晓东. 2019. 青藏高原中东部地壳和上地幔顶部P波层析成像[J]. 地震学报,41(4):411–424. doi: 10.11939/jass.20190003 Zhang G M,Li X B,Zheng C,Song X D. 2019. Crustal and uppermost mantle velocity structure beneath the central eastern Tibetan Plateau from P-wave tomography[J]. Acta Seismologica Sinica,41(4):411–424 (in Chinese).
张衡,赵俊猛,徐强. 2011. 西藏东部地区层析成像及东南部裂谷成因讨论[J]. 科学通报,56(27):2328–2334. Zhang H,Zhao J M,Xu Q. 2011. Seismic P-wave tomography in eastern Tibet:Formation of the rifts[J]. Chinese Science Bulletin,56(23):2450–2455. doi: 10.1007/s11434-011-4577-x
张红亮. 2010. 藏北水热活动的地质背景研究[D]. 北京: 中国地质大学（北京）: 13–17. Zhang H L. 2010. Study on the Geological Background of Hydrothermal Activities in North Tibet[D]. Beijing: China University of Geosciences （Beijing）: 13–17 （in Chinese）.
张锦玲,朱新运,马起杨. 2019. 宁夏地区Lg波衰减及场地响应特征[J]. 地震学报,41(4):425–434. doi: 10.11939/jass.20180134 Zhang J L,Zhu X Y,Ma Q Y. 2019. Lg-wave attenuation and site response in Ningxia region[J]. Acta Seismologica Sinica,41(4):425–434 (in Chinese).
赵连锋,谢小碧,王卫民,姚振兴. 2018. 中国东北和朝鲜半岛地区地壳Lg波宽频带衰减模型[J]. 地球物理学报,61(3):856–871. doi: 10.6038/cjg2018L0394 Zhao L F,Xie X B,Wang W M,Yao Z X. 2018. A broadband crustal Lg wave attenuation model in Northeast China and the Korean Peninsula[J]. Chinese Journal of Geophysics,61(3):856–871 (in Chinese).
赵文津,吴珍汉,史大年,熊嘉育,薛光琦,宿和平,胡道功,叶培盛. 2008. 国际合作INDEPTH项目横穿青藏高原的深部探测与综合研究[J]. 地球学报,29(3):328–342. doi: 10.3321/j.issn:1006-3021.2008.03.007 Zhao W J,Wu Z H,Shi D N,Xiong J Y,Xue G Q,Su H P,Hu D G,Ye P S. 2008. Comprehensive deep profiling of Tibetan Plateau in the INDEPTH project[J]. Acta Geoscientica Sinica,29(3):328–342 (in Chinese).
周蕙兰. 1990. 地球内部物理[M]. 北京: 地震出版社: 218–237. Zhou H L. 1990. Physics of the Earth’s Interior[M]. Beijing: Seismological Press: 218–237 （in Chinese）.
周连庆,赵翠萍,修济刚,陈章立. 2008a. 川滇地区Lg波Q值层析成像[J]. 地球物理学报,51(6):1745–1752. Zhou L Q,Zhao C P,Xiu J G,Chen Z L. 2008a. Tomography of Q Lg in Sichuan-Yunnan zone[J]. Chinese Journal of Geophysics,51(6):1745–1752 (in Chinese).
周连庆,赵翠萍,修济刚,陈章立,郑斯华. 2008b. 利用天然地震研究地壳Q值的方法和进展[J]. 国际地震动态,(2):1–11. Zhou L Q,Zhao C P,Xiu J G,Chen Z L,Zheng S H. 2008b. Methods and developments of research on crustal Q value by using earthquakes[J]. Recent Developments in World Seismology,(2):1–11 (in Chinese).
朱新运. 2016. 华北盆地Lg波衰减及台站场地响应特征[J]. 地球科学,41(12):2109–2117. Zhu X Y. 2016. Characteristics of Lg wave attenuation and site response in North China basin[J]. Earth Science,41(12):2109–2117 (in Chinese).
Al-Damegh K,Sandvol E,Al-Lazki A,Barazangi M. 2004. Regional seismic wave propagation (Lg and Sn) and Pn attenuation in the Arabian Plate and surrounding regions[J]. Geophys J Int,157(2):775–795. doi: 10.1111/j.1365-246X.2004.02246.x
Armijo R,Tapponnier P,Han T. 1989. Late Cenozoic right-lateral strike-slip faulting in southern Tibet[J]. J Seismol Res,94(B3):2787–2838.
Bao X Y,Sandvol E,Ni J,Hearn T,Chen Y J,Shen Y. 2011a. High resolution regional seismic attenuation tomography in eastern Tibetan Plateau and adjacent regions[J]. Geophys Res Lett,38(16):L16304.
Bao X Y,Sandvol E,Zor E,Sakin S,Mohamad R,Gök R,Mellors R,Godoladze T,Yetirmishli G,Türkelli N. 2011b. Pg attenuation tomography within the northern Middle East[J]. Bull Seismol Soc Am,101(4):1496–1506. doi: 10.1785/0120100316
Brown L D,Zhao W,Nelson K D,Hauck M,Alsdorf D,Ross A,Cogan M,Clark M,Liu X,Che J. 1996. Bright spots,structure,and magmatism in southern Tibet from INDEPTH seismic reflection profiling[J]. Science,274(5293):1688–1690. doi: 10.1126/science.274.5293.1688
Campillo M,Plantet J L,Bouchon M. 1985. Frequency-dependent attenuation in the crust beneath central France from Lg waves:Data analysis and numerical modeling[J]. Bull Seismol Soc Am,75(5):1395–1411.
Chen Q Z,Freymueller J T,Wang Q,Yang Z Q,Xu C J,Liu J N. 2004. A deforming block model for the present‐day tectonics of Tibet[J]. J Geophys Res,109(B1):B01403.
England P C,Houseman G A. 1988. The mechanics of the Tibetan Plateau[J]. Phil Trans R Soc Lond A,326(1589):301–320. doi: 10.1098/rsta.1988.0089
Fan G W,Lay T. 2002. Characteristics of Lg attenuation in the Tibetan Plateau[J]. J Geophys Res,107(B10):2256.
Fan G W,Lay T. 2003. Strong Lg attenuation in the northern and eastern Tibetan Plateau measured by a two-station/two-event stacking method[J]. Geophys Res Lett,30(10):1530.
Ford S R,Dreger D S,Mayeda K,Walter W R,Malagnini L,Phillips W S. 2008. Regional attenuation in northern California:A comparison of five 1D Q methods[J]. Bull Seismol Soc Am,98(4):2033–2046. doi: 10.1785/0120070218
Gök R,Sandvol E,Türkelli N,Seber D,Barazangi M. 2003. Sn attenuation in the Anatolian and Iranian Plateau and surrounding regions[J]. Geophys Res Lett,30(24):8042.
Herrmann R B. 1980. Q estimates using the coda of local earthquake[J]. Bull Seismol Soc Am,70(4):447–468.
Hu S B,He L J,Wang J Y. 2000. Heat flow in the continental area of China:A new data set[J]. Earth Planet Sci Lett,179(2):407–419. doi: 10.1016/S0012-821X(00)00126-6
Li S H,Unsworth M J,Booker J R,Wei W B,Tan H D,Jones A G. 2003. Partial melt or aqueous fluid in the mid-crust of southern Tibet?Constraints from INDEPTH magnetotelluric data[J]. Geophys J Int,153(2):289–304. doi: 10.1046/j.1365-246X.2003.01850.x
McNamara D T,Owens T J,Walter W R. 1996. Propagation characteristics of Lg across the Tibetan Plateau[J]. Bull Seismol Soc Am,86(2):457–469.
Mitchell B J. 1980. Frequency dependence of shear wave internal friction in the continental crust of eastern North America[J]. J Geophys Res,85(B10):5212–5218. doi: 10.1029/JB085iB10p05212
Myers S C,Beck S,Zandt G,Wallace T. 1998. Lithospheric-scale structure across the Bolivian Andes from tomographic images of velocity and attenuation for P and S waves[J]. J Geophys Res,103(B9):21233–21252. doi: 10.1029/98JB00956
Nelson K D,Zhao W J,Brown L D,Kuo J,Che J K,Liu X W,Klemperer S L,Makovsky Y,Meissner R,Mechie J,Kind R,Wenzel F,Ni J,Nabelek J,Leshou C,Tan H D,Wei W B,Jones A G,Booker J,Unsworth M,Kidd W S F,Hauck M,Alsdorf D,Ross A,Cogan M,Wu C D,Sandvol E,Edwards M. 1996. Partially molten middle crust beneath southern Tibet:Synthesis of project INDEPTH results[J]. Science,274(5293):1684–1688. doi: 10.1126/science.274.5293.1684
Ni J,Barazangi M. 1983. High-frequency seismic wave propagation beneath the Indian shield,Himalayan arc,Tibetan Plateau and surrounding regions:High uppermost mantle velocities and efficient Sn propagation beneath Tibet[J]. Geophys J R astr Soc,72(3):665–689. doi: 10.1111/j.1365-246X.1983.tb02826.x
Nuttli O W. 1980. The excitation and attenuation of seismic crustal phases in Iran[J]. Bull Seismol Soc Am,70(2):469–485.
Paige C C,Saunders M A. 1982. LSQR:An algorithm for sparse linear equations and sparse lest squares[J]. ACM Trans Math Software,8(1):43–71. doi: 10.1145/355984.355989
Phillips W S,Hartse H E,Taylor S R,Randall G E. 2000. 1 Hz Lg Q tomography in central Asia[J]. Geophys Res Lett,27(20):3425–3428. doi: 10.1029/2000GL011482
Press F,Ewing M. 1952. Two slow surface waves across North America[J]. Bull Seismol Soc Am,42(3):219–228. doi: 10.1785/BSSA0420030219
Reese C C,Rapine R R,Ni J F. 1999. Lateral variation of Pn and Lg attenuation at the CDSN station LSA[J]. Bull Seismol Soc Am,89(1):325–330. doi: 10.1785/BSSA0890010325
Ringdal F,Marshall P D,Alewine R W. 1992. Seismic yield determination of Soviet underground nuclear explosions at the Shagan River test site[J]. Geophys J Int,109(1):65–77. doi: 10.1111/j.1365-246X.1992.tb00079.x
Ruzaikin A I,Nersesov I L,Khalturin V I,Molnar P. 1977. Propagation of Lg and lateral variations in crustal structure in Asia[J]. J Geophys Res,82(2):307–316. doi: 10.1029/JB082i002p00307
Sato R. 1967. Attenuation of seismic waves[J]. J Phys Earth,15(2):32–61. doi: 10.4294/jpe1952.15.32
Shin T C,Herrmann R B. 1987. Lg attenuation and source studies using 1982 Miramichi data[J]. Bull Seismol Soc Am,77(2):384–397.
Singh C,Shekar M,Singh A,Chadha R K. 2012. Seismic attenuation characteristics along the Hi-CLIMB profile in Tibet from Lg Q inversion[J]. Bull Seismol Soc Am,102(2):783–789. doi: 10.1785/0120110145
Singh S,Herrmann R B. 1983. Regionalization of crustal coda Q in the continental United States[J]. J Geophys Res,88(B1):527–538. doi: 10.1029/JB088iB01p00527
Tapponnier P,Mercier J L,Armijo R,Han T L,Zhou J. 1981. Field evidence for active normal faulting in Tibet[J]. Nature,294(5840):410–414. doi: 10.1038/294410a0
Tapponnier P,Peltzer G,Le Dain A Y,Armijo R,Cobbold P R. 1982. Propagating extrusion tectonics in Asia:New insights from simple experiments with plasticine[J]. Geology,10(12):611–616. doi: 10.1130/0091-7613(1982)10<611:PETIAN>2.0.CO;2
Tapponnier P,Xu Z Q,Roger F,Meyer B,Arnaud N,Wittlinger G,Yang J S. 2001. Oblique stepwise rise and growth of the Tibet Plateau[J]. Science,294(5547):1671–1677. doi: 10.1126/science.105978
USGS. 2019. Search earthquake catalog[EB/OL]. [2019-09-01]. https://earthquake.usgs.gov/earthquakes/search.
Wei W B,Unsworth M,Jones A,Booker J,Tan H D,Nelson D,Chen L S,Li S H,Solon K,Bedrosian P,Jin S,Deng M,Ledo J,Kay D,Roberts B. 2001. Detection of widespread fluids in the Tibetan crust by magnetotelluric studies[J]. Science,292(5517):716–719. doi: 10.1126/science.1010580
Xie J. 2002. Lg Q in the eastern Tibetan Plateau[J]. Bull Seismol Soc Am,92(2):871–876. doi: 10.1785/0120010154
Xie J,Mitchell B J. 1990a. Attenuation of multiphase surface waves in the Basin and Range Province,part I:Lg and Lg coda[J]. Geophys J Int,102(1):121–137. doi: 10.1111/j.1365-246X.1990.tb00535.x
Xie J,Mitchell B J. 1990b. A back-projection method for imaging large-scale lateral variations of Lg coda Q with application to continental Africa[J]. Geophys J Int,100(1):161–181.
Xie J,Gok R,Ni J,Aoki Y. 2004. Lateral variations of crustal seismic attenuation along the INDEPTH profiles in Tibet from Lg Q inversion[J]. J Geophys Res,109(B10):B10308.
Xie J K. 1998. Spectral inversion using Lg from earthquakes:Improvement of the method with applications to the 1995,western Texas earthquake sequence[J]. Bull Seismol Soc Am,88(6):1525–1537.
Yang X N. 2002. A numerical investigation of Lg geometrical spreading[J]. Bull Seismol Soc Am,92(8):3067–3079. doi: 10.1785/0120020046
Zhao L F,Xie X B,Wang W M,Zhang J H,Yao Z X. 2010. Seismic Lg-wave Q tomography in and around Northeast China[J]. J Geophys Res,115(B8):B08307.
Zhao L F,Xie X B,He J K,Tian X B,Yao Z X. 2013a. Crustal flow pattern beneath the Tibetan Plateau constrained by regional Lg-wave Q tomography[J]. Earth Planet Sci Lett,383:113–122. doi: 10.1016/j.jpgl.2013.09.038
Zhao L F,Xie X B,Wang W M,Zhang J H,Yao Z X. 2013b. Crustal Lg attenuation within the North China Craton and its surrounding regions[J]. Geophys J Int,195(1):513–531. doi: 10.1093/gji/ggt235
Zhao W L,Morgan W J. 1987. Injection of Indian crust into Tibetan lower crust:A two-dimensional finite element model study[J]. Tectonics,6(4):489–504. doi: 10.1029/TC006i004p00489
-
期刊类型引用(3)
1. 李宏男，郭鑫，张皓，田威. 强风和地震同时发生及对结构的作用研究进展. 中国科学:技术科学. 2023(12): 2005-2020 . 
百度学术
2. 赵家庆，陈界宏，林立青. 2016年台风“尼伯特”激发微地动的频率特征分析. 中国地震. 2021(02): 523-531 . 百度学术
3. 杨锦玲，杨小林，李军，韦进. 东南沿海区域重力微震信号特征分析. 地球物理学进展. 2021(06): 2362-2370 . 百度学术
其他类型引用(0)
下载:
计量
- 文章访问数: 1002
- HTML全文浏览量: 402
- PDF下载量: 85
- 被引次数: 3
