Strain segmentation and lateral membrane deformation rate of the subducted Ryukyu slab

Abstract Focal mechanisms of intermediate-depth earthquakes within the subducted Ryukyu slab indicate a systematic pattern of down-dip extension in the northern part under Kyushu. These mechanisms switch rapidly around the Tokara channel to down-dip compression in the southern part toward Taiwan. This intriguing pattern of strain segmentation, together with the differences between the slab geometry of the northern and southern parts, has raised the question of whether stress segmentation, as defined possibly by a tear fault, might exist in the slab. However the Ryukyu trench has a concave oceanward shape in the northern segment while it is convex in the southern oceanward part toward Taiwan. The inflection zone is located around the Tokara channel. The concept of Gaussian curvature of a curved surface suggests that the along-arc variation of the geometric configuration of a subducted slab is related to the shape of the trench. This is in order to accommodate the lateral membrane deformation of the slab as the oceanic lithosphere subducts from a spherical shell to the geometry delineated by the Wadati-Benioff zone. The membrane deformation regime of the subducted Ryukyu slab and its relation with the trench geometry was examined by assuming that the subduction be modelled by the flow field of a thin viscous sheet. A projection operator was utilized to compute the membrane strain-rate tensor of an arbitrary non-Euclidean surface. Numerical experiments indicated that the northern segment of the slab was dominated by lateral compression and down-dip extension and the southern part by lateral extension and down-dip compression. This transition is sharply located near the Tokara channel. These patterns were compatible with what had been observed from studies of focal mechanisms, suggesting that the strain segmentation might be controlled, at least in part, by the lateral membrane deformation within the slab due to the trench shape in this subduction zone. The slab geometry was predicted by minimizing the integrated total dissipation power; this revealed distinct features that were consistent with observations. This implies that the slab geometry may also be affected by the membrane deformation in a systematic fashion.     

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.     

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

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

Seismicity pattern: an indicator of source region of volcanism at convergent plate margins

The results of detailed investigation into the geometry of distribution of earthquakes around and below the volcanoes Korovin, Cleveland, Makushin, Yake-Dake, Oshima, Lewotobi, Fuego, Sangay, Nisyros and Montagne Pelée at convergent plate margins are presented. The ISC hypocentral determinations for the period 1964-1999, based on data of global seismic network and relocated by Engdahl, van der Hilst and Buland, have been used.The aim of this study has been to contribute to the solution of the problem of location of source regions of primary magma for calc-alkaline volcanoes spatially and genetically related to the process of subduction. Several specific features of seismicity pattern were revealed in this context. (i) A clear occurrence of the intermediate-depth aseismic gap (IDAG) in the Wadati-Benioff zone (WBZ) below all investigated active volcanoes. We interpret this part of the subducted slab, which does not contain any teleseismically recorded earthquake with magnitude greater than 4.0, as a partially melted domain of oceanic lithosphere and as a possible source of primary magma for calc-alkaline volcanoes. (ii) A set of earthquakes in the shape of a seismically active column (SAC) seems to exists in the continental wedge below volcanoes Korovin, Makushin and Sangay. The seismically active columns probably reach from the Earth surface down to the aseismic gap in the Wadati-Benioff zone. This points to the possibility that the upper mantle overlying the subducted slab does not contain large melted domains, displays an intense fracturing and is not likely to represent the site of magma generation. (iii) In the continental wedge below the volcanoes Cleveland, Fuego, Nisyros, Yake-Dake, Oshima and Lewotobi, shallow seismicity occurs down to the depth of 50 km. The domain without any earthquakes between the shallow seismically active column and the aseismic gap in the Wadati-Benioff zone in the depth range of 50-100 km does not exclude the melting of the mantle also above the slab. (iv) Any earthquake does not exist in the lithospheric wedge below the volcano Montagne Pelée. The source of primary magma could be located in the subducted slab as well as in the overlying mantle wedge. (v) Frequent aftershock sequences accompanying stronger earthquakes in the seismically active columns indicate high fracturing of the wedge below active volcanoes. (vi) The elongated shape of clusters of epicentres of earthquakes of seismically active columns, as well as stable parameters of the available fault plane solutions, seem to reflect the existence of dominant deeply rooted fracture zones below volcanoes. These facts also favour the location of primary magma in the subducting slab rather than in the overlying wedge.We suppose that melts advancing from the slab toward the Earth surface may trigger the observed earthquakes in the continental wedge that is critically pre-stressed by the process of subduction. However, for definitive conclusions it will be necessary to explain the occurrence of earthquake clusters below some volcanoes and the lack of seismicity below others, taking into account the uncertainty of focal depth determination from global seismological data in some regions.     

