Viscous mantle tangential stresses in a numerical model of a supercontinental cycle

The values of maximum tangential stresses, as well as their spatial orientations, are calculated at each point of the mantle. These calculations are performed for successive stages of the supercontinental cycle in terms of a numerical 2-D model of assembling and diverging continents driven by mantle flows; in turn, the flows themselves are formed under the thermal and mechanical action of continents. Zones of maximum tangential stresses, within which their values are equal to or exceed 30 MPa (300 bar), are shown to be located in places where ascending or descending mantle flows approach the horizontal mantle boundaries. These zones have the shape of a bowl with the broad side toward the boundary, and their maximum tangential stresses are oriented at an angle of about 45° to the horizontal. The dimensions of these zones amount approximately to 300 km laterally and 300 km in depth. Their subhorizontal and subvertical tangential stresses are close to zero. It is also shown that the stress values vary significantly during the evolution of the suboceanic mantle under a new opening ocean (after the breakup of a supercontinent). Thus, the maximum tangential stresses in this vast mantle region at the stage of fast ocean opening are, on average, more than two times higher than at the subsequent slow stage. Note that, at all stages, the axes of the maximum tangential stresses are oriented nearly horizontally or vertically almost throughout the aforementioned area.     

Horizontal stresses in the mantle and in the moving continent for the model of two-dimensional convection with varying viscosity

Numerical experiments on studying the spatial fields and evolution of viscous overlithostatic horizontal stresses and pressure in the mantle and in the moving continent are carried out. The continent moves consistently with time-dependent forces, which act from the viscous mantle. By introducing the varying viscosity, we gain the possibility for taking into account the oceanic lithosphere and the difference between the viscosity of the upper and the lower mantle in the context of a purely viscous model. The typical overlithostatic horizontal stresses in the main part of the mantle are ±(7–9) MPa (70–90 bar); in the highly viscous regions and, particularly, in the subduction zones they are at least three times larger. The descending mantle flows in the depth interval from approximately 50 km to about 300 km are more sharply pronounced in the pressure field than in the field of horizontal stresses. At the considered stages of motion and in different parts, the continent is characterized by the following typical values of stresses: the overlithostatic pressure ranges from ?5 to +15 MPa; the horizontal overlithostatic tensile stress amounts up to ?4MPa (?40 bar); and the compressive stress in case of the overriding of the subduction zone attains +35 MPa (350 bar).     

Implications of melting for stabilisation of the lithosphere and heat loss in the Archaean

The increased depth and volume of melting induced in a higher temperature Archaean mantle controls the stability of the lithosphere, heat loss rates and the thickness of the oceanic crust. The relationship between density distributions in oceanic lithosphere and the depth of melting at spreading centres is investigated by calculating the mineral proportions and densities of residual mantle depleted by extraction of melt fractions. The density changes related to compositional gradients are comparable to those produced by thermal effects for lithosphere formed from a mantle which is 200°C or more hotter than modern upper mantle. If Archaean continental crust formed initially above oceanic lithosphere, the compositional density gradients may be sufficient to preserve a thick Archaean continental lithosphere within which the Archaean age diamonds are preserved. The amount of heat advected by melts at mid-ocean ridges today is small but heat advected by melting becomes proportionally more important as higher mantle temperatures lead to a greater volume of melt and as the rate of production of oceanic plates increases. Archaean tectonics could have been dominated by spreading rates 2–3 times greater than now and with mantle temperatures between ca. 1600°C and 1800°C at the depth of the solidus. Mid-ocean ridge melting would produce a relatively thick but light refractory lithosphere on which continents could form, protected from copious volcanism and high mantle temperatures.     

The role of serpentine in preferential craton formation in the late Archean by lithosphere underthrusting

