首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 687 毫秒
1.
Thermal convection in the mantle is caused by the heat transported upwards from the core and by the heat produced by the internal radioactive sources. According to the data on the heat transfer by the mantle plumes and geochemical evidence, only 20% of the total heat of the Earth is supplied to the mantle from the core, whereas most of the heat is generated by the internal sources. Along with the models that correctly allow for the internal heat sources, there are also many publications (including monographs) on the models of mantle convection that completely ignore the internal heating or the heat flux from below. In this study, we analyze to what extent these approximations could be correct. The analytical distributions of temperature and heat flux in the case of internal heating without convection and the results of the numerical modeling for convection with different intensity are presented. It is shown that the structure of thermal convection is governed by the distribution of the heat flux in the mantle but not by the heat balance, as it is typically implicitly assumed in most works. Heat production by the internal sources causes the growth of the heat flux as a function of radius. However, in the spherical mantle of the Earth, the heat flux decreases with radius due to the geometry. It turned out that with the parameters of the present Earth, both these effects compensate each other to a considerable extent, and the resulting heat flux turns out to be nearly constant as a function of radius. Since the structure of the convective flows in the mantle is determined by the distributions of heat flux and total heat flux, in the Cartesian models of the mantle convection the effective contribution of internal heating is small, and ignoring the heat flux from the core significantly distorts the structure of the convective currents and temperature distributions in the mantle.  相似文献   

2.
Traditional models for the heat loss in oceanic and continental regions are combined into a regionalized model for the thermal evolution of the Earth. The need for regionalization is obvious when one considers that the mantle loses 3 to 4 times as much heat per unit area in oceanic regions than in continental areas. The present-day rate of heat loss together with a geochemical estimate of the concentration of heat-producing elements in the Earth fixes the response time of the thermally convecting mantle. The response time in turn can be used to select the most reasonable representation for mantle convection in terms of the sensitivity of viscosity on temperature and layering versus mantle-wide circulation. Present geochemical estimates of the bulk composition of the Earth are most easily reconciled with the observed heat flow if the mantle is layered and its rheology is slightly less temperature dependent than generally assumed. The layered system can produce sufficiently high temperatures to explain the high-magnesian komatiites of the Archean. One difficulty with the models is that they predict widespread melting at shallow depth in the early stages of Earth history but do not address how such melting affects and alters the heat transfer mechanisms.  相似文献   

3.
The evolution of the Earth is characterized by irreversible processes: radioactive decay of the major heat-producing elements, thermal convection and chemical segregation. The prevailing heating from within and the temperature dependence of the viscosity are essential for thermal convection. In the present paper, the chemical and thermal evolution of the mantle and the generation of the continent material are represented by a two-dimensional and finite-difference Boussinesq convection model. We have incorporated the above-mentioned principal features in this model, a geochemical paper by Hofmann (1988, Earth Planet. Sci. Lett., 90: 297–314) constituting our starting point for the distribution of the radionuclides. The concentration of the radionuclides and the viscosity are functions of the location and time developing according to our system of differential equations. Although the real Earth is a much more complex system, we have dared to make a comparison with observed geophysical and geological data; we obtain a depleted upper mantle and acceptable values for the heat flow on the surface of the Earth as well as for the distribution of temperature, viscosity and of the velocity of creep in the mantle. The ups and downs of the convective vigour of the model roughly resemble the supercontinental cycles, the world-wide distribution of mineral dates in time, the sea-level variations and the variations of a number of geochemical parameters.  相似文献   

4.
We investigate the thermal and degassing history of the Earth with the help of a parameterized mantle convection model including the volatile exchange between mantle and surface reservoirs. The weakening of mantle silicates by dissolved volatiles is described by a functional relationship between creep rate and water fugacity. We use flow law parameters of diffusion creep in olivine under dry and wet conditions. The mantle degassing rate is considered as directly proportional to the seafloor spreading rate, which is also dependent on the mantle heat flow and the continental area. To calculate the spreading rate, we assume three different continental growth models: constant growth, delayed growth, and the one proposed by Reymer and Schubert (1984, Tectonics, 3: 63–77). The rate of regassing also depends on the seafloor spreading rate, as well as on other factors. Both mechanisms (degassing and regassing) are coupled self-consistently with the help of a parameterized convection model under implementation of a temperature and volatile-content dependent mantle viscosity. We calculate time series for the Earth's evolution over 4.6 Gyr for the average mantle temperature, the mantle heat flow, the mantle viscosity, the Rayleigh number, the Urey ratio, the volatile loss, and the seafloor spreading rate. In those numerical simulations with continental growth from the beginning and a high initial average mantle temperature water is outgassed rapidly. In the delayed continental growth model there is a very early outgassing event and the delayed continental growth has no remarkable influence on the thermal and outgassing history. A similar situation is found for the linear continental growth model but not for the Reymer and Schubert (1984) model.  相似文献   

