基于震源机制解的鄂霍次克微板块东部俯冲带地区构造应力场反演

基于国际地震中心(ISC)提供的1970年1月～2016年12月期间的地震震源机制解,对鄂霍次克微板块东部俯冲带地区进行了应力张量反演,得到了日本海沟、千岛海沟和勘察加海沟3个俯冲带区域的构造应力场特征。研究结果显示:①海沟地区浅部区域(h100km)的水平主压应力轴与西北太平洋板块的俯冲方位一致,与海沟走向近似垂直,其洋壳一侧以拉张型应力状态为主,而陆壳一侧则以挤压型应力为主,且在弧后区域均存在拉张的应力状态;80～200km深度范围区域表现出双地震带"Ⅰ"型构造应力场特征。②日本海沟带由于俯冲角相对较小(相比于千岛海沟和勘察加海沟),水平方向沿NWW向延伸更远,大洋板块与上覆板块之间耦合更加强烈,逆冲型地震发生数量最多。③对于深部区域(h300km),千岛地区应力场表现出非均匀性特征,可能是由地幔阻力导致的;而勘察加地区应力场表现出拉张型,可能是因为俯冲板片的拉伸拖曳作用更强。     

Metamorphism,fluid behavior and magmatism in oceanic subduction zones

Based on the updated results of experimental petrology and phase equilibria modelling and combined with the available thermal structure models of subduction zones, this paper presents an overview on the dehydration and melting of basic,sedimentary and ultrabasic rocks that occur in the different stages during oceanic subduction processes and their influences on magmatism above subduction zones. During the subduction at the forearc depth of <90–100 km, the basic and ultrabasic rocks from most oceanic slabs can release very small amounts of water, and significant dehydration may occur in the slab superficial sediments. Strong dehydration occurs in both basic and ultrabasic rocks during subduction at the subarc depth of 90–200 km. For example, more than 90% water in basic rocks is released by the successive dehydration of chlorite, glaucophane, talc and lawsonite in the subarc depths. This is diversely in contrast to the previous results from synthetic experiments. Ultrabasic rocks may undergo strong dehydration through antigorite, chlorite and phase 10 ? at the subarc depth of 120–220 km. However,sediments can contribute minor fluids at the subarc depth, one main hydrous mineral in which is phengite(muscovite). It can stabilize to ~300 km depth and transform into K-hollandite. After phengite breaks down, there will be no significant fluid release from oceanic slab until it is subducted to the mantle transition zone. In a few hot subduction zones, partial melting(especially flux melting) can occur in both sediments and basic rocks, generating hydrous granitic melts or supercritical fluids, and in carbonates-bearing sediments potassic carbonatite melts can be generated. In a few cold subduction zones, phase A occurs in ultrabasic rocks, which can bring water deep into the transition zone. The subducted rocks, especially the sediments, contain large quantities of incompatible minor and trace elements carried through fluids to greatly influence the geochemical compositions of the magma in subduction zones. As the geothermal gradients of subduction zones cannot cross the solidi of carbonated eclogite and peridotite during the subarc subduction stage, the carbonate minerals in them can be carried into the deep mantle.Carbonated eclogite can melt to generate alkali-rich carbonatite melts at >400 km depth, while carbonated peridotite will not melt in the mantle transition zone below a subduction zone.     

2015年3月新不列颠MS7.4地震震源及邻区构造应力场特征

2015年3月30日至5月15日,巴布亚新几内亚-新不列颠地区发生了一系列地震.为研究该地区的构造应力环境及孕震背景,本文基于Global CMT目录,对新不列颠区域浅部进行构造应力场反演,拟得到高精度的应力图像.反演结果显示：（1）沿着南、北俾斯麦块体边界的区域构造应力场呈走滑体系,最大主压应力轴方位呈SWW-NEE向.（2）所罗门海的NW和NE走向的海沟处于压缩状态,所罗门海块体向新不列颠和所罗门群岛俯冲的板块弯曲部分是局部拉张.（3）受俯冲带的北向推挤,南俾斯麦板块顺时针旋转的挤压,太平洋板块向西部运动汇聚作用,新不列颠岛东北部与新爱尔兰岛南部交汇区域呈现明显非均匀应力状态.（4）此次地震序列的大多数走滑型和逆冲型地震,可能是所罗门海块体俯冲运动,和南俾斯麦块体与太平洋板块的近EW向挤压作用共同引发.     

