首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
A temperature and pressure increase in the mantle causes phase transitions and related density changes in its material. The transition boundary in the pressure-temperature phase diagram is determined by the curve of phase equilibrium with the slope γ = dp/dT. If the slope is nonzero, a phase transition in hot ascending and cold descending mantle flows occurs at different depths and, therefore, either enhances (γ > 0) or slows down convection (γ < 0). The mantle material has a multicomponent composition. Therefore, phase transitions in the mantle are distributed over an interval of pressures and depths. In this interval, the concentration of one phase smoothly decreases and the concentration of the other increases. The widths of phase transition zones in the Earth’s mantle vary from 3 km for the endothermic transition in olivine at a depth of 660 km to 500 km for the exothermic transition in perovskite, and the high-to-low spin change in the atomic state of iron takes place at a depth of about 1500 km. This work presents results of calculations demonstrating the convection effect of phase transitions as a function of the transition zone width. Transitions of both types with different slopes of the phase curve and different intensities of mantle convection are examined. For the first time, the convection enhancement and an increase in the mass transfer across the phase boundary are quantitatively investigated in the presence of an exothermic phase transition as a function of the slope of the phase curve. The mixing of material under conditions of partially layered convection is examined with the help of markers.  相似文献   

2.
The strong control that the endothermic phase change from spinel to perovskite and magnesiowüstite at a depth of 660 km has on mantle convection is discussed. The phase transition determines the morphology and length scales of upflow and downflow structures and, through retardation of sinking slabs, can cause an avalanche phenomenon involving rapid flushing of cold upper mantle material down to the base of the lower mantle. The phase change significantly heats plumes that rise from the lower mantle and penetrate into the upper mantle. The exothermic phase change from olivine to spinel at a depth of 400 km in the mantle mitigates the effects of the dynamically and thermally dominant endothermic phase transition.  相似文献   

3.
Geochemical data indicate that two major reservoirs 1–2 Ga in age are present in the mantle. The upper mantle, feeding mid-ocean ridges, is depleted in chemical elements carried away into the continental crust. The lower mantle, feeding hotspot plumes, is close in composition to primordial matter. The 660-km depth of an endothermic phase transition in olivine has been considered over the last two decades as a possible boundary between the reservoirs. In this period, many models of mantle convection were constructed that used values of the phase transition parameters, which led to temporal (up to 1 Gyr long) convection layerings and periodic avalanche-induced mantle intermixing events. However, laboratory measurements with new improved instrumentation give other values of the phase transition parameters that require a revision of the majority of the existence of large-scale avalanches in the Earth’s history becomes disputable. The paper is devoted to comprehensive study of the phase transition effect on the structure of mantle flows with different values of phase transition parameters and Rayleigh numbers; in particular, the mass transfer through the phase boundary is calculated for different regimes of steady-state convection.  相似文献   

4.
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.  相似文献   

5.
The thermochemical model of the authors is shown to be naturally related to the general theory of V.P. Myasnikov. A heterogeneous modification of this homogeneous theory is described in light of the present ideas on the differentiation of the mantle substance at the boundary with the core and its eclogitization during submersion from the outer boundary and at the endothermic phase transition at a depth of 670 km. The Earth’ls evolution from an initial hot state is numerically modeled. The evolution is shown to start with an abrupt mantle overturn followed by a long period of steady evolution. Global mantle overturns recur a few times, gradually weaken, and are transformed into regional avalanches. The spatial configuration of overturns is represented by a predominant funnel-shaped sink and a few (three to five) ascending superplumes, which convincingly explains the causes of the formation of supercontinents, the opening of oceans, and the observed asymmetry of the planet. The times of overturns remarkably correlate with geological data on the existence of supercontinents. The processes of core growth, mantle cooling, and crust formation exhibit a clearly expressed stepwise behavior. The supplementation of the endothermic phase transition by chemical transformations favors the overcoming of the phase barrier between the upper and lower mantle, enhances the nonlinearity of mantle convection, and imparts a heterocyclic pattern to the process of evolution. It is shown that the lower mantle plume of chemical origin is fragmented by the phase transition into parts that, interacting with the thermal convection, generate a system of upper mantle plumes. This modeling provides an explanation of the coeval systems of oceanic plateaus and continental traps observed on the surface.  相似文献   

