Tomographic structure of East Asia:Ⅱ.Stagnant slab above 660 km discontinuity and its geodynamic implications

P-wave arrival times of both regional and teleseismic earthquakes were inverted to obtain mantle structures of East Asia.No fast(slab) velocity anomalies was not find beneath the 660-km discontinuity through tomography besides a stagnant slab within the transition zone.Slow P-wave velocity anomalies are present at depths of 100-250 km below the active volcanic arc and East Asia.The western end of the flat stagnant slab is about 1 500 km west to active trench and may also be correlated with prominent surface topographic break in eastern China.We suggested that active mantle convection might be operating within this horizontally expanded "mantle wedge" above both the active subducting slabs and the stagnant flat slabs beneath much of the North China plain.Both the widespread Cenozoic volcanism and associated extensional basins in East Asia could be the manifestation of this vigorous upper mantle convection.Cold or thermal anomalies associated with the stagnant slabs above the 660-km discontinuity have not only caused a broad depression of the boundary due to its negative Clapeyron slope but also effectively shielded the asthenosphere and continental lithosphere above from any possible influence of mantle plumes in the lower mantle.     

The structure of thermal convection and geotherms beneath a continent and ocean

The formation of the thermal cross section of the lithosphere and mantle upon the interaction between the mantle convection and the immobile continent surrounded by the oceanic lithosphere is studied by numerical modeling. The convective temperature and velocity fields and then the averaged geotherms for subcontinental and suboceanic regions up to the boundary with the core are calculated from the solution of convection equations with a jump in viscosity in the continental zone. Using the experimental data on the solidus temperature in the rocks of the upper mantle, the average thickness of the continental and oceanic lithosphere is estimated at 190 and 30 km, respectively. The effect of a hot spot formed in the subcontinental upper mantle at a depth of 250–500 km, which has not been previously noted, is revealed. Although the temperature in this zone is typically assumed to be close to adiabatic, the calculations show that it is actually higher than adiabatic by up to 200°C. The physical mechanism responsible for this effect is associated with the accumulation of convective heat beneath the thermally insulating layer of the continental lithosphere. The revealed anomalies can be important in studying the phase and mineral transformations at the base of the lithosphere and in the regional geodynamical reconstructions.     

Large aspect ratio cells in two-dimensional thermal convection

Numerical experiments have been carried out on two-dimensional thermal convection, in a Boussinesq fluid with infinite Prandtl number, at high Rayleigh numbers. With stress free boundary conditions and fixed heat flux on upper and lower boundaries, convection cells develop with aspect ratios (width/depth) λ? 5, if heat is supplied either entirely from within or entirely from below the fluid layer. The preferred aspect ratio is affected by the lateral boundary conditions. If the temperature, rather than the heat flux, is fixed on the upper boundary the cells haveλ ≈ 1. At Rayleigh numbers of 2.4 × 10⁵ and greater, small sinking sheets are superimposed on the large aspect ratio cells, though they do not disrupt the circulation. Similar two-scale flows have been proposed for convection in the earth's mantle. The existence of two scales of flow in two-dimensional numerical experiments when the viscosity is constant will allow a variety of geophysically important effects to be investigated.     

Flow visualization studies of two-phase thermal convection in a porous layer

Summary Two-phase thermal convection has been studied in a porous layer heated from below. A water saturated porous layer was heated so that boiling occurred on the lower boundary. In order to observe flow patterns one lateral dimension of the apparatus was made small. At moderate heat fluxes a water zone overlay a two-phase, steam-water zone. Water velocities and streamlines were obtained as well as the location of the two-phase zone for several heat fluxes. Within the water zone heat transfer took place due to both conduction and convection. In the two-phase zone heat transfer took place due to counterpercolation of steam and water.     

Phase transformations in the earth's mantle and subducting slabs: Implications for their compositions, seismic velocity and density structures and dynamics

