1982年日本茨城地震序列表现出的局部板块俯冲动态过程

板块构造的运动和动力学是地学研究的前沿课题．板块边缘强震是板块运动的直接结果．分析1982年日本茨城地震序列发现，震源断层面存在分阶段沿走向扩展、倾向调整的特殊图象．结合震源机制的走向、倾向以及滑动方向等资料，确认属于日本海沟下方的俯冲带逆冲地震．文中讨论了地震序列过程的成因，认为采用介于年尺度的空间大地测量结果与秒尺度的强震矩张量分析之间的地震序列时空分布研究，表明太平洋板块向欧亚板块局部下方俯冲的动态过程．根据地震平均错距和板块运动速度，推测近年当地可能再次发生强震．      

Large earthquake doublets and fault plane heterogeneity in the northern Solomon Islands subduction zone

In the Solomon Islands and New Britain subduction zones, the largest earthquakes commonly occur as pairs with small separation in time, space and magnitude. This doublet behavior has been attributed to a pattern of fault plane heterogeneity consisting of closely spaced asperities such that the failure of one asperity triggers slip in adjacent asperities. We analyzed body waves of the January 31, 1974,M _w =7.3, February 1, 1974,M _w =7.4, July 20, 1975 (1437)M _w =7.6 and July 20, 1975 (1945),M _w =7.3 doublet events using an iterative, multiple station inversion technique to determine the spatio-temporal distribution of seismic moment release associated with these events. Although the 1974 doublet has smaller body wave moments than the 1975 events, their source histories are more complicated, lasting over 40 seconds and consisting of several subevents located near the epicentral regions. The second 1975 event is well modeled by a simple point source initiating at a depth of 15 km and rupturing an approximate 20 km region about the epicenter. The source history of the first 1975 event reveals a westerly propagating rupture, extending about 50 km from its hypocenter at a depth of 25 km. The asperities of the 1975 events are of comparable size and do not overlap one another, consistent with the asperity triggering hypothesis. The relatively large source areas and small seismic moments of the 1974 doublet events indicate failure of weaker portions of the fault plane in their epicentral regions. Variations in the roughness of the bathymetry of the subducting plate, accompanying subduction of the Woodlark Rise, may be responsible for changes in the mechanical properties of the plate interface.To understand how variations in fault plane coupling and strength affect the interplate seismicity pattern, we relocated 85 underthrusting earthquakes in the northern Solomon Islands Are since 1964. Relatively few smaller magnitude underthrusting events overlap the Solomon Islands doublet asperity regions, where fault coupling and strength are inferred to be the greatest. However, these asperity regions have been the sites of several previous earthquakes withM _s 7.0. The source regions of the 1974 doublet events, which we infer to be mechanically weak, contain many smaller magnitude events but have not generated any otherM _s 7.0 earthquakes in the historic past. The central portion of the northern Solomon Islands Arc between the two largest doublet events in 1971 (studied in detail bySchwartz et al., 1989a) and 1975 contains the greatest number of smaller magnitude underthrusting earthquakes. The location of this small region sandwiched between two strongly coupled portions of the plate interface suggest that it may be the site of the next large northern Solomon Islands earthquake. However, this region has experienced no known earthquakes withM _s 7.0 and may represent a relatively aseismic portion of the subduction zone.     

The 1986 Kermadec earthquake and its relation to plate segmentation