6.
The new post-perovskite phase near the core-mantle boundary has important ramifications on lower mantle dynamics. We have investigated the dynamical impact arising from the interaction of temperature- and depth-dependent viscosity with radiative thermal conductivity, up to a lateral viscosity contrast of 104, on both the ascending and descending flows in the presence of both the endothermic phase change at 670 km depth and an exothermic post-perovskite transition at 2650 km depth. The phase boundaries are approximated as localized zones. We have employed a two-dimensional Cartesian model, using a box with an aspect-ratio of 10, within the framework of the extended Boussinesq approximation. Our results for temperature- and depth-dependent viscosity corroborate the previous results for depth-dependent viscosity in that a sufficiently strong radiative thermal conductivity plays an important role for sustaining superplumes in the lower mantle, once the post-perovskite phase change is brought into play. This aspect is especially emphasized, when the radiative thermal conductivity is restricted only to the post-perovskite phase. These results revealed a greater degree of asymmetry is produced in the vertical flow structures of the mantle by the phase transitions. Mass and heat transfer between the upper and lower mantle will deviate substantially from the traditional whole-mantle convection model. Streamlines revealed that an overall complete communication between the top and lower mantle is difficult to be achieved.  相似文献   

7.
The temperature gradient in the lower mantle is fundamental in prescribing many transport properties, such as the viscosity, thermal conductivity and electrical conductivity. The adiabatic temperature gradient is commonly employed for estimating these transport properties in the lower mantle. We have carried out a series of high-resolution 3-D anelastic compressible convections in a spherical shell with the PREM seismic model as the background density and bulk modulus and the thermal expansivity decreasing with depth. Our purpose was to assess how close under realistic conditions the horizontally averaged thermal gradient would lie to the adiabatic gradient derived from the convection model. These models all have an endothermic phase change at 660 km depth with a Clapeyron slope of around −3 MPa K−1, uniform internal heating and a viscosity increase of 30 across the phase transition. The global Rayleigh number for basal heating is around 2×106, while an internal heating Rayleigh number as high as 108 has been employed. The pattern of convection is generally partially layered with a jump of the geotherm across the phase change of at most 300 K. In all thermally equilibrated situations the geothermal gradients in the lower mantle are small, around 0.1 K km−1, and are subadiabatic. Such a low gradient would produce a high peak in the lower-mantle viscosity, if the temperature is substituted into a recently proposed rheological law in the lower mantle. Although the endothermic phase transition may only cause partial layering in the present-day mantle, its presence can exert a profound influence on the state of adiabaticity over the entire mantle.  相似文献   

8.
The process of multiple self-nucleation and ascent of mantle plumes is studied in the numerical models of thermal convection. The plumes are observed even in the simplest isoviscous models of thermal convection that leave aside the more complex rheology of the material, thermochemical effects, phase transformations, etc., which, although controlling the features of plumes, are not necessary for their formation. The origin of plumes is mainly due to the instability of the mantle flows at highly intense (low-viscous) thermal convection. At high viscosity, convective flows form regular cells. As viscosity decreases, the ascending and descending flows become narrower and unsteady. At a further decrease in viscosity, the ascending plumes assume a mushroom-like shape and occasionally change their position in the mantle. The lifetime of each flow can attain 100 Ma. Using markers allows visualizing the evolution of the shape of the mantle plumes.  相似文献   

9.
A two-dimensional numerical convection model in cartesian geometry is used to study the influence of trench migration on the ability of subducted slabs to penetrate an endothermic phase boundary at 660 km depth. The transient subduction history of an oceanic plate is modelled by imposing plate and trench motion at the surface. The viscosity depends on temperature and depth. A variety of styles of slab behaviour is found, depending predominantly on the trench velocity. When trench retreat is faster than 2–4 cm/a, the descending slab flattens above the phase boundary. At slower rates it penetrates straight into the lower mantle, although flattening in the transition zone may occur later, leading to a complex slab morphology. The slab can buckle, independent of whether it penetrates or not, especially when there is a localised increase in viscosity at the phase boundary. Flattened slabs are only temporarily arrested in the transition zone and sink ultimately into the lower mantle. The results offer a framework for understanding the variety in slab geometry revealed by seismic tomography.  相似文献   