日本海沟俯冲带MW9.0地震震源区应力场演化分析

于2011年3月11日发生在日本东北部的M_W9.0级逆冲型板间地震是日本有地震记录以来震级最大的一次地震.本研究基于NIED F-net矩张量解目录中的震源机制解,选取两个长轴相互垂直的矩形区域进行应力场2D反演,获取了日本海沟俯冲带地区应力场的空间及时间分布图像.结果表明：主震前,俯冲带地区应力状态在空间上大体趋于一致,即应力轴（P轴、σ₁轴及S_Hmax轴）系统性地倾向板块汇聚方向,P轴、σ₁轴倾角整体偏缓（<30°）,且远离震源区及日本海沟东侧区域内的应力轴倾角普遍大于主震震源区内应力轴倾角;主震前,受2003年5月26日在宫城县北部发生的M_W7.0地震影响,位于M_W9.0地震震源区西北侧的应力场出现明显扰动,σ₁轴倾向顺时针偏转150°~180°,并于之后大体恢复至震前状态,同期其他地区没有明显变化,这种情况可能和主震断层局部（深部）的前兆性滑动有关;主震后,距离震源区较远处应力场变化不大,主震震源区内应力场发生显著改变,P轴及σ₁轴均以大角度（>60°）倾伏于板块汇聚方向,S_Hmax轴顺时针偏转60°~90°且在日本海沟附近普遍平行于海沟轴.这项研究以时空图像的方式展示了大地震前应力场变化的特点,反映了大地震孕震过程中构造与地震的相互作用,对于理解大地震孕震过程有重要意义.     

华北北部地区地壳三维P波速度结构的双差地震层析成像

利用区域固定台站和华北科学探测台阵记录的10 461个近震事件的183 909个Pg波绝对走时和495 753个相对走时数据,采用双差地震层析成像获得华北北部(37.5°～41.5°N,111.5°～119.5°E)范围内的地壳三维P波速度结构模型。结果表明:研究区内各主要构造单元具有明显不同的速度结构特征,速度异常的走向与区域构造的走向一致,浅层速度图像很好地反映了地表地质和岩性的变化;重定位后的大部分地震集中在0～20km的深度上,主要位于低速区的内部或高速和低速区的交界部位;三河—平谷和唐山地震震源区中、下地壳的低速异常可能是流体的显示。结合前人成果和本文模型所揭示的深、浅结构,我们认为太平洋板块在中国东部之下的俯冲和滞留引起板块脱水、软流圈物质上涌等一系列过程,软流圈热物质到达上地幔顶部并沿超壳断裂上侵进入地壳,致使上地幔顶部和下地壳中的含水矿物发生脱水作用产生流体,流体继续上移造成中、上地壳发震层的弱化,从而导致大地震的发生。因此华北北部地区的强震活动,以及地壳结构的非均匀性应是与太平洋板块俯冲、滞留引起的深部过程密切相关的。     

考虑俯冲带潜在震源区的地震危险性算法研究及在中国海域的应用

本文研究了俯冲带潜在震源区离散化方法及考虑俯冲带高震级地震震源破裂面和震源深度的场点地震动计算方法,推导了俯冲带潜在震源区地震危险性计算公式,并使用中国海域及邻区地震危险性模型进行地震危险性试算.结果表明,本文建立的考虑俯冲带潜在震源区的地震危险性算法能够实现场点地震危险性计算时对俯冲带高震级地震震源破裂面和震源深度的...     

Cymrite as mineral clathrate: An overlooked redox insensitive transporter of nitrogen in the mantle