Abstract Phase transformations in model mantle compositions and those in subducting slabs have been reviewed to a depth of 800 km on the basis of recent high-pressure experimental data. Seismic velocity and density profiles in these compositions have also been calculated using these and other mineral physics data. The nature of the seismic velocity and density profiles calculated for a pyrolite composition was found to generally agree with those determined by seismic observations (e.g. PREM). The locations of the seismic discontinuities at 400 and 670 km correspond almost exactly to the depths where the transformations of the olivine component to denser phases take place. Moreover, the steep gradients in the seismic velocity/density profiles observed between these depths are qualitatively consistent with those expected from the successive transformations in the complementary pyroxene-garnet component in the pyrolite composition. Further, the calculated seismic velocity and density values agree well with those observed in the upper mantle and mantle transition region within the uncertainties attached to these calculations and observations. Pyrolite or peridotite compositions are thus most likely to represent the composition of the mantle above 670 km depth, although some degrees of chemical heterogeneity may exist in the transition region. The observed sharp discontinuous increases of seismic velocities and density at this depth may be attributed either to the phase transformation to a perovskite-bearing assemblage in pyrolite or to chemical composition changes. Density profiles in subducted slabs have been calculated along adequate geotherms assuming that the slabs are composed of the former oceanic crust underlain by a thicker harzburgitic layer. It is shown that the former oceanic crust is substantially less dense than the surrounding pyrolite mantle at depths below 670 km, while it is denser than pyrolite in the upper mantle and the transition region. The subducted former oceanic crust may be trapped in this region, forming a geochemically enriched layer at the upper mantle-lower mantle boundary. Thick and cool slabs may penetrate into the lower mantle, but the chemically derived buoyancy may result in strong deformation and formation of megalith structures around the 670 km seismic discontinuity. These structures are consistent with those detected by recent seismic tomography studies for subduction zones.     

On the geodynamic setting of kimberlite genesis

The emplacement of kimberlites in the North American and African continents since the early Palaeozoic appears to have occurred during periods of relatively slow motion of these continents. The distribution of kimberlites in time may reflect the global pattern of convection, which forces individual plates to move faster or slower at different times. Two-dimensional numerical experiments on a convecting layer with a moving upper boundary show two different regimes: in the first, when the upper boundary velocity is high, heat is transferred by the large-scale circulation and in the second, when the upper boundary velocity is lower, heat is predominantly transferred by thermal plumes rising from the lower boundary layer. For a reasonable mantle solidus, this second regime can give rise to partial melting beneath the moving plate, far from the plate boundaries. The transition between these modes takes place over a small range of plate velocities; for a Rayleigh number of 10⁶ it occurs around 20 mm yr^?1. We suggest that the generation of kimberlite magmas may result from thermal plumes incident on the base of a slowly moving plate.     

Long wavelength thermal convection between non-conducting boundaries

Non-linear Rayleigh-Bénard convection in a fluid layer is considered as a model of convection in the Earth's upper mantle. Previous studies have shown that when the temperature is held fixed at one of the boundaries of the layer, convection takes place in cells of width of the order of the layer depth or less. We investigate the effects of a different thermal boundary condition, in which the flux of heat is held fixed on both layer boundaries; then if this flux is just greater than that required for the onset of convection, motion takes place on horizontal scales much greater than the layer depth. An analytical treatment of the equations, based on an expansion in the depth-to-width ratio of the cells, shows that cells of a definite horizontal scale are the fastest growing according to linearised theory, but that these cells are unstable to ones of larger wavelength than themselves. Thus the dominant wavelength lengthens with time. The results hold whether the heat flux is generated internally of comes from beneath the layer. These results produce flow patterns similar to those found when the heat flux is much greater than the critical value. The results have important consequences for the understanding of mantle convection.     

Geophysical observations pertaining to solid-state convection in the terrestrial planets