10.
The upper mantle flow beneath the North China Platform   总被引：2，自引：0，他引：2  
In this paper we establish an upper mantle convection model which is constrained by regional isostatic gravity anomalies. Comparing the computed convection patterns with the tectonic features of the North China Platform we find that there are two positive anomaly centers connected with upward flows. These anomalies belong to the tectonic units of the Shan-Xi geoanticline and the Lu-Xi geoanticline. The centers of downward flows are connected with the tectonic units of the Liao-Ji geosyncline. It is reasonable to suggest that the upward mantle flows push the lithosphere upward and generate the observed positive isostatic gravity anomaly. The downward mantle flows pull the lithosphere down and generate the negative anomaly. However, the use of simple analysis makes it difficult to explain the complex lithospheric dynamics of this region. In order to understand lithospheric structures and tectonic features we must investigate the mechanical properties of the lithosphere and the relationship between the lithosphere and the mantle. These problems are discussed in the last section of this paper.  相似文献   

11.
This paper presents a study on the effects of phase transitions on the mantle convection of Venus in a three-dimensional(3D)spherical shell domain.Our model includes strong depth-and temperature-dependent viscosity and exothermic phase change from olivine to spinel as well as endothermic phase change from spinel to perovskite.From extensive numerical simulations of the effects of Rayleigh number(Ra),and the Clapeyron slopes and depths of phase changes,we found the following:(1)The endothermic phase change prevents mass flow through the interface.Increasing the absolute value of the Clapeyron slopes decreases radial mass flux and normalized radial mass flux at the endothermic phase boundary,and decreases the number of mantle plumes.In other words,mass flow through the phase boundary decreases.The inhibition influence of phase changes increases,as do convective wavelengths.(2)Increasing Ra also increases the convective wavelength and decreases the number of mantle plumes,but it has less influence on the mass exchange.As Ra increases,the convective vigor increases along with the radial mass flux and the mass flow through the phase boundary;however,the normalized mass flux through the phase boundary varies little with Ra,which is different from the conclusion that increasing Ra will greatly increase the inhibition of mass flow through the phase boundary based on two-dimensional(2D)modeling.(3)Increasing the depth of endothermic phase change will slightly decrease the number of mantle plumes,but has little effect on the mass flow through the phase boundary.Consistent with previous studies,our results show that the phase change from spinel to perovskite could inhibit the mass flow through the phase boundary,but they also show that the buildup of hot materials under the endothermic phase boundary in the 3D model could not be so large as to cause strong episodic overturns of mantle materials,which is quite different from previous 2D studies.Our results suggest that it is difficult for phase changes to cause significant magmatism on Venus;in other words,phase changes may not be the primary cause of catastrophic resurfacing on Venus.  相似文献   

12.
The Runcorn stress equations and 2–30° harmonic coefficients of the geopotential have been applied to determine the mantle convection pattern beneath China. The pattern is compared with geophysical and geological observations and it is found that the directional change belts of mantle flows coincide with the major fault belts between tectonic units of China. The stress field generated by mantle flows, except in the Tian Shan region, also coincide with the stress field of recent tectonic movement in China. The Tarim and Junggar basins are formed by tensional stresses due to divergent mantle convection currents under northwest China. The formation of the Qinhai-Xizang (Tibet) plateau is due mainly to the compression of the Tarim block and Indian plate, caused by convergent mantle convection currents. The shear-fault belts in central China (100–105°E) are generated by the running change belt of mantle flows, a well-known N-S seismic zone. In eastern China, tensional faults, grabens, lake and sea depressions are related to the eastward displacement of continental lithosphere exerted by eastward dispersal mantle flows under this region.This paper provides new material for further study of the force source mechanism of recent tectonic movement from the viewpoint of mantle convection currents.  相似文献   