We study the phase relations and mineral chemistry in the systems muscovite–NH₃–N₂-H₂O and eclogite + muscovite–NH₃–N₂-H₂O at 6.3–7.8 GPa, 1000–1200 °C, and oxygen fugacity (fO₂) from ∼IW (Fe–FeO) to ∼ NNO (Ni–NiO) equilibria. The quenched H₂O-bearing fluids differ in nitrogen speciation from NH₃-rich to N₂–rich, and the respective N₂/(NH₃+N₂) ratio varies from <0.1 to ∼ 1. N-bearing K-cymrite is obtained in association with a kyanite-garnet-jadeite ± muscovite ± corundum assemblage in the muscovite–NH₃–N₂-H₂O system and coexists with pyrope-almandine garnet and omphacite in the eclogite + muscovite–NH₃–N₂-H₂O system. The presence of an N-bearing fluid in the studied systems stabilizes the K-cymrite structure. Muscovite does not convert to K-cymrite in the absence of NH₃–N₂-bearing fluid up to 7.8 GPa and 1070–1120 °C. According to FTIR and Raman spectroscopy, K-cymrite in equilibrium with an N-rich fluid can capture a huge amount of nitrogen in cages of its framework, mainly as N₂ molecules at fO₂ ∼NNO and predominantly as NH₃ molecules at fO₂ ∼IW. The storage capacity of K-cymrite with respect to nitrogen increases from 2.9 to 6.3 wt% with increase of fO₂. FTIR spectroscopy of muscovite equilibrated with K-cymrite shows that the clathrate mechanism of nitrogen entrapment by aluminosilicates (as neutral N₂ and NH₃ molecules) is much more efficient than the K⁺ → (NH₄)⁺ substitution. The structure of N-bearing K-cymrite (K,(NH₄⁺))[AlSi₃O₈]·(N₂,NH₃,H₂O) determined using X-ray single-crystal diffraction is very similar to that of H₂O-bearing K- and Ba-cymrites. It includes aluminosilicate layers consisting of double six-membered tetrahedral rings and cation sites statistically occupied with K⁺, Ba²⁺ and (NH₄)⁺ on the six-fold symmetry axis in interlayer space. The N₂ and NH₃ molecules are located near the cage centers and, unlike H₂O molecules, are included in the coordination environment of the cations. Our study confirms that NH₃- and N₂-rich K-cymrite may be stable in metapelites and can act as a redox insensitive carrier of nitrogen to >250 km mantle depths in downgoing slabs. The stability field of N-rich K-cymrite in the presence of an N₂–H₂O–NH₃-bearing fluid is inferred to be P ≥ 4 GPa in metasediments rich in K-feldspar and P ≥ 6 GPa in those containing phengite. As the slab material sinks deeper than 250–300 km where N-bearing K-cymrite may lose stability, the releasing nitrogen may migrate to metal-saturated mantle and become stored there in γ−Fe, Fe₃C, metal melt, or even iron nitride phases.     

The potential influence of subduction zone polarity on overriding plate deformation, trench migration and slab dip angle

A geodynamic model exists, the westward lithospheric drift model, in which the variety of overriding plate deformation, trench migration and slab dip angles is explained by the polarity of subduction zones. The model predicts overriding plate extension, a fixed trench and a steep slab dip for westward-dipping subduction zones (e.g. Mariana) and predicts overriding plate shortening, oceanward trench retreat and a gentle slab dip for east to northeastward-dipping subduction zones (e.g. Chile). This paper investigates these predictions quantitatively with a global subduction zone analysis. The results show overriding plate extension for all dip directions (azimuth α = − 180° to 180°) and overriding plate shortening for dip directions with α = − 90° to 110°. The wide scatter in data negate any obvious trend and only local mean values in overriding plate deformation rate indicate that overriding plate extension is somewhat more prevalent for west-dipping slabs. West-dipping subduction zones are never fixed, irrespective of the choice of reference frame, while east to northeast-dipping subduction zones are both retreating and advancing in five out of seven global reference frames. In addition, westward-dipping subduction zones have a range in trench-migration velocities that is twice the magnitude of that for east to northeastward-dipping slabs. Finally, there is no recognizable correlation between slab dip direction and slab dip angle. East to northeast-dipping slabs (α = 30° to 120°) have shallow (0–125 km) slab dip angles in the range 10–60° and deep (125–670 km) slab dip angles in the range 40–82°, while west-dipping slabs (α = − 60° to − 120°) have shallow slab dip angles in the range 19–50° and deep slab dip angles in the range 25–86°. Local mean deep slab dip angles are nearly identical for east and west-dipping slabs, while local mean shallow slab dip angles are lower by only 4.7–8.1° for east to northeast-dipping slabs. It is thus concluded that overall, there is no observational basis to support the three predictions made by the westward drift model, and for some sub-predictions the observational basis is very weak at most. Alternative models, which incorporate and underline the importance of slab buoyancy-driven trench migration, slab width and overriding plate motion, are better candidates to explain the complexity of subduction zones, including the variety in trench-migration velocities, overriding plate deformation and slab dip angles.     

