Crustal structure and seismic activity in the Hanzhong basin and its adjacent areas
-
摘要: 为了进一步探查汉中盆地的深部动力学机制和孕震构造特征，本文基于背景噪声成像、多频接收函数和面波联合反演以及莫霍面Ps震相时深转换方法反演了汉中盆地及其邻区的地壳S波速度和厚度，并进一步对比分析了研究区深部结构与地震活动性之间的关系。结果表明：汉中盆地不同区域的浅表沉积厚度和速度存在差异；部分区域莫霍面处的速度变化平缓，Ps震相与P震相的振幅比＜0.2；汉中盆地内部鲜有地震发生，其周边10 km范围内地震分布主要受到断层控制；4—16 km震源深度上下界面大致对应于低速体底层和高速体顶层。本文获得的非均匀分布的沉积厚度、渐变的壳幔过渡带结构与汉中盆地长期处于秦岭构造带、大巴山褶皱带以及青藏地块交界区的三联点构造位置密切相关。Abstract: In order to explore the deep dynamic mechanism and seismogenic structural characteristics in the Hanzhong basin, this paper used the ambient noise tomography, joint inversion of the multi-frequency receiver function and surface wave, and the time-depth conversion method of Ps phase from the Moho to image the crustal S-wave velocity and thickness in the Hanzhong basin and its adjacent areas. Then, we further analyzed the relationship between the crustal structure and seismic activity in the studied region. The results show that the sedimentary thickness and velocity of shallow surface are different in different areas of the Hanzhong basin; in some areas of the studied region Moho surface varies gently with a Ps/P ratio less than 0.2. There are few earthquakes in the Hanzhong basin, and the earthquake distribution within 10 km around the basin is mainly controlled by faults. The upper and lower interfaces of the focal depth （4−16 km） correspond to the bottom of the low-velocity region and the top of the high-velocity region, respectively. The uneven sedimentary thickness and gradual crust-mantle transitional zone obtained in this paper are closely related to the triple junction tectonic position of the Hanzhong basin among the Qinling tectonic belt, the Dabashan folded belt and the Tibetan Plateau. The crustal structure and seismicity characteristics obtained in this paper can provide useful information for the study on the deep dynamic mechanism and seismogenic environment in the Hanzhong basin.
-
Keywords:
- Hanzhong basin /
- crustal structure /
- seismic activity /
- ambient noise tomography /
- receiver function
-
引言
我国大陆地区自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 汉中盆地及其邻区地质构造背景（邓起东等，2002） （a）以及2009年以来的地震和台站分布（b）
图中红色圆圈为4次M5—5½ 历史地震。F1：勉略断裂；F2：汉中北缘断裂；F3：青川断裂；F4：茶坝—林庵寺断裂；F5：梁山南缘断裂；F6：汉中南缘断裂，下同
Figure 1. Regional geological tectonic settings （Deng et al，2003） （a） and distribution of earthquakes since 2009 and stations （b） in the Hanzhong basin and its adjacent regions
The red circles represent the four historical earthquakes with M5−5½ . F1：Mianlüe fault；F2：Hanzhong north edge fault；F3：Qingchuan fault；F4：Chaba-Lin’ansi fault；F5：Liangshan south edge fault；F6：Hanzhong south edge fault；the same below
![]()
图 2 典型台站对垂直分量波形之间的互相关函数（a），邻域算法反演结果（b），群速度频散曲线（c，黄色线条所示）及其拟合结果（d）
Figure 2. Vertical-component cross-correlation functions （a），neighborhood algorithm inversion results （b），group velocity dispersion measurement denoted by the yellow curve （c） and the fitting of group velocity dispersion curves （d） for typical station-pairs
![]()
图 3 汉中盆地及其邻区浅层频散曲线路径覆盖（a）和不同深度h的S波速度分布（b−f）
Figure 3. Path coverage of dispersion curves （a） and distribution of S-wave velocity at different depths h （b−f） in the Hanzhong basin and its surrounding areas
（a） h＝0.3 km；（b） h＝0.7 km；（c） h＝1.1 km；（d） h＝1.7 km；（e） h＝2.5 km；（f） h＝3.9 km
![]()
图 4 多频接收函数及面波频散联合反演示例图
图（a），（b）和（c）分别为典型台站NSHT反演所得S波速度以及接收函数和频散曲线的拟合情况，其中图（b）中F为拟合系数，G为高斯系数，D为平均震中距；图（d）为研究区6个宽频带地震台的联合反演结果
Figure 4. The schematic diagrams for the joint inversion of multi-frequency receiver function and surface wave dispersion
Figs. （a），（b） and （c） show the inverted S-wave velocity，fitting for the receiver functions and dispersion curves for a typical station NSHT. In Fig. （b），F，G and D are the fitting coefficient，Gaussian factor and average epicentral distance，respec-tively. Fig. （d） shows the inversion results for six stations equipped with broadband seismometers in the studied region
![]()
图 5 16个台站远震接收函数叠加后的前8 s波形（a）和两个典型台站S22 （b）和LUYA （c）叠加前的接收函数波形
Figure 5. The stacked teleseismic receiver functions from −1 s to 8 s for 16 stations （a）and original receiver functions for two typical stations S22 （b） 和LUYA （c）
![]()
图 6 汉中盆地及其邻区地形（a）和地壳厚度（b）
图（b）中黄色三角形旁边的数字标明该处的地壳厚度值，单位为km
Figure 6. Topography （a） and crustal thickness （b） in Hanzhong basin and its adjacent areas
In Fig.（b），the numbers next to the yellow triangles （stations） show the crustal thickness values with unit of km beneath the stations
![]()
图 7 2009年1月1日至2018年9月1日汉中盆地及其周边10 km范围内的地震分布（a）及最大震级统计分析（b）
Figure 7. Distribution of the earthquakes in the Hanzhong basin and its surrounding areas （＜10 km） from 1 January 2009 to 1 September 2018 （a） and statistic on corresponding maximum magnitude （b）
-
邓起东,张培震,冉勇康,杨晓平,闵伟,楚全芝. 2002. 中国活动构造基本特征[J]. 中国科学:D辑,32(12):1020–1030. Deng Q D,Zhang P Z,Ran Y K,Yang X P,Min W,Chu Q Z. 2003. Basic characteristics of active tectonics of China[J]. Science in China:Series D,46(4):356–372.
甘卫军,沈正康,张培震,任金卫,万永革,周德敏. 2004. 青藏高原地壳水平差异运动的GPS观测研究[J]. 大地测量与地球动力学,24(1):29–35. Gan W J,Shen Z K,Zhang P Z,Ren J W,Wan Y G,Zhou D M. 2004. Horizontal crustal movement of Tibetan Plateau from GPS measurements[J]. Journal of Geodesy and Geodynamics,24(1):29–35 (in Chinese).
国家测震台网数据备份中心. 2007. 国家测震台网地震波形数据[DB/OL]. 中国地震局地球物理研究所. http://www.seisdmc.ac.cn. doi: 10.11998/SeisDmc/SN. Data Management Centre of China National Seismic Network. 2007. Waveform data of China National Seismic Network[DB/OL]. Institute of Geophysics, China Earthquake Administration. http://www.seisdmc.ac.cn. doi:10.11998/SeisDmc/SN （in Chinese）.
国家地震局震害防御司. 1995. 中国历史强震目录（公元前23世纪—公元1911年）[M]. 北京: 地震出版社: 1−514. Department of Earthquake Disaster Prevention, State Seismological Bureau. 1995. Catalog of Chinese History Strong Earthquakes （From 23rd Century BC to 1911 AD）[M]. Beijing: Seismological Press: 1−514 （in Chinese）.
胡家富,朱雄关,夏静瑜,陈赟. 2005. 利用面波和接收函数联合反演滇西地区壳幔速度结构[J]. 地球物理学报,48(5):1069–1076. doi: 10.3321/j.issn:0001-5733.2005.05.013 Hu J F,Zhu X G,Xia J Y,Chen Y. 2005. Using surface wave and receiver function to jointly inverse the crust-mantle velocity structure in the west Yunnan area[J]. Chinese Journal of Geophysics,48(5):1069–1076 (in Chinese). doi: 10.1002/cjg2.750
李晓妮,冯希杰,任隽,李高阳,李苗,王夫运,姬计法. 2013. 陕南汉中盆地西部梁山南缘断裂隐伏段的活动性鉴定[J]. 地震学报,35(4):534–542. doi: 10.3969/j.issn.0253-3782.2013.04.008 Li X N,Feng X J,Ren J,Li G Y,Li M,Wang F Y,Ji J F. 2013. Activity identification of the buried segment of Liangshan south margin fault in the west of Hanzhong basin[J]. Acta Seismologica Sinica,35(4):534–542 (in Chinese).
刘启元,李昱,陈九辉,van der Hilst R D,郭飚,王峻,齐少华,李顺成. 2010. 基于贝叶斯理论的接收函数与环境噪声联合反演[J]. 地球物理学报,53(11):2603–2612. Liu Q Y,Li Y,Chen J H,van der Hilst R D,Guo B,Wang J,Qi S H,Li S C. 2010. Joint inversion of receiver function and ambient noise based on Bayesian theory[J]. Chinese Journal of Geophysics,53(11):2603–2612 (in Chinese).
唐明帅,葛粲,郑勇,王海涛. 2013. 短周期地震仪接收函数的可行性分析:以新疆和田地震台阵为例[J]. 地球物理学报,56(8):2670–2680. doi: 10.6038/cjg20130816 Tang M S,Ge C,Zheng Y,Wang H T. 2013. Feasibility analysis of short-period seismograph receiver function:An example of Hotan Seismic Array,Xinjiang[J]. Chinese Journal of Geophysics,56(8):2670–2680 (in Chinese).
王明明. 2013. 汉中盆地发育机制及构造演化研究[D]. 北京: 中国地震局地质研究所: 4−8. Wang M M. 2013. A Study on Developmental Mechanism and Tectonic Evolution of the Hanzhong Basin[D]. Beijing: Institute of Geology, China Earthquake Administrator: 4−8 （in Chinese）.
徐杰,计凤桔,周本刚. 2012. 有关我国新构造运动起始时间的探讨[J]. 地学前缘,19(5):284–292. Xu J,Ji F J,Zhou B G. 2012. On the lower chronological boundary of the Neotectonic period in China[J]. Earth Science Frontiers,19(5):284–292 (in Chinese).
张国伟,程顺有,郭安林,董云鹏,赖绍聪,姚安平. 2004. 秦岭—大别中央造山系南缘勉略古缝合带的再认识:兼论中国大陆主体的拼合[J]. 地质通报,23(9/10):846–853. Zhang G W,Cheng S Y,Guo A L,Dong Y P,Lai S C,Yao A P. 2004. Mianlüe paleo-suture on the southern margin of the Central Orogenic System in Qinling-Dabie:With a discussion of the assembly of the main part of the continent of China[J]. Geological Bulletin of China,23(9/10):846–853 (in Chinese).
郑秀芬,欧阳飚,张东宁,姚志祥,梁建宏,郑洁. 2009. " 国家数字测震台网数据备份中心”技术系统建设及其对汶川大地震研究的数据支撑[J]. 地球物理学报,52(5):1412–1417. doi: 10.3969/j.issn.0001-5733.2009.05.031 Zheng X F,Ouyang B,Zhang D N,Yao Z X,Liang J H,Zheng J. 2009. Technical system construction of Data Backup Centre for China Seismograph Network and the data support to researches on the Wenchuan earthquake[J]. Chinese Journal of Geophysics,52(5):1412–1417 (in Chinese). doi: 10.3969/j.issn.0001-5733.2009.05.031
Aki K, Richards P G. 1980. Quantitative Seismology: Theory and Methods[M]. San Francisco: W. H. Freeman and Company: 133−155.
Bao X W,Song X D,Xu M J,Wang L S,Sun X X,Mi N,Yu D Y,Li H. 2013. Crust and upper mantle structure of the North China Craton and the NE Tibetan Plateau and its tectonic implications[J]. Earth Planet Sci Lett,369/370:129–137. doi: 10.1016/j.jpgl.2013.03.015
Bensen G D,Ritzwoller M H,Barmin M P,Levshin A L,Lin F,Moschetti M P,Shapiro N M,Yang Y. 2007. Processing seismic ambient noise data to obtain reliable broad-band surface wave dispersion measurements[J]. Geophys J Int,169(3):1239–1260. doi: 10.1111/gji.2007.169.issue-3
Campillo M, Roux P. 2014. Crust and lithospheric structure: Seismic imaging and monitoring with ambient noise corre-lations[G]//Treatise on Geophysics （2nd ed）. Amsterdam: Elsevier: 391–417.
Chang S J,Baag C E. 2005. Crustal structure in southern Korea from joint analysis of teleseismic receiver functions and surface-wave dispersion[J]. Bull Seismol Soc Am,95(4):1516–1534. doi: 10.1785/0120040080
Chen Y L,Niu F L,Liu R F,Huang Z B,Tkalčić H,Sun L,Chan W. 2010. Crustal structure beneath China from receiver function analysis[J]. J Geophys Res,115(B3):B03307. doi: 10.1029/2009JB006386
Ditmar P G, Yanovskaya T B. 1987. A generalization of the Backus-Gilbert method for estimation of lateral variations of surface wave velocity[J]. Izv Akad Nauk SSSIt, Fiz Zemli, 6: 30−60 （in Russian）.
He R Z,Shang X F,Yu C Q,Zhang H J,van der Hilst R D. 2014. A unified map of Moho depth and vP/vS ratio of continental China by receiver function analysis[J]. Geophys J Int,199(3):1910–1918. doi: 10.1093/gji/ggu365
Herrmann R B. 2013. Computer programs in seismology:An evolving tool for instruction and research[J]. Seismol Res Lett,84(6):1081–1088. doi: 10.1785/0220110096
Julià J,Ammon C J,Herrmann R B,Correig A M. 2000. Joint inversion of receiver function and surface wave dispersion observations[J]. Geophys J Int,143(1):99–112. doi: 10.1046/j.1365-246x.2000.00217.x
Laske G, Masters G, Ma Z T, Pasyanos M. 2013. Update on CRUST1.0: A 1-degree global model of Earth’s crust[C]//EGU General Assembly Conference Abstract, Vol. 15. Vienna, Austria: EGU: 2658.
Li Z W,Ni S D,Zhang B L,Bao F,Zhang S Q,Deng Y,Yuen D A. 2016. Shallow magma chamber under the Wudalianchi Volcanic Field unveiled by seismic imaging with dense array[J]. Geophys Res Lett,43(10):4954–4961. doi: 10.1002/grl.v43.10
Lin A M,Rao G,Yan B. 2014. Structural analysis of the right-lateral strike-slip Qingchuan fault,northeastern segment of the Longmen Shan thrust belt,central China[J]. J Struct Geol,68:227–244. doi: 10.1016/j.jsg.2014.09.014
Ling Y,Chen L,Wei Z G,Jiang M M,Wang X. 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[J]. Science China Earth Science,60(10):1752–1768. doi: 10.1007/s11430-017-9092-8
Liu Z,Tian X B,Gao R,Wang G C,Wu Z B,Zhou B B,Tan P,Nie S T,Yu G P,Zhu G H,Xu X. 2017. New images of the crustal structure beneath eastern Tibet from a high-density seismic array[J]. Earth Planet Sci Lett,480:33–41. doi: 10.1016/j.jpgl.2017.09.048
Okay A I,Kaşlılar-Özcan A,Imren C,Boztepe-Güney A,Demirbağ E,Kuşçu İ. 2000. Active faults and evolving strike-slip basins in the Marmara Sea,northwest Turkey:A multichannel seismic reflection study[J]. Tectonophysics,321(2):189–218. doi: 10.1016/S0040-1951(00)00046-9
Sambridge M. 1999. Geophysical inversion with a neighbourhood algorithm:I. Searching a parameter space[J]. Geophys J Int,138(2):479–494. doi: 10.1046/j.1365-246X.1999.00876.x
Shen W S,Ritzwoller M H,Kang D,Kim Y H,Lin F C,Ning J Y,Wang W T,Zheng Y,Zhou L Q. 2016. A seismic reference model for the crust and uppermost mantle beneath China from surface wave dispersion[J]. Geophys J Int,206(2):954–979. doi: 10.1093/gji/ggw175
Wang C Y,Sandvol E,Zhu L P,Lou H,Yao Z X,Luo X H. 2014. Lateral variation of crustal structure in the Ordos block and surrounding regions,North China,and its tectonic implications[J]. Earth Planet Sci Lett,387:198–211. doi: 10.1016/j.jpgl.2013.11.033
Wang W L,Wu J P,Fang L H,Lai G J,Cai Y. 2017. Sedimentary and crustal thicknesses and Poisson’s ratios for the NE Tibetan Plateau and its adjacent regions based on dense seismic arrays[J]. Earth Planet Sci Lett,462:76–85. doi: 10.1016/j.jpgl.2016.12.040
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[J]. Geophys J Int,186(2):385–389. doi: 10.1111/j.1365-246X.2011.05089.x
Wei Z G,Chen L,Wang B Y. 2013. Regional variations in crustal thickness and vP/vS ratio beneath the central-western North China Craton and adjacent regions[J]. Geol J,48(5):531–542. doi: 10.1002/gj.v48.5
Wei Z G,Chu R S,Chen L. 2015. Regional differences in crustal structure of the North China Craton from receiver functions[J]. Science China Earth Sciences,58(12):2200–2210. doi: 10.1007/s11430-015-5162-y
Wei Z G,Chen L,Li Z W,Ling Y,Li J. 2016. Regional variation in Moho depth and Poisson’s ratio beneath eastern China and its tectonic implications[J]. J Asian Earth Sci,115:308–320. doi: 10.1016/j.jseaes.2015.10.010
Xu X W,Wen X Z,Yu G H,Chen G H,Klinger Y,Hubbard J,Shaw J. 2009. Coseismic reverse- and oblique-slip surface faulting generated by the 2008 MW7.9 Wenchuan earthquake,China[J]. Geology,37(6):515–518. doi: 10.1130/G25462A.1
Yao H J,van der Hilst R D,de Hoop M V. 2006. Surface-wave array tomography in SE Tibet from ambient seismic noise and two-station analysis:I. Phase velocity maps[J]. Geophys J Int,166(2):732–744. doi: 10.1111/gji.2006.166.issue-2
Zhang P Z,Wen X Z,Shen Z K,Chen J H. 2010. Oblique,high-angle,listric-reverse faulting and associated development of strain:The Wenchuan earthquake of May 12,2008,Sichuan,China[J]. Ann Rev Earth Planet Sci,38:353–382. doi: 10.1146/annurev-earth-040809-152602
-
期刊类型引用(1)
1. 侯争，郭增长，王小旗，杜久升. 川滇地壳形变信息的GNSS时空探测及分析. 测绘科学. 2022(06): 30-37+64 . 
百度学术
其他类型引用(4)
下载:
计量
- 文章访问数: 1454
- HTML全文浏览量: 721
- PDF下载量: 102
- 被引次数: 5