Cratons form the cores of continents and were formed within a narrow window of time (2.5–3.2 Gy ago), the majority having remained stable ever since. Petrologic evidence suggests that the thick mantle roots underlying cratons were built by underthrusting of oceanic and arc lithosphere, but paradoxically this requires that the building blocks of cratons are weak even though cratons must have been strong subsequent to formation. Here, we propose that one form of thickening could be facilitated by thrusting of oceanic lithospheres along weak shear zones, generated in the serpentinized upper part of the oceanic lithosphere (crust + mantle) due to hydrothermal interaction with seawater. Conductive heating of the shear zones eventually causes serpentine breakdown at ~ 600 °C, shutting down the shear zone and culminating in craton formation. However, if shear zones are too thin, serpentine breakdown and healing of the shear zone occurs too soon and underthrusting does not occur. If shear zones are too thick, serpentine breakdown takes too long so healing and lithospheric thickening is not favored. Shear zone thicknesses of ~ 18 km are found to be favorable for craton formation. Because the maximal depth of seawater-induced serpentinization into the lithosphere is limited by the depth of the isotherm for serpentine breakdown, shear zone thicknesses should have increased with time as the Earth's heat flux and depth to the serpentine breakdown isotherm decreased and increased, respectively, with time. We thus suggest that the greater representation of cratons in the late Archean might not necessarily be explained by preferential recycling in the early Archean but may simply reflect preferential craton formation in the late Archean. That is, our model predicts that the early Archean was too hot, the Phanerozoic too cold, and the late Archean just right for making cratons.     

Numerical models of subduction of the oceanic crust with basaltic plateaus

Light continents and islands characterized by a crustal thickness of more than 30 km float over a convective mantle, while the thin basaltic oceanic crust sinks completely in subduction zones. The normal oceanic crust is 7 km thick. However, anomalously thick basaltic plateaus forming as a result of emplacement of mantle plumes into moving oceanic lithospheric plates are also pulled into the mantle. One of the largest basaltic plateaus is the Ontong Java plateau on the Pacific plate, which arose during the intrusion of a giant superplume into the plate ～100 Myr ago. Notwithstanding its large thickness (averaging ～30 km), the Ontong Java plateau is still experiencing slow subduction. On the basis of numerical modeling, the paper analyzes the oceanic crust subduction process as a function of the mantle convection vigorousness and the density, thickness, viscosity, and shape of the crust. Even a simplified model of thermocompositional convection in the upper mantle is capable of explaining the observed facts indicating that the oceanic crust and sediments are pulled into the mantle and the continental crust is floating on the mantle.     

The Role of History-Dependent Rheology in Plate Boundary Lubrication for Generating One-Sided Subduction

We have developed a two-dimensional dynamical model of asymmetric subduction integrated into the mantle convection without imposed plate velocities. In this model we consider that weak oceanic crust behaves as a lubricator on the thrust fault at the plate boundary. We introduce a rheological layer that depends on the history of the past fracture to simulate the effect of the oceanic crust. The thickness of this layer is set to be as thin as the Earth's oceanic crust. To treat 1-kilometer scale structure at the plate boundary in the 1000-kilometer scale mantle convection calculation, we introduce a new numerical method to solve the hydrodynamic equations using a couple of uniform and nonuniform grids of control volumes. Using our developed models, we have systematically investigated effects of basic rheological parameters that determine the deformation strength of the lithosphere and the oceanic crust on the development of the subducted slab, with a focus on the plate motion controlling mechanism. In our model the plate subduction is produced when the friction coefficient (0.004–0.008) of the modeled oceanic crust and the maximum strength (400 MPa) of the lithosphere are in plausible range inferred from the observations on the plate driving forces and the plate deformation, and the rheology experiments. In this range of the plate strength, yielding induces the plate bending. In this case the speed of plate motion is controlled more by viscosity layering of the underlying mantle than by the plate strength. To examine the setting of the overriding plate, we also consider the two end-member cases in which the overriding plate is fixed or freely-movable. In the case of the freely-movable overriding plate, the trench motion considerably changes the dip angle of the deep slab. Especially in the case with a shallow-angle plate boundary, retrograde slab motion occurs to generate a shallow-angle deep slab.     

The relationship between glaciation and volcanism: Numerical simulation and the Holocene volcanism in Kamchatka