5.
The theory of three-dimensional and finite-amplitude convection in a viscous spherical shell with temperature and pressure dependent physical parameters is developed on the basis of a modified Boussinesq fluid assumption. The lateral dependences of the variables are resolved through their spherical harmonic representations, whereas their radial and time dependences are determined by numerical procedures. The theory is then applied to produce thermal evolution models for Venus. The emphasis is on illustrating the effects of certain physical parameters on the thermal evolution rather than proposing a specific thermal history for the planet. The main conclusions achieved in this paper are (1) a significant portion of the present temperature in the mantle and heat flux at the surface of Venus is probably owing to the decay of a high temperature established in the planet at the completion of its core formation, (2) the effective Rayleigh number of the mantle is so high that even the lower order modes of convection cool the planet sufficiently and maintain an almost adiabatic temperature gradient in the convecting region and high temperature gradients in the thermal boundary layers, (3) the convection is oscillatory with avalanche type properties which induces oscillatory features to the surface heat flux and the thickness of the crustal layer, and (4) a planetary model with a recycling crust cools much faster than those with a permanently buoyant crust.The models presented in this paper suggest that Venus has been highly convective during its history until ～ 0.5 Ga ago. The vigorous convection was bringing hot and fresh material from the deep interior to the surface and dragging down the crustal slags, floating on the surface, in to the mantle. The rate of cooling of the planet was so high that its core has solidified. In the last 0.5 Ga the vigour of convection diminished considerably and the crustal slags developed into a global and permanently buoyant crustal layer. The tectonic style on Venus has, consequently, changed from the recycling of crustal plates to hot spot volcanics. At the present time the planet is completely solid, except in the upper part of its mantle where partial melting may occur.  相似文献   

6.
Calculations of the radial distribution of the energy released in core formation indicate that the cores of all the terrestrial planets may be expected to receive a disproportionate share of the gravitational energy released. Since the model of the process used in these calculations favors transfer of energy to the mantle, it is likely that other reasonable models of the process will result in even more energy being deposited in the cores of the early planets. The calculations also show that it is necessary for a certain amount of core phase to separate and accumulate, before the energy released by gravitational settling is sufficient to supply the latent heat of fusion of the core phase. The amount of melting required to reach this point varies according to the total mass of the planet, and mass fraction of core, but is not particularly great (<5% for the Earth to ～ 37% for the Moon). In the case of the Moon, this amount of segregation, although large, amounts to a surface layer about 260 km thick, similar to the proposed depth of early wholesale melting. Core separation in terrestrial planets appears to be a self-sustaining process even for fairly small bodies, provided that a small amount of a dense potential core phase is present. Although it seems likely to occur rapidly (within 106–107 years) even for small (Moon-size) bodies, detailed kinetic models will be necessary to specify the time scale.  相似文献   

7.
A mechanism for the production of a chemical change in the mantle, from primordial silicate compositions above the 650 km discontinuity to differentiated compositions below, is reviewed. Some consequences of this are the stabilization of two layer convection with a temperature contrast between the anhydrous mantle solidus and the geotherm which, at 650 km depth, is lower than any other location in the mantle. With thermal contributions from the concentration of the heat producing elements U, Th and K below the 650 km mantle boundary layer and the higher geotherms in the past, widespread or catastrophic melting may have taken place at this location. An episodic breach of this boundary layer by extensive heat and mass transport may have periodically destroyed any simple two-layer convection geometry in the mantle. Such episodic injections of mass and energy into the upper mantle from below may have been the mechanism responsible for episodes of enhanced surface tectonism and thermal activity which appear to be recorded in apparent polar wandering paths and radiometric ages of continental rocks.  相似文献   