To evaluate the tectonic significance of the October 20, 1986 Kermadec earthquake (M _w =7.7), we performed a comprehensive analysis of source parameters using surface waves, body waves, and relocated aftershocks. Amplitude and phase spectra from up to 93 Rayleigh waves were inverted for centroid time, depth, and moment tensor in a two-step algorithm. In some of the inversions, the time function was parameterized to include information from the body-wave time function. The resulting source parameters were stable with respect to variations in the velocity and attenuation models assumed, the parameterization of the time function, and the set of Rayleigh waves included. The surface wave focal mechanism derived (=275°, =61°, =156°) is an oblique-compressional mechanism that is not easy to interpret in terms of subduction tectonics. A seismic moment of 4.5×10²⁰ N-m, a centroid depth of 45±5 km, and a centroid time of 13±3 s were obtained. Directivity was not resolvable from the surface waves. The short source duration is in significant contrast to many large earthquakes.We performed a simultaneous inversion ofP andSH body waves for focal mechanism and time function. The focal mechanism agreed roughly with the surface wave mechanism. Multiple focal mechanisms remain a possibility, but could not be resolved. The body waves indicate a short duration of slip (15 to 20 s), with secondary moment release 60s later. Seismically radiated energy was computed from the body-wave source spectrum. The stress drop computed from the seismic energy is about 30 bars. Sixty aftershocks that occurred within three months of the mainshock were relocated using the method of Joint Hypocentral Determination (JHD). Most of the aftershocks have underthrusting focal mechanisms and appear to represent triggered slip on the main thrust interface. The depth, relatively high stress drop, short duration of slip, and paucity of true aftershocks are consistent with intraplate faulting within the downgoing plate. Although it is not clear on which nodal plane slip occurred, several factors favor the roughly E-W trending plane. The event occurred near a major segmentation in the downgoing plate at depth, near a bend in the trench, and near a right-lateral offset of the volcanic are by 80 km along an E-W direction. Also, all events in the region from 1977 to 1991 with CMT focal mechanisms similar to that of the Mainshock occurred near the mainshock epicenter, rather than forming an elongate zone parallel to the trench as did the aftershock activity. We interpret this event as part of the process of segmentation or tearing of the subducting slab. This segmentation appears to be related to the subduction of the Louisville Ridge, which may act as an obstacle to subduction through its buoyancy.     

西北太平洋和汤加俯冲地区深震特殊聚簇的特性及成因探究

对比研究了具有不同热参数、不同俯冲形态的西北太平洋俯冲地区和汤加俯冲地区的深震特殊聚簇的地震学特性和成因.利用单键群方法探测到两个特殊的深震聚簇G1N和G1T.聚簇G1N位于地震空区下方,具有极低的b值（~0.54）,完全不同于具有高b值（~1.04）的汤加俯冲地区聚簇G1T.通过对聚簇地区板块形态、地震主应力轴、地震深度分布特征的分析,以及和汤加典型的板片折曲处地震活动性的对比,我们认为深震聚簇G1N附近的板块表现出板片折曲的特征,板块俯冲到地幔过渡带底部受到下地幔的黏性阻力,板片局部向上凸起发生折曲,产生局部的拉张应力,叠加在俯冲造成的压缩构造背景上,应力状态发生改变,从而影响该深震聚簇的地震活动性.汤加地区G1T聚簇深震的成因则完全不同,没有体现出板片折曲、应力变化的特征;相反,这些深震发生在较冷的Vitiza-Fiji俯冲板块上,该板块在5~8 Ma年前先行俯冲到G1T区域并与Tonga板块发生拆离,G1T聚簇深震就发生在这些温度依然很低、滞留于~500 km深度处的高速板片残留体上.

     

What controls the lateral variation of large earthquake occurrence along the Japan Trench?

Abstract The lateral (along trench axis) variation in the mode of large earthquake occurrence near the northern Japan Trench is explained by the variation in surface roughness of the subducting plate. The surface roughness of the ocean bottom near the trench is well correlated with the large-earthquake occurrence. The region where the ocean bottom is smooth is correlated with'typical'large underthrust earthquakes (e.g. the 1968 Tokachioki event) in the deeper part of the seismogenic plate interface, and there are no earthquakes in the shallow part (aseismic zone). The region where the ocean bottom is rough (well-developed horst and graben structure) is correlated with large normal faulting earthquakes (e.g. the 1933 Sanriku event) in the outer-rise region, and large tsunami earthquakes (e.g. the 1896 Sanriku event) in the shallow region of the plate interface zone. In the smooth surface region, the coherent metamorphosed sediments form a homogeneous, large and strong contact zone between the plates. The rupture of this large strong contact causes great under-thrust earthquakes. In the rough surface region, large outer-rise earthquakes enhance the well-developed horst and grabens. As these structure are subducted with sediments in the graben part, the horsts create enough contact with the overriding block to cause an earthquake in the shallow part of the interface zone, and this earthquake is likely to be a tsunami earthquake. When the horst and graben structure is further subducted, many small strong contacts between the plates are formed, and they can cause only small underthrust earthquakes.     

