Scattered PKKP: Further evidence for scattering at a rough core-mantle boundary

PKKP signals from Novaya Zemlya recorded at LASA at distances around 60° show consistent anomalies in both slowness and azimuth. The observed anomaly suggests that the signal is a BC branch arrival, although the arrival time corresponds to the DF branch. The BC branch, however, does not extend back to this distance. The azimuth of approach is in the range 229–245°, instead of the expected 186°. These anomalies are associated only with PKKP; analysis of the core phases PKiKP and P′P′ (BC) from the same events show that they arrived at LASA with the appropriate slownesses and azimuths.The PKKP signals can be interpreted as “scattered” PKKP; the scattering occurs on underside reflection at the core-mantle boundary and is probably caused by topographic irregularities on the boundary itself. The calculated scattering region has a surface projection at about 60°S, 134°E, which is outside the diametral plane through source and receiver, and about 21° from the expected PKKP reflection point at 76°S, 95°E.Both the “direct” and “scattered” arms of the PKKP signal have a PK path close to that of the “C” end of the BC branch. The unexpectedly large amplitude of the arrival suggests that there may be a focusing of energy at C, which would indicate a change in velocity gradient just above the inner core boundary. The observations nevertheless require, on the scattering interpretation, lateral variations in the topography of the core-mantle boundary and a region of relatively large topography responsible for the anomalous PKKP observations.     

On the theory of core-mantle coupling

This article commences by surveying the basic dynamics of Earth's core and their impact on various mechanisms of core-mantle coupling. The physics governing core convection and magnetic field production in the Earth is briefly reviewed. Convection is taken to be a small perturbation from a hydrostatic, “adiabatic reference state” of uniform composition and specific entropy, in which thermodynamic variables depend only on the gravitational potential. The four principal processes coupling the rotation of the mantle to the rotations of the inner and outer cores are analyzed: viscosity, topography, gravity and magnetic field. The gravitational potential of density anomalies in the mantle and inner core creates density differences in the fluid core that greatly exceed those associated with convection. The implications of the resulting “adiabatic torques” on topographic and gravitational coupling are considered. A new approach to the gravitational interaction between the inner core and the mantle, and the associated gravitational oscillations, is presented. Magnetic coupling through torsional waves is studied. A fresh analysis of torsional waves identifies new terms previously overlooked. The magnetic boundary layer on the core-mantle boundary is studied and shown to attenuate the waves significantly. It also hosts relatively high speed flows that influence the angular momentum budget. The magnetic coupling of the solid core to fluid in the tangent cylinder is investigated. Four technical appendices derive, and present solutions of, the torsional wave equation, analyze the associated magnetic boundary layers at the top and bottom of the fluid core, and consider gravitational and magnetic coupling from a more general standpoint. A fifth presents a simple model of the adiabatic reference state.     

The interaction of elastic waves with a solid-liquid interface,with applications to the core-mantle boundary

The paper discusses basic ideas and principles underlying methods, which have proved useful in the interpretation of diffraction and scattering phenomena by a smooth or slightly rough solid-liquid interface. Generally, the wave interaction may be formulated as an excitation problem; it involves (1) finding an equivalent dislocation or source distribution on the interface, and (2) evaluating the excited wave field. These steps are taken through perturbation theory and/or adopting the appropriate generalization of ray theory. In this context an explicit form of Green's function is also given. The methods have been applied to the core-mantle boundary, with a view toward the interpretation of recent data: (a) diffracted P and S waves around the core (in particular, their attenuation); (b) scattered short-period core phases (in particular, precursors to PKP and PKKP). Other types of wave interaction, and implications for models of the core-mantle boundary structure, are briefly mentioned.Paper presented at the EGS/ESC workshop on Generation and propagation of seismic waves in Neustadt, Federal Republic of Germany, August 1978.     

帕米尔-兴都库什地区中源地震的空间分布和震源机制解特征