We present a broad-based review of the observational evidence that pertains to or otherwise implies solid-state convection to be occurring (or have occurred) in the interiors of the terrestrial planets.For the Earth, the motion of the plates is prima facie evidence of large-scale mantle convection. Provided we understand upper-mantle thermal conductivity correctly, heat flow beneath the old ocean basins may be too high to be transported conductively from the upper mantle through the base of the lithosphere and therefore convection on a second smaller scale might be operative. The horizontal scale of plate dimensions implies, due to typical cell aspect ratios observed in convection, that the motion extends to the core-mantle boundary. Improved global data coverage and viscoelastic modeling of isostatic rebound due to Pleistocene deglaciation imply a uniform mantle viscosity, and thus indicate that whole-mantle convection could exist. Additionally, there is some seismic evidence of lithospheric penetration to depths deeper than 700 km. We discuss some salient features and assumption boundedness of arguments for convection confined to the upper mantle and for convection which acts throughout the mantle since the vertical length scale has a profound effect upon the relevance of geophysical observations. The horizontal form of mantle convection may be fully three-dimensional with complex planform and, therefore, searching for correlative gravity patterns in the ocean basins may not be useful without additional geophysical constraints. Many long-wavelength gravity anomalies may arise from beneath the lithosphere and must be supported dynamically, although thermal convection is not a unique explanation. Topography is an additional geophysical constraint, but for wavelengths greater than a few hundred kilometers, a general lack of correlation exists between oceanic residual gravity and topography, except at specific locations such as Hawaii. Theoretical calculations predict a complex relationship between these two observational types. Oceanic gravity data alone shows no regular planform and there is no correlation with any small-scale convective pattern predicted by laboratory experiments.All of the observational evidence argues against Martian plate tectonics occurring now or over much of the history of this planet, but lack of plate tectonics is not an argument against interior convection. The Tharsis uplift on Mars may have resulted from convective processes in the mantle, and the present-day gravity anomaly associated with Tharsis must be supported by the finite strength of the lithosphere or by mantle convection. Stresses imparted by the present topographic load would be greater than a kilobar, in excess of long-term finite strength. Observed fracture patterns are probably a direct result of this load, and the key question concerns the level of resultant strain relief. The global topographic and geomorphic dichotomy between the northern and southern hemisphere required a solid-state flow process to create the accompanying center-of-figure to center-of-mass offset.Lunar heat flow values, in analogy with oceanic heat flow on the Earth, strongly imply a convective mechanism of heat transport in the interior which, based on seismic Q values, is limited to the lower mantle. The presence of moonquakes in this region does not preclude solid-state convective processes. Lunar conductivity profiles provide no information on convection because of the difficulty in conductivity modeling, uniqueness of models, and the uncertainty in the conductivity-temperature relationship. The excess oblateness of the lunar figure over the hydrostatic value does not require convective support; in fact, such a mechanism is unlikely.The presence of a dipole magnetic field on Mercury does not provide a constraint on mantle convection unless its existence can be inextricably linked to a molten core. The non-hydrostatic shape of the equatorial figure, required for the observed $\text{3}\text{2}$ resonance between Mercury's rotational and orbital periods, is most likely related to surface processes, as opposed to convection. The $\text{3n}\text{2}$ resonance implies escape from a 2n resonance and, therefore, is related to the question of a molten core. Further dynamical data is needed to constrain interior models.Interpretation of limited radar imagery for the surface of Venus is enigmatic in terms of plate tectonics and therefore interior convection. Linear tensional and possibly compressional features are observed, but there are also crustal regions which appear to show large impact structures and are thus geologically old and may not have been recycled.     

华南陆缘高热流区的壳幔温度结构与动力学背景

华南陆缘是我国重要的矿产、地热资源区.晚中生代以来,在太平洋板块西向俯冲,地幔热对流活动共同作用下,该区出现多期岩浆-热事件和大规模爆发式成矿作用.在前人研究基础上,本文利用地表热流观测资料、地震剪切波资料、重力位球谐系数,计算了壳-幔温度结构,分析了动力学背景.计算结果表明：华南陆缘东南沿海地带,地壳10 km以浅温度达200℃以上,居里点温度475℃,莫霍面平均温度550℃.地壳浅层较热,花岗岩中放射性元素衰变放热是地壳浅层地下水热活动的重要热源,但地壳总体温度不高,为"冷壳热幔"型热结构.地幔中,90 km深度,温度950~1250℃;120 km深度,温度1050~1400℃;150 km深度,温度1200~1450℃;220 km深度,温度1500~1700℃."热"岩石圈底界深度在110~150 km之间,西深东浅.岩石圈内,地幔应力场为挤压-伸展相间格局;岩石圈之下,地幔应力场为一个以南昌为中心、长轴NE-SW向的椭圆.分析认为,晚中生代以来,太平洋板块的西向俯冲,导致华南陆缘在区域性SE向地幔对流背景上叠加局域性不稳定热扰动,在175~85Ma期间,上地幔物质向上流动,形成不同的岩浆活动高峰期.同时,岩石圈地幔受俯冲洋壳流体的影响,含水量高,黏度小,在地幔流切向应力场作用下,岩石圈底界由西向东"波浪"状减薄.现今岩石圈之下仍具备地幔小尺度热对流温度条件,但除地表浅层外,地壳整体温度不高,岩石圈构造稳定.     