上地幔俯冲板块的动力学过程:数值模拟

大洋板块俯冲到地幔转换带,进而可形成不同的形态:板块可以停滞在660km不连续面,抑或穿过地幔转换带进入下地幔.这些不同的俯冲模式可进一步影响到海沟的运动.为更好地理解上地幔中俯冲板片的变形行为以及俯冲过程与海沟运动之间的关系,本文通过建立一系列高精度二维热-力学自由俯冲的数值模型,揭示了俯冲板块在上地幔中的变形方式及其与地幔转换带之间的相互作用过程.模拟结果显示,在俯冲板块与地幔转换带的相互作用过程中,其动力学过程可以分为以海沟后撤主导、海沟前进主导以及稳定型海沟等三种主要动力学类型.对于年龄较老,厚度较大的俯冲板块容易形成海沟后撤型俯冲,俯冲板块停滞在660km不连续面.相反,年龄较小,塑性强度较小的板块容易形成海沟前进型俯冲,俯冲板块穿越660km不连续面.     

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.     

汤加—克马德克俯冲带现今非均匀应力场特征及其动力学意义

汤加—克马德克俯冲带是太平洋板块向澳大利亚板块俯冲碰撞的动力作用区,是全球俯冲带动力学研究的热点区域.本研究基于EHB地震目录,对汤加—克马德克俯冲带(18.5°S—28.5°S)区域进行平面拟合,得到该范围内俯冲带走向约为196°,倾角约为48°;利用该俯冲带研究区域内Global CMT目录,对不同位置、不同深度进行区域应力张量反演,得到汤加—克马德克俯冲带研究区内精细的应力图像.结果显示:(1)俯冲带浅部(60~300km)应力结构非均匀特征明显,主应力轴倾伏角变化多样,并且最大主压应力轴方位在24°S左右发生明显偏转,我们推测这可能与洋底构造路易斯维尔海链俯冲有关;(2)中部(300~500km)最大主压、主张应力轴由北向南逐渐发生偏转,这可能与由北向南流动的地幔流对俯冲板片产生推挤作用有关,并且这种推挤作用向南逐渐减弱;(3)深部(500~700km)最大主压应力轴沿俯冲方向分布;(4)本文的结果还发现了主俯冲带深部西侧"偏移"板片与主俯冲带应力结构不同,表明"偏移"板片与主俯冲带是分离的.     

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.     

基于重力数据特征的东北地区深部动力机制解析

中国东北地区处于古亚洲构造域、蒙古—鄂霍茨克构造域和环太平洋构造域叠加作用最为显著的地区,是地学研究的热点区域.为了探析欧亚大陆下西太平洋板片的俯冲形态以及揭示该区域深部地球动力学机制,利用卫星重力数据通过预处理共轭梯度快速密度反演算法获得了包含东北地区、华北部分地区及日本海海域在内的研究区域上地幔三维密度结构,结合天然地震三维层析成像结果共同揭示太平洋板片的俯冲形态和深部动力机制.俯冲的太平洋板片在日本海沟处呈高密度异常,向西横向持续扩张,深度方向上有逐渐增加趋势.不连续的高密度体俯冲至地幔转换带(410～660km)后继续水平西向俯冲,部分滞留板片可能进入下地幔;在大兴安岭断裂带下面转换带中同样发现水平分布的高密度体,推断是大兴安岭断裂带下方地幔岩石圈拆沉的结果,横向不均匀分布的俯冲板片边缘已抵至大兴安岭造山带附近,这对于研究东北地区深部动力学机制具有重要的意义.     