利用国际地震中心（ISC）提供的1964至2003年的高精度地震资料,给出了帕米尔－兴都库什地区中源地震带更完整、更明确的几何形态.兴都库什中源地震带（H带）和帕米尔中源地震带（P带）的倾向和最大深度沿走向有变化,可以进一步划分为HW段（H带西段）、HE段（H带东段）、PSW段（P带西南段）、PM段（P带中段）和PNE段（P带东北段）.H带在170～190 km 深度附近存在地震空区,其下方地震带的倾角明显大于上方,接近垂直.同时,空区下方的地震带沿东西方向成连续的倒“V"字形,两个分支的交角接近垂直.西段分支属于HW段,较浅,没有双层结构;东段分支属于HE段,较深,有双层结构.中源地震的震源机制解规律明显.兴都库什地区主要以近垂直的T轴和近水平的垂直于地震带走向的P轴为特征.帕米尔地区的震源机制解由西南到东北逐渐由P轴近水平并沿地震带走向,转变为B轴近水平并沿走向,直至T轴近水平并沿走向.同时,P轴方向由西向东逐渐转为恒定地与地震带的走向相垂直.最后,我们讨论了这一地区地震活动可能的动力学成因.     

Global P-wave tomography: On the effect of various mantle and core phases

In this work, many global tomographic inversions and resolution tests are carried out to investigate the influence of various mantle and core phase data from the International Seismological Center (ISC) data set on the determination of 3D velocity structure of the Earth's interior. Our results show that, when only the direct P data are used, the resolution is good for most of the mantle except for the oceanic regions down to about 1000 km depth and for most of the D″ layer, and PP rays can provide a better constraint on the structure down to the middle mantle, in particular for the upper mantle under the oceans. PcP can enhance the ray sampling of the middle and lower mantle around the Pacific rim and Europe, while Pdiff can help improve the spatial resolution in the lowermost mantle. The outer core phases (PKP, PKiKP and PKKP) can improve the resolution in the lowermost mantle of the southern hemisphere and under oceanic regions. When finer blocks or grid nodes are adopted to determine a high-resolution model, pP data are very useful for improving the upper mantle structure. The resulting model inferred from all phases not only displays the general features contained in the previous global tomographic models, but also reveals some new features. For example, the image of the Hawaiian mantle plume is improved notably over the previous studies. It is imaged as a continuous low velocity anomaly beneath the Hawaiian hotspot from the core-mantle boundary (CMB) to the surface, implying that the Hawaiian mantle plume indeed originates from the CMB. Low-velocity anomalies along some mid-oceanic ridges extend down to about 600 km depth. Our results suggested that later seismic phases are of great importance in better understanding the structure and dynamics of the Earth's interior.     

Core mantle boundary topography from short period PcP, PKP, and PKKP data

We use a total of 839,369 PcP, PKPab, PKPbc, PKPdf, PKKPab, and PKKPbc residual travel times from [Bull. Seism. Soc. Am. 88 (1998) 722] grouped in 29,837 summary rays to constrain lateral variation in the depth to the core-mantle boundary (CMB). We assumed a homogeneous outer core, and the data were corrected for mantle structure and inner-core anisotropy. Inversions of separate data sets yield amplitude variations of up to 5 km for PcP, PKPab, PKPbc, and PKKP and 13 km for PKPdf. This is larger than the CMB undulations inferred in geodetic studies and, moreover, the PcP results are not readily consistent with the inferences from PKP and PKKP. Although the source-receiver ambiguity for the core-refracted phases can explain some of it, this discrepancy suggest that the travel-time residuals cannot be explained by topography alone. The wavespeed perturbations in the tomographic model used for the mantle corrections might be too small to fully account for the trade off between volumetric heterogeneity and CMB topography. In a second experiment we therefore re-applied corrections for mantle structure outside a basal 290 km-thick layer and inverted all data jointly for both CMB topography and volumetric heterogeneity within this layer. The resultant CMB model can explain PcP, PKP, and PKKP residuals and has approximately 0.2 km excess core ellipticity, which is in good agreement with inferences from free core nutation observations. Joint inversion yields a peak-to-peak amplitude of CMB topography of about 3 km, and the inversion yields velocity variations of ±5% in the basal layer. The latter suggests a strong trade-off between topography and volumetric heterogeneity, but uncertainty analyses suggest that the variation in core radius can be resolved. The spherical averages of all inverted topographic models suggest that the data are best fit if the actual CMB radius is 1.5 km less than in the Earth reference model used (i.e. the average outer core radius would be 3478 km).     