This study is concerned with numerical simulation of the strain due to glaciation and glacial melting, when a magma zone (a layer containing inclusions of magma and magma cumulates) is present at the crust–mantle boundary. According to analytical solutions of this problem that involves viscous relaxation of an uncompensated depression at the place of the molten glacier, the depth to the zone of increased shear stresses beneath the depression is proportional to its width, while the relaxation duration is proportional to viscosity of the lithosphere and is a few thousand years. These fundamental estimates are corroborated by our numerical simulation. According to it, the magma zone at the Moho boundary shields the zone of increased shear stresses, limiting it from below. The maximum values (12–25 MPa) with glacial thickness 500–1000 m are reached at the top of this layer of low viscosity. The directions of maximum compression (s₁) as calculated for the time after the melting indicate that the magma that rises along dikes is displaced from the center of the magma lens toward its periphery. It is found that glacial unloading makes the dipping faults in the crust above the low-viscosity layer attractors for the rising magma. Glacial unloading accelerates, by factors of a few times, the magma generation in the mantle that occurs following the mechanism of adiabatic decompression, as well as facilitating the accumulation of mantle fluids in the zone of increased shear stresses at the boundary of the low viscosity layer. The magma traverses this deep fluid collector and increases the intensity and explosivity of eruptions at the beginning of an interglacial period. Our numerical simulation results are in general agreement with published data on Early Holocene volcanic eruptions that occurred after the second phase of the Late Pleistocene glaciation in Kamchatka.     

Seismic tomography and continental drift

Based on data of seismic tomography, the structure of the mantle flows of the contemporary Earth and the continental drift are calculated. Results of calculation of the contemporary motion of continents and their future drift for 150 Myr are presented. The present-day positions of six continents and the nine largest islands are taken as an initial state. The contemporary temperature distribution in the mantle is calculated according to the data of seismic tomography. The 3-D distribution of seismic wave velocities is converted into the density distribution and then into the temperature distribution. The Stokes equation is numerically solved for flows in a viscous mantle with floating continents for the given initial temperature distribution. In this way, the velocities of convective flows are determined in the entire present-day mantle and the surface distribution for the Earth’s heat flux is obtained. The reliability of the calculated flows in the mantle is estimated by the comparison of the calculated velocities of the contemporary continents and oceanic lithosphere with data of satellite measurements. Further, evolutionary equations of convection with floating continents were numerically solved. The calculated structure of mantle flows, temperature distribution, and position of continents are presented for a time moment 150 Myr in the future. The resulting successive changes in the position of continents in time show how islands (in particular, Japan and Indonesia) will be attached to continents and how continents will converge, exhibiting a tendency toward the formation of a new supercontinent in the southern hemisphere of the Earth.     

Isotopic inferences on the chemical structure of the mantle

Variations in the isotopic composition of rocks derived from the upper mantle can be used to infer the chemical history and structure of the Earth's interior. The most prominent material in the upper mantle is the source of mid-ocean ridge basalts (MORB). The MORB source is characterized by a general depletion in incompatible elements caused by the extraction of the continental crust from the mantle. At least three other isotopically distinct components are recognized in the suboceanic mantle. All three could be generated by the recycling of near surface materials (oceanic crust, pelagic sediments, continental lithospheric mantle) into the mantle by subduction. Therefore, the isotope data do not require a compositionally layered mantle, but neither do they deny the existence of such layering. Correlations between the volumetric output of plume volcanism with the reversal frequency of the Earth's magnetic field, and between the geographic distribution of isotopic variability in oceanic volcanism with seismic tomography suggest input of deep mantle material to surface volcanism in the form of deep mantle plumes. Volcanism on the continents shows a much wider range in isotopic composition than does oceanic volcanism. The extreme isotopic compositions observed for some continental magmas and mantle xenoliths indicate long-term (up to 3.3 Gyr) preservation of compositionally distinct material in thick (>200 km) sections of continental lithospheric mantle.     

俯冲带的后撤与弧后扩张