Low Q zone at the base of the mantle: Evidence for lower mantle convection?

Temperatures in the lower mantle of the Earth are estimated from the observed Q distribution. A thermal boundary-layer where temperatures rise rapidly is found at the base of the mantle, corresponding to the low Q zone described by Anderson and Hart (1978a,b). The existence of this thermal boundary layer indicates that the lower mantle participates in convection, and also that some of the energy driving the convection is coming from the core.     

A theoretical study of hydrothermal convection and the origin of the ophiolitic sulphide ore deposits of Cyprus

Hydrothermal circulation of seawater has been suggested as a mass transport mechanism for the formation of sulphide ore deposits in the ophiolitic rocks of Cyprus. Since ophiolitic sequences are generally regarded as fragments of oceanic crust and upper mantle, hydrothermal circulation of a form inferred from geological observations on Cyprus may be analogous to that thought to occur in oceanic crust at spreading ridges. The hypothesis that ore deposits were formed in ascending plumes of hot, buoyant fluid is examined by considering thermal convection in a permeable medium. To match the inferred pattern of circulation, finite amplitude convection in a cylindrical geometry is studied using finite difference approximations. These results combined with available geological and geochemical data are applied to understand better the physical controls on mineralisation.A simple model for the formation of the hydrothermal ore deposits of Cyprus is discussed. The model is semi-quantitatively reasonable in terms of vertical fluid flow rate, thermal structure, permeability and basal heat flow, and predicts volumes of maximum mineralisation similar to those observed. Three factors are identified which were important in confining mineralisation to a small volume immediately beneath the sea water/rock boundary: (1) hot fluid was confined to a narrow core zone of a rising plume, (2) the upward fluid flux was greatest in this same core zone, and (3) significant temperature decrease occurred within a thin surface boundary layer.     

Numerical model for the generation of the ensemble of lithospheric plates and their penetration through the 660-km boundary

In the kinematic theory of lithospheric plate tectonics, the position and parameters of the plates are predetermined in the initial and boundary conditions. However, in the self-consistent dynamical theory, the properties of the oceanic plates (just as the structure of the mantle convection) should automatically result from the solution of differential equations for energy, mass, and momentum transfer in viscous fluid. Here, the viscosity of the mantle material as a function of temperature, pressure, shear stress, and chemical composition should be taken from the data of laboratory experiments. The aim of this study is to reproduce the generation of the ensemble of the lithospheric plates and to trace their behavior inside the mantle by numerically solving the convection equations with minimum a priori data. The models demonstrate how the rigid lithosphere can break up into the separate plates that dive into the mantle, how the sizes and the number of the plates change during the evolution of the convection, and how the ridges and subduction zones may migrate in this case. The models also demonstrate how the plates may bend and break up when passing the depth boundary of 660 km and how the plates and plumes may affect the structure of the convection. In contrast to the models of convection without lithospheric plates or regional models, the structure of the mantle flows is for the first time calculated in the entire mantle with quite a few plates. This model shows that the mantle material is transported to the mid-oceanic ridges by asthenospheric flows induced by the subducting plates rather than by the main vertical ascending flows rising from the lower mantle.     

流变边界层及其对华北克拉通热/地震岩石圈底界差异的意义

二维热传导/对流数值模型显示,纯传导的固体岩石圈与纯对流的流体软流圈之间存在一过渡层,即流变边界层,其间传导与对流共同作用来传递热量.流变边界层厚度主要由软流圈黏性系数（η）控制,而受固体岩石圈厚度及热状态影响很小.随着η从1×10²¹Pa·s降低至1×10¹⁹ Pa·s,流变边界层也随之减薄,流变边界层的厚度与lg（η）成正比. 流变边界层的存在是造成热岩石圈与地震岩石圈厚度差异的重要因素. 全球典型克拉通岩石圈的对比结果表明,地震岩石圈厚度普遍大于热岩石圈厚度,二者的差异多数在70~90 km,很好地验证了流变边界层的存在. 研究发现二者的差异在华北克拉通自西向东逐渐减小：由西部鄂尔多斯的约80 km减少至渤海湾盆地的约20 km. 反映出华北克拉通岩石圈下部流变边界层厚度自西向东减薄,意味着软流圈黏性系数自西向东逐渐降低.这可能与中生代太平洋俯冲脱水形成的低黏大地幔楔有关,从一侧面印证了太平洋俯冲对华北克拉通破坏的影响.     

