Seismic evidence for a metastable olivine wedge in the subducting Pacific slab under Japan Sea

We apply a forward-modeling approach to high-quality arrival time data from 23 deep earthquakes greater than 400 km depth to investigate the detailed structure of the subducting Pacific slab beneath the Japan Sea. Our results show that a finger-like anomaly exists within the subducting Pacific slab below 400 km depth, which has a P-wave velocity 5% lower than the surrounding slab velocity (or 3% lower than that of the normal mantle), suggesting the existence of a metastable olivine wedge (MOW) in the slab. The MOW top and bottom depths are 400 and 560 km, respectively. The MOW is estimated to be about 50 km wide at 400 km depth and close to the slab upper boundary. At 560 km depth the MOW is located at about 25 km below the slab upper boundary. Most of the deep earthquakes are located in the MOW. Our results favor transformational faulting as the mechanism for deep earthquakes.     

Subduction zone earthquakes and stress in slabs

Summary The pattern of seismicity as a function of depth in the world, and the orientation of stress axes of deep and intermediate earthquakes, are explained using viscous fluid models of subducting slabs, with a barrier in the mantle at 670 km. 670 km is the depth of a seismic discontinuity, and also the depth below which earthquakes do not occur. The barrier in the models can be a viscosity increase of an order of magnitude or more, or a chemical discontinuity where vertical velocity is zero. LongN versus depth, whereN is the number of earthquakes, shows (1) a linear decrease to about 250–300 km depth, (2) a minimum near that depth, and (3) an increase thereafter. Stress magnitude in a subducting slab versus depth, for a wide variety of models, shows the same pattern. Since there is some experimental evidence thatN is proportional toe , where is a constant and is the stress magnitude, the agreement is encouraging. In addition, the models predict down-dip compression in the slab at depths below 400 km. This has been observed in earlier studies of earthquake stress axes, and we have confirmed it via a survey of events occurring since 1977 which have been analysed by moment tensor inversion. At intermediate depths, the models predict an approximate but not precise state of down-dip tension when the slab is dipping. The observations do not show an unambiguous state of down-dip tension at intermediate depths, but in the majority of regions the state of stress is decidedly closer to down-dip tension than it is to down-dip compression. Chemical discontinuities above 670 km, or phase transitions with an elevation of the boundary in the slab, predict, when incorporated into the models, stress peaks which are not mirrored in the profile of seismicity versus depth. Models with an asthenosphere and mesosphere of appropriate viscosity can not only explain the state of stress observed in double Benioff zones, but also yield stress magnitude profiles consistent with observed seismicity. Models where a nonlinear rheology is used are qualitatively consistent with the linear models.     

2-D viscoelastic FEM simulation on stress state in the deep part of a subducted slab