西太平洋地壳年龄较老,因而岩石层较冷和比重较大,俯冲带的角度也较大,活动和成熟的弧后盆地则较多;条件与之相反的东太平洋弧后盆地则较少.本文探讨这种相关关系的力学成因,计算了俯冲板块诱生的弧后上涌地幔流动.计算表明,俯冲角度大及存在后撤俯冲时,有利于在弧后地区产生明显的上涌地幔流,这种深部热物质的上涌会导致弧后扩张.反之,年龄较轻的海洋地块较热和较轻,俯冲角度一般也较小,不易诱生上涌地幔物质流动和弧后扩张.大陆地壳密度小于地幔物质,大陆碰撞区就更不具备弧后扩张的条件.     

海洋板块俯冲作用下上覆大陆岩石层减薄机制的动力学模拟

几乎所有大陆岩石层的减薄现象,可能都与海洋板块的俯冲作用相关,但是两者之间的内在联系迄今仍不十分明确,为此,我们设计了一系列包含洋-陆俯冲系统的二维数值模型,来探讨海洋板块的俯冲作用对上覆大陆岩石层变形行为的影响,尤其对大陆岩石层减薄效应的制约.模型结果表明,海洋板块俯冲过程中的地幔楔熔体对大陆岩石层地幔的热侵蚀以及由熔体上升所诱发的地幔局部对流的强烈扰动会导致上覆大陆岩石层的减薄效应.这种效应不仅表现在横向上的向陆内蔓延,还表现在垂向上的向浅部发展.且多类动力学参数都能制约大陆岩石层的减薄效应.具体地,随着汇聚速率和洋壳厚度的增加,上覆大陆岩石层在横向上的减薄范围越大,在垂向上的减薄程度也越深;而随着俯冲海洋板块年龄的增加,上覆大陆岩石层在横向上的减薄范围增大,但在垂向上的减薄程度会减小;随着上覆大陆岩石层厚度的增加,其横向减薄范围会减小,但在垂向上的减薄程度会加深.本文研究成果能为揭示华北克拉通减薄/破坏的动力学过程提供一定的理论参考依据.     

大陆岩石圈、地幔底部异常体与地幔对流相互作用的数值模拟

在地球表层存在着占地表面积约30%的具有低固有密度、高黏度的大陆岩石圈.由于其特殊的物理化学性质,大陆岩石圈通常不直接参与下方的地幔对流,但其与地幔对流格局有着重要的相互影响.大量研究显示,在中太平洋和非洲的下地幔底部,存在着两块占核幔边界（CMB）面积约20%的高密度热化学异常体（由于其剪切波速度较低,常称作低剪切波速度省（LSVPs））.LSVPs的演化既受地幔对流的影响,同时也影响地幔物质运动的格局和动力学过程.本文系统研究了存在大陆岩石圈,下地幔LSVPs的地幔对流模型.模拟结果显示：（1）当大陆体积较小时,其边缘常伴随着俯冲,大陆区域地幔常处于下涌状态,其上地幔温度较低,大陆岩石圈在水平方向处于压应力状态.随着大陆体积的增大,大陆边缘的俯冲逐渐减弱,大陆区域地幔由下涌转为上涌,其上地幔温度较高,大陆岩石圈水平方向处于拉应力状态.（2） 岩石圈与软流圈边界（LAB）在大陆下方较深,温度较低;在海洋区域较浅,温度较高.随着大陆体积的增大,陆洋之间LAB深度、温度的差异逐渐减小.（3）大陆区域地幔底部LSVPs物质的丰度与大陆的体积呈正相关.当大陆体积较小时,大陆下方的LSVPs丰度比海洋区域少.随着大陆体积的增大,大陆下方LSVPs的丰度逐渐增大.（4）海洋地区地表热流高,且随时间波动大,大陆地区地表热流低,随时间波动较小;LSVPs区域的核幔边界热流低.     

利用横波分裂到时差确定北美夏洛特皇后群岛岩石圈各向异性分层参数