Core cooling by subsolidus mantle convection

Although vigorous mantle convection early in the thermal history of the Earth is shown to be capable of removing several times the latent heat content of the core, we are able to construct a thermal evolution model of the Earth in which the core does not solidify. The large amount of energy removed from the model Earth's core by mantle convection is supplied by the internal energy of the core which is assumed to cool from an initial high temperature given by the silicate melting temperature at the core-mantle boundary. For the smaller terrestrial planets, the iron and silicate melting temperatures at the core-mantle boundaries are more comparable than for the Earth, and the cores of these planets may not possess enough internal energy to prevent core solidification by mantle convection. Our models incorporate temperature-dependent mantle viscosity and radiogenic heat sources in the mantle. The Earth models are constrained by the present surface heat flux and mantle viscosity. Internal heat sources produce only about 55% of the Earth model's present surface heat flow.     

Heat flow and depth versus age for the Mesozoic northwest Atlantic Ocean: results from the Sohm abyssal plain and implications for the Bermuda Rise

Heat flow in the Sohm abyssal plain is measured to be 53 mW/m² at an age of 163 Ma. This is 25% higher than predicted by conductive cooling models, even though the sediment-corrected basement depth of 6.5 km at this location is normal for its age. An analysis of existing heat flow, depth and geoid anomalies in the northwest Atlantic shows that there is little correlation between heat flow and depth throughout the entire region. Depth and geoid are clearly related to the Bermuda swell while the associated heat flow anomaly, once adjusted for variations with age, is limited to 5 mW/m² and only decays to the south. This means that the Bermuda swell is probably not caused by extensive thermal reheating within the lithosphere, but instead by dynamic uplift at its lower boundary due to the convective upwelling of a mantle plume. The regionally high heat flow in the northwest Atlantic may be a thermal remanent of previous plumes which passed beneath this region early in its history. Therefore, depth and heat flow anomalies from this region cannot be used to provide constraints on steady-state parameters of the lithosphere, such as the presence or absence of a long-term boundary layer at its base.     

基于波形拟合的中国东海地区410 km间断面附近速度结构研究

410 km间断面是地幔转换带的顶界面,对其速度结构和起伏形态开展地震学探测有助于认识地球内部物质组成和相关的地球动力学过程.本文选取了由中国数字地震台网记录到的位于琉球俯冲区的一个中源地震P波宽频带波形资料,利用三重震相波形拟合研究了中国东海地区410 km间断面附近的精细速度结构.结果表明：中国东海地区下方410 km间断面整体表现为一尖锐的速度界面且有8~15 km的小幅抬升;该间断面之上存在52~62 km厚的低速层,其P波速度降低0.5%~1.6%;440 km深度以下存在1.0%~3.0%的P波高速异常.结合前人在该地区的层析成像结果,我们推测该高速异常体与西太平洋俯冲板片在中国东海地区地幔转换带内的滞留有关;板片内水相E分解使得转换带内水含量增加,这引发了410 km间断面的抬升;410 km间断面之上的低速层应与含水矿物脱水导致的地幔橄榄岩部分熔融有关.     

Three-dimensional numerical simulations of crustal deformation and subcontinental mantle convection