Introduction It is found that there are some relationships between the thermal structures of subduction zones and the deep seismicity, while the mechanism relates the thermal structure and the deep seismicity is still unsure (Helffrich, Brodholt, 1991; Furukawa, 1994; Kirby, et al, 1996). From 1980s, geoscientists have constituted a series of numerical simulations on the stress states of subduction slabs. Based on the kinetic computation of Sung and Burns (1976a, b), Goto, et al (1983, 1987…     

俯冲带深部应力场的二维粘弹性有限元数值模拟

基于粘弹性平面应变有限元数值模拟方法，研究了俯冲带深部的应力场特征．当根据Karato等的研究结果给出俯冲板块的粘度结构时，橄榄石——尖晶石相变界面以下中心低粘区两侧出现应力集中区，而且其主压应力方向沿俯冲方向走向．其特征与已有的深源地震主压应力方向沿俯冲方向及深源地震有两个条带的地震观测结果相吻合．相变过渡区因矿物更小的颗粒粒度而导致的更小的等效粘度对俯冲带深部应力场的影响不大；存在和不存在亚稳态橄榄石楔的情况下，橄榄石——尖晶石相变过渡区附近都有最大剪应力的最大值出现，而其主压应力方向有垂直于橄榄石——尖晶石相变过渡区走向的趋势，尚与地震观测资料不符．      

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

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

On the origin of El Chichón volcano and subduction of Tehuantepec Ridge: A geodynamical perspective

The origin of El Chichón volcano is poorly understood, and we attempt in this study to demonstrate that the Tehuantepec Ridge (TR), a major tectonic discontinuity on the Cocos plate, plays a key role in determining the location of the volcano by enhancing the slab dehydration budget beneath it. Using marine magnetic anomalies we show that the upper mantle beneath TR undergoes strong serpentinization, carrying significant amounts of water into subduction. Another key aspect of the magnetic anomaly over southern Mexico is a long-wavelength (∼ 150 km) high amplitude (∼ 500 nT) magnetic anomaly located between the trench and the coast. Using a 2D joint magnetic-gravity forward model, constrained by the subduction P–T structure, slab geometry and seismicity, we find a highly magnetic and low-density source located at 40–80 km depth that we interpret as a partially serpentinized mantle wedge formed by fluids expelled from the subducting Cocos plate. Using phase diagrams for sediments, basalt and peridotite, and the thermal structure of the subduction zone beneath El Chichón we find that ∼ 40% of sediments and basalt dehydrate at depths corresponding with the location of the serpentinized mantle wedge, whereas the serpentinized root beneath TR strongly dehydrates (∼90%) at depths of 180-200 km comparable with the slab depths beneath El Chichón (200-220 km). We conclude that this strong deserpentinization pulse of mantle lithosphere beneath TR at great depths is responsible for the unusual location, singularity and, probably, the geochemically distinct signature (adakitic-like) of El Chichón volcano.     

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.     

亚稳态橄榄石相变与深源地震研究进展

由于岩石圈俯冲板块内部温度低, 在400km左右的深度橄榄石不发生相变, 低压相的橄榄石便以亚稳态的形式进入更深的深度。 高温高压实验研究表明, 亚稳态橄榄石相变机制以相变形成的超塑性透镜状尖晶石反裂纹破裂为深源地震的触发机制, 能很好的解释深源地震的特征。 亚稳态橄榄石在俯冲带中存在的范围,是制约亚稳态相变机制的重要因素。 然而亚稳态橄榄石在俯冲带存在的范围尚存争议, 不同模型给出结果相差较大, 文中给出一些解决这个问题的建议。     

Time dependence of negative buoyancy and the subduction of continental lithosphere

The time evolution of negative buoyancy of a subducting slab is modelled from the beginning of subduction under various kinematic conditions (dip angle and subduction velocity). The calculations take into account the thermal and density effects of the variations of the thermophysical parameters with temperature and pressure, and of phase transitions. The magnitude of the negative buoyancy increases during subduction of oceanic lithosphere, up to values in the (2–4) × 10¹³ N m⁻¹ range when the tip of the slab reaches a depth of 600–700 km. If continental material arrives at the trench and is subducted, the downward buoyancy decreases by an amount proportional to the volume of the subducted continental crust. Assuming that subduction stops when the buoyancy becomes zero, and that delamination of the continental crust or slab breakoff do not occur, the maximum downdip length of the subductable continental crust is estimated as a function of the dip angle, subduction velocity and geometry of the margin. In most cases, subduction of continental material down to depths of 100–250 km is possible, and continental subduction can continue for times up to 10–15 Ma if the velocity is low. These estimates are not significantly affected by the hypothetical occurrence of a metastable olivine wedge within the slab, and could be lower bounds if the lower continental crust is mafic and transforms to eclogite.     

The thermal structure of subduction zones constrained by seismic imaging: Implications for slab dehydration and wedge flow

Thermal models of subduction zones often base their slab–wedge geometry from seismicity at mantle depths and, consequently, cannot be used to evaluate the relationship between seismicity and structure. Here, high-resolution seismic observations from the recent Broadband Experiment Across the Alaska Range (BEAAR) constrain, in a rare instance, the subducting slab geometry and mantle wedge temperature independent of seismicity. Receiver functions reveal that the subducting crust descends less steeply than the Wadati-Benioff Zone. Attenuation tomography of the mantle wedge reveals a high Q and presumably cold region where the slab is less than 80 km deep. To understand these two observations, we generate thermal models that use the improved wedge geometry from receiver functions and that incorporate temperature- and strain-rate-dependent olivine rheology. These calculations show that seismicity within the subducting crust falls in a narrow belt of pressure–temperature conditions, illuminating an effective Clapeyron slope of 0.1 K/MPa at temperatures of 450–750 °C. These conditions typify the breakdown of high-pressure hydrous minerals such as lawsonite and suggest that a single set of dehydration reactions may trigger intermediate-depth seismicity. The models also require that the upper, cold nose of the mantle wedge be isolated from the main flow in the mantle wedge in order to sustain the cold temperatures inferred from the Q tomography. Possibly, sufficient mechanical decoupling occurs at the top of the downgoing slab along a localized shear zone to 80 km depth, considerably deeper than inferred from thrust zone seismicity.     

亚稳态橄榄石对俯冲带负浮力的影响及其动力学意义

利用有限元方法,计算了不同俯冲速度及热传导系数俯冲岩石层的负浮力.在h-lt;400km和h-gt;740km的深度范围,低温高密俯冲带的负浮力随深度单调增加.因为尖晶石相到后尖晶石相的相变有负的克兰帕龙斜率,俯冲带冷的物质在660km间断面以下不到100km的深度范围内仍以低密度的低压物质相存在.所以在该深度范围,俯冲带受到了周围地幔阻止其插入下地幔的浮力作用.在400km-lt;h-lt;660km深度范围,由于受橄榄石相变的影响,不同计算模型负浮力随深度的变化有明显的不同.对于可能的相变动力学模型,亚稳态橄榄石的存在使负浮力随深度的增加值减小.其作用是不利于俯冲带直接穿透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.     

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.     

Hugoniot equation of state of olivine and its geodynamic implications

Large olivine samples were hot-pressed synthesized for shock wave experiments. The shock wave experiments were carried out at pressure range between 11 and 42 GPa. Shock data on olivine sample yielded a linear relationship between shock wave velocity D and particle velocity u described by D=3.56(?0.13)+2.57(?0.12)u. The shock temperature is determined by an energy relationship which is approximately 790°C at pressure 28 GPa. Due to low temperature and short experimental duration, we suggest that no phase change occurred in our sample below 30 GPa and olivine persisted well beyond its equilibrium boundary in metastable phase. The densities of metastable olivine are in agreement with the results of static compression. At the depth shallower than 410 km, the densities of metastable olivine are higher than those of the PREM model, facilitating cold slab to sink into the mantle transition zone. However, in entire mantle transition zone, the shock densities are lower than those of the PREM model, hampering cold slab to flow across the "660 km" phase boundary.     

Crustal structure beneath Aso Caldera,Southwest Japan,as derived from receiver function analysis

Aso Volcano experienced a huge pyroclastic eruption 90 thousand years ago, and formed a large caldera (18 km × 25 km). In order to test the hypothesis of a magma body in the mid and lower crust that has been suggested geophysically and geochemically, we investigated seismic velocity discontinuities and velocity structure beneath Aso Caldera using receiver functions and a genetic algorithm inversion. We confirm the existence of the Moho at depths between 30 km and 35 km and a large velocity anomaly should exist in the deep portion of the crust beneath Aso Caldera, from imaging of receiver functions observed only at stations outside the caldera. As a result of a more detailed examination with GA inversion, a low velocity layer is detected at depths between 10 km and 24 km beneath the western part of the caldera. S-wave velocity of the layer is estimated to be 2.0–2.4 km/s. We estimate that the low velocity layer contains at most 15% melt or 30% aqueous fluid. The layer exists near the Conrad and at the same depths as the swarm of the low frequency earthquakes and a compressional and dilatational deformation source which are expected to be caused by fluid movement beneath the middle-eastern part of the caldera. Fluid contained in the layer might be related with huge pyroclastic eruptions of Aso Volcano.     

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.     

Seismicity of the Prince William Sound Region for over Thirty Years Following the 1964 Great Alaskan Earthquake

—In this paper we present results of body wave-form modeling of 19 earthquakes (generally m _b ≥5.7) occurring from 1964 to 1983 in the vicinity and down-dip of the large asperity within the Prince William Sound region that ruptured in 1964. These data are supplemented with source parameters from studies of more recent (post-1980) events. Our results suggest that moderate earthquakes which occurred in the region between 1964 and 1984 were predominantly located in the vicinity of the Prince William Sound asperity and could be assigned to two groups. The first group consists of events occurring above the plate interface within Prince William Sound along reverse faults or low angle thrusts. The second group occurs at 35 to 60 km depth in the region north of Prince William Sound, and represents normal to normal-oblique faulting within the subducted Pacific crust or upper mantle. These earthquakes occur below the northern edge of the 1964 asperity in a region where the subducting plate undergoes a rapid change in strike and dip. A third group of events occurs in Cook Inlet well down-dip of the 1964 asperity and below the plate interface. These events exhibit a variety of mechanisms and many at depths of 50 to 70 km may be associated with complexities in the shape of the downgoing slab. Most of the Cook Inlet events occurred after 1984, whereas a few events of similar magnitude have occurred in the vicinity of the Prince William Sound asperity since 1984.     

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.     

一组可用于探测亚稳态橄榄石存在情况的震相

基于简化的波速结构,利用二维伪谱法研究了俯冲带不同位置震源激发的地震波的传播特征.发现当俯冲板块中存在6%的低速异常时,直达P波和直达S波之间存在波阵面近平面的易于识别的P-S转换震相,其传播方向不因震源位置偏离低速层的中心而改变.对于低速层的边界存在±2 km的随机起伏、低速层为楔形、周围介质为水平分层模型这三种情况,这组震相仍然清晰可辨.对于150 km长度的橄榄石楔,模拟结果显示仍然能在一定震中距范围内接收到这一震相.因此,可以尝试用这组震相来探测俯冲带亚稳态橄榄石的存在状态.