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
-
引言
我国大陆地区自1976年唐山MS7.8地震发生后，经历了12年无MS≥7.0地震的平静期。从1988年起进入了MS≥7.0地震的活跃期，到2017年底已发生14次MS≥7.0地震。这些大地震均发生在青藏高原地区地壳构造块体的边界断层上，表明现阶段青藏高原地壳构造块体系统中，在块体边界断层的多个部位上，应变已积累到或接近地下岩石的破裂强度水平。特别是，从1996年至今，我国大陆地区所有MS≥7.0地震，包括2次MS≥8.0地震（2001年昆仑山口西MS8.1地震、2008年汶川MS8.0地震）和7次MS≥7.0地震（1996年喀喇昆仑山MS7.1地震、1997年玛尼MS7.9地震、2008年于田MS7.0地震、2010年玉树MS7.3地震、2013年芦山 MS7.0地震、2014年于田MS7.0地震、2017年九寨沟MS7.0地震），均发生在巴颜喀拉块体的边界断层上。
图1是我国大陆地区构造块体及1988年以来MS≥6.9地震震中分布图，其中，块体的划分引自张培震等（2003），震中分布图引自本研究团队于2009年发表的论文（陈祖安等, 2009），但增补了自2009年以来发生的MS≥7.0地震。由图1可见，青藏高原地区各构造块体的边界断层上大震活动从未有过的活跃，显示了该地区的大震活动与构造块体运动、变形的密切关联。
![]()
图 1 我国大陆地区构造块体及1988年以来MS≥6.9地震震中分布图Figure 1. Distribution of tectonic blocks and MS≥6.9 earthquakes since 19882008年汶川MS8.0地震发生前，在巴颜喀拉地块的西、南、北边界断层上，分别发生了多次MS≥7.0地震，唯独在其东边界断层龙门山断裂带上未发生MS≥7.0地震（陈祖安等，2009）。2001年发生在巴颜喀拉地块北边界断层东昆仑断裂带上的昆仑山口西MS8.1地震，其震源机制是左旋走滑错动，破裂过程成像显示该地震的破裂总体上自西向东扩展（许力生，陈运泰，2004）。近东—西走向的昆仑山口西地震发震断层南侧向东的错动，必然会推动巴颜喀拉地块进一步东移，从而加速龙门山断裂带的应变积累。从大震活动与构造块体运动变形相互影响的角度分析，青藏高原构造块体边界断层的某些地段，特别是龙门山断裂带，的中长期大震活动危险性应予以特别关注。自汶川大地震前两年起，作者以国家自然科学基金项目“玛尼、昆仑山大震发生对青藏高原构造块体稳定性影响的数值模拟研究”为题，开始对青藏高原构造块体系统的稳定性、各构造块体边界断层上失稳危险度的分布，以及1997年玛尼地震和2001年昆仑山口西地震的发生对构造块体系统稳定性的影响等方面的问题进行了研究，获得了一些有意义的结果（陈祖安等，2008）。
2008年汶川大地震发生近5年之后，2013年4月20日在龙门山断裂带西南段上汶川发震断层的西南端外侧发生了芦山地震。由地震波资料反演结果可知，这两次大震均为全球内陆地区罕见的低倾角逆冲型大地震（陈运泰等，2013）。有关这两次大震的孕震机理及震前中长期地震危险性，是值得深入研究的问题。近年来，本文第一作者及其研究团队成员运用三维流变非连续变形法（描述地壳块体边界）与有限元法（描述地壳块体内部）相结合的数值模拟方法，研究了青藏高原及其邻区（东侧四川盆地和鄂尔多斯地块）内共12个构造块体的运动、变形、相互作用及其与近几十年发生于该区的大地震之间的关系（陈祖安等，2008，2009，2011；Chen et al，2011 ；林邦慧等，2012，2014，2016）。本文将在系统总结上述研究工作的基础上，进一步研究汶川地震和芦山地震的孕震机理以及大震前的中长期地震危险性。
2. 方法与模型
印度板块向北挤压造成青藏高原的隆起并产生厚约70 km的地壳。由于地壳内温度上升，下地壳层速度低且高导层发育，介质呈较强黏塑性。在印度板块向北的推挤下，下地壳的黏塑性流动驱动着上地壳脆性层运动与变形，形成了该地区复杂的孕震环境。考虑到上述因素，宜采用描述块体边界相互作用十分有效的非连续变形法（discontinuous deformation analysis，简写为DDA）与描述块体内相互作用的有限元法（finite element method，简写为FEM）相结合的方法研究这一问题。非连续变形分析法是一个分析块体系统中块体相互作用的有效工具（石根华，1997），可以用来计算块体系统在动力及静力作用下块体的滑移和不连续变形。对于各块体，允许有滑移、压性或剪切变形；对于整个块体系统，允许块体边界断层有滑动。迄今，已有一些研究人员使用该方法进行了多项研究，并取得重要成果（Cai et al，2000 ；白武明等，2003），这些结果证明了该方法适合用于模拟大震的破裂过程和大震活动与构造块体相互作用的关系。
陈祖安等（2011）在石根华（1997）非连续变形法和二维情形程序的基础上，编制了三维流变DDA方法与FEM方法相结合的程序“三维DDA＋FEM程序”。该程序的关键之处在于将块体之间的相互作用转变为它们之间的接触问题（Perić，Owen，1992）；同时假设位移很小，块体之间的接触点已知，详细情况请参见陈祖安等（2011），这里从略。本文将这个程序用于汶川、芦山地震孕震机理等研究。
研究区包括青藏高原及其东侧的四川盆地和鄂尔多斯地区（图2a）。参考张培震等（2003）对我国大陆地区构造块体的划分结果，将研究区划分为塔里木、柴达木、阿拉善、祁连、昆仑、羌塘、拉萨、川滇、滇西南、鄂尔多斯、四川盆地及华南等12个地块。地块边界断层包括阿尔金、祁连山、西秦岭、东昆仑、西昆仑、六盘山、贺兰山、玛尼—玉树、喀喇昆仑、喜马拉雅、龙门山、鲜水河、三江、红河、安宁河—小江、澜沧江、秦岭大别山、汾渭和阴山等断裂带。
将上述区域垂直深度达120 km的空间划分为5层，在青藏高原东缘以西地区划分上地壳20 km，中地壳20 km，下地壳30 km及上地幔。受资料限制，为简化起见，除龙门山断裂带以外，其它边界断层均为直立。在青藏高原东缘及以东地区，考虑了龙门山断裂带东西两侧地势、地壳厚度和分层的明显变化（滕吉文等，2008）：地势从四川盆地海拔500 m向西到龙门山地区陡升至4 000 m；地壳厚度在四川盆地为45 km，向西逐渐增厚，在龙门山断裂带东侧地区为60 km，并逐渐向西增厚，其西侧的青藏高原地区地壳厚度为70 km。
我们参考了陈运泰研究团队有关汶川地震破裂过程的反演结果（陈运泰等，2008，2013；张勇等，2008），对本研究团队2009年论文中所使用的汶川地震的发震断层模型（陈祖安等，2009）进行了改进。考虑到龙门山断裂带上汶川大震的发震断层由3个相连接的子事件组成，倾角由西南向东北扩展且逐渐增大，南段倾角在中地壳为40°，在上地壳中为60°，北段的倾角接近直立。参考万永革等（2009）一文中的龙门山断裂带三维模型，改进并构制了龙门山断裂带上汶川发震断层的三维分布模型（图2b）。
![]()
图 2 研究区构造块体模型图（a）和龙门山断裂带上汶川MS8.0地震发震断层的三维示意图（b）（引自陈祖安等，2011）Figure 2. Tectonic block system model used in this study （a） and 3D schematic diagram of the seismogenic faults of the 2008 Wenchuan MS8.0 earthquake on the Longmenshan fault zone （b） （after Chen et al，2011 ）3. 研究区块体边界的失稳危险度、应变率强度及应变能密度分布
This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
$G = \frac{\tau }{{{\sigma _{\rm{n}}}{\rm{tan}}\phi }}, $
(1) This page contains the following errors:
error on line 1 at column 1: Start tag expected, '<' not foundBelow is a rendering of the page up to the first error.
我们模拟计算了研究区现今构造块体边界断层上的失稳危险度分布（林邦慧等，2012，2014）。计算结果表明，上、中地壳构造块体边界断层上的失稳危险度高的地段多数与近30年来发生MS≥7.0地震的区域基本一致，包括1990年共和MS6.9地震、1997年玛尼MS7.9地震、2001年昆仑山口西MS8.1地震、2008年汶川MS8.0地震、2013年芦山MS7.0地震、2014年于田MS7.3地震和2017年九寨沟MS7.0地震的发震断层。此外，龙门山断裂带（包括汶川地震和芦山地震的发震断层）除北西端外均为失稳危险度最高值（图4）。本文参考陈运泰等（2008，2013）和张勇等（2008）关于汶川地震破裂过程的反演结果提出的汶川发震断层的三维模型的模拟结果表明，边界断层特别是龙门山断裂带的失稳危险度分布明显地比2009年用较简单的汶川发震断层模型（陈祖安等，2009）结果更符合近30年MS≥7.0地震的地震活动性。
![]()
图 3 由数值模拟计算得到的研究区地表速度场（a）和上地壳层震源机制（b）分布图（引自陈祖安等，2009）Figure 3. Distributions of the ground surface velocity calculated （a） and the focal mechanism of upper crust （b） in the study area （after Chen et al，2009）![]()
图 4 研究区上地壳（a）和中地壳（b）构造块体边界断层上的失稳危险度分布1. 1990年共和MS6.9地震；2. 1997年玛尼MS7.9地震；3. 2001年昆仑山口西MS8.1地震；4. 2008年汶川MS8.0地震；5. 2013年芦山MS7.0地震；6. 2014年于田MS7.3地震；7. 2017年九寨沟MS7.0地震Figure 4. Distribution of instability risk factor of boundary faults between tectonic blocks of upper crust （a） and lower crust （b） in the studied area1. MS6.9 Gonghe earthquake in 1990；2. MS7.9 Mani earthquake in 1997；3. MS8.1 western Kunlun Mountain Pass earthquake in 2001；4. MS8.0 Wenchuan earthquake in 2008；5. MS7.0 Lushan earthquake in 2013；6. MS7.3 Yutian earthquake in 2014；7. MS7.0 Jiuzhaigou earthquake in 2017为了进一步研究汶川大震的孕震机理，我们还计算了研究区上、中地壳的应变率强度分布，结果如图5所示。计算结果表明，青藏高原绝大部分地区均为高剪切应变率地区，其东缘基本上为应变率强度的急剧变化带，特别是龙门山断裂带东西两侧应变率强度的变化最为急剧，其西侧的应变率强度约是东侧的4倍。我们还计算了研究区应变能密度分布（陈祖安等，2009），结果如图6所示，可以看出：龙门山断裂带及邻近地区为宽度与龙门山断裂带宽度基本相同、走向一致的高应变能密度带，其东西两侧则是应变能密度较低的地区，且断裂带东侧应变能密度较西侧衰减得快。
![]()
图 5 研究区上地壳（上）和中地壳（下）应变率强度分布图（左）及其等值线图（右）（引自陈祖安等，2009）Figure 5. Distribution of strain-rate intensity （left） and contours of strain-rate （right） of the upper crust layer （upper） and the middle crust layer （lower） in the studied area （after Chen et al，2009）![]()
图 6 研究区上地壳（上）、中地壳（下）应变能密度分布图（左）及其等值线图（右）（引自陈祖安等，2009）Figure 6. Distribution of strain energy-density （left） and contours of strain energy-density （right） of the upper crust layer （upper） and the middle crust layer （lower） in the studied area （afterChen et al，2009 ）4. 2001年昆仑山口西地震和2008年汶川地震的影响
4.1 2001年昆仑山口西MS8.1地震的影响
为了探讨2001年昆仑山口西MS8.1地震对2008年汶川MS8.0地震产生的影响，我们通过数值计算模拟了昆仑山口西地震的破裂过程，研究了该地震对研究区各块体边界断层库仑破裂应力分布的影响（陈祖安等，2011）。图7给出了昆仑山口西地震引起的研究区各边界断层上库仑破裂应力变化的分布，可见，该地震的发生除了引起其发震断层两端部库仑破裂应力增加外，对阿尔金断裂带、玛尼—玉树断裂带中段和西南地段，以及从龙门山断裂带的北端到南端，包括汶川地震和芦山地震的发震断层，库仑破裂应力均有较大的增加，增幅约为0.016 MPa。汶川地震和芦山地震的发震断层是同时具有高失稳危险度（林邦慧等，2012）且大震引起库仑破裂应力增加这两个特点的地带。计算得出，这两条断层上约0.016 MPa的库仑破裂应力增加相当于龙门山断裂带约两年的应力积累，也就是说，昆仑山口西地震的发生使汶川、芦山地震的发震断层失稳破裂提前了大约两年。
![]()
图 7 昆仑山口西MS8.1地震对青藏高原及邻近地区上地壳（a）和中地壳（b）各块体边界断层上库仑破裂应力的影响（引自陈祖安等，2011）Figure 7. The influence of the 2001 western Kunlun Mountain Pass MS8.1 earthquake on the Coulomb failure stress of the boundary faults of the tectonic blocks of upper crust （a） and middle crust （b） in the Qinghai-Tibet Plateau and its vicinity （after Chen et al，2011）4.2 汶川地震的影响及2013年芦山MS7.0地震的成功预测
汶川地震的发生对青藏高原及邻近地区各边界断层地震危险性影响是本文特别予以关注的问题。汶川地震发生后，我们根据上述汶川大震破裂过程由3个倾角不同的子事件组成的研究结果，提出了汶川地震的发震断层模型（图2b），数值模拟了汶川地震的破裂过程，研究了大震引起各边界断层应力状态变化的特征（陈祖安等，2011），并在其文中的图6加上2013年MS7.0芦山地震和2017年MS7.0九寨沟地震的位置，如图8所示。结果显示，在上地壳，汶川发震断层南端外侧到龙门山断裂带南端、汶川发震断层北端外侧、鲜水河断裂带南段、东昆仑断裂带东南段等地段的库仑破裂应力增加幅度最大，达到约0.2 MPa，而2013年芦山地震的发震断层正处于汶川发震断层南端外侧，即龙门山断裂带西南段的中间部位，龙门山断裂带西南段是满足失稳危险度大（林邦慧等，2012）且经历昆仑山口大震和汶川大震两次库仑破裂应力增幅较大这两个条件的地区，是大震发生的危险区，林邦慧等（2012）据此对这次地震成功地进行了中长期地震预测。
![]()
图 8 汶川地震的发生对青藏高原及邻近地区上地壳（a）和中地壳（b）各块体边界断层库仑破裂应力的影响（修改自陈祖安等，2011）Figure 8. The influence of the 2008 Wenchuan MS8.0 earthquake on the Coulomb failure stress of the boundary faults of the tectonic blocks of upper crust （a） and middle crust （b） in the Qinghai-Tibet Plateau and its vicinity （modified from Chen et al，2011 ）5. 讨论与结论
本文概要叙述了作者运用三维流变非连续变形分析法与有限元法相结合的方法研究青藏高原及其东侧邻区构造块体的运动、变形、相互作用及其与近30年来发生于该区的大地震活动关系以及汶川、芦山大震的孕震机理，得出：
1） 失稳危险度高的地段与近30年来发生的MS≥7.0地震所在区域基本一致，其中龙门山断裂带是力学上失稳危险度最高的地段，也是应变率强度的急剧变化带，西侧应变率强度约是东侧应变率强度的4倍，而且断裂带东侧的应变率强度等值线衰减得比西侧快。
2） 龙门山断裂带在上、中地壳中均位于与龙门山断裂带宽度相同、走向一致的高应变能密度带中，在上地壳的这个带的东西两侧则是应变能密度较低的地区，而在中地壳其强度在断裂带东侧逐渐向东衰减；西侧应变能密度高，而东侧应变能密度较低。这表明，汶川、芦山大震前，其发震断层已积累了大量的应变能，失稳危险度高，处于力学上的不稳定状态。
3） 2001年昆仑山口西MS8.1地震的发生使得低倾角的汶川发震断层和包含芦山发震断层在内的龙门山断裂带西南段的库仑破裂应力增大约0.016 MPa，相当于龙门山断裂带约两年的应力积累，也就是说，昆仑山口西MS8.1地震的发生使汶川、芦山地震发震断层的失稳破裂大约提前了两年。
4） 2008年汶川MS8.0地震的发生也使芦山大震的发震断层和龙门山断裂带西南段的库仑破裂应力在昆仑山口西大震时增大，之后再次增大，应变能积累增强，对已处于失稳危险度较高状态的芦山大震发震断层的提前失稳起到了促进作用，据此成功地作了中长期地震预测（陈祖安等，2011；林邦慧等，2012）。陈运泰等（2013）的研究结果表明，龙门山断裂带西南段上，位于芦山大震发震断层南端至龙门山断裂带南端之间，以及芦山发震断层北端至汶川发震断层南端之间存在两个有发生大震危险的地震矩释放“亏空区”，长度分别约为70 km和30 km。在上述研究结果中，这两个区正好是满足失稳危险度大和库仑破裂应力多次增大这两个条件的地区。从青藏高原地区大震活动与构造块体运动变形的相互影响角度来看，汶川、芦山大震后龙门山断裂带上这些相关地段仍存在明显的中长期地震危险性。
5） 2008年汶川地震的发生引起东昆仑断裂带东南段库仑破裂应力的增幅也最大，表明汶川地震对于2017年九寨沟大震发震断层的提前失稳起到了促进作用（陈祖安等，2011；林邦慧等， 2012，2016）。
作者自2001年昆仑山口西MS8.1地震发生以来，便开始关注巴颜喀拉块体边界断层上发生的大震活动。2006年启动的国家自然科学基金项目“玛尼、昆仑山大震发生对青藏高原构造块体稳定性影响的数值模拟研究”，使得我们对各构造块体边界断层上失稳危险度分布以及玛尼、昆仑山口西地震的发生对构造块体系统稳定性的影响等问题获得了一些新的认识。从2008年汶川地震至今近十年来，我们对发生在巴颜喀拉块体边界断层上的一系列地震（包括1990年共和MS6.9地震、1997年玛尼MS7.9地震、2001年昆仑山口西MS8.1地震、2008年汶川MS8.0地震、2013年芦山MS7.0地震、2014年于田MS7.3地震以及2017年九寨沟MS7.0地震）的研究结果表明了这些大地震的发生与高失稳危险度、高应变能密度分布以及高应变率强度分布的急剧变化密切关联，这些新的认识对于中长期地震危险性的研究和地震预测预报研究是有益的参考。
作者对陈运泰院士以及两位匿名审稿专家对本文的修改建议表示衷心的感谢。
-
![]()
图 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
-
期刊类型引用(1)
1. 侯争，郭增长，王小旗，杜久升. 川滇地壳形变信息的GNSS时空探测及分析. 测绘科学. 2022(06): 30-37+64 . 
百度学术
其他类型引用(4)
下载:
计量
- 文章访问数: 1022
- HTML全文浏览量: 413
- PDF下载量: 86
- 被引次数: 5