Seismicity in subduction zones from local and regional network observations

The data provided by local and regional seismograph networks are essential for the solution of many problems of subduction-zone seismology. The capabilities of such networks are limited by the instrumentation currently in common use and the unfavorable source-station geometry often imposed by the regional geography. Nevertheless, important contributions have come from the data gathered in many of the earth's subduction zones. The accuracy of hypocenter locations based on regional data is affected by the complex velocity structures characteristic of subduction zones, but the problems are now well-understood. Examples of numerous studies of the spatial configurations of the seismicity in subduction zones and consequent interpretations of seismogenesis and subduction processes are reviewed. Studies of the distributions of earthquakes in time and with magnitude, for events down to the microearthquake level, have the potential for clarifying the earthquake-generating processes and, possibly, a basis for earthquake prediction. Other uses of local and regional network data have been for investigations of coda-Q and the identification of asperities in subduction zones.     

2016年青海门源M_S 6.4地震特征与地质构造

2016年青海门源M_S 6.4地震具有显著特征,如:余震强度低、能量释放水平低、震害较轻等。震区地震构造表现出一定规律性,如:冷龙岭断裂呈周期性破裂,发震断裂含多支相互斜交的分支断裂等。虽然地震前小震活动未能提供有价值的预测信息,但震中落入2015年甘肃省年度危险区。本文对此次地震的震源机制解、地震序列衰减、震害特点及冷龙岭断裂带构造进行分析,给出地震序列属性、发震断层及错动动力源,提出地震并未发生在冷龙岭主断裂的证据,从而为震区及邻区地震活动状态与孕震机制判定提供参考,并为后续震情判定提供震例与数据。     

2003年宁河地震序列特征的分析

2003年4月天津宁河地震序列的发生, 使天津市大部分地区明显有感, 震中烈度Ⅴ度。 分析结果表明, 宁河地震序列从空间、 时间、 强度分布等方面均具有震群性质的基本特征, 但余震贫乏、 序列的频度和强度衰减快, 不同于传统意义的震群; 宁河地震序列发生的位置靠近唐山余震区的西南端, 又具有相对独立的空间特征; 3次4级左右地震的震源机制解不一致, 可能与震源区构造复杂有关, 是震源区调整让位的一种表现形式, 说明震源区应力仍在调整中。     

1996年四川田湾地震序列揭示的水库诱震结构

１９９６年３月２８日田湾地震序列是四川二滩水库地震遥测台网自１９９２年１月正式运行以来,发生在库区的规模最大的一次地震序列,主震震级ＭＬ＝４．０,持续时间近２月,后续余震３ ３个,其中震级在３．０以上的余震３个,因该序列地震位于二滩台网监控范围以外,为改善观测台站的分布条件,本文综合使用了二滩水库地震遥测台网,西昌地震遥测台网和四川区域台网约２５个台站的观测资料,利用工业爆破和天然地震资料,计算耻     

日本俯冲带地震发生过程的数值模拟研究

了解地震发生的动力学机制是研究地震发震原因的关键,而数值模拟的方法是高速、有效的手段.2011年3月11日日本东北部宫城县发生9.0级大地震,文中以该次大地震所在的日本俯冲带为研究对象,通过使用黏弹性有限元数值模拟,并引用接触对,建立了研究区二维数值模型,模拟俯冲带与上覆板片之间的滑动、黏滞到再滑动的过程,亦即断层失稳发生地震的过程.模拟结果显示,随着太平洋板块不断俯冲,在俯冲带上自发出现了断层闭锁、解锁到再闭锁的黏滑过程,且这种过程呈现一定的准周期性,大事件主要集中分布在20~30 km的深度范围内.根据俯冲带可能在俯冲过程中角度的变化,建立了不同的模型,进行模拟对比研究,结果表明,俯冲带的几何形态,以及俯冲角度变化所在的不同深度,对模拟的结果有不同的影响.     