The amplitude ratioPcP/P and the core-mantle boundary

Summary Using the Haskell matrix formulation, theoretical reflection coefficient curves have been calculated for a multi-layered core-mantle boundary for comparison with observational data. Two cases are considered, first when the shear velocity in the core is equal to zero and second when the core has a finite rigidity. If the velocity contrast is large between the imbedded layer and the mantle, the reflection coefficient curves for the multi-layered medium are irregular in shape as compared to those for two half-spaces, representing the core and the mantle, respectively. The reflection coefficient curves show an oscillatory character if the imbedded layer is thick and has a high velocity contrast.The observational data consist of short-period vertical-component seismograph records ofP andPcP from nuclear explosions in the Aleutian chain, Nevada, Novaya Zemlya, Kazakh and Sahara. Attenuation and geometrical spreading are taken into consideration. Four different models for the quality factorQ are applied to the observational data. The data are found to be much affected by theQ-model used for the corrections.Based on proposedQ-values, a model for the core-mantle boundary is found, characterized by two low-velocity layers at the bottom of the mantle. The thicknesses are 16.10 km (outer layer) and 19.96 km (inner layer), the compressional wave velocities 12.17 km/sec and 10.94 km/sec and the shear wave velocities are 6.29 km/sec and 5.33 km/sec, respectively. A better fit to this model is found when in addition the shear velocity in the outer core is 2.20 km/sec and the density ratio at the core-mantle boundary is 1.07. In other words, the observations favour a layer of finite rigidity in the outer core rather than a fluid one.     

Dynamos with weakly convecting outer layers: implications for core-mantle boundary interaction

Convection in the Earth's core is driven much harder at the bottom than the top. This is partly because the adiabatic gradient steepens towards the top, partly because the spherical geometry means the area involved increases towards the top, and partly because compositional convection is driven by light material released at the lower boundary and remixed uniformly throughout the outer core, providing a volumetric sink of buoyancy. We have therefore investigated dynamo action of thermal convection in a Boussinesq fluid contained within a rotating spherical shell driven by a combination of bottom and internal heating or cooling. We first apply a homogeneous temperature on the outer boundary in order to explore the effects of heat sinks on dynamo action; we then impose an inhomogeneous temperature proportional to a single spherical harmonic Y ₂² in order to explore core-mantle interactions. With homogeneous boundary conditions and moderate Rayleigh numbers, a heat sink reduces the generated magnetic field appreciably; the magnetic Reynolds number remains high because the dominant toroidal component of flow is not reduced significantly. The dipolar structure of the field becomes more pronounced as found by other authors. Increasing the Rayleigh number yields a regime in which convection inside the tangent cylinder is strongly affected by the magnetic field. With inhomogeneous boundary conditions, a heat sink promotes boundary effects and locking of the magnetic field to boundary anomalies. We show that boundary locking is inhibited by advection of heat in the outer regions. With uniform heating, the boundary effects are only significant at low Rayleigh numbers, when dynamo action is only possible for artificially low magnetic diffusivity. With heat sinks, the boundary effects remain significant at higher Rayleigh numbers provided the convection remains weak or the fluid is stably stratified at the top. Dynamo action is driven by vigorous convection at depth while boundary thermal anomalies dominate in the upper regions. This is a likely regime for the Earth's core.     

帕米尔—兴都库什地区构造应力场反演及拆离板片应力形因子特征研究

从Global CMT目录搜集了1976年1月至2016年6月之间的震源深度大于70km的255个震源机制解,用阻尼应力反演方法,分70～160km和170～310km两个深度,计算了帕米尔—兴都库什地区的构造应力场;同时以10km为间隔计算了兴都库什地区深度介于70～310km之间的应力形因子.得到以下初步结论:兴都库什板片向下俯冲和帕米尔地区断裂带的横向拉张,可能是导致应力场不同的原因.兴都库什俯冲带与帕米尔俯冲带碰撞,导致交汇地区(37°N—37.5°N)的应力场参数突变.兴都库什俯冲板片受到深部温度、压力等因素,出现薄弱面进而形成拆离板片.其脱离了主俯冲板片的束缚后,重力的上下拉张作用导致空区附近张轴倾伏角接近90°,拆离板片俯冲至上地幔不连续面,导致板片部分熔融进而应力形因子随着深度变小.而拆离板片受到地幔挤压其内部发生破碎,其压应力轴由西部的NS到东部NW-SE方向偏转及纵向张应力轴倾伏角变小.     