本文收集到北美夏洛特皇后群岛VIB和DIB台站2005年1月1日至2014年12月31日2700余条近震三分量横波观测记录.经严格筛选提取到55条35°横波窗内的快慢波到时差值.通过建立到时差与多界面横波分裂路径关系,实现了对横波不同深度分裂界面的可靠识别和深度确定.依据对横波到时差及归一化到时差的分析结果,确认地处俯冲带的夏洛特皇后群岛的地震各向异性分别存在于该地区的上层陆壳、中层洋壳俯冲带及俯冲带底部或下层岩石圈顶部附近.发现除地幔岩石圈外,上两个分层的各向异性强度由北向南减弱,同时向西运动的洋壳俯冲带向南加深.不仅如此,归一化到时差结果显示该地区2012年M7.8强震前后各向异性强度减弱,但在3个各向异性层存在较大差异.无论在分层界面较浅的北部格雷厄姆岛附近还是在分层界面较深的南部莫尔兹比岛附近,强震后最上层的陆壳和最下层地幔岩石圈的各向异性强度均没有变化.M7.8后各向异性减弱全部发生在中层洋壳俯冲带内.并且,距离强震震中越近各向异性强度减弱越大,在一定程度上揭示出强震前后应力变化的空间分布特征.表明该地区洋壳俯冲是引起M7.8强震和现今构造运动的主要动力源.     

The structure of 2D mantle convection and stress fields: Effects of viscosity distribution

Spatial fields of temperature, velocity, overlithostatic pressure, and horizontal stresses in the Earth’s mantle are studied in two-dimensional (2D) numerical Cartesian models of mantle convection with variable viscosity. The calculations are carried out for three different patterns of the viscosity distribution in the mantle: (a) an isoviscous model, (b) a four-layer viscosity model, and (c) a temperature- and pressure-dependent viscosity model. The pattern of flows, the stresses, and the surface heat flow are strongly controlled by the viscosity distribution. This is connected with the formation of a cold highly viscous layer on the surface, which is analogous to the oceanic lithosphere and impedes the heat transfer. For the Rayleigh number Ra = 10⁷, the Nusselt number, which characterizes the heat transfer, is Nu = 34, 28, and 15 in models with constant, four-layered, and p, T-dependent viscosity, respectively. In all three models, the values of overlithostatic pressure and horizontal stresses σ _xx in a vast central region of the mantle, which occupies the bulk of the entire volume of the computation domain, are approximately similar, varying within ±5 MPa (±50 bar). This follows from the fact that the dimensionless mantle viscosity averaged over volume is almost similar in all these models. In the case of temperature- and pressure-dependent viscosity, the overlithostatic pressure and stress σ _xx fields exhibit much stronger concentration towards the horizontal boundaries of the computation domain compared to the isoviscous model. This effect occurs because the upwellings and downwellings in a highly viscous region experience strong variations in both amplitude and direction of flow velocity near the horizontal boundaries. In the models considered with the parameters used, the stresses in the upper and lower mantle are approximately identical, that is, there is no denser concentration of stresses in the upper or lower mantle. In contrast to the overlithostatic pressure field, the fields of horizontal stresses σ _xx in all models do not exhibit deep roots of highly viscous downwelling flows.     

Numerical modeling of mantle diapirism as a cause of intracontinental rifting

The numerical model of mantle diapirism and active rifting is developed. The model describes the possibility of extension and thinning of the Earth’s crust under the action of a local 100-km long heat source in the sublithospheric mantle, which causes melting and rising of the magmatic diapir through the cratonic lithosphere. The model combines the mechanisms of the uplifting of the anomalously hot material due to its gravitational instability, underplating of magma beneath the continental crust, and its extension by the forces of the convective flows at the base of the plate. The obtained results shed light on some geological features of the joint formation of the large Vilyui igneous province and Vilyui sedimentary basin.     

Rheology of the mantle and tectonics of the oceanic lithospheric plates

The modern concepts of the rheology of viscous mantle and brittle lithosphere, as well as the results of the numerical experiments on the processes in a heated layer with a viscosity dependent on pressure, temperature, and shear stress, are reviewed. These dependences are inferred from the laboratory studies of olivine and measurements of postglacial rebound (glacial isostatic adjustment) and geoid anomalies. The numerical solution of classical conservation equations for mass, heat, and momentum shows that thermal convection with a highly viscous rigid lithosphere develops in the layer with the parameters of the mantle with the considered rheology under a temperature difference of 3500 K, without any special additional conditions due to the self-organization of the material. If the viscosity parameters of the lithosphere correspond to dry olivine, the lithosphere remains monolithic (unbroken). At a lower strength (probably due to the effects of water), the lithosphere splits into a set of separate rigid plates divided by the ridges and subduction zones. The plates submerge into the mantle, and their material is involved in the convective circulation. The results of the numerical experiment may serve as direct empirical evidence to validate the basic concepts of the theory of plate tectonics; these experiments also reveal some new features of the mantle convection. The probable structure of the flows in the upper and lower mantle (including the asthenosphere), which shows the primary role of the lithospheric plates, is demonstrated for the first time.     

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.     