The effect of the descending lithosphere beneath the Tonga island arc on P-wave travel-time residuals at the Warramunga Seismic Array

P-wave travel-time residuals at the Warramunga Seismic Array (WRA) in the Northern Territory, Australia, have been studied from 49 earthquakes with epicenters south of 19°S in the Fiji-Tonga region. Focal depths are between 42 and 679 km as determined from pP-P. Using the Jeffreys-Bullen and the Herrin travel-time tables the epicentral parameters have been redetermined by considering only “normal” seismic stations in the location procedure. These are those stations where P-wave travel times are probably not affected by lateral heterogeneities caused by the lithosphere descending beneath the Tonga trench. Epicenters of deep earthquakes below 300 km have been relocated by using stations at Δ > 25° only. Epicenters from shallower-depth earthquakes have been recalculated without using stations between 35 < Δ < 75° epicentral distance. In both cases focal depths were determined from pP-P times. The resulting pattern of P-residuals at WRA does not show any significant change with depth below 350 km. The residuals become more negative for shallower earthquakes above about 250 km. P-waves to WRA are advanced by approximately 2 s compared with those from deep earthquakes. The results do not essentially differ for the two different travel-time tables used. The observations can be interpreted by P-wave velocities that are higher in the sinking slab down to 350–400 km by 5±2% than in both the Jeffreys-Bullen and Herrin models. Without considering possible elevations of phase boundaries this estimate yields a temperature contrast of 1000±450°C between slab and normal mantle material in this depth range.     

Gravity and seismicity analyses of the El Chichon volcano,Chiapas, Mexico

Subduction induced compressional stress, lateral displacement faults probably linked to the NW extension of the Polochic-Motagua fault system, and local magma emplacement create a complex geological setting for the El Chichon volcano.Gravity analysis reveals that the principal structures in the area correspond to synclines and anticlines. Downward continuation of the gravity field suggests the presence of magma emplacement SE of the volcano. This magmatic emplacement appears to be the principal cause of the truncation of the syncline in which the volcano lies.Analyses of earthquakes occurring after eruptions show that the vertical distribution of events is confined to a roughly cylindrical zone which correlates well with the region of magmatic emplacement inferred from the gravity analysis.Even though more than 350 km separate the volcano from the Middle American trench, the compressional stress regime reported for the area may correlate with subduction stress (N30°E). However, the strike direction for the principal faults in the zone, N60°W, permits us to postulate that the NW extension of the Polochic-Motagua fault system is implicated in the local tectonics of the El Chichon area.     

Marine seismological observations in the Central Kurile segment before catastrophic earthquakes on November, 2006 (<Emphasis Type="Italic">M</Emphasis> = 8.3) and January, 2007 (<Emphasis Type="Italic">M</Emphasis> = 8.1)