姚安地震序列与大姚地震序列震源参数的对比研究

2000年1月15日在云南省姚安县发生了M_S6.5地震, 2003年7月21日和10月16日在云南省大姚县发生了M_S6.2、 M_S6.1地震。 这两个地震序列均发生在滇中块体, 震中位置相距42.5 km, 构造、 应力背景相似, 但序列类型不同, 时间间隔短。 为深入研究姚安地震和大姚地震序列的特征, 根据昆明区域数字台网记录的数字地震波资料, 利用遗传算法计算了2000年姚安地震和2003年大姚地震序列的震源参数。 研究结果表明, 姚安地震序列和大姚地震序列的地震矩和近震震级有很好的线性关系, 它们的地震矩M₀均在10¹²~10¹⁴ N·m, 地震矩与震源半径之间呈线性关系。 应力降随震级变化不明显。 拐角频率和地震矩之间有明显的依赖关系, 用最小二乘法拟合拐角频率和地震矩的关系式, 利用此关系式可计算出给定地震矩的拐角频率估计值f_a, 进而可分析实际计算的拐角频率与估计值之差(f_c-f_a)的时间变化曲线。 姚安地震序列和大姚地震序列实际计算的拐角频率与估计值之差(f_c-f_a)及应力降随时间的变化表现出了不同的特征: 姚安地震序列的应力降和拐角频率差值在4.3级和4.6级强余震前均有一个上升-下降过程, 反映余震区应力场有一个增高-下降过程; 大姚地震序列应力降和拐角频率差值在M_S≥4.6强余震发生之前均会出现高值异常, 但在M_S4.6～4.7地震之前出现的高值异常持续时间短, 在M_S6.1地震前出现的高值异常持续较长。 文中结合两个地震的发震构造, 序列类型等特征进行了分析讨论。     

Rejuvenated extension of the Philippine Sea plate and its effect on subduction dynamics in the Nankai Trough

The Nankai Trough, Japan, is a subduction zone characterized by the recurrence of disastrous earthquakes and tsunamis. Slow earthquakes and associated tremor also occur intermittently and locally in the Nankai Trough and the causal relationship between slow earthquakes and large earthquakes is important to understanding subduction zone dynamics. The Nankai Trough off Muroto, Shikoku Island, near the southeast margin of the rupture segment of the 1946 Nankai earthquake, is one of three regions where slow earthquakes and tremor cluster in the Nankai Trough. On the Philippine Sea plate, the rifting of the central domain of the Shikoku Basin was aborted at ~15 Ma and underthrust the Nankai forearc off Muroto. Here, the Tosa-Bae seamount and other high-relief features, which are northern extension of the Kinan Seamount chain, have collided with and indented the forearc wedge. In this study, we analyzed seismic reflection profiles around the deformation front of accretionary wedge and stratigraphically correlated them to drilling sites off Muroto. Our results show that the previously aborted horst-and-graben structures, which were formed around the spreading center of the Shikoku Basin at ~15 Ma, were rejuvenated locally at ~6 Ma and more regionally at ~3.3 Ma and have remained active since. The reactivated normal faulting has enhanced seafloor roughness and appears to affect the locations of slow earthquakes and tremors. Rejuvenated normal faulting is not limited to areas near the Nankai Trough, and extends more than 200 km into the Shikoku Basin to the south. This extension might be due to extensional forces applied to the Philippine Sea plate, which appear to be driven by slab-pull in the Ryukyu and Philippine trenches along the western margin of the Philippine Sea plate.     

Characteristics of earthquake sequences in the Kuril Islands

IntroductionEarthquakesshowanon-homogeneousdistributionontheEarth'ssurface;inmostregions,therearenoorfeweanhquakestooccur.Mosteanhquakesconcentrateinsomenarrowandlonghugeseismicbelts,forexample,thecircum-Pacificseismicbelt,Alpine-Himalnyanseismicbeltandmid-oceanridgeseismicbelt.Themid-oceanridgeseismicbelttraversesthePacific,Indian,AtlanticandArcticOceanstoformaunifiedandinterconnectedmid-oceanridgeseismicbelt;itisalsoconnectedwiththecircum-Pacific,Alpine-Himalayanseismicbelts.Therefore,t…     