13.
In the PREM seismic model, the boundary between the upper and the lower mantle is accepted at a depth of 670 km, where seismic velocities and density increase. However, until recently there was an obvious inconsistency in this model. The density increases abruptly, and the velocities, in addition to the jumps, have also the subsequent zones of increased gradient. The discontinuity between the upper and the lower mantle is related to the transition of olivine from the ringwoodite phase into the mixture of perovskite and magnesiowustite. However, in the pyrolyte model, the transition zone of the upper mantle consists not wholly of olivine, but partly of olivine (60%) and partly of garnet (40%). The latest data of the garnet measurement at high pressures show that it also experiences phase transition, being converted into magnesium perovskite with the impurity of calcium perovskite. In contrast to the sharp transition in olivine (within a depth interval of only 5 km), the transition in garnet is spread over the interval of depths of 660–710 km. In the widely used PREM and AK135 models, this additional transition corresponds to the zone of the increased gradient in seismic velocities, while in the density distribution it is included in the sharp transition of ringwoodite. Thus, the mineralogy data indicate the need for correction of the PREM and AK135 seismic models: the density jump at a depth of 660 km should be reduced by approximately a factor of two, and a subjacent layer with the increased density gradient should be added at the depth interval of 660–710 km. The phase transition in olivine hampers the mantle flows, although in garnet it accelerates them. Therefore, with an allowance for the smaller jump in density, the decelerating effect of the subducting plates, caused by the phase transition in olivine, decreases, and, furthermore, the effect of their acceleration, caused by the phase transition in garnet, is added. The decrease in the density jump by almost a factor of two will lead to essential changes in the results of the majority of recent works addressing the assessment of the deceleration of convection at the upper/lower mantle discontinuity on the basis of the PREM model.  相似文献   

14.
—Recently a high-resolution tomographic model, the P1200, based on P-wave travel times was developed, which allowed for detailed imaging of the top 1200 km of the mantle. This model was used in diverse ways to study mantle viscosity structure and geodynamical processes. In the spatial domain there are lateral variations in the transition zone, suggesting interaction between the lower-mantle plumes and the region from 600 km to 1000 km. Some examples shown here include the continental region underneath Manchuria, Ukraine and South Africa, where horizontal structures lie above or below the 660 km discontinuity. The blockage of upwelling is observed under central Africa and the interaction between the upwelling and the transition zone under the slow Icelandic region appears to be complex. An expansion of the aspherical seismic velocities has been taken out to spherical harmonics of degree 60. For degrees exceeding around 10, the spectra at various depths decay with a power-law like dependence on the degree, with the logarithmic slopes in the asymptotic portion of the spectra containing values between 2 and 2.6. These spectral results may suggest the time-dependent nature of mantle convection. Details of the viscosity structure in the top 1200 km of the mantle have been inferred both from global and regional geoid data and from the high-resolution tomographic model. We have considered only the intermediate degrees (l = 12–25) in the nonlinear inversion with a genetic algorithm approach. Several families of acceptable viscosity profiles are found for both oceanic and global data. The families of solutions for the two data sets have different characteristics. Most of the solutions asociated with the global geoid data show the presence of asthenosphere below the lithosphere. In other families a low viscosity zone between 400 and 600 km depth is found to lie atop a viscosity jump. Other families evidence a viscosity decrease across the 660 km discontinuity. Solutions from oceanic geoid show basically two low viscosity zones one lying right below the lithosphere; the other right under 660-km depth. All of these results bespeak clearly the plausible existence of strong vertical viscosity stratification in the top 1000 km of the mantle. The presence of the second asthenosphere may have important dynamical ramifications on issues pertaining to layered mantle convection. Numerical modelling of mantle convection with two phase transitions and a realistic temperature- and pressure-dependent viscosity demonstrates that a low viscosity region under the endothermic phase transition can indeed be generated self-consistently in time-dependent situations involving a partially layered configuration in an axisymmetric spherical-shell model.  相似文献   

15.
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.  相似文献   

16.
Gravitational field models derived from satellite tracking and surface gravity data have been used to derive the forces in the earth's mantle under Asia. Based on studies of tectonic forces from these models, a subcrustal stress field under China has been obtained. The stresses are due to mantle convection. According to the stress patterns, the east and west China blocks and five seismic zones are identified. The tensional stresses exerted by the upwelling mantle convection flows under the crust of Tibet seem to be related to the Tibetan uplift. The compressional orogenic region from the southern tip of Lake Baikal, through Tien Shan, Hindu Kush and the Himalayas to northern Burma appears to be connected with the downwelling mantle convection flows. It is found that the directions of the subcrustal stresses under China are disposed perpendicularly to the major fault systems and seismic belts. The results of stress calculations show that the crust of north China should be in compression and that stresses within it should be sufficient to form the Shansi Graben and Linfen Basin Systems and fracture the lithosphere. This gives a possible explanation of why strong earthquakes occurred in north China which is an isolated narrow region of highest seismicity far from plate boundaries. The tensional stress fields, caused by the upwelling mantle convection flows, are found to be regions of structural kinship characterized by major concentrations of mineral and metal deposits in China.  相似文献   