8.
Formation,history and energetics of cores in the terrestrial planets   总被引：1，自引：0，他引：1  
The cores of the terrestrial planets Earth, Moon, Mercury, Venus and Mars differ substantially in size and in history. Though no planet other than the Earth has a conclusively demonstrated core, the probable cores in Mercury and Mars and Earth's core show a decrease in relative core size with solar distance. The Moon does not fit this sequence and Venus may not. Core formation must have been early (prior to ～4 · 109 yr. ago) in the Earth, by virtue of the existence of ancient rock units and ancient paleomagnetism and from UPb partitioning arguments, and in Mercury, because the consequences of core infall would have included extensional tectonic features which are not observed even on Mercury's oldest terrain. If a small core exists in the Moon, still an open question, completion of core formation may be placed several hundred million years after the end of heavy bombardment on tectonic and thermal grounds. Core formation time on Mars is loosely constrained, but may have been substantially later than for the other terrestrial planets. The magnitude and extent of early heating to drive global differentiation appear to have decreased with distance from the sun for at least the smaller bodies Mercury, Moon and Mars.Energy sources to maintain a molten state and to fuel convection and magnetic dynamos in the cores of the terrestrial planets include principally gravitational energy, heat of fusion, and long-lived radioactivity. The gravitational energy of core infall is quantifiable and substantial for all bodies but the Moon, but was likely spent too early in the history of most planets to prove a significant residual heat source to drive a present dynamo. The energy from inner core freezing in the Earth and in Mercury is at best marginally able to match even the conductive heat loss along an outer core adiabat. Radioactive decay in the core offers an attractive but unproven energy source to maintain core convection.  相似文献   

9.
Using density–pressure relationships for mantle silicate and core alloy closely matching PREM we have constructed six models of the Earth in different evolutionary states. Gravitational energies and elastic strain energies are calculated for models with homogeneous composition, separated mantle and liquid core, separated inner and outer cores with the inner core either liquid or solid and models with increased densities, representing cooling of either the mantle or core. In this way we have isolated the gravitational energy released by each of several evolutionary processes and subtracted the consequent increase in strain energy to obtain the net energy released as heat or geodynamo power. Radiogenic heat (∼7.8×1030 J) is found to contribute only about 25% of the total heat budget, the balance originating as residual gravitational energy from the original accretion and from core separation (14×1030 J). The total energy of compositional convection, driven by inner core formation, is 3.68×1028 J and this is the most important (or even the only) energy source for the dynamo for the most recent 2 billion years. It appears unlikely that the inner core existed much before that time. The total net (gravitational minus strain) energy released in the core by the process of inner core formation, 11.92×1028 J, is not much less than the thermal energy released in this process, 15.1×1028 J. In the mantle the net (gravitational minus strain) energy released by thermal contraction is about 20% of the heat release. All of the numerical results are presented in a manner that allows simple rescaling to any revised density estimates.  相似文献   

10.
Plate tectonics on the Earth is a surface manifestation of convection within the Earth’s mantle, a subject which is as yet improperly understood, and it has motivated the study of various forms of buoyancy-driven thermal convection. The early success of the high Rayleigh number constant viscosity theory was later tempered by the absence of plate motion when the viscosity is more realistically strongly temperature dependent, and the process of subduction represents a continuing principal conundrum in the application of convection theory to the Earth. A similar problem appears to arise if the equally strong pressure dependence of viscosity is considered, since the classical isothermal core convection theory would then imply a strongly variable viscosity in the convective core, which is inconsistent with results from post-glacial rebound studies. In this paper we address the problem of determining the asymptotic structure of high Rayleigh number convection when the viscosity is strongly temperature and pressure dependent, i.e. thermobaroviscous. By a method akin to lid-stripping, we are able to extend numerical computations to extremely high viscosity contrasts, and we show that the convective cells take the form of narrow, vertically-oriented fingers. We are then able to determine the asymptotic structure of the solution, and it agrees well with the numerical results. Beneath a stagnant lid, there is a vigorous convection in the upper part of the cell, and a more sluggish, higher viscosity flow in the lower part of the cell. We then offer some comments on the possible meaning and interpretation of these results for planetary mantle convection.  相似文献   

11.
全球地表热流的产生与分布   总被引：7，自引：1，他引：7       下载免费PDF全文
全球地表热流是反映地球内部热与动力学过程的一种主要能流.本文在三维球坐标框架下，就几个不同的粘度模型分别研究地幔内部密度异常(基于全球地震层析结果)以及板块运动激发的地幔流动的热效应及其对于观测地表热流产生和分布特征的贡献.由于地幔动力系统具有较高的Pe数，可以期望由板块运动激发的地幔流动将强烈地扰动地幔内部初始传导状态下的温度场以及地表热的热流分布.结果表明，与地幔内部密度异常产生的热效应相比，运动的板块及其激发的地幔流动在全球地表观测热流的产生和分布特征上起着更为重要的作用.观测到的大洋中脊处的高热流在很大程度上可以归因于板块激发的地幔流动的热效应.计算的平均温度剖面较好地揭示了岩石圈和D″层的温度特征，即温度随深度的剧烈变化，这与我们目前通过其他手段对岩石圈和D″层的温度结构了解是一致的.一个下地幔粘度比上地幔高出30倍的粘度结构(文中使用的粘度模型2)较之其余模型的拟合程度似乎更好.  相似文献   

