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

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

琉球海沟的构造和运动特征

琉球海沟是欧亚板块与菲律宾海板块之间的边界，海沟西坡是大陆性质的琉球岛弧，东坡是大洋性质的菲律宾海地壳．海底地震反射探测和地震震源定位表明，菲律宾海板块沿海沟向琉球岛弧下俯冲，俯冲角度与深度沿海沟走向变化．有证据显示，由于俯冲板前缘的横向移动，海沟和岛弧正朝大洋方向后退，弧后盆地-冲绳海槽发生拉张变形．最近一次的海沟后退与冲绳海槽扩张可能是从上新世末(2Ma前)开始的，岛弧的后退移动和弧后拉张在南部与海沟走向垂直，在中部和北部与海沟走向斜交，总体上向南的运动分量占优势、与海沟后退相关的弧后拉张集中在冲绳海槽，没有证据表明对其西侧的中国东海陆架盆地产生影响、海沟后退的原因可能与俯冲板的动力不平衡以及它与周围地幔的相互作用有关．     

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

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

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.     

Crustal structure and local seismicity in Colombia

Using P-wave travel time data from local seismicity, the crustal structure ofthe central and southern part of Colombia was determined. A very stableand narrow range of possible velocity models for the region was obtainedusing travel time inversion. This range of models was tested with earthquakelocations to select the best velocity model. The 1D velocity modelproposed has five layers over a halfspace, with interfaces at depths of 4,25, 32, 40 and 100 km and P-wave velocities of 4.8, 6.6, 7.0, 8.0, 8.1and 8.2 km/sec, respectively. According to this model the Moho lies at32 km depth on average. For P-waves, the station corrections range from–0.62 to 0.44 sec and for S-wave they range from –1.17 to 0.62 sec.These low variations in station residuals indicate small lateral velocitychanges and therefore the velocity model found should be well suited forearthquake locations and future starting model for 3D tomography studies.Using this new velocity model, the local earthquakes were relocated. Theshallow seismicity, < 30 km, clearly shows the borders betweentectonic plates and also the main fault systems in the region. The deepseismicity, > 80 km, shows two subduction zones in the country: theCauca subduction zone with a strike of N120^°E, dip of 35^°and thickness of 35 km, and the Bucaramanga subduction zone which has,for the northern part, a strike of N103^°E, dip of 27^° andthickness undetermined and, for the southern part, a strike ofN115^°E, dip of 40^° and thickness of 20 km. Based ondifferences of thickness of brittle crust in the subducted slab and spatialdistribution of the seismicity, the Cauca and Bucaramanga subduction zonesseem to represent independent processes. The Cauca subduction seems tobe connected to the process of the Nazca plate being subducted under theNorth Andes Block. In the Bucaramanga subduction zone, the transitionbetween southern and northern parts and changes in geometry of the slabseem to be gradual and there is no evidence of a tear in the slab, howeverthe local seismicity does not allow us to determine which plate or plates arebeing subducted. The Bucaramanga nest appears to be included into thesubducted slab.     

Two-phase opening of Andaman Sea: a new seismotectonic insight

High-resolution reconstruction of Benioff zone depth–dip angle trajectory for Burma–Java subduction margin between 2° and 17°N Lat. reveals two major episodes of plate geometry change expressed as abrupt deviation in subduction angle. Estimation of effective rate of subduction in different time slices (and then length of subducted slab) allowed drawing of isochrones in Ma interval through these trajectories for the time period 5–12 Ma. With these isochrones, the deformation events on the subducting Indian plate are constrained in time as of 4–5 and 11 Ma old. This well-constrained time connotation offered scope for the correlation of slab deformation events with the well-established two-phase opening history of the Andaman Sea. While the 11 Ma event recorded from southern part of the study area is correlated with early stretching and rifting phase, the 4–5 Ma event is interpreted as major forcing behind the spreading phase of the Andaman Sea. Systematic spatio-temporal evaluation of Indian plate obliquity on the Andaman Sea evolution shows its definite control on the early rifting phase, initiated towards south near northwest Sumatra. The much young spreading phase recorded towards north of 7° Lat. is possibly the result of late Miocene–Pliocene trench retreat and follow-up transcurrent movement (along Sagaing and Sumatran fault system) with NW–SE pull-apart extension.Nonconformity between plate shape and subduction margin geometry is interpreted as the causative force behind Mid-Miocene intraplate extension and tearing. Enhanced stretching in the overriding plate consequently caused active forearc subsidence, recorded all along this plate margin. Initial phase of the Andaman Sea opening presumably remains concealed in this early–middle Miocene forearc subsidence history. The late Miocene–Pliocene pull-apart opening and spreading was possibly initiated near the western part of the Mergui–Sumatra region and propagated northward in subsequent period. A temporary halt in rifting at this pull-apart stage and northeastward veering of the Andaman Sea Ridge (ASR) are related with uplifting of oceanic crust in post-middle Miocene time in form of Alcock and Sewell seamounts, lying symmetrically north and south of this spreading ridge.     