甘肃北山明舒井岩体形成年龄、地球化学特征及其地质意义

北山明舒井岩体为一由辉长岩、闪长岩和花岗岩构成的岩浆杂岩体。LA-ICP-MS锆石U-Pb定年分别获得辉长岩442.9±1.2Ma、闪长岩437.0±2.3Ma、花岗岩435.9±2.2Ma的形成年龄。其中,辉长岩为富钙、高铁和镁的钙碱性岩,相对富集LREE和LILE、贫化HFSE,亏损Nb、Ta、Ti、Zr和P,Pb正异常,具中等偏高的(87Sr/86Sr)i(0.705892)、高放射性成因Pb,低εNd(t)(0.1),锆石εHf(t)变化大(-7.1~+8.3)。闪长岩和花岗岩富钠、低钙、铁和镁,为准铝质钙碱性岩;它们富集LREE和LILE、亏损HFSE,Nb、Ta及Sr、P、Ti明显负异常;也显示中等偏高的(87Sr/86Sr)i(0.705951~0.706196)和高放射性成因Pb、低εNd(t)(-0.2~-0.1),但锆石εHf(t)明显偏正(-1.1~+12.0)。不同岩类的地球化学特征及其相对一致的Sr-Nd-Pb同位素及明显不同的锆石Hf同位素,指示它们形成于板块消减带构造环境,为源自受流体交代和陆壳物质改造的幔源岩浆与新元古代陆壳物质部分熔融的岩浆混合和分离结晶演化的产物,代表北山南部早古生代晚奥陶世-早志留世晚期与大洋俯冲消减作用有关的构造岩浆事件。     

中国大陆显生宙俯冲型、碰撞型和复合型片麻岩穹窿(群)

片麻岩穹窿(gneiss dome)是中下地壳热动力学过程产生的、与岩浆作用(或混合岩化作用)密切相关的穹窿状构造。片麻岩穹窿大部分是地壳深层次变动的产物,在世界范围内几乎出露在所有的折返造山带中,反映了所在地区地壳的大幅度抬升。片麻岩穹窿核部主要是无或弱岩浆组构的花岗岩体和高级变质岩(例如混合岩),边部是具有岩浆组构的花岗片麻岩,幔部由来自地壳深部的高级片岩和片麻岩组成。片麻岩穹窿的形成经历从垂直上升的地壳流导致的岩浆上涌的挤压收缩机制到岩浆体侵位的顶部伸展机制的转化过程。根据片麻岩穹窿的岩石组合、组构特征、成因机制和大地构造背景以及片麻岩穹窿与地壳流关系的分析,结合中国大陆典型片麻岩穹窿的研究,提出中国大陆显生宙的片麻岩穹窿和片麻岩穹窿群可以划分为与大洋岩石圈板片俯冲增生与随后的折返造山相关的"俯冲型"片麻岩穹窿(群),如秦岭片麻岩穹窿;与陆陆碰撞折返造山有关的"碰撞型"片麻岩穹窿(群),如北喜马拉雅拉轨岗日片麻岩穹窿(群)和松潘甘孜雅江片麻岩穹窿(群);与俯冲和碰撞的叠合作用有关的"复合式"片麻岩穹窿(群),如帕米尔空喀山片麻岩穹窿和东冈底斯林芝片麻岩穹窿(群)。