共查询到20条相似文献，搜索用时 796 毫秒
1.
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. 相似文献
2.
Hans-Peter Plag Bjrn Engen Thomas A. Clark John J. Degnan Bernd Richter 《Journal of Geodynamics》1998,25(3-4)
The WEGENER† activities related to the study of post-glacial rebound are presented together with a review of the present state-of-the-art in this study field. Post-glacial rebound research is an unique tool for studying the viscoelastic behaviour of the Earth's mantle on time scales of thousands of years. The viscosity structure of the Earth's mantle determined from an inversion of observations of glacially induced deformations is a basic requirement for modelling long-term phenomena such as the convection in the Earth's mantle, and for better understanding unsolved questions such as the nature of the mantle discontinuities or the vertical scale of convection.First, an introduction to the scientific background is given, and the three principal ingredients for post-glacial rebound studies, namely the ice model, the Earth model, and the observations are briefly considered. For the ice model, the uncertainties due to a trade-off between ice model and Earth rheology are outlined. The different approaches used to model the Earth and its deformations in post-glacial rebound studies are discussed emphasising the preliminary nature of the derived rheologies and depth dependencies. The observations, in particular the relative sea-level changes and three-dimensional surface deformations, are described with special emphasis on observational gaps. Based on the discussion of the ingredients, an outline of the future developments in post-glacial rebound research is attempted with particular emphasis on the Earth model and the theory of deformations.For several decades extreme efforts have been made to precisely monitor the land uplift in Scandinavia. However, for the height component the existing data still are associated with large uncertainties while reliable data on the horizontal component are practically nil. The ongoing long-term (longer than ten years) spacegeodetic measurements are likely to provide the three-dimensional deformations with the spatial resolution and accuracy required in order to substantially contribute to post-glacial rebound studies. Thus, present-day three-dimensional deformations of the Earth's surface beneath and around the former ice sheets as a constraint for the mantle rheology and viscosity structure will increasingly become important as they become known from space-geodetic measurements with high spatial resolution and an accuracy approaching the mm/a-level. 相似文献
3.
This study considers two-dimensional mantle flow beneath a rigid lithosphere. The lithosphere which forms the upper boundary of a convecting region moves with a prescribed uniform horizontal velocity, and thickens with distance from the accreting plate boundary as it cools. Beneath the lithosphere, the mantle deforms viscously by diffusion creep and is heated radiogenically from within. Solutions for thermal convection beneath the lithosphere are obtained by finite-difference methods. Two important conclusions have resulted from this study: (1) convective patterns of large aspect ratio are stable beneath a rigid moving lithosphere; (2) even for a lithosphere velocity as small as 3 cm/yr. and a Rayleigh number as large as 106, mantle circulation with large aspect ratio is driven dominantly by the motion of the lithosphere rather than by temperature gradients within the flow. Gravity, topography and heat flow are determined and implications for convection in the upper mantle are discussed. 相似文献
4.
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 resonance between Mercury's rotational and orbital periods, is most likely related to surface processes, as opposed to convection. The 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. 相似文献
5.
In our previous works, based on numerical models, it was shown that under certain conditions a hot material can rise in portions in the tails of thermal mantle plumes. The spectrum of these pulsations can correspond to the observed spectra of catastrophic hotspot eruptions. Since most of the existing numerical models of thermal convection for the mantle of the present Earth do not reveal these pulsations, in this work, we analyze the physical cause and initiation conditions of pulsations of thermal plumes. The results of a numerical solution of the thermal convection equations for a material with varying parameters in the extended Boussinesq approximation are presented. It is shown how the structure of the convection is transformed with the increase of convection intensity. At the Rayleigh numbers Ra > 106, convection becomes unsteady, and the configuration of the ascending and descending flows changes. The new flow emerging at the mantle bottom acquires a mushroom shape with a head and a tail. After the rise of the plume’s head to the surface, the tail remains in the mantle in the form of a quasi-stationary hot steam. It turns out that at Ra ~ 5 × 107, the thermal mantle plume becomes pulsating and its tail is in fact a heated channel through which the hot material rises in successive portions. At the Rayleigh numbers Ra > 5 × 108, the tail of the thermal plume breaks and the plume becomes a regular conveyor of separate ascending portions of the hot material, which are referred to as thermals. Thus, thermal convection with pulsating plumes takes place at the transitional stage from the regime of quasi-stationary plumes to the regime of thermals. 相似文献
6.
A wide class of equations is defined for a high pressure and subcritical temperature range of a fluid state whose thermodynamic properties enable the construction of a polytropic model of the mantle. A variant of deep convection equations of the Ogura and Phillips type is substantiated in terms of the polytropic mantle model. The proposed system of the deep convection equations includes fluctuation of the generalized potential temperature, has a quasi-incompressible form, and is transformed into Mihaljan’s system of shallow convection equations with a decrease in the layer depth. This circumstance is of great importance because it validates the use of the same dimensionless parameters as in the shallow convection model. The advantage of the proposed variant of the deep convection equations is its complete conservatism, which allows one to gain constraints on the efficiency of energy conversion in deep mantle processes and the thermal energy power expended on the generation rate of the convection kinetic energy and associated processes. This power is shown to be of the order of half the geothermal flux measured on the Earth’s surface. 相似文献
7.
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. 相似文献
8.
采用较为符合实际岩石层变形的非线性幂指数本构关系,基于ANSYS有限元平台, 模拟了近20万年来中国大陆地区地表运动及演化过程,探讨了印度板块挤压作用和地幔对流拖曳力各自对于中国大陆地区地表形变运动格局的影响.模拟结果与观测数据的比较表明：在印度板块的挤压和地幔拖曳力联合作用下,中国及东亚大陆岩石层运动形变模式能够和现代GPS观测有较好的吻合; 印度大陆和欧亚大陆的碰撞以及印度大陆的持续向北推进、挤压所产生的应力环境,一直主导了以青藏高原为核心的我国西部地域岩石圈构造、运动和演化,但其影响随着远离青藏高原地区而逐渐变小;地幔对流产生的作用于岩石层底部的拖曳力是中国大陆（特别是远离碰撞带）岩石层运动构造变形的重要驱动力.然而在构造复杂和东部靠近太平洋板块的区域,模型预测结果和GPS观测还存在一定的差距,这说明在未来的中国大陆岩石层变形运动的数值模拟中,应当采用更为复杂的构造模型和驱动力因素. 相似文献
9.
A model is proposed for the origin of hot spots that depends on the existence of major-element heterogeneities in the mantle. Generation of basaltic crust at spreading centers produces a layer of residual peridotite ～20–25 km thick directly beneath the crust which is depleted in Fe/Mg, TiO2, CaO, Al2O3, Na2O and K2O, and which has a slightly lower density than undepleted peridotite beneath it. Upon recycling of this depleted peridotite back into the deep mantle at subduction zones, it becomes gravitationally unstable, and tends to rise as diapirs through undepleted peridotite. For a density contrast of 0.05 g cm?3, a diapir 60 km in diameter would rise at roughly 8 cm y?1, and could transport enough heat to the base of the lithosphere to cause melting and volcanism at the surface. Hot spots are thus viewed as a passive consequence of mantle convection and fractionation at spreading centers rather than a plate-driving force.It is suggested that depleted diapirs exist with varying amounts of depletion, diameters, upward velocities and source volumes. Such variations could explain the occurrence of hot spots with widely varying lifetimes and rates of lava production. For highly depleted diapirs with very low Fe/Mg, the diapir would act as a heat source and the asthenosphere and lower lithosphere drifting across the diapir would serve as the source region of magmas erupted at the surface. For mildly depleted diapirs with Fe/Mg only slightly less than in normal undepleted mantle, the diapir could provide not only the source of heat but also most or all of the source material for the erupted magmas. The model is consistent with isotopic data that require two separate and ancient source regions for mid-ocean ridge and oceanic island basalts. The source for mid-ocean ridge basalts is considered to be material upwelling at spreading centers from the deep mantle. This material forms the oceanic lithosphere. Oceanic island basalts are considered to be derived from varying mixtures of sublithospheric and lower lithospheric material and the rising diapir itself. 相似文献
10.
V. P. Trubitsyn 《Izvestiya Physics of the Solid Earth》2008,44(11):857-872
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. 相似文献
11.
The geological-geophysical and petrological-geochemical studies of the Earth’s crust and upper mantle are combined to estimate the state of the lithosphere at the junction zone of Tarim and Tien Shan. The laboratory measurements of electric conductivity in the rocks sampled from the upper mantle and lower crust considered against the geoelectrical and thermal models revealed lherzolite, granulite, and eclogite massifs in the deep section of the Tarim and Tien Shan junction zone. The experimental results suggest that the crustal thickness in the southern Tien Shan attained 35–40 km 70 Ma ago. 相似文献
12.
13.
Based on the recent observations about the movement and rheological structure of the lithosphere and deformation pattern of the crust, we developed a three-dimensional finite element model for the northeastern margin of the Tibetan Plateau. The model considered the impacts of both external and internal conditions, including mantle convection, gravitational potential energy and block interactions. We compared the simulated surface movement rates to the observed GPS velocities, and the results revealed that crustal movement gradually decreased toward the edge of the plateau. The factors controlling this pattern are the interactions of adjacent blocks, gravitational potential energy of the plateau, and also mantle convection as well. Additionally, according to the observation that there was an apparent difference between the horizontal movement rate of the lithosphere and convective velocity of the underlying mantle, and also based on the results of seismic anisotropy studies that suggest different strengths and deformation regimes of the lithosphere in different tectonic blocks, we proposed that the impact of mantle convection on the lithosphere may have varied in space, and introduced a parameter named mantle convection intensity factor in numerical simulations. Our simulation results show consistent surface movement rates with GPS observations, which further supports the viewpoint of seismic anisotropy studies, i.e., the degree of coupling between the crust and mantle varies significantly among different blocks. 相似文献
14.
对于大地测量应用来说,目前IERS机构在定义地球参考系时推荐采用岩石圈无整体旋转(No-Net-Rotation-NNR)约束条件,然而对于地球物理应用来说,相对于NNR参考基准的绝对板块运动数据可能会对地幔对流等研究结果产生误导.考虑到热点的运动,提出建立平均热点(MHS-Medial HotSpot)参考基准的方法,给出建立该基准的约束准则,分别以地学模型NNR-NUVEL1A和实测模型ITRF2005VEL为基础,建立了平均热点参考基准MHS-NUVEL1A和MHS-ITRF2005,并与其它基于热点的绝对板块运动模型进行了比较和分析;讨论了岩石圈的西向漂移,给出了岩石圈相对于下地幔整体旋转的更精确的定量估计,即基于实测的热点参考架MHS-ITRF2005和地学模型NNR-NUVEL1A之间的整体旋转为0.26°/Ma,旋转极在(50°S, 62°E),这与由板块的受力模型给出的岩石圈的整体旋转的旋转极很接近,旋转速率大致快了10%. 相似文献
15.
16.
利用地热学、流变学和重力学方法，计算了南海岩石层温度结构、流变特征及地幔对流格局.南海莫霍面温度在600-1000℃之间.岩石层底界面温度在1150-1300℃之间，有效粘滞系数为1020-1021Pa·s，与冰期回弹资料确定的地幔粘度吻合，表明南海深部具备产生地幔热对流的物理条件.研究认为地幔物质由北西向南东方向的运移以及印澳-欧亚板块的碰撞，导致南海北部大陆边缘向洋扩张、离散和断裂解体.在向洋离散过程中，陆-洋岩石层底部地幔局部对流使中央海盆扩张和北部陆缘发生差异性块断运动. 相似文献
17.
In this work we have modeled the thermal structure of the lithosphere of the Spanish Central System and the Tajo Basin, and their implications for lithospheric strength. For this, we have used refined heat-producing elements (HPE) values to obtain new estimates of heat production rates in the Spanish Central System and Tajo Basin areas, which have been used joined to the relation between topography and thermal structure of the lithosphere to calculate the best-fit surface heat flows in the study area. Moreover, we have implemented a temperature-dependent thermal conductivity (appropriate for olivine) for the lithospheric mantle to improve the calculations of temperature profiles in the mantle. The geotherms so obtained, together with the implementation of a new rheological law for the upper lithospheric mantle, have been used to calculate refined estimations of the strength and effective elastic thickness of the lithosphere. We have obtained surface heat flow values of 84 mW m−2 and ∼82 mW m−2 for the Spanish Central System and the Tajo Basin, respectively. The thermal state of the lithosphere affects mantle temperatures, and hence may be playing an important role in the uplift and maintenance of the Spanish Central System. 相似文献
18.
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. 相似文献
19.
S. M. Molodenskii 《Izvestiya Physics of the Solid Earth》2017,53(2):173-180
The question of ambiguity in the solution of the inverse problem for determining the Brünt-Väisäla frequency in the Earth’s mantle from the entire set of the up-to-date data on seismicity, free oscillations, and forced nutations of the Earth, as well as the data on the Earth’s total mass and total moment of inertia, is considered. Based on the results of a series of numerical experiments, the band of admissible distributions of the Brünt-Väisäla frequency and mantle density with depth is calculated. This estimate is used for investigating the convective and gravitational stability of the different regions of the mantle against relatively small adiabatic and nonadiabatic perturbations. The generalization of the known Rayleigh criterion of convective stability of homogeneous and a nonself-gravitating incompressible viscous fluid for the case of a compressible self-gravitating fluid is given. A system of the ordinary eight-order differential equations with complex coefficients and homogeneous boundary conditions, whose eigenvalues determine the transition from the stable state to instability, is obtained. Examples of the numerical determination of these eignevalues are presented. For interpreting the data about the band of the admissible distributions of the Brünt-Väisäla frequency with depth, the notion of the effective bulk modulus of the medium at different depths is introduced. This quantity governs the depth changes in temperature in a convecting mantle and allows us to make a conclusion about the role of heat conduction and the radial heterogeneity of the mantle composition without imposing any constraints on the convection mechanism. It is shown that within the present-day observation errors in the frequencies of the Earth’s free oscillations, the simplest reasonable model is that in which the ratio of the effective bulk modulus to its adiabatic value in the lower and middle mantle is 1.043 ± 0.05. The closeness of this value to unity indicates that convection in the lower and middle mantle is fairly close to adiabatic. At the same time, when the analysis only relies on seismic data and on the information about the periods of the free oscillations of the Earth, there is a significant uncertainty in the models of the effective bulk modulus distribution in the upper mantle and crust. This uncertainty precludes us from making purely empirically derived conclusions that reliably and unambiguously describe the role of the effects of heat conduction and radially heterogeneous material composition in the convection in the upper mantle. 相似文献
20.
V. P. Trubitsyn 《Izvestiya Physics of the Solid Earth》2012,48(6):467-485
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. 相似文献