Topography on the D′′ region from analysis of a thin dense layer beneath a convecting cell

In this study, we examine the development of topography on a thin dense layer at the base of the lower mantle. The effect of the convecting mantle above is represented as a traction acting on the upper surface of the layer. Topography on the layer boundaries is predicted by a balance of dynamic flow stress and external traction. The nature of boundary topography depends on the magnitude of the driving tractions and the density variation within the layer. If we assume that the layer density is greatest beneath areas of mantle downwelling and decreases to a minimum beneath areas of mantle upwelling (the layer is thermally coupled to the convection in the overlying mantle) then its upper boundary develops a cusp-like peak beneath the upwelling mantle. The height of this peak is potentially much greater than the layer thickness. If, however, the layers are effectively coupled by viscous shear then internal density gradients of the opposite sign may be established. In this case, we observe solutions where the layer is completely swept away beneath areas of mantle downwelling leaving steep-sided ‘islands’ of dense material. This mechanism therefore provides a possible explanation for steep-sided anomalously slow regions at the base of the mantle observed by seismic methods (e.g. beneath south Africa) or for discrete ultralow velocity zones detected at the core-mantle boundary beneath locations of surface hotspots. The magnitude of the upper boundary driving tractions compared to the density gradient within the layer is the key parameter that determines the nature of flow in, and consequently boundary topography of, the layer. The deflection of the core-mantle boundary is small compared with that of the top of the dense layer, but a change in sign of the ratio of these deflections is observed as the magnitude of the driving tractions changes relative to the magnitude of the internal density gradient. We compare seismic measurements of core-mantle boundary topography and D′′ topography with the predictions of this model in an attempt to constrain model parameters, but no clear correlation seems to exist between D′′ thickness and CMB topography.     

印度板块与欧亚板块在兴都库什-帕米尔地区相互俯冲的动力作用分析

兴都库什-帕米尔地区是印度板块与欧亚板块相互碰撞的强烈变形区域,也是中国大陆与周边板块动力传递的关键部位,明确该地区两大板块俯冲接触的几何形态和动力作用对研究区域动力环境具有实际意义.本文首先基于Hayes等在2009和2010年提出的Slab1.0的研究思路,利用地震参数准定量地给出两大板块在兴都库什-帕米尔地区碰撞接触的几何形态.结果表明,印度板块在兴都库什地区呈现自南往北的俯冲;欧亚板块在帕米尔地区呈现由北往南的俯冲;同时在兴都库什和帕米尔之间存在俯冲交汇区,在该区印度板块以北西方向、欧亚板块以南东方向相互俯冲.其次基于哈佛大学提供的震源机制解,对不同接触部位进行了应力张量反演,结果显示在兴都库什俯冲区域主要表现为逆冲性质,帕米尔弧西段主要表现为走滑性质,且均具有较好的一致性;而在俯冲交汇区域,走滑、逆冲性质并存,表现为震源机制一致性紊乱.结合两大板块接触的几何形态和区域应力场反演结果,认为印度板块在兴都库什地区主动往北俯冲,而欧亚板块在帕米尔地区被动往南东-南向俯冲,形成两大板块的相互俯冲.本文从几何形态和应力场反演分析两大板块在兴都库什-帕米尔地区碰撞的动力作用方式,可为该区域地球动力学相关研究提供基础资料.     

Intermediate-depth earthquakes beneath the Pamir Hindu Kush Region: Evidence for collision between two opposite subduction zones

We employed a double-difference algorithm (hypoDD) to relocate earthquakes within the region bounded by 66°E-78°E and 32°N-42°N in the period of 1964-2003 reported by the International Seismological Center (ISC). The improved hypocentral locations delineate a double-layered Wadati-Benioff zone in the eastern Hindu Kush intermediate seismic belt. Based on this feature and other evidences, we propose that the intermediate-depth earthquakes beneath the Pamir-Hindu Kush region may occur in two collided subduction zones with opposite dip directions.     

Evidences for a low-velocity core-mantle transition zone