Oblique subduction, collision of microcontinents and subduction of oceanic ridge: Their implications on the Cretaceous tectonics of Japan

Abstract The Izanagi plate subducted rapidly and obliquely under the accretionary terrane of Japan in the Cretaceous before 85 Ma. A chain of microcontinents collided with it at about 140 Ma. In southwest Japan the major part of it subducted thereafter, but in northeast Japan it accreted and the trench jumped oceanward, resulting in a curved volcanic front. The oblique subduction and the underplated microcon-tinent caused uplifting of high-pressure (high-P) metamorphic rocks and large scale crustal shortening in southwest Japan. The oblique subduction caused left-lateral faulting and ductile shearing in northeast Japan. The arc sliver crossed over the high-temperature (high-T) zone of arc magmatism, resulting in a wide high-T metamorphosed belt. At about 85 Ma, the subduction mode changed from oblique to normal and the tectonic mode changed drastically. Just after this the Kula/Pacific ridge subducted and the subduction rate of the Pacific plate decreased gradually, causing the intrusion of huge amounts of granite magma and the eruption of acidic volcanics from large cauldrons. The oblique subduction of the Pacific plate resumed at 53 Ma and the left-lateral faults were reactivated.     

海洋板块俯冲作用下上覆大陆岩石层减薄机制的动力学模拟

几乎所有大陆岩石层的减薄现象,可能都与海洋板块的俯冲作用相关,但是两者之间的内在联系迄今仍不十分明确,为此,我们设计了一系列包含洋-陆俯冲系统的二维数值模型,来探讨海洋板块的俯冲作用对上覆大陆岩石层变形行为的影响,尤其对大陆岩石层减薄效应的制约.模型结果表明,海洋板块俯冲过程中的地幔楔熔体对大陆岩石层地幔的热侵蚀以及由熔体上升所诱发的地幔局部对流的强烈扰动会导致上覆大陆岩石层的减薄效应.这种效应不仅表现在横向上的向陆内蔓延,还表现在垂向上的向浅部发展.且多类动力学参数都能制约大陆岩石层的减薄效应.具体地,随着汇聚速率和洋壳厚度的增加,上覆大陆岩石层在横向上的减薄范围越大,在垂向上的减薄程度也越深;而随着俯冲海洋板块年龄的增加,上覆大陆岩石层在横向上的减薄范围增大,但在垂向上的减薄程度会减小;随着上覆大陆岩石层厚度的增加,其横向减薄范围会减小,但在垂向上的减薄程度会加深.本文研究成果能为揭示华北克拉通减薄/破坏的动力学过程提供一定的理论参考依据.

     

Temporal stress variations along a seismogenic megasplay fault in the subduction zone: An example from the Nobeoka Thrust,southwestern Japan

The Nobeoka Thrust, an ancient megasplay fault in the Shimanto Belt, southwestern Japan, contains fault rocks from the seismogenic zone, providing an accessible analog of active megasplay faults in deep subduction settings. In this study, the paleostress along the Nobeoka Thrust was analyzed using multiple inversion techniques, including k‐means clustering of fault datasets acquired from drillcores that intersected the thrust. The six resultant stress orientation clusters can be divided into two general groups: stress solutions with north–south‐trending σ₁ axes, and those with east–west‐trending σ₁ axes. These groups are characterized by the temporal changes for the orientations of the σ₁ and σ₃ principal stress axes that involve alternation between horizontal and vertical. The findings are probably due to a change in stress state before and after earthquakes that occurred on the fault; similar changes have been observed in active tectonic settings, such as the 2011 Tohoku‐Oki earthquake (Japan).     

地震发生过程的有限单元法模拟——以苏门答腊俯冲带上的大地震为例

