西北太平洋俯冲带日本本州至中国东北段应力场反演

本研究基于Global CMT提供的1196个1976年11月—2017年1月MW4.6地震矩心矩张量解,对西北太平洋俯冲带日本本州至中国东北段的应力场进行反演计算,得到了从浅表到深部俯冲带应力状态的完整分布.结果显示:俯冲带浅表陆壳一侧应力场呈现水平挤压、垂向拉伸状态,洋壳一侧的应力状态则相反,即近水平拉张、近垂向压缩.沿着俯冲板片向下,应力主轴逐渐向俯冲板片轮廓靠拢,其中位于双地震层(120km深度附近)之上的部分,主张应力轴沿俯冲板片轮廓展布而又比其更为陡倾;双地震层内的应力模式同典型I型双层地震带内的应力模式一致,即上层沿俯冲板片轮廓压缩、下层沿俯冲板片轮廓拉伸;双地震层之下,应力模式逐步转变为主压应力轴平行于俯冲板片轮廓.通观所研究的整个俯冲系统,水平面内主压和主张应力轴基本保持了与西北太平洋板片俯冲方向上的一致性,同经典俯冲板片的应力导管模型所预言的俯冲带应力模式相符;而主张应力轴在俯冲板片表面之下的中源地震深度范围内转向海沟走向,或许同研究区域横跨日本海沟与千岛海沟结合带,改变的浅部海沟形态致使完整俯冲板片下部产生横向变形有关.     

Where and why do large shallow intraslab earthquakes occur?

We try to find how often, and in what regions large earthquakes (M≥7.0) occur within the shallow portion (20-60 km depth) of a subducting slab. Searching for events in published individual studies and the Harvard University centroid moment tensor catalogue, we find twenty such events in E. Hokkaido, Kyushu-SW, Japan, S. Mariana, Manila, Sumatra, Vanuatu, N. Chile, C. Peru, El Salvador, Mexico, N. Cascadia and Alaska. Slab stresses revealed from the mechanism solutions of these large intraslab events and nearby smaller events are almost always down-dip tensional. Except for E. Hokkaido, Manila, and Sumatra, the upper plate shows horizontal stress gradient in the arc-perpendicular direction. We infer that shear tractions are operating at the base of the upper plate in this direction to produce the observed gradient and compression in the outer fore-arc, balancing the down-dip tensional stress of the slab. This tectonic situation in the subduction zone might be realized as part of the convection system with some conditions, as shown by previous numerical simulations.     

Relocation of shallow earthquakes in southern Alaska using Joint Hypocenter Determination method

A subset of 2660 shallow earthquakes (0–50 km) that occurred from 1988 to 1996 in south central Alaska between 155°W and 145°W and 59°N and 63°N was relocated using the joint hypocenter determination (JHD) method. Both P- and S-wave observations recorded by the regional seismic network were used. Events were relocated in twenty different groups based on their geographic location and depth using two velocity models. As a result of the relocation, the majority of the hypocenters shifted downward, while the epicenter locations did not change significantly. The distribution of the shallow subduction zone earthquakes indicates the existence of two seismically independent blocks, with one block occupying the northeastern part and the other occupying the central and western parts of the study area. The boundary between the blocks is marked by a 15 to 20 km wide seismicity gap to the southeast of 149.5°W and 62°N. The analysis of the fault plane solutions for shallow subduction zone earthquakes shows that an overwhelming majority of the solutions represent normal, oblique-normal or strike-slip faulting with predominant WNW-ESE orientation of T-axes. This indicates a down-dip extensional regime for the subducting slab at shallow depths. Very few earthquakes yielded fault plane solutions consistent with thrusting on a contact zone between the overriding and subducting plates. This result may be an indication that currently either the strain energy is not released at the contact zone or it is associated with aseismic motion.     

Stresses along the metastable wedge of olivine in a subducting slab: possible explanation for the Tonga double seismic layer