Using the spectral ratios $\text{P}\text{PcP}\text{,}\text{ScS}{}_{\mathrm{n+1}}\text{ScS}{}_{\mathrm{n}}\text{,}\text{sScS}{}_{\mathrm{n+1}}\text{sScS}{}_{\mathrm{n}}\text{and}\text{SKS}\text{ScS}$, models for the core-mantle boundary are found. The models have close similarity with each other, implying an irregular surface with lateral variation in the core-mantle properties. The models are characterized by two to four low-velocity, high-density layers imbedded between the mantle and the core half space. The velocities of the imbedded layers decrease towards the core boundary with a lower bound of 9.3 km/sec for the compressional wave and 3.5 km/sec for the shear wave. The models fitted to the empirical data support the hypothesis of a finite rigid outer core with a higher bound for the shear velocity of 1.4 km/sec. Based on this finite rigidity in the outer core and a layered core-mantle transition zone, the value of Q for the whole mantle is 2,000. For the outer core Q ranges from 100–1,000 , which may indicate that it is chemically zoned.     

Seismotectonics of the Himalaya and its vicinity from centroid-moment tensor (CMT) solution of earthquakes

The centroid-moment tensor solutions of more than 300 earthquakes that occurred in the Himalayas and its vicinity regions during the period of 1977–1996 are examined. The resultant seismic moment tensor components of these earthquakes are estimated. The Burmese arc region shows prominent east–west compression and north–south extension with very little vertical extension. Northeast India and Pamir–Hindu Kush regions show prominent vertical extension and east–west compression. The Indian plate is subducting eastward beneath the northeast India and Burmese arc regions. The overriding Burmese arc has overthrust horizontally with the underthrusting Indian plate at a depth of 20–80 km and below 80 km depth, it has merged with the Indian plate making “Y” shape structure and as a result the aseismic zone has been formed in the region lying between 26°N–28°N and 91.5°E–94°E at a depth of 10–50 km. Similarly, the Indian plate is underthrusting in the western side beneath the Pamir–Hindu Kush region and the overriding Eurasian plate has overthrust it to form a “Y” shape structure at a depth of 10–40 km and below 60 km depth, it has merged with the Indian plate and both the plates are subducting below 60–260 km depth. Further south, the overriding Eurasian plate has come in contact with the Indian plate at a depth of 20–60 km beneath northwest India and Pakistan regions with left lateral strike slip motion.     

On the formation of a stably stratified layer near the core-mantle boundary

Within the framework of a model of liquid immiscibility in the outer core, we calculate a stably stratified layer about 11 km thick near the core-mantle boundary and discuss its reflection and scattering properties for seismic waves.     

Research on the scatter features of the PKiKP/PcP amplitude ratio and the inner core boundary density contrasts beneath Northeast Asia

Quantifying the density contrasts of the Earth's inner core boundary(ICB) is crucial to understand core-mantle coupling and the generation of the geodynamo. The PKiKP/PcP amplitude ratio is commonly used to obtain the density contrast at the ICB, but its applications are limited by scattered observed data. In this study, we selected the PKiKP and PcP phases reflected at the same region of inner-core and core-mantle boundaries beneath Northeast Asia from different earthquakes for the first time, and the observations suggested that the PKiKP/PcP amplitude ratio is widely scattered. We also compared the PKiKP and PcP amplitudes, which demonstrated that the scatter cannot be attributed only to ICB anomalies but might also arise from raypath differences and heterogeneities throughout the crust and mantle. By fitting the observed PKiKP/PcP amplitude ratio, we obtained a density contrast of approximately 0.65 g cm~(-3) and a compressional velocity contrast of approximately 0.87 km s~(-1) at the ICB beneath Northeast Asia. The larger contrast values indicate the possible occurrence of local crystallization occurring at the inner core surface.     

Crustal thickness variations in the Zagros continental collision zone (Iran) from joint inversion of receiver functions and surface wave dispersion