3-D simulations of mantle convection allowing for continental crust are explored to study the effects of crustal thickening on lithosphere stability and of continents on large-scale mantle flow. Simulations begin with a crustal layer within the upper thermal boundary layer of a mantle convection roll in a 1 × 1 × 1 Cartesian domain. Convective stresses cause crust to thicken above a sheet-like mantle downwelling. For mild convective vigor an initial crustal thickness variation is required to induce 3-D lithospheric instability below the zone of crustal convergence. The amplitude of the required variation decreases with increasing convective vigor. Morphologically, instability is manifest in formation of drip-like thermals that exist within the large-scale roll associated with initial crustal thickening. A strong surface signature of the drips is their ability to cause deviations from local Airy compensation of topography. After the initial thickening phase, the crustal accumulation that forms serves as a model analog to a continent. Its presence leads to mantle flow patterns distinctly different from the steady-state roll that results in its absence. Large lateral thermal gradients are generated at its edge allowing this region to be the initiation site for continued small-scale thermal instabilities. Eventually these instabilities induce a restructuring of large-scale mantle flow, with the roll pattern being replaced by a square cell. Although preliminary and idealized, the simulations do show the fluid dynamical plausibility behind the idea that significant mantle variations can be generated along the strike of a largely 2-D mountain chain by the formation of the chain itself. The ability of a model continent to cause a change in fundamental convective planform also suggests that the effects of continental crust on mantle convection may be low-order despite the seemingly trivial volume of crust relative to mantle.     

Upper mantle heterogeneity in the northern part of Eurasia

A two dimensional velocity model of the upper mantle has been compiled from a long-range seismic profile crossing the West Siberian young plate and the old Siberian platform. It revealed considerable horizontal and vertical heterogeneity of the mantle. A sharp seismic boundary at a depth of 400 km outlines the high-velocity gradient transition zone, its base lying at a depth of 650 km. Several layers with different velocities, velocity gradients and wave attenuation are distinguished in the upper mantle. They likewise differ in their inner structure. For instance, the uppermost 50–70 km of the mantle are divided into blocks with velocities from 7.9–8.1 to 8.4–8.6 km s^?1.Comparison of the travel-time curves for the Siberian long-range profile with those compiled from seismological data for Europe distinguished large-scale upper mantle inhomogeneities of the Eurasian continent and allowed for the correlation of tectonic features and geophysical fields. The velocity heterogeneity of the uppermost 50–100 km of the mantle correlates with the platform age and heat flow, i.e., the young plates of Western Europe and Western Siberia have slightly lower velocities and higher heat flows than the ancient East European and Siberian platforms. At greater depths (150–250 km) the upper mantle velocities increase from the ocean to the inner parts of the continent. The structure of the transition zone differs significantly beneath Western Europe and the other parts of Eurasia. The sharp boundary at a depth of 400 km, traced throughout the whole continent as the boundary reflecting intensive waves, transforms beneath Western Europe into a gradient zone. This transition zone feature correlates with positions of the North Atlantic-west Europe geoid and heat-flow anomalies.     

Rheology of the lower mantle

The rheology of the lower mantle of the Earth is examined from the viewpoint of solid state physics. Recent developments in high-pressure research suggest that the lower mantle contains a considerable amount of (Mg, Fe)O with Fe/Mg + Fe = 0.2–0.3. The pressure and temperature dependences of diffusion in (Mg, Fe)O are estimated by the theory of diffusion in ionic solids. Of the materials composing the lower mantle, (Mg, Fe)O may be the “softest”, and therefore the rheology of the lower mantle may be that of (Mg, Fe)O, unless the framework effect is important.Temperatures in the lower mantle are inferred from the depths of phase transitions and the melting temperatures of the core materials. A thermal boundary layer at the base of the mantle is suggested. The physical mechanisms of creep are examined based on a grain size-stress relation and non-Newtonian flow is shown to be the dominant flow mechanism in the Earth's mantle.The effective viscosity for the temperature models, with and without the thermal boundary layer, is calculated for constant stress and constant strain rate (with depth). For constant strain rate, which may be appropriate for discussing the mechanics of descending slabs, the increase in effective viscosity with depth is smaller than for the constant-stress case, which may be appropriate for discussing the flow induced by the surface motion of plates.The relatively small depth gradient of viscosity, for constant strain rate, suggests that the lower mantle could also participate in convection. The effective viscosity increases with depth, however, by at least 10² to 10³ from the top to the bottom of the lower mantle, for a reasonable range of activation volumes and temperatures. There will be a low-viscosity layer at the base of the mantle, in contrast to the high-viscosity layer at the top of the mantle (plates), if a thermal boundary layer is present. The constant Newtonian viscosity inferred from rebound data may be an apparent feature resulting from the difference in deformation mechanisms between isostatic rebound and large-scale flow.