A numerical calculation of the stresses associated with changes in volume during phase transitions of olivine in a descending slab results in a double layer of high shear stress along the metastable olivine wedge in the depth range 350-460 km. Stress in the upper layer is in-plane tensional and stress in the lower layer is down-dip compressional. The modeled stress field agrees with observations of stress in the Tonga double seismic zone. High shear stress also exists in the slab at depths below the metastable wedge. This stress distribution involves down-dip compression and trench-parallel tension, which agrees with about half of the focal mechanisms in the Tonga slab at depths of 460-690 km. The model supports the idea that at least two possible stress release mechanisms for deep earthquakes may act in the Tonga subducting slab. One, transformational faulting, is restricted to the metastable wedge while the other one acts below the metastable wedge.     

琉球岛弧地区的地震分布、Benioff带形态及应力状态h

利用ISC及中国台网的资料,研究了琉球岛弧及冲绳海槽的地震分布及震源机制解,讨论了Benioff带的形态及应力状态.mb4.0的地震主要分布于琉球海沟西侧的弧形带,并形成明显的Benioff带.吐噶喇海峡以北俯冲带弯曲明显,深部倾角大,约92,70km以下张应力轴沿俯冲方向;吐噶喇海峡以南,俯冲带较平直,深部倾角较小,约55,压应力轴基本沿俯冲方向。冲绳海槽内处于NNW向近水平的拉伸,华北应力场与之类似,没有受到菲律宾海块挤压作用的影响.      

帕米尔-兴都库什地区地震空间分布特征及应力场特征

本文研究了兴都库什及帕米尔地区地震的空间分布.发现h<70km的地震分布广泛,h≥100km的地震形成-S形的倾斜中源地震带.在71.5°E以西,中源地震带倾向接近正北,倾角随深度变化,在深部接近垂直,且倾角自西向东逐渐变陡,在71.5°E以东,倾向逐渐由东南变为正南. 分析了121个m_b≥5.0地震的机制解.浅源地震机制解的P轴大多位于NS和NNW-SSE方向,且多近水平,反映此区受到NS或NNW-SSE方向挤压.各剖面应力轴分布规律性强,在150km以下,总的趋势是机制解的T轴接近于倾斜的中源地震带的下倾方向,而P轴倾角较小且垂直于倾斜的中源地震带的走向.     

Active deformation in the Vrancea region,Rumania

Recent and historical seismicity as well as reliable fault plane solutions are used to study the active deformation caused by the occurrence of intermediate depth (60–170 km) earthquakes of the Vrancea region, Rumania. In this area, located in the southeastern part of the Carpathian arc, the westward subduction of the Carpathian trench has terminated, leaving continental lithosphere, at present, at the arc. The principalT axis of the intermediate depth events trends N159°E and has a plunge of 74°, which is the same as the dip of the subducted plate. TheP axis has a trend of 314° and a shallow plunge of 15°. The analysis of the moment tensor of six focal mechanisms showed that the dominant mode of deformation of the subducted lithosphere is a down-dip extension at a rate of about 2 cm/yr, based on seismicity data.     

The geometry of the Burmese-Andaman subducting lithosphere