12.
Thermal history of Mimas and Enceladus is investigated from the beginning of accretion to 400 Myr. The numerical model of convection combined with the parameterized theory is used. The following heat sources are included: short lived and long lived radioactive isotopes, accretion, serpentinization, and phase changes. The heat transfer processes are: conduction, solid state convection, and liquid state convection. We find that temperature of Mimas’ interior was significantly lower than that of Enceladus. If Mimas accreted 1.8 Myr after CAI then the internal melting and differentiation did not occur at all. Comparison of thermal models of Mimas and Enceladus indicates that conditions favorable for the start of tidal heating lasted for a short time (~107 yr) in Mimas and for ~108 yr in Enceladus. This could explain the Mimas—Enceladus paradox. In fact, in view of the chronology based on cometary impact rate, one cannot discard a possibility that also Mimas was for some time active and it has the interior differentiated on porous core and icy mantle.  相似文献   

13.
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 102 to 103 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.  相似文献   

14.
The case is presented that the efficiency of variable viscosity convection in the Earth's mantle to remove heat may depend only very weakly on the internal viscosity or temperature. An extensive numerical study of the heat transport by 2-D steady state convection with free boundaries and temperature dependent viscosity was carried out. The range of Rayleigh numbers (Ra) is 104?107 and the viscosity contrast goes up to 250000. Although an absolute or relative maximum of the Nusselt number (Nu) is obtained at long wavelength in a certain parameter range, at sufficiently high Rayleigh number optimal heat transport is achieved by an aspect ratio close to or below one. The results for convection in a square box are presented in several ways. With the viscosity ratio fixed and the Rayleigh number defined with the viscosity at the mean of top and bottom temperature the increase of Nu with Ra is characterized by a logarithmic gradient β = ?ln(Nu)/? ln(Ra) in the range of 0.23–0.36, similar to constant viscosity convection. More appropriate for a cooling planetary body is a parameterization where the Rayleigh number is defined with the viscosity at the actual average temperature and the surface viscosity is fixed rather than the viscosity ratio. Now the logarithmic gradient β falls below 0.10 when the viscosity ratio exceeds 250, and the velocity of the surface layer becomes almost independent of Ra. In an end-member model for the Earth's thermal evolution it is assumed that the Nusselt number becomes virtually constant at high Rayleigh number. In the context of whole mantle convection this would imply that the present thermal state is still affected by the initial temperature, that only 25–50% of the present-day heat loss is balanced by radiogenic heat production, and the plate velocities were about the same during most of the Earth's history.  相似文献   

15.
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 32 resonance between Mercury's rotational and orbital periods, is most likely related to surface processes, as opposed to convection. The 3n2 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.  相似文献   

16.
地幔对流的数值模拟及其与表面观测的关系   总被引：9，自引：3，他引：9       下载免费PDF全文
本文从基本的热对流方程出发,并结合地幔对流特点,特别考虑到自重及非线性影响,探讨地幔对流及其与表面观测的关系,发展了相应的数值方法.结果表明,计算得到的长波大地水准面、地表地形、板块速度场水平散度与观测值符合程度较好.上、下地幔的非绝热温度异常与由地震层析得到的地震波速异常显示一定的相关性.地幔内部的流动呈现复杂形态,反映了高瑞利数对流的特征.  相似文献   

17.
张晨  张双喜 《地震学报》2014,36(5):872-882
热传导系数和热膨胀系数是影响板块俯冲动力学过程的两个重要参数. 由于地球介质的不均匀性,热系数也会随深度发生变化.然而,这种变化在地球动力学模拟研究中往往被忽略.本文针对随温度变化的热传导系数和热膨胀系数, 模拟板块俯冲的动力学过程,分析热系数、黏度对板块俯冲形态的影响及其对应的地幔对流特征.结果表明,依温度变化的热传导系数和热膨胀系数会影响地幔温度及黏度分布,进而改变板块的俯冲角度;黏度是控制板块俯冲动力学演化过程的重要因素;地幔对流受黏度结构的影响,呈现分层对流及局部多个对流环等多种不同形态的对流场特征.   相似文献   