Joint bulk-sound and shear tomography for Western Pacific subduction zones

Detailed regional body wave tomographic inversion of the Western Pacific region has been performed using P and S travel times from common sources and receivers, with a joint inversion in terms of bulk-sound and shear wave-speed variations in the mantle. This technique allows the separation of the influence of bulk and shear moduli, and hence a more direct comparison with mineral physics information. The study region is parameterized with cells of side 0.5° to 2° and 19 layers to a depth of 1500 km, while the rest of the mantle was parameterized with 5×5° cells with 16 layers between the surface and the core–mantle boundary. A simultaneous inversion is made for regional and global structures to minimize the influence of surrounding structures on the regional image. A nested iterative inversion scheme is employed with local linearization and three-dimensional ray tracing through the successive model updates. The results of the regional tomographic inversion reveal the penetration of a subducted slab below the 660 km discontinuity at the Kurile–Kamchatka trench, while flattening of slabs above this depth is observed in the Japan and Izu–Bonin subduction zones on both the bulk-sound and shear wave-speed images. The penetration of a subducted slab down to a depth of at least 1200 km is seen below the southern part of the Bonin trench, Mariana, Philippine, and Java subduction zones. Fast shear wave-speed perturbations associated with the subducted slabs, down to the 410 km transition zone, are larger than the comparable bulk-sound perturbations for all these subduction zones except the Philippines. The bulk-sound signature for the subducted slab is more pronounced than for shear in the Philippines, Talaud, New Guinea, Solomon, and Tonga subduction zones, where penetration of the slab into the middle mantle is observed. Variation in the amplitude ratio between bulk-sound and shear wave-speed anomalies correlates well with the subduction parameters of the descending slab. Slabs younger than 90 Ma at the trench show bulk-sound dominance in the upper mantle, while older slabs have a stronger shear wave-speed signature. Spreading of the fast shear wave-speed zone between 800 and 1000 km is observed in the areas of deep subducted slab penetration, but has no comparable expression in the bulk-sound images. This high-velocity feature may reflect physical or chemical disequilibria introduced to the lower mantle by subducted slabs.     

The influence of trench migration on slab penetration into the lower mantle

A two-dimensional numerical convection model in cartesian geometry is used to study the influence of trench migration on the ability of subducted slabs to penetrate an endothermic phase boundary at 660 km depth. The transient subduction history of an oceanic plate is modelled by imposing plate and trench motion at the surface. The viscosity depends on temperature and depth. A variety of styles of slab behaviour is found, depending predominantly on the trench velocity. When trench retreat is faster than 2–4 cm/a, the descending slab flattens above the phase boundary. At slower rates it penetrates straight into the lower mantle, although flattening in the transition zone may occur later, leading to a complex slab morphology. The slab can buckle, independent of whether it penetrates or not, especially when there is a localised increase in viscosity at the phase boundary. Flattened slabs are only temporarily arrested in the transition zone and sink ultimately into the lower mantle. The results offer a framework for understanding the variety in slab geometry revealed by seismic tomography.     

Deep geodynamics of far field interconti- nental back-arc extension: Formation of Cenozoic volcanoes in northeastern China

Introduction Northeastem China has the most strong Cenozoic volcanism in China (Liu, 1999), where dis-tributes more than 500 Cenozoic volcanoes, including sleeping volcanoes of Tianchi Lake (Celes-tial Pond) of Changbai Mountain, and Wudalianchi (Five linked Lakes) (LIU, 1999). Vo lcano ofTianchi Lake of Changbai Mountain consists of basaltic rocks of shield-forming stage andtrachytes and pantellerites in cone-forming stage. It is suggested by study of REE, incompatibleelements a…     

Deep geodynamics of far field intercontinental back-arc extension: Formation of Cenozoic volcanoes in northeastern China