地震预报由目前的经验预报走向物理预报,数值模拟地震过程是其中的关键.文中应用统一的数学公式表述了速率相关的摩擦接触中黏着（sticking）和滑移(sliding)这两种不同的运动状态;有限元计算中采用静力显示的时间积分方法,基于R最小策略,控制时间步长以保持力学状态变化稳定,从而保证有限元计算过程平稳、收敛.以2004年发生过Mm=9.3特大地震的苏门答腊俯冲带为例,模拟了俯冲带上俯冲板片与上伏板块之间的闭锁、解锁、滑动到再闭锁这一准周期性过程,即地震的孕育、发生过程.计算结果表明,俯冲带上具有较大尺度、介质均匀、摩擦系数相同的区域是产生大规模、大幅度整体突然滑动（即大地震）的条件;模拟的苏门答腊俯冲带上的大地震在时间上有准周期性,空间上有迁移特征,破裂由深部向浅部进行;此外,俯冲带的几何特征对大地震的震源位置有很大的影响.     

Rupture model of the great AD 365 Crete earthquake in the southwestern part of the Hellenic Arc

An M8.3 earthquake struck the southwestern part of the Hellenic Arc, near the Island of Crete, in AD 365, generating a tsunami that affected almost the entire eastern Mediterranean region. Taking into account that the time history of seismicity in this region is fairly complete for such earthquakes in the historical catalog, which can be dated as back as the 5th century B.C., there is no indication that this segment of plate boundary has been fully ruptured again. The seismic hazard associated with this part of the Hellenic Arc necessitates the evaluation of the rupture characteristics of this great event. The constraint of the faulting geometry was initially achieved by using information from seismicity, and the focal mechanisms of earthquakes that occurred during the instrumental period. A rupture model for this great earthquake is constructed by assuming an elastic medium and calculating the theoretical surface displacements for various fault models that are matched with the observed surface deformation gleaned from historical reports. The resulted fault model concerns thrust faulting with a rupture length of 160 km and a seismic moment of 5.7 × 10²⁸ dyn·cm, an average slip of 8.9 m and a corresponding moment magnitude equal to 8.4, in excellent agreement with the macroseismic estimation. The absence of such events recurrence is an indication of the lack of complete seismic coupling that is common in subduction zones, which is in accordance with the back arc spreading of the Aegean microplate and with previous results showing low coupling for extensional strain of the upper plate.     

1999年岫岩地震序列尾波Qc的变化过程

采用Sato单次散射模型，利用近场台站岫岩台和营口台波形资料，计算了1999年岫岩地震序列s尾波品质因子Q，得到了3个分量10个频率点上Q随时间进程的变化。将序列分为3个时段：前震时段、强余震活动时段、序列结束时段，分别计算出每个时段内Qr与频率的关系Q(f)，最后计算了整个序列的Qr(f)。结果表明，Q在大震前后有较明显的变化，基本趋势是震前Q增高，震后降低。另外，在不同的台站和不同的分量上Q变化有不同的反映灵敏度。比较该方法与Aki单次散射模型方法在计算结果上的差别，认为Sato方法在本的数据处理上更为合理。     

Reconciling strong slab pull and weak plate bending: The plate motion constraint on the strength of mantle slabs

Although subducting slabs undergo a bending deformation that resists tectonic plate motions, the magnitude of this resistance is not known because of poor constraints on slab strength. However, because slab bending slows the relatively rapid motions of oceanic plates, observed plate motions constrain the importance of bending. We estimated the slab pull force and the bending resistance globally for 207 subduction zone transects using new measurements of the bending curvature determined from slab seismicity. Predicting plate motions using a global mantle flow model, we constrain the viscosity of the bending slab to be at most ~ 300 times more viscous than the upper mantle; stronger slabs are intolerably slowed by the bending deformation. Weaker slabs, however, cannot transmit a pull force sufficient to explain rapid trenchward plate motions unless slabs stretch faster than seismically observed rates of ~ 10^− 15 s^− 1. The constrained bending viscosity (~ 2 × 10²³ Pa s) is larger than previous estimates that yielded similar or larger bending resistance (here ~ 25% of forces). This apparent discrepancy occurs because slabs bend more gently than previously thought, with an average radius of curvature of 390 km that permits subduction of strong slabs. This gentle bending may ultimately permit plate tectonics on Earth.