18.
Recognition that the cooling of the core is accomplished by conduction of heat into a thermal boundary layer (D″) at the base of the mantle, partly decouples calculations of the thermal histories of the core and mantle. Both are controlled by the temperature-dependent rheology of the mantle, but in different ways. Thermal parameters of the Earth are more tightly constrained than hitherto by demanding that they satisfy both core and mantle histories. We require evolution from an early state, in which the temperatures of the top of the core and the base of the mantle were both very close to the mantle solidus, to the present state in which a temperature increment, estimated to be ～ 800 K, has developed across D″. The thermal history is not very dependent upon the assumption of Newtonian or non-Newtonian mantle rheology. The thermal boundary layer at the base of the mantle (i.e., D″) developed within the first few hundred million years and the temperature increment across it is still increasing slowly. In our preferred model the present temperature at the top of the core is 3800 K and the mantle temperature, extrapolated to the core boundary without the thermal boundary layer, is 3000 K. The mantle solidus is 3860 K. These temperatures could be varied within quite wide limits without seriously affecting our conclusions. Core gravitational energy release is found to have been remarkably constant at ～ 3 × 1011 W. nearly 20% of the core heat flux, for the past 3 × 109 y, although the total terrestrial heat flux has decreased by a factor of 2 or 3 in that time. This gravitational energy can power the “chemical” dynamo in spite of a core heat flux that is less than that required by conduction down an adiabatic gradient in the outer core; part of the gravitational energy is used to redistribute the excess heat back into the core, leaving 1.8 × 1011 W to drive the dynamo. At no time was the dynamo thermally driven and the present radioactive heating in the core is negligibly small. The dynamo can persist indefinitely into the future; available power 1010 y from now is estimated to be 0.3 × 1011 W if linear mantle rheology is assumed or more if mantle rheology is non-linear. The assumption that the gravitational constant decreases with time imposes an implausible rate of decrease in dynamo energy. With conventional thermodynamics it also requires radiogenic heating of the mantle considerably in excess of the likely content of radioactive elements.  相似文献   

19.
The mantle convection model with phase transitions, non-Newtonian viscosity, and internal heat sources is calculated for two-dimensional (2D) Cartesian geometry. The temperature dependence of viscosity is described by the Arrhenius law with a viscosity step of 50 at the boundary between the upper and lower mantle. The viscosity in the model ranges within 4.5 orders of magnitude. The use of the non-Newtonian rheology enabled us to model the processes of softening in the zone of bending and subduction of the oceanic plates. The yield stress in the model is assumed to be 50 MPa. Based on the obtained model, the structure of the mantle flows and the spatial fields of the stresses σxz and σxx in the Earth’s mantle are studied. The model demonstrates a stepwise migration of the subduction zones and reveals the sharp changes in the stress fields depending on the stage of the slab detachment. In contrast to the previous model (Bobrov and Baranov, 2014), the self-consistent appearance of the rigid moving lithospheric plates on the surface is observed. Here, the intense flows in the upper mantle cause the drift and bending of the top segments of the slabs and the displacement of the plumes. It is established that when the upwelling plume intersects the boundary between the lower and upper mantle, it assumes a characteristic two-level structure: in the upper mantle, the ascending jet of the mantle material gets thinner, whereas its velocity increases. This effect is caused by the jump in the viscosity at the boundary and is enhanced by the effect of the endothermic phase boundary which impedes the penetration of the plume material from the lower mantle to the upper mantle. The values and distribution of the shear stresses σxz and superlithostatic horizontal stresses σxx are calculated. In the model area of the subducting slabs the stresses are 60–80 MPa, which is by about an order of magnitude higher than in the other mantle regions. The character of the stress fields in the transition region of the phase boundaries and viscosity step by the plumes and slabs is analyzed. It is established that the viscosity step and endothermic phase boundary at a depth of 660 km induce heterogeneities in the stress fields at the upper/lower mantle boundary. With the assumed model parameters, the exothermic phase transition at 410 km barely affects the stress fields. The slab regions manifest themselves in the stress fields much stronger than the plume regions. This numerically demonstrates that it is the slabs, not the plumes that are the main drivers of the convection. The plumes partly drive the convection and are partly passively involved into the convection stirred by the sinking slabs.  相似文献   

20.
A new approach to analytical and numerical study of the process of the post-glacial uplifting of the Earth’s surface was proposed within the framework of a viscous model. Displacement of the Earth’s surface is considered as the motion of the density boundary due to chemico-density convection. It is shown that the incorporation of the non-Newtonian rheology at observed velocities of post-glacial uplifts requires an obligatory presence of faults in the lithosphere and gives rise to quasi-uniform motion of the mantle material, whose viscosity under the lithosphere is, on the average, sufficiently small and amounts to ～1019 Pa. The study of the stability of the constructed model of the post-glacial uplift considered as the chemico-density convection relative to the thermal convection shows that the velocity of thermal convection developing in the presence of a quasiuniform mantle flow related to the post-glacial recovery is ～1 m/yr.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号