The results of detailed seismological observations with bottom seismographs in the Central Kurile segment in August-September, 2006 are discussed. The system of six bottom seismographs was placed on the island slope of the Kurile deep-sea trench southeast of Urup Island and southwest of the Bussol Strait. Over 230 earthquakes with M _LH = 0.5–5.5 were registered in the area with a radius of 150 km around the center of the observation system at depths up to 300 km during 16 days. Records of 80 earthquakes with hypocenters in the earth crust (h = 0–30 km) beneath the island slope of the Kurile deep-sea trench were first obtained by bottom seismographs. These data are inconsistent with previous concepts of aseismicity of this zone. The discovery of the unique morphological structure of the Benioff zone beneath the central Kurile Arc represents the most important result of detailed seismological observations. The zone consists of an inner seismoactive subzone, which is located beneath the island slope of the arc at depths of 15–210 km, being characterized by an angle of incline of 50° under the latter and crosses the ocean bottom approximately 80 km away from the trench axis, and outer low-activity subzone. The latter is traceable beyond the trench almost parallel to the inner zone beginning from a depth of 50 km below the sea bottom up to a depth of approximately 300 km. Due to the slightly lower incline (∼45°) of the outer subzone, both subzones gradually converge downward. The integral thickness of the Benioff zone varies from 150 km in its upper part to 125 km at depths of 210–260 km. The medium sandwiched between these subzones is practically aseismic. The reality of this defined structure is confirmed by the distribution of aftershocks of the earthquake that occurred on November 15, 2006 (M = 8.3). These seismic events served as foreshocks for the subsequent strong earthquake of January 13, 2007 (M = 8.1) with the hypocenter located beyond the trench under the ocean bottom. Such a structure of this zone within the central Kurile Arc segment is unique, having no analogues either in the flanks of the Kurile-Kamchatka Arc or other arcs. The results of detailed seismological observations obtained two months before the first of the catastrophic Central Kurile earthquakes appeared to be typical for the period of foreshocks (the lower seismic activity of the Simushir block, which hosted the hypocenter of the earthquake that occurred on November 15, 2006, particularly at depths of 0–50 km, the gentler incline of the recurrence plot, and other features).     

Landward flow in the upper mantle: effects of the heat sink and viscous coupling of the sinking slab

Regional studies landward of the trench show that gravity anomalies are low compared with those expected from the cold sinking slab and that residual depth anomalies are strongly negative at distances up to 1200 km from the trench. It is shown that flow models taking simultaneously account of the viscous coupling and of the cooling effect of the sinking slab are able to explain these facts: the viscous coupling being related only to observed data down to 400 km from the trench, while the cooling effect mainly explains data at larger distances from the subduction zone.     

汕头－吕宋岛岩石圈速度结构与震源构造──台湾浅滩7．3级地震构造之一

汕头－吕宋岛岩石圈速度结构剖面，划分出华南陆缘古生代陆壳、陆架区晚古生代－中生代陆壳、陆坡带中生代－早第三纪过渡壳、新生代南海海盆洋壳及吕宋岛中生代－新生代岛弧陆壳与东吕宋海槽洋壳等地壳构造组分，并确定了上述地壳构造之间的边界断裂构造及其性质。结合地震震源分布及机制，初步确定了华南陆架盆岭构造带北、南两侧地震构造的控震构造与发震构造性质及其震源力学特征；１）指出１９９４年９月１６日台湾浅滩７．３级地震属于板缘壳幔地震及造成一千公里有感范围的原因；２）马尼拉海沟的海底地堑构造与南海海盆岩石圈地幔上隆是马尼拉海沟俯冲带震源显示正断层性质的原因，且为被动的或转换俯冲带；３）东吕宋海槽仍属于菲律宾海俯冲带性质；吕宋岛东西两侧俯冲带岩石圈板片震源深度的准三层分布，可能表明俯冲带岩石圈板片存在相应的低速滑移层。     

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.     

Thermal Structure of the Ionian Slab