The gross seismotectonic features for the Burmese-Andaman arc system which defines the northeast margin of the Indian plate are rather well known but variations in the subduction zone geometry along and across the arc and fault pattern within the subducting Indian plate have not been studied. Present workaims to study these by using seismicity data whose results are presented in the form of: (a) Lithospheric across-the-arc sections at about every 100–120 km (approximately one degree latitude apart) covering the 3500 km longBurmese-Andaman arc system, (b) a structure contour map showing the depth tothe top surface of the seismically active lithosphere and (c) interpretationof focal mechanism solutions for 148 Benioff zone earthquakes. Both penetrationdepth and the dip of the Benioff zone vary considerably along the arc in correspondence to the curvature of the fold-thrust belt which varies from concave to convex in different sectors of the arc. Several extensive `Hinge faults' that abut at high angles to the arc orientation, are inferred from aninterpretation of the structure contour map. Active nature of the hinge faultsis established in several areas by their association with earthquakes andcorroborated through fault plane solutions. At shallow level of the Benioffzone along these faults, focal mechanism solutions display left lateral strikeslip movement while at deeper levels reverse fault solutions are common.     

日本海-鄂霍次克海下俯冲带的应力状态及其深部形变

通过地震分布及地震机制解所反映的日本海-鄂霍次克海俯冲带的形态及应力状态,研究了俯冲带深部形变及650km间断面的穿透问题.日本海Benioff带较直,连续性较好;鄂霍次克海Benioff带弯度稍大,220-320km深度之间地震很少.两俯冲带在浅部及深部地震密集,100-200km深度之间有双地震层.应力状态随深度变化,200km深度以下P,T轴方向相对集中,P轴接近俯冲方向,在约100-200km深度附近,P,T轴均接近俯冲方向.观测和理论地震图拟合分析表明,地震断层面走向接近俯冲带走向,断裂的结果使俯冲带在深部倾角变小.     

Lateral structure of the subducting pacific plate beneath the Hokkaido corner from intermediate and deep earthquakes

We present a study of the lateral structure and mode of deformation in the transition between the Kuril and Honshu subduction zones. We begin by examining the source characteristics of the January 19, 1969, intermediate depth earthquake north of Hokkaido in the framework of slab-tearing, which for the December 6, 1978 event has been well documented by previous studies. We use a least-squares body wave inversion technique, and find that its focal mechanism is comparable to the 1978 event. To understand the cause of these earthquakes, which in the case of the 1978 event occurred on a vertical tear fault but does not represent hinge faulting, we examine the available International Seismological Centre [ISC] hypocenters and Harvard centroid-moment tensor [CMT] solutions to determine the state of stress, and lateral structure and segmentation in the Kuril and northern Honshu slabs. These data are evaluated in the framework of two models. Model (A) requires the subducting slab at the Hokkaido corner to maintain surface area. Model (B) requires slab subduction to be dominated by gravity, with material subducting in the down-dip direction. The distribution of ICS hypocenters shows a gap in deep seismicity down-dip of the Hokkaido corner, supporting model (B). From the CMT data set we find that three types of earthquake focal mechanisms occur. The first (type A) represents dip-slip mechanisms consistent with down-dip tension or compression in the slab in a direction normal to the strike of the trench. These events occur throughout the Honshu and Kuril slabs with focal mechanisms beneath Hokkaido showing NNW plungingP andT axes consistent with the local slab geometry. The second (type B) occurs primarily at depths over 300 km in the southern part of the Kuril slab with a few events in the northern end of the Honshu deep seismicity. These earthquakes have focal mechanisms with P axes oriented roughly E-W, highly oblique to the direction of compression found in the type A events, with which they are spatially interspersed. The third (type C) group of earthquakes are those events which do not fit in either of the first two groups and consist of either strike-slip focal mechanisms, such as the tearing events, or oddly oriented focal mechanisms. Examination of the stress axes orientations for these three types reveals that the compressional axes of the type C events are consistent with those of type B. The slab tearing events are just differential motion reflecting the E-W compressive states of stress which is responsible for the type B family of events. There is no need to invoke down-dip extension which does not fit the slab geometry. We conclude that these two states of stress can be explained as follows: 1) The type A events and the seismicity distribution support model (B). 2) The type B and C events upport model (A). The solution is that the slab subducts according to model (B), but the flow in the mantle maintains a different trajectory, possibly induced by the plate motions, which produces the second state of E-W compressive stress.     

缅甸山弧地区Benioff带的形态及其应力状态

研究了缅甸山弧附近的中源地震分布,发现h>70km的地震主要分布在20°N-27°N之间,形成明显的条带分布,24°N以南走向近南北,24°N以北走向逐渐接近NE;通过垂直剖面的研究,发现缅甸山弧下的Benioff带形态是变化的,在南北两端倾角较小,且较为平直,延伸深度浅,小于100km;在地震带的中间部分,Benioff带的倾角逐渐加大,且倾角随深度加深而增大,延伸深度可达180km。在一些剖面上出现双地震层,一般出现在45-100km的深度范围内,两层间的距离从10-25km不等;在同一剖面上,两层间在浅部间距大,在深部间距小。研究了沉降带上的应力状态,发现沉降带上P轴的优势方向位于NE-SW,T轴分布较分散;P、T轴随深度没有明显变化;在上地震层中,T轴明显接近于Benioff带的倾向;通过地壳内及沉降带上地震机制解的对比,发现前者的优势方向相对于后者逆时针旋转了一定角度。     

Space Variations of Stress Along the Tyrrhenian Wadati-Benioff Zone

—Gephart and Forsyth’s (1984) algorithm for stress inversion of earthquake fault-plane solutions has been applied to a set of ninety intermediate and deep events occurring in the southern Tyrrhenian region between 1976 and 1995. P- and S-wave data from local seismic networks in southern Italy, the Italian National Network and international bulletins, have been used for hypocenter and focal mechanism computations. Stress inversion runs performed after accurate selection and weighting of fault-plane solutions have allowed us to identify stress space variations at a higher level of detail than available from all previous investigations carried out in the study area. The maximum compressive stress has been shown to follow the depth-decreasing dip of the Wadati-Benioff zone, along the entire zone from a depth of 90 km, to the depth of the deepest events (about 500 km). Variations to such a stress pattern have been found, possibly related to mantle dynamics and the complex composition of the subducting structure. The diffused state of down-dip compression suggests that the Tyrrhenian subduction has already evolved to the point where the lower end of the slab has reached high-strength mantle materials, the load of the excess mass is entirely supported from below and most of the subducted slab is under compression. In agreement with the lack of large, shallow thrusting events in the immersion zone, the findings of the present study appear to agree well with geodynamic models assuming a passive subduction process with eastward roll-back of the Ionian lithosphere in the study area. In this context, the depth-decrease of the slab dip may also find a reasonable explanation.     

永清MS4.3地震和廊坊MS3.0地震的发震构造研究

基于首都圈数字地震台网的宽频带资料,首先采用CAP方法确定了永清M_S4.3地震和廊坊M_S3.0地震的震源机制解:永清地震节面Ⅰ的走向、倾角和滑动角分别为52°,62°和?140°,节面Ⅱ的走向、倾角和滑动角分别为300°,55°和?35°;廊坊地震节面I的走向、倾角和滑动角分别为48°,57°和?147°,节面Ⅱ的走向、倾角和滑动角分别为299°,63°和?38°。两次地震的震源机制解较为一致,推测它们可能具有相同的发震断层。利用近震转换波获得两次地震的震源深度,分别为19 km和13 km。利用双差法对两次地震的主余震进行重新定位,结果显示:两个地震序列的震中均呈NE向分布,余震震源深度均浅于主震震源深度,震源深度分别集中在17—20 km和12—13 km范围内,两个序列的短轴剖面揭示了震源分布均呈现倾向SE,倾角陡立的特点。将地震序列的分布与震源机制解的结果进行对比,认为两个序列的水平展布方向与其对应的主震震源机制解中节面Ⅰ的走向比较接近,深度分布的高倾角特征也与节面Ⅰ比较相似,因此认为发震断层面均为节面Ⅰ。通过将震源机制解中节面Ⅰ的参数和地震序列的分布与区域活动断层的产状性质进行比较,取得了一些关于发震构造和地震成因的重要认识:① 永清M_S4.3地震和廊坊M_S3.0地震的发震构造不是上地壳的先存正断裂?河西务断裂,不排除与中下地壳的新生构造或深大断裂有关;② 永清、廊坊地震发生在13—19 km深度上,结合地壳结构、断裂构造以及区域流变结构等资料,推测该深度范围可能是廊固凹陷的壳内脆性?韧性转换区域,是地震孕育和发生的有利构造部位。      

The seismicity and tectonic stress field characteristics of the Longmenshan fault zone before the Wenchuan MS8.0 earthquake

The seismicity of Longmenshan fault zone and its vicinities before the 12 May 2008 Wenchuan M_S8.0 earthquake is studied. Based on the digital seismic waveform data observed from regional seismic networks and mobile stations, the focal mechanism solutions are determined. Our analysis results show that the seismicities of Longmenshan fault zone before the 12 May 2008 Wenchuan earthquake were in stable state. No obvious phenomena of seismic activity intensifying appeared. According to focal mechanism solutions of some small earthquakes before the 12 May 2008 Wenchuan earthquake, the direction of principal compressive stress P-axis is WNW-ESE. The two hypocenter fault planes are NE-striking and NW-striking. The plane of NE direction is among N50°?70°E, the dip angles of fault planes are 60°?70° and it is very steep. The faultings of most earthquakes are dominantly characterized by dip-slip reverse and small part of faultings present strike-slip. The azimuths of principal compressive stress, the strikes of source fault planes and the dislocation types calculated from some small earthquakes before the 12 May 2008 Wenchuan earthquake are in accordance with that of the main shock. The average stress field of micro-rupture along the Longmenshan fault zone before the great earthquake is also consistent with that calculated from main shock. Zipingpu dam is located in the east side 20 km from the initial rupture area of the 12 May 2008 Wenchuan earthquake. The activity increment of small earthquakes in the Zipingpu dam is in the period of water discharging. The source parameter results of the small earthquakes which occurred near the initial rupture area of the 12 May 2008 Wenchuan earthquake indicate that the focal depths are 5 to 14 km and the source parameters are identical with that of earthquake.     

Depth distribution of stresses in the Hokkaido Wadati-Benioff zone as deduced by inversion of earthquake focal mechanisms

This study is based on the detailed geometry of the Hokkaido Wadati-Benioff zone and the paleosubduction zone as delineated by Hanus and Vanek (1984). The used data includes 217 CMT Harvard solutions for earthquakes, which belong to the Wadati-Benioff zone and 13 for the paleosubduction zone. The inverse technique by Gephart and Forsyth (1984) was incorporated for determining the best fit principal stress directions σ₁, σ₂, σ₃ and the ratio (R=σ₂−σ₁/σ₃−σ₁) for 20 km depth intervals in the Wadati-Benioff zone and for the paleosubduction zone considered as a single body. In almost all the considered depth layers, the maximum compressive stress σ₁ is normal to the strike of the slab and dips less than 25°, indicating the NW-SE convergence between the Pacific and Eurasian lithospheric plates. Exceptions are in the depth layer 81–120 km, the paleosubduction zone with steeply dipping along-strike σ₁, and the lower part of the subduction zone (161–220 km) where σ₁ is almost horizontal and of E trend. The minimum compressive stress σ₃ is mostly along-strike and of a different dip with the exception of the 21–60 km layer wher they are down-dipping. The results obtained for the depth ranges 0–20 km, 81–100 km, 121–160 km, and the paleosubduction zone indicate heterogeneous stress fields. These results show that the slab pull and the mantle resistance, acting on the slab edge, are not the main forces which control the contemporary plate tectonics in the Hokkaido region. Along-strike compression at depths 81–120 km and along-strike extension at 0–20 and 61–220 km are involved in the slab dynamics. These can be related to horizontal bending of the subducting Pacific plate.     

An unusual zone of seismic coupling in the Bonin arc: The 1972 Hachijo-Oki earthquakes and related seismicity

The 1972 February and December Hachijo-Oki earthquakes (M _s=7.3 and 7.4), in the northernmost part of the Izu-Bonin subduction zone, are the only major events (M _s>7.0) in the Bonin arc for the past 80 years. Relocation of the hypocenters, using one smaller event having a wellconstrained focal depth as a master event, shows that the depth of the February event is 10 km shallower than that of the December event. We have determined the rupture process for both events by minimizing the error in waveform between observed and synthetic seismograms. Although the number of available stations are limited, the depth range of the major energy release for the December event extends deeper than for the February one. The rupture propagated up-dip for both events. It is likely that the rupture zone of the two events overlapped, and that the December event ruptured the deeper part. This suggestion is consistent with the observation that the aftershock zones of both events overlap with that of the December event shifted landward. The waveforms of the December event have a smaller high frequency component than those of the February event, suggesting that the stress at the thrust zone became more uniform or reduced after the February event.No thrust type smaller event occurred near the rupture zone. Instead, theP-axes of smaller events are parallel to the dip of the slab and theirT-axes dip to the southwest. Focal depths of these events estimated byP-wave forward modeling are generally between 40–50 km and located beneath the thrust zone. We thus interpret them as the events within the Pacific slab near the zone ruptured by the two major events. The stress concentration around the rupture zone of the major events is suggested to have triggered these slab events. After the occurrence of the large events, the slab events are concentrated near the deeper portion of the rupture zone. These events may have been caused by the loading of the down-dip compressional stress near the down-dip end of the rupture zone due to the rupture. The occurrence of the doublet of large earthquakes and a number of down-dip compressional events beneath their rupture zones in a shallow portion of the subducting slab indicates an unusual zone of seismic coupling in the Bonin arc, most of which is seismically quiescent.     

汶川余震震源机制变化的原因

2008年5月12日,汶川M8.0主震是一个主压应力轴NW-SE的逆倾滑型的地震,而之后的余震震源机制解有的与主震震源机制解一致,有的发生了明显变化,由南至北逐步变成走滑型地震.主震和同主震震源机制解一致的部分余震,在构造应力场直接作用下,龙门山推复体向四川地块逆冲,致使在逆断层的上下盘之间的断层面上产生粘滑而发生逆倾...     

2018年新疆伽师MS5.5地震的发震构造初探

运用CAP方法反演2018年9月4日新疆伽师M_S5.5地震及M_S≥3.0余震的震源机制解,计算得出伽师M_S5.5地震的震源机制解为:节面Ⅰ:走向48°,倾角83°,滑动角3°;节面Ⅱ:走向318°,倾角87°,滑动角173°;主压应力P轴方位角为3°,倾角为3°,主张应力T轴方位角273°,倾角为7°;矩震级为MW5.3。使用双差定位法对主震及余震共计129个M_S≥1.5地震进行重新定位,并对震源机制解和重定位结果进行综合分析,发现此次重定位地震结果与CAP方法反演结果的展布方向一致,地震集中分布在NEE向,因此认为节面I是此次地震的主破裂面;重定位后NS、EW和UD方向的平均相对误差分别为0.25、0.23及0.09 km,平均走时残差为0.026 s,震源深度集中分布在5~15 km。此次地震及其余震附近地表无明显的断层出露,所以初步判定2018年新疆伽师M_S5.5地震可能受控于柯坪断裂带附近的隐伏断裂。     

2021年吴忠—灵武ML3.6震群重新定位及震源机制研究

2021年7月18日—8月7日,宁夏吴忠—灵武地区发生M_L3.6显著震群活动。本文利用多阶段定位方法对该震群进行了重新定位,并根据gCAP方法反演了2021年7月20日灵武M_L3.6地震的震源机制及震源矩心深度,采用Snoke方法计算了震群中3次M_L3.0以上地震的震源机制,测定了同一地震多个震源机制的中心解。结果表明,该震群中最大的地震即7月20日02时40分M_L3.6地震的震源机制为节面Ⅰ走向289°,倾角72°,滑动角?22°,节面Ⅱ走向26°,倾角69°,滑动角?161°,震源矩心深度为12 km,初始破裂深度为12.5 km;7月20日03时15分M_L3.2地震的震源机制为节面Ⅰ走向290°,倾角82°,滑动角?2°,节面Ⅱ走向20°,倾角88°,滑动角?172°,初始破裂深度为11.9 km;7月21日04时55分M_L3.1地震的震源机制为节面Ⅰ走向285°,倾角53°,滑动角2°,节面Ⅱ走向194°,倾角88°,滑动角143°,初始破裂深度为11.6 km,这些地震震源机制的主压应力轴主要为NE向。该震群序列的震源深度主要相对集中在7—15 km之间,其中M_L3.0以上地震的震源深度主要介于11—13 km,震源深度剖面显示震群相对集中的区域由深到浅大体呈现近似于陡立的展布。本文进一步研究发现区域应力场在灵武M_L3.6地震震源机制NNE向节面产生的相对剪应力为0.393,而在NWW向节面产生的相对剪应力为0.945。结合地质构造和已有断层资料初步分析认为,若NNE向的崇兴隐伏断裂为灵武M_L3.6地震的发震断层,则表明崇兴断裂可能是一条断裂薄弱带,地震破裂方式主要为右旋走滑;若NWW向的未知隐伏断裂为发震断层,则表明NWW向断裂可能为该地震在区域应力场下的剪应力相对最大释放节面,其破裂方式为左旋走滑。