中国及其邻区地球三维结构初始模型的建立

对人工地震测深及天然地震面波体波三维层折反演数据进行统一处理，建立了中国及其邻区地球三维结构初始模型．此模型图像表明，中国及其邻区地球各圈层横向变化明显．岩石圈及软流圈内速度分布主要反映这一区域自古生代以来板块及地块拼合模式．各主要板块或地块（塔里木、扬子、中朝、青藏、哈萨克斯坦、印度、印度支那）岩石圈增厚或有很深的地慢根，板块或地块间的造山带岩石圈减薄，软流圈速度降低。下地幔底部及核幔边界D″层出现高速异常，表明古太平洋及古特提斯洋俯冲板块因重力坍塌已进入地球深层，形成亚洲超级下降地幔柱。这一下降地幔柱引起地球表层物质向中亚、东亚地区集中，印度半岛、青藏高原、新疆、蒙古至贝加尔一带，成为全球岩石圈最大的汇聚场所．     

论青藏高原及邻区板片构造的一个新模式

本文首先论述了板块学说提出的过程和存在的一些不足与疑问,特别是该学说将Holmes（1948）的地幔热对流说作为驱使岩石圈板块运动的动力机制.而后又以青藏高原及邻区为例,根据区域地质、蛇绿岩和地质构造研究的成果,特别是地震测深研究的成果,详细地论证了本区不存在有大洋中脊扩张成为大洋盆地的新大洋和大洋板块简单的B型俯冲模式,但存在有海底扩张的陆间海和海洋地壳板片（蛇绿岩构造岩片）的仰冲以及大陆岩石圈板片复杂的A型俯冲新模式.新模式不是以地幔对流运动,而是以扩张分离A型俯冲的大陆岩石圈板片与软流圈之间的水平剪切相对运动机制作为它的躯动力.     

上通流对大陆岩石圈地幔-地壳热结构模式的潜在影响

本文以多孔介质中大尺度传热问题为基础,结合热平衡理论分析与数值计算,探讨了上通流对大陆岩石圈地幔-地壳热结构模式的潜在影响.根据大陆岩石圈中孔隙波传热概念模型的初步理论分析结果,指出了采用理论分析和数值模拟相结合的方法在研究大陆岩石圈地幔-地壳热结构模式时的重要性.理论分析方法可用来确定岩石圈尺度范围内大陆岩石圈的厚度和大陆地壳相关的边界条件,从而为地壳范围内数值模型的建立提供一些重要信息.数值模拟方法可以用来模拟地壳尺度范围内地壳的详细结构和复杂几何形状.如果地壳内的热分布是所考虑的主要因素,采用具有地壳尺度的合理数值模型可以有效减少计算机工作量.利用理论分析方法求出的岩石圈尺度范围内大陆岩石圈厚度与地幔传导热流之间关系的理论解,不仅可以用来验证模拟大陆岩石圈内传热问题所采用的数值方法, 而且可以用来初步研究大陆岩石圈内热分布的基本规律,为研究岩石圈地幔热事件中大陆岩石圈热减薄过程提供相应的边界条件.本文从理论分析的观点初步探讨了中国大陆不同构造背景下大陆岩石圈的热结构模式,其结果与从地球物理和地质资料中获得的大陆岩石圈热结构模式十分吻合.研究结果表明由大陆岩石圈中孔隙波传播所导致的上通流是影响大陆岩石圈地幔-地壳热结构模式及大陆岩石圈地幔与地壳之间物质和能量交换的可能机制之一.