There are three cases of variation of trench location possible to occur during subduction: trench fixed, trench advancing, and trench retreating. Retreat of trench may lead to back-arc extension. The Pacific plate subducts at low angle beneath the Eurasia plate, tomographic results indicate that the subducted Pacific slab does not penetrate the 670 km discontinuity, instead, it is lying flat above the interface. The flattening occurred about 28 Ma ago. Geodynamic computation suggests: when the frontier of the subducted slab reaches the phase boundary of lower and upper mantle, it may be hindered and turn flat lying above the boundary, facilitates the retreat of trench and back-arc extension. Volcanism in northeastern China is likely a product of such retreat of subduction, far field back-arc extension, and melting due to reduce of pressure while mantle upwelling. Foundation item: National Natural Science Foundation of China (40234042 and 40174027).     

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.     

Cenozoic volcanism, stress gradient and back-arc opening in the North Island, New Zealand: Origin of Taupo-Rotorua Depression

Abstract Extensive subduction-related and intraplate volcanism characterize Cenozoic magmatism in the North Is., New Zealand. Volcanics in the central North Is., predominantly intermediate to felsic, form above the dipping seismic zone and show tectonic/geochemical features common to magmatism in most subduction zones. Basaltic volcanism in Northland, the northern part of the North Is., has chemical characteristics typical of intraplate magmatism and may be caused by the upwelling of asthenospheric materials from deeper parts of the mantle. The rifting just behind the present volcanic front (the Taupo-Rotorua Depression), which follows the trench ward migration of the volcanic front and the gradual steepening of the subducted slab, is also a feature of the North Is. A possible mechanism for the back-arc rifting in the area is injection of asthenospheric materials into the mantle wedge; this asthenospheric flow results from the mantle upwelling beneath Northland and pushes both the rigid fore-arc mantle wedge and the subducted slab trenchwards. This mechanism is also consistent with the stress fields in the North Is.: dilatation in Northland, northwest-southeast tension in the Taupo-Rotorua Depression, and the northeast-southwest compression in the fore-arc region.     

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.     

Effects of relative plate motion on the deep structure and penetration depth of slabs below the Izu-Bonin and Mariana island arcs

An increasing number of seismological studies indicate that slabs of subducted lithosphere penetrate the Earth's lower mantle below some island arcs but are deflected, or, rather, laid down, in the transition zone below others. Recent numerical simulations of mantle flow also advocate a hybrid form of mantle convection, with intermittent layering. We present a multi-disciplinary analysis of slab morphology and mantle dynamics in which we account explicitly for the history of subduction below specific island arcs in an attempt to understand what controls lateral variations in slab morphology and penetration depth. Central in our discussion are the Izu-Bonin and Mariana subduction zones. We argue that the differences in the tectonic evolution of these subduction zones—in particular the amount and rate of trench migration—can explain why the slab of subducted oceanic lithosphere seems to be (at least temporarily) stagnant in the Earth's transition zone below the Izu-Bonin arc but penetrates into the lower mantle below the Mariana arc. We briefly speculate on the applicability of our model of the temporal and spatial evolution of slab morphology to other subduction zones. Although further investigation is necessary, our tentative model shows the potential for interpreting seismic images of slab structure by accounting for the plate-tectonic history of the subduction zones involved. We therefore hope that the ideas outlined here will stimulate and direct new research initiatives.     

智利三联点南部扩张洋脊俯冲区域岩石层热结构的数值模拟

智利三联点(CTJ)位于纳兹卡板块、南极洲板块与南美板块的交界处,由南极洲—纳兹卡板块之间的智利洋脊俯冲到智利海沟而形成.巴塔哥尼亚板片窗的发展是智利洋脊长期扩张俯冲的结果之一.随着纳兹卡板块的不断东向俯冲,纳兹卡板块范围逐渐变小,CTJ同时向北移动.本文采用数值模拟方法,建立了关于洋脊海沟碰撞的简单二维模型,来研究智利三联点南部扩张洋脊俯冲区域岩石层的热结构.模拟结果表明,洋脊的位置、板块相对汇聚速度及上覆大陆板块的存在均对俯冲区域海洋板块的温度结构有着很大影响,并且大陆板块下方海洋板块温度变化最大的位置距洋脊的水平距离与洋脊到板片窗范围的水平距离两者之间具有较好的一致性.同时,当存在两两板块间的相对汇聚时,洋脊右侧大陆板块下表面的温度升高,俯冲带内海洋板块温度接近于地幔温度.纳兹卡板块以7.8 cm·a~(-1)的速度急速俯冲于南美板块之下的过程中,同时伴随着智利洋脊的持续扩张俯冲,在智利三联点南部,南美板块之下纳兹卡板块的温度因而可以更快地达到地幔软流层的约1300℃温度,并最终消亡于地幔之中.     

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.     