We propose a thermal model of the subducting Ionian microplate. The slab sinks in an isothermal mantle, and for the boundary conditions we take into account the relation between the maximum depth of seismicity and the thermal parameter L_th of the slab, which is a product of the age of the subducted lithosphere and the vertical component of the convergence rate. The surface heat-flux dataset of the Ionian Sea is reviewed, and a convective geotherm is calculated in its undeformed part for a surface heat flux of 42 mW m^–2, an adiabatic gradient of 0.6 mK m^–1, a mantle kinematic viscosity of 10¹⁷ m² s^–1 and an asthenosphere potential temperature of 1300°C. The calculated temperature-depth distribution compared to the mantle melting temperature indicates the decoupling limit between lithosphere and asthenosphere occurs at a depth of 105 km and a temperature of 1260°C. A 70–km thick mechanical boundary layer is found. By considering that the maximum depth of the seismic events within the slab is 600 km, a L_th of 4725 km is inferred. For a subduction rate equal to the spreading rate, the corresponding assimilation and cooling times of the microplate are about 7 and 90 Myr, respectively. The thermal model assumes that the mantle flow above the slab is parallel and equal to the subducting plate velocity of 6 cm yr^–1, and ignores the heat conduction down the slab dip. The critical temperature, above which the subduced lithosphere cannot sustain the stress necessary to produce seismicity, is determined from the thermal conditions governing the rheology of the plate. The minimum potential temperature at the depth of the deepest earthquake in the slab is 730°C.     

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.     

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.     

马尼拉俯冲带的地震层析成像研究

基于国际地震中心的P波走时数据和层析成像反演方法,获得了具有较高分辨率的马尼拉俯冲带的深部速度模型.结果表明,（1）高速的南海俯冲板片沿马尼拉俯冲带的俯冲形态随纬度发生变化,在14°N和16°N之间,板片俯冲角度较大,俯冲深度可达400~500 km,在17°N附近,俯冲板片角度和深度较南部变小,而在18°N附近,俯冲板片以近垂直角度俯冲到地幔转换带;（2）17°N和18°N之间俯冲角度的变化意味着南海板片发生了撕裂;（3）在14°N附近,南海板片由300 km以上的近垂直俯冲转为200~300 km深度的近水平展布,与震源分布存在较大的差异,表明南海板片发生了撕裂,并且导致410 km间断面抬升.根据成像结果计算的不同位置南海板片的俯冲长度和时间表明,南海板片俯冲之前的面积为现今面积的两倍,14°N最先开始发生俯冲,并由南向北扩展.     

Regional focal mechanisms for earthquakes in the Aegean area

The distribution of the focal mechanisms of the shallow and intermediate depth (h>40 km) earthquakes of the Aegean and the surrounding area is discussed. The data consist of all events of the period 1963–1986 for the shallow, and 1961–1985 for the intermediate depth earthquakes, withM _s 5.5. For this purpose, all published fault plane solutions for each event have been collected, reproduced, carefully checked and if possible improved accordingly. The distribution of the focal mechanisms of the earthquakes in the Aegean declares the existence of thrust faulting following the coastline of southern Yugoslavia, Albania and western Greece extending up to the island of Cephalonia. This zone of compression is due to the collision between two continental lithospheres (Apulian-Eurasian). The subduction of the African lithosphere under the Aegean results in the occurrence of thrust faulting along the convex side of the Hellenic arc. These two zones of compression are connected via strike-slip faulting observed at the area of Cephalonia island. TheP axis along the convex side of the arc keeps approximately the same strike throughout the arc (210° NNE-SSW) and plunges with a mean angle of 24° to southwest. The broad mainland of Greece as well as western Turkey are dominated by normal faulting with theT axis striking almost NS (with a trend of 174° for Greece and 180° for western Turkey). The intermediate depth seismicity is distributed into two segments of the Benioff zone. In the shallower part of the Benioff zone, which is found directly beneath the inner slope of the sedimentary arc of the Hellenic arc, earthquakes with depths in the range 40–100 km are distributed. The dip angle of the Benioff zone in this area is found equal to 23°. This part of the Benioff zone is coupled with the seismic zone of shallow earthquakes along the arc and it is here that the greatest earthquakes have been observed (M _s 8.0). The deeper part (inner) of the Benioff zone, where the earthquakes with depths in the range 100–180 km are distributed, dips with a mean angle of 38° below the volcanic arc of southern Aegean.