17.
假设地震层析成像提供的地震波速异常对应于上地幔物质的密度异常分布,而该密度异常直接源于上地幔热对流相应的温度扰动. 在给定边界条件下,利用三维傅里叶变换,在波数域内求解控制流体行为的运动方程和连续性方程,得到上地幔小尺度对流流场. 利用密度异常驱动上地幔小尺度对流的数学 物理模型,采用胥颐、刘福田等提供的地震层析成像数据计算得到了我国西北及周边地区上地幔对流模式. 结果表明,对流流场的顶部在岩石圈较薄的盆地区域呈现上升发散流动特征,如塔里木盆地、柴达木盆地、哈萨克斯坦块体及准噶尔盆地;岩石圈较厚的山脉则对应了会聚下降的流动特征,如天山山脉、昆仑山山脉和祁连山山脉. 同时,塔里木盆地处于拉张状态,驱动其上地幔物质南下向青藏高原北部西昆仑运动,以及北上向天山下部流动,这可能是天山隆升的原因之一.   相似文献   

18.
金星表层年龄和构造活动特点表明其岩石层在最近的地质历史时期经历过广泛的更新.这种全球性的表层改造与其内部热演化历史进程密切相关.如果金星存在相变形成的上、下地幔,依据现今所了解的金星物理性质和参量化的热对流理论,并且考虑金星地幔相变边界层状态对对流的控制作用,我们计算了金星热演化历史.结果表明,金星的热演化历史是一种非单调的冷却过程,在这种非单调的热演化历史进程中,金星地幔会出现大体等周期的翻转.由于参数选取的不同,翻转时金星上、下地幔的温差随时间可能出现稳定变化、逐渐加强、逐渐减弱三种不同演化模式,目前尚不能确定实际金星热演化历史究竟是哪一种模式.金星地幔相变边界层的穿透对流可能是推动其表层岩石层全球性更新的关键,导致其表层火山活动和地表构造以大致500 Ma时间间隔更新和重造.  相似文献   

19.
Introduction The northwest of China includes Tarim, Junggar and Qaidam basins, and Kunlun, Tianshan, Altun and Qilian mountains, as well as the north part of the Tibetan Plateau. For a long time, the study of lithosphere structures and dynamics in this area has been a popular topic in geoscience, and has yielded many results. For example, TANG (1994) and LI et al (1998) suppose that the tectonic structure of Tarim Basin is various with geological periods, which changes many times betwee…  相似文献   

20.
Some consequences arising from the superposition of flows of two different kinds or scales in a non-Newtonian mantle are discussed and applied to the cases mantle convection plus postglacial rebound flow as well as small- plus large-scale mantle convection. If the two flow types have similar magnitude, the apparent rheology of both flows becomes anisotropic and the apparent viscosity for one flow depends on the geometry of the other. If one flow has a magnitude significantly larger than the other, the apparent viscosity for the weak flow is linear but develops direction-dependent variations about a factorn (n being the power exponent of the rheology). For the rebound flow lateral variations of the apparent viscosity about at least 3 are predicted and changes in the flow geometry and relaxation time are possible. On the other hand, rebound flow may weaken the apparent viscosity for convection. Secondary convection under moving plates may be influenced by the apparent anisotropic rheology. Other mechanisms leading to viscous anisotropy during shearing may increase this effect. A linear stability analysis for the onset of convection with anisotropic linear rheology shows that the critical Rayleigh number decreases and the aspect ratio of the movement cells increases for decreasing horizontal shear viscosity (normal viscosity held constant). Applied to the mantle, this model weakens the preference of convection rolls along the direction of plate motion. Under slowly moving plates, rolls perpendicular to the plate motion seem to have a slight preference. These results could be useful for resolving the question of Newtonian versus non-Newtonian or isotropic versus anisotropic mantle rheology.  相似文献   

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

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