A present-day tectonic stress map for eastern Asia region

Introduction Tectonically the eastern Asia refers to the region bounded by the following three active tec-tonic zones: in the east the western Pacific subduction zone, including Japan trench, Ryukyu trench and Philippine trench; in the southwest the Himalaya continental collision zone and the Burma-arc-Java-trench subduction zone; in the northwest the Tianshan-Baikal continental defor-mation zone (Figure 6). In the world the eastern Asia is one of the regions with the strongest pre-sent-da…     

对日本俯冲带与IBM俯冲带俯冲特征的地球物理研究:来自重力与震源分布数据的启示

日本俯冲带与IBM俯冲带位于太平洋板块、菲律宾海板块和欧亚板块三者的交汇地带,是典型的"俯冲工厂"地区,具有重要的研究意义.本文利用震源分布资料与卫星重力数据对日本俯冲带与IBM俯冲带进行了研究.通过空间重力异常反映了俯冲带地区的区域构造形态,在此基础上基于艾利模式计算了均衡异常以反映地壳均衡特征.利用震源分布资料,分别从垂直俯冲带走向与沿俯冲带走向划定了横截剖面(cross-sections)进行了地震提取,讨论了俯冲带地区的Wadati-Benioff带形态特征,并借助于俯冲带地震等深线图直观描述了俯冲带的俯冲形态.在日本俯冲带与伊豆—小笠原俯冲带各选取了一条典型剖面进行了重力2.5D反演,研究了俯冲带地区的壳幔结构特征.研究结果表明,九州—帕劳海脊与IBM岛弧在均衡异常上存在差异,前者已逐渐趋向于地壳均衡.IBM的Wadati-Benioff带存在明显的南北差异,反映出伊豆—小笠原俯冲板片停留在了660km转换带中,而马里亚纳俯冲板片很可能垂直穿过了这一转换带,造成这种南北差异的原因与板块相对运动、岩石圈黏性和年龄差异以及俯冲板片的重力效应等因素有关.在IBM的中部和南部存在板片撕裂现象.日本俯冲带的俯冲洋壳密度随俯冲深度变化较小,洋幔存在一定程度的蛇纹岩化,地幔楔蛇纹岩化作用不典型,海沟处有一范围较小的含水畸变带;伊豆—小笠原俯冲带俯冲洋壳密度随深度增大而明显增大,洋幔蛇纹岩化程度较日本俯冲带低,地幔楔蛇纹岩化作用强烈,板块交汇处存在明显的蛇纹岩底辟.日本俯冲带与IBM俯冲带一线自北向南板片俯冲变陡,两侧板块耦合度降低,与俯冲带两侧的板块运动速率差异有关.     

南海东部现时地壳运动、震源机制及晚中新世以来的板块相互作用

南海东部的板块汇聚带是了解南海和菲律宾海晚中新世以来构造演化的一个重要窗口.针对这一区域地壳运动的研究,获得了该区内微块体的现时地壳运动特征及其动力机制,在此基础上结合俯冲板片形态和震源机制资料提出了该区晚中新世以来的板块汇聚作用特征.研究发现:晚中新世晚期,菲律宾海板块西缘在南、北部的西向运动均受到限制的情况下,位于中间区域的吕宋岛北部则由于其西侧相对自由而继续往NW方向运动.该过程中,菲律宾大断裂等汇聚带内部的走滑断裂对于协调不同块体之间的地壳运动速度差异有着重要作用.受南、北阻挡的影响,中间部分的西向运动速率呈现出中间大(吕宋岛北部)两端小的特征.因而北吕宋西侧的马尼拉海沟也以相对于南部更快的速率不断向NW迁移.然而俯冲的南海岩石圈受其下方SE向地幔流的影响,未能发生相应的俯冲板片后撤,而是在两板块之间的直接接触面形成强烈推挤并发生反向弯曲.结合这一板块作用特征认为,马尼拉海沟现今的构造形迹是在上述背景下北段多次向NW方向变迁形成的,而双火山弧的形成则主要是由菲律宾海板块在吕宋岛弧南、北部的西向运动速率差异引起的.