Variations in crustal thickness in the Zagros determined by joint inversion of P wave receiver functions (RFs) and Rayleigh wave group and phase velocity dispersion. The time domain iterative deconvolution procedure was employed to compute RFs from teleseismic recordings at seven broadband stations of INSN network. Rayleigh wave phase velocity dispersion curves were estimated employing two-station method. Fundamental mode Rayleigh wave group velocities for each station is taken from a regional scale surface wave tomographic imaging. The main variations in crustal thickness that we observe are between stations located in the Zagros fold and thrust belt with those located in the Sanandaj–Sirjan zone (SSZ) and Urumieh–Dokhtar magmatic assemblage (UDMA). Our results indicate that the average crustal thickness beneath the Zagros Mountain Range varies from ～46 km in Western and Central Zagros beneath SHGR and GHIR up to ～50 km beneath BNDS located in easternmost of the Zagros. Toward NE, we observe an increase in Moho depth where it reaches ～58 km beneath SNGE located in the SSZ. Average crustal thickness also varies beneath the UDMA from ～50 km in western parts below ASAO to ～58 in central parts below NASN. The observed variation along the SSZ and UDMA may be associated to ongoing slab steepening or break off in the NW Zagros, comparing under thrusting of the Arabian plate beneath Central Zagros. The results show that in Central Iran, the crustal thickness decrease again to ～47 km below KRBR. There is not a significant crustal thickness difference along the Zagros fold and thrust belt. We found the same crystalline crust of ～34 km thick beneath the different parts of the Zagros fold and thrust belt. The similarity of crustal structure suggests that the crust of the Zagros fold and thrust belt was uniform before subsidence and deposition of the sediments. Our results confirm that the shortening of the western and eastern parts of the Zagros basement is small and has only started recently.     

On the Mechanism of the Secular Tidal Acceleration of the Solid Inner Core and the Viscosity of the Liquid Core

As is known, the secular deceleration of the Earth's diurnal rotation is explained mainly by the tidal friction in the ocean. Below we consider this mechanism in some detail, taking into account also elastic deformations of the mantle under the action of ocean loading and the interaction between the tide-generating body, ocean tidal wave, liquid outer core, and solid inner core. It is shown that elastic displacements of the core-mantle boundary under the action of ocean loading are of about the same amplitude and phase as the elastic loading displacements of the Earth's outer surface. As a result, side by side with the mechanism of secular deceleration of diurnal rotation of the mantle, there are also (1) the opposite mechanism of secular acceleration of diurnal rotation of the outer liquid core and of the solid inner core and (2) the mechanism of excitation of differential rotation in the liquid core. Taking these effects into account, we compare theoretical and modern observed data on the eastward drift of the solid inner core. It is shown that the best agreement may be obtained if the turbulent viscosity of the liquid core is about 2 × 10 ³ Poise     

Research on the seismogenic environment for deep earthquake and the cause of earthquakes in Xinjiang and its adjacent areas

ＲｅｓｅａｒｃｈｏｎｔｈｅｓｅｉｓｍｏｇｅｎｉｃｅｎｖｉｒｏｎｍｅｎｔｆｏｒｄｅｅｐｅａｒｔｈｑｕａｋｅａｎｄｔｈｅｃａｕｓｅｏｆｅａｒｔｈｑｕａｋｅｓｉｎＸｉｎｊｉａｎｇａｎｄｉｔｓａｄｊａｃｅｎｔａｒｅａｓＱｉａｎｇＬＩ（李强）；Ｒｕｉ－Ｆｅｎｇ...     

中国东部地区的壳-幔过渡带结构

莫霍面是地壳和上地幔的边界,但莫霍面并不是一个简单的"面",而是一个反映地壳和地幔物质交换、相互作用等动力学意义的"过渡带".本文综合深地震反射、宽角地震折射和高温高压岩石物理实验结果,确定壳-幔过渡带的地震P波速度变化范围为6.8~7.5 km·s^-1.在克拉通等构造活动稳定地区壳-幔过渡带内的速度梯度强且壳-幔过渡带厚度薄,而在造山带等构造活动区域壳-幔过渡带内的速度梯度弱且壳-幔过渡带厚.中国东部地区的壳-幔过渡带的平均厚度约为5~10 km,在四川盆地下方最薄（<5 km）,而在华北克拉通中部造山带下方的壳-幔过渡带最厚（~30 km）.综合地球化学结果,华北中部巨厚壳-幔过渡带主要是幔源岩浆的底侵作用和堆晶作用而形成.