2. 中国海洋大学海洋地球科学学院,青岛市松岭路238号,266100;
3. 中国地震局地质研究所,北京市华严里甲1号,100029
日本-千岛海沟转折处俯冲带位于西北太平洋,包括北日本海沟俯冲带和南千岛海沟俯冲带,由“沟-弧-盆”体系构成(图 1)。东北北海道火山的火山爆发指数(voclanic explosivity index,VEI)≤3,而南北海道火山曾在1640~1739年发生4次VEI=5的喷发,是日本火山喷发最强烈的地区[1]。该地区地震活动频繁,北海道之下的震源深度可达300 km [2]。
晚元古代(750~600 Ma BP)扬子克拉通发生大陆裂陷作用导致古日本诞生,之后西南日本因太平洋超级地幔柱的裂陷作用而逐渐形成被动大陆边缘[3-5]。古生代(540~250 Ma BP)扬子克拉通从被动大陆边缘转变为主动大陆边缘(500 Ma BP)。中生代(250~65 Ma BP)中-朝板块和扬子板块碰撞(250~210 Ma BP),鄂霍次克板块形成(90 Ma BP),古千岛之下开始发生俯冲作用而逐渐形成火山弧(90 Ma BP)。新生代(65~0 Ma BP)古伊豆-博宁弧形成(40 Ma BP),千岛、日本等边缘海盆地发育(30~15 Ma BP),日本海的扩张使得日本分离成为岛弧(28~15 Ma BP)[6]。
前人应用地震层析成像技术对该区域进行了一些研究。刘鑫[2]应用台网外层析成像法反演得到日本海沟俯冲带三维地震波衰减结构、南千岛海沟俯冲带三维地震波速度结构及P波各向异性参数。Miler等[7]通过反演北海道转角处的三维P波、S波速度结构,认为北海道岛之下太平洋俯冲板块是连续且弯曲的。但目前缺乏应用地震层析成像对日本-千岛海沟转折处俯冲带分带及分析断裂特征的研究。本文利用P波走时数据对日本-千岛海沟转折处俯冲带进行震源重定位和速度结构反演,结合前人分带和断裂特征资料,对研究区俯冲带进行分带,并划分断裂级别,为全面、深入理解西北太平洋俯冲带结构和构造特征提供地震层析成像依据和研究成果。
1 数据和方法 1.1 数据P波走时目录数据来自国际地震中心(International Seismological Centre,ISC),它经过了EHB方法校正[17],数据详情见表 1、图 2。需要说明的是,设置反演网格的时候设置了一个包含整个研究区在内的边界,因此反演过程中图 2所示网格范围外的地震事件以及台站也有使用。
初始模型的建立参考了IASP91模型(源于地震合作研究所,Incorporated Research Institutions for Seismology,IRIS)。该模型认为上地壳厚度为20 km,P波速度为5.8 km/s;下地壳厚度为15 km,P波速度为6.5 km/s。因此,康拉德面深度为20 km,莫霍面深度则为35 km。在俯冲地区引入太平洋俯冲板块岩石圈的上下界面,并设置太平洋俯冲板块岩石圈内的初始P波速度比相同深度处的壳幔高4%[18-19]。引入高速异常体之后能有效减小走时残差均方根。
1.3 反演方法和分辨率分析使用Zhang[20]基于双差定位方法研发的层析成像反演软件TomoDD,其能够同时解决定位和速度结构联合反演问题。正演计算采用伪弯曲算法[21],反演计算采用最小二乘正交分解法[22]。
根据研究区地形地貌、海沟形态、震源和台站位置等特点,设置反演网格(表 2)。
分辨率分析采用射线密度法(图 2)和偏微分加权数法(derivative weight sum,DWS)[23-24],射线密度越大,分辨率和可信度越高,一般情况下DWS≥100时分辨率好[25]。这两种方法结果均表明,反演得到的速度结果分辨率和可信度较高,且海沟以西、80 km深度范围内的反演结果质量最佳。
2 结果及其初步解释经多次震源重定位及速度反演迭代,研究区的P波走时残差均方根大幅下降(表 3),将第16次反演迭代结果作为最终结果。据美国地质勘探局(United States Geological Survey,USGS)提供的0~700 km太平洋板块俯冲深度等深图、反演得到的震源重定位及P波速度结构结果、前人研究成果[2]等,确定太平洋俯冲板块上下界面,最终得到图 3,并进行如下讨论。
1) 俯冲倾角。由于本文测线近垂直于海沟走向,太平洋板块俯冲方向也近垂直于海沟走向,又由于图 3中横纵向比例尺为1 :1,所以根据图 3中太平洋板块的俯冲形态,可以测出太平洋板块的视俯冲倾角,其变化范围为18°~34°。太平洋板块初始俯冲倾角较平缓,至火山弧之下变陡。由南至北不同测线的平均视俯冲倾角依次为:18°(t1剖面)、21°(t2剖面)、19°(t3剖面)、19°(t4剖面)、19°(t5剖面)、20°(t6剖面)、21°(t7剖面)、20°(t8剖面)、28°(t9剖面)、30°(t10剖面)、24°(t11剖面)、30°(t12剖面)、31°(t13剖面)、31°(t14剖面)、28°(t15剖面)、29°(t16剖面)、31°(t17剖面)、33°(t18剖面)、31°(t19剖面)、33°(t20剖面)、34°(t21剖面)。整体呈现小(t1~t8剖面)-大(t9~t21剖面)的变化趋势,分别对应北日本海沟俯冲带、南千岛海沟俯冲带。
2) 地震分布。某些测线剖面中(如t11剖面)可观察到双地震带。
3) 高速异常体。高速异常体成片、较连续地分布在俯冲板块之内,能够较好地反映太平洋板块的俯冲形态和位置。t22~t36剖面中,高速异常体对应的太平洋俯冲板块在剖面西部-中部弯曲上拱,而在剖面东部弯曲下凹,这与海沟转折处复杂的应力环境有关。
4) 低速异常体。t22~t37剖面中,上覆板块地壳中的低速异常体水平展布,连续性好。由低速异常体形态及位置推测,研究区地壳发生局部熔融,如t5剖面。t7、t8剖面中,在约50 km深度处,低速异常体碎片落在了太平洋俯冲板块上界面,通过观察其上方低速异常体的形态,推断该低速异常体对应的地质体很可能与上方大片的地壳曾经是一个整体,后由于某种原因向下脱落[26]。
5) 与前人地震层析成像剖面对比分析。如图 4,本文的t7、t28测线剖面的俯冲倾角、地震分布、高速异常体和低速异常体的位置及形态、太平洋俯冲板块上下界面位置及形态与前人研究结果(c1、c2测线剖面)基本吻合,表明本文反演结果可信。与前人相比,本文反演结果的细节更多(如能看出局部地壳脱落到太平洋俯冲板块上界面),因此能够更好地反映研究区的结构和构造特征。
6) 断裂几何学特征。地形的变化能指示断层、火山等地质活动。双差层析成像能对断层、火山等构造单元附近的地震进行精确定位并进行速度反演,反过来,反演得到的震源位置及速度结构能够指示断层、火山的特征。一般地,在海拔骤降、地震发生、P波速度变化的地方发育断裂。如图 5,由前人断裂资料可确定各测线剖面上各断裂的位置、倾向、性质,在此基础上,根据地形、震源重定位及速度反演结果,对研究区内各剖面上断裂的几何学特征进行精细刻画(重点确定各断裂的倾角、切割深度),最终得到典型测线的断裂剖面图(图 6)。
P波速度为7 km/s的等值线能一定程度上指示莫霍面[20]。据图 6,地形与Vp=7 km/s的等值线深度大致呈“镜面对称”的关系,具体表现在:在火山弧(日本北海道岛弧和千岛岛弧)及弧前(如增生楔),海拔较高,Vp=7 km/s的等值线深度较大;在海沟(日本海沟和千岛海沟)附近,水深达全区最大,Vp=7 km/s的等值线深度减小。
研究区内发育具有右旋走滑性质的Tartary-Ishikari断裂(TIF)、中央萨哈林-北海道断裂(CHSF)、Abashiri断裂(AF)、北海道断裂(R3~R8)以及逆断层(如R1~R2、R9~R10)。断裂呈NNW或N走向,倾向NEE或E,中等-较大倾角,纵向延伸距离变化大,发育深度可至莫霍面或太平洋俯冲板块上边界,平面形态为弯曲形,剖面形态为铲状,断裂性质为逆断型,大型断裂还具有右旋走滑的性质(如TIF、CHSF、AF)。
该地区最引人注目的是TIF、CHSF、AF及R3~R8断裂组合成的断层组合样式,它们呈阶梯状分布,均向东逆冲,从西向东倾角增大、上盘上升,在地壳的底部可能发育滑脱层,与上覆断裂相连,特殊的断裂组合与研究区内特殊的应力挤压状态有关:既受到东南部太平洋板块NW向俯冲汇聚作用,又受到东北部千岛弧前碎片带SW向挤压碰撞作用。
3 讨论 3.1 俯冲带分带前人对海沟转折处俯冲带分带及不同带之间是否存在撕裂等问题早有探讨[27-29],但就分带数量、边界线位置、撕裂等问题上存在诸多争议。
本次研究结果表明,一方面,北日本海沟俯冲带及南千岛海沟俯冲带在板块俯冲倾角、地震分布、高速异常体、低速异常体等方面存在差异,据此对研究区进行分带(图 7);另一方面,200 km深度范围内日本-千岛海沟转折处俯冲带之下不存在撕裂,为韧性变形的边界,能够看到明显的挤压变形导致的太平洋俯冲板块弯曲上拱和下凹现象(图 3)。
对北日本海沟俯冲带和南千岛海沟俯冲带的特征进行整理归纳,结果见表 4。
日本海沟和千岛海沟控制了西北太平洋边缘俯冲带生长演化等重要过程,而断裂活动是海沟俯冲带内强烈应力作用的主要表现形式之一,研究区内长期而复杂的构造演化导致不同分布范围、规模、形态、性质断裂的发育。对断裂级别进行划分有助于更好地理解研究区的构造特征。
断裂分级主要考虑两个原则:断裂的空间范围(平面延伸长度及切割深度)及其贯穿构造单元的范围(即具体贯穿的构造单元个数);断裂对构造单元发育演化的贡献大小。综合考虑分级原则,对日本-千岛海沟转折处俯冲带进行断裂分级(图 8)。
一级断裂包括:Tartary断裂(TF)、Tartary-Tanakura断裂(TTF)、Tartary-Ishikari断裂(TIF)、东日本海断裂(EJSF)、中央萨哈林-北海道断裂(CHSF)、Abashiri断裂(AF)、千岛盆地南部边界(KBB2)。它们的空间延伸长度大,贯穿多个构造单元,是大型构造单元的边界。EJSF、TF、TTF共同构成了欧亚板块与鄂霍次克板块(北美板块)的边界;TIF、CHSF、AF组成的右旋走滑剪切带是控制其西侧日本盆地及东侧千岛盆地发育的重要边界;KBB2是鄂霍次克海的南边界。
二级断裂包括:Hatagawa构造线(HTL)、北海道断裂(R3~R8断裂)。它们空间延伸长度中等或较大,贯穿整个构造单元,是构造单元内部次级单元的重要边界。HTL将日本本州岛分为中部日本和东北日本;R3~R8断裂是中央北海道-萨哈林岛剪切带的组成部分,也是日高山脉等北海道山系的边界。
三级及以上断裂包括逆断性质的R1~R2断裂、R9~R10断裂。它们空间延伸长度小,产状、性质多样化,对局部小构造起控制和调节作用。
4 结语本文获得了日本-千岛海沟转折处俯冲带之下200 km深度范围内的震源重定位及P波速度结构图件结果,并将日本-千岛海沟转折处俯冲带划分为北日本海沟俯冲带及南千岛海沟俯冲带。虽然北日本海沟俯冲带和南千岛海沟俯冲带在海沟走向、太平洋板块俯冲倾角、地震分布、速度异常体等方面存在差异,但两者之下200 km深度范围内的太平洋俯冲板块是连续的,且在日本海沟俯冲带之下弯曲上拱,在千岛海沟俯冲带之下弯曲下凹。
通过整合前人及本文的研究成果,综合分析断裂的走向、倾向、倾角、性质、平面展布长度、切割深度及对相关构造单元的控制作用等特征,制定统一的划分标准,将研究区内的断裂带大致分为三级。一级断裂包括Tartary断裂(TF)、Tartary-Tanakura断裂(TTF)、Tartary-Ishikari断裂(TIF)、东日本海断裂(EJSF)、中央萨哈林-北海道断裂(CHSF)、Abashiri断裂(AF)、千岛盆地南部边界(KBB2);二级断裂包括Hatagawa构造线(HTL)、北海道断裂(R3~R8断裂);三级及以上断裂包括:逆断性质的R1~R2断裂、R9~R10断裂。
[1] |
洪汉净. 火山预测与预警[M]. 北京: 地震出版社, 2013 (Hong Hanjing. Prediction and Early-Warning of Volcano[M]. Beijing: Seismological Press, 2013)
(1) |
[2] |
刘鑫. 西北太平洋俯冲带构造特征及其对弧前大地震成因的影响[D]. 青岛: 中国海洋大学, 2013 (Liu Xin. Structual Heterogeneity of Northwest-Pacific Subduction Zones and Its Implications for Interplate Megathrust Earthquakes [D]. Qingdao: Ocean University of China, 2013)
(5) |
[3] |
Li Z X, Zhang L, Powell C M A. South China in Rodinia: Part of the Missing Link between Australia East Antarctica and Laurentia?[J]. Geology, 1995, 23(5): 407-411 DOI:10.1130/0091-7613(1995)023<0407:SCIRPO>2.3.CO;2
(1) |
[4] |
Maruyama S, Isozaki Y, Kimura G, et al. Paleogeographic Maps of the Japanese Islands: Plate Tectonic Synthesis from 750 Ma to the Present[J]. Island Arc, 1997, 6(1): 121-142 DOI:10.1111/iar.1997.6.issue-1
(0) |
[5] |
Maruyama S, Liou J G, Zhang R. Tectonic Evolution of the Ultrahigh-Pressure(UHP) and High-Pressure(HP) Metamorphic Belts from Central China[J]. Island Arc, 1994, 3(2): 112-121 DOI:10.1111/iar.1994.3.issue-2
(1) |
[6] |
Isozaki Y. Studies on Orogeny Based on Plate Tectonics in Japan and New Geotectonic Subdivision of the Japanese Islands[J]. Journal of Geography, 1991, 100(5): 697-761 DOI:10.5026/jgeography.100.5_697
(1) |
[7] |
Miller M S, Kennett B L N, Gorbatov A. Morphology of the Distorted Subducted Pacific Slab beneath the Hokkaido Corner, Japan[J]. Physics of the Earth & Planetary Interiors, 2006, 156(1-2): 1-11
(1) |
[8] |
Fournier M, Jolivet L, Huchon P, et al. Neogene Strike-slip Faulting in Sakhalin and the Japan Sea Opening[J]. Journal of Geophysical Research, 1994, 99(B2): 2 701-2 725 DOI:10.1029/93JB02026
(2) |
[9] |
Gaku K, Kensaku T. Collision, Rotation, and Back-Arc Spreading in the Region of the Okhotsk and Japan Seas[J]. Tectonics, 1986, 5(3): 389-401 DOI:10.1029/TC005i003p00389
(0) |
[10] |
Ienaga M, Mcneill L C, Mikada H, et al. Borehole Image Analysis of the Nankai Accretionary Wedge, ODP Leg 196: Structural and Stress Studies[J]. Tectonophysics, 2006, 426(1-2): 207-220 DOI:10.1016/j.tecto.2006.02.018
(0) |
[11] |
Isozaki Y, Ehiro M, Nakahata H, et al. Cambrian Plutonism in Northeast Japan and Its Significance for the Earliest Arc-Trench System of Proto-Japan: New U-Pb Zircon Ages of the Oldest Granitoids in the Kitakami and Ou Mountains[J]. Journal of Asian Earth Sciences, 2015, 108: 136-149 DOI:10.1016/j.jseaes.2015.04.024
(0) |
[12] |
Lallemand S, Jolivet L. Japan Sea: A Pull-Apart Basin?[J]. Earth & Planetary Science Letters, 1986, 76(3-4): 375-389
(0) |
[13] |
Martin A K. Double Saloon Door Tectonics in the Japan Sea, Fossa Magna, and the Japanese Island Arc[J]. Tectonophysics, 2011, 498(1-4): 45-65 DOI:10.1016/j.tecto.2010.11.016
(0) |
[14] |
Naruse H. Cretaceous to Paleocene Depositional History of North-Pacific Subduction Zone: Reconstruction from the Nemuro Group, Eastern Hokkaido, Northern Japan[J]. Cretaceous Research, 2003, 24(1): 55-71 DOI:10.1016/S0195-6671(03)00024-7
(0) |
[15] |
Tujino T, Murakami F. Changes in Planview Geometry and Erosion-Deposition Patterns along a Channel-Levee Complex in the Southern Corner of the Kurile Basin, Okhotsk Sea[J]. Marine Geology, 2006, 230(1-2): 127-146 DOI:10.1016/j.margeo.2006.04.007
(0) |
[16] |
Wakita K. Origin of Chaotically Mixed Rockbodies in the Early Jurassic to Early Cretaceous Sedimentary Complex of the Mino Terrane, Central Japan[J]. Bulletin of the Geological Survey of Japan, 1988, 39: 675-757
(2) |
[17] |
Engdahl E R, Hilst R V D, Buland R. Global Teleseismic Earthquake Location with Improved Travel Times and Procedures for Depth Determination[J]. Bulletin of the Seismological Society of America, 1998, 88(3): 722-743
(1) |
[18] |
Zhao D, Hasegawa A, Horiuchi S. Tomographic Imaging of P and S Wave Velocity Structure beneath Northeastern Japan[J]. Journal of Geophysical Research Solid Earth, 1992, 97(B13): 19 909-19 928 DOI:10.1029/92JB00603
(1) |
[19] |
Zhao D, Yanada T, Hasegawa A, et al. Imaging the Subducting Slabs and Mantle Upwelling under the Japan Islands[J]. Geophysical Journal International, 2012, 190(2): 816-828 DOI:10.1111/gji.2012.190.issue-2
(1) |
[20] |
Zhang H H. Double-Difference Seismic Tomography Method and Its Applications[D]. Madiso: University of Wisconsin, 2003
(2) |
[21] |
Thurber C H. Earthquake Locations and Three-Dimensional Crustal Structure in the Coyote Lake Area, Central California[J]. Journal of Geophysical Research Atmospheres, 1983, 88(B10): 8 226-8 236 DOI:10.1029/JB088iB10p08226
(1) |
[22] |
Paige C C. Algorithm 583: LSQR: Sparse Linear Equations and Least Squares Problems[J]. Acm Transactions on Mathematical Software, 1982, 8(2): 195-209 DOI:10.1145/355993.356000
(1) |
[23] |
Martin R, Donna E, Graham S. A Three Dimensional Image of Shallow Subduction: Crustal Structure of the Raukumara Peninsula, New Zeland[J]. Geophysical Journal International, 1999, 137(3): 873-890 DOI:10.1046/j.1365-246x.1999.00842.x
(1) |
[24] |
Toomey D R, Foulger G R. Tomographic Inversion of Local Earthquake Data from the Hengill-Grensdalur Central Volcano Complex, Iceland[J]. Journal of Geophysical Research Atmospheres, 1989, 94(B12): 17 497-17 510 DOI:10.1029/JB094iB12p17497
(1) |
[25] |
Scarfì L, Giampiccolo E, Musumeci C, et al. New Insights on 3D Crustal Structure in Southeastern Sicily(Italy) and Tectonic Implications from an Adaptive Mesh Seismic Tomography[J]. Physics of the Earth & Planetary Interiors, 2007, 161(15): 74-85
(1) |
[26] |
Kita S, Okada T, Hasegawa A, et al. Anomalous Deepening of a Seismic Belt in the Upper-Plane of the Double Seismic Zone in the Pacific Slab beneath the Hokkaido Corner: Possible Evidence for Thermal Shielding Caused by Subducted Forearc Crust Materials[J]. Earth and Planetary Science, 2010, 290(3-4): 415-426 DOI:10.1016/j.epsl.2009.12.038
(1) |
[27] |
任建业. 海洋底构造导论[M]. 武汉: 中国地质大学出版社, 2008 (Ren Jianye. An Introduction to Ocean Floor Tectonics[M]. Wuhan: China University of Geosciences Press, 2008)
(1) |
[28] |
Burbach V N, Frohlich C. Intermediate and Deep Seismicity and Lateral Structure of Subducted Lithosphere in the Circum-Pacific Region[J]. Reviews of Geophysics, 1986, 24(4): 833-874 DOI:10.1029/RG024i004p00833
(0) |
[29] |
Wada I, He J, Hasegawa A, et al. Mantle Wedge Flow Pattern and Thermal Structure in Northeast Japan: Effects of Oblique Subduction and 3-D Slab Geometry[J]. Earth & Planetary Science Letters, 2015, 426: 76-88
(1) |
2. College of Marine Geo-Sciences, Ocean University of China, 238 Songling Road, Qingdao 266100, China;
3. Institute of Geology, CEA, A1 Huayanli, Beijing 100029, China