首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Heat and mass transfer processes in the conduit of a thermochemical plume located beneath an oceanic plate far from a mid-ocean ridge (MOR) proceed under conditions of horizontal convective flows penetrating the plume conduit. In the region of a mantle flow approaching the plume conduit (in the frontal part of the conduit), the mantle material heats and melts. The melt moves through the plume conduit at the average velocity of flow v and is crystallized on the opposite side of the conduit (in the frontal part of the conduit). The heat and the chemical dope transferred by the conduit to the mantle flow are carried away by crystallized mantle material at the velocity v.The local coefficients of heat transfer at the boundary of the plume conduit are theoretically determined and the balance of heat fluxes through the side of the plume conduit per linear meter of the conduit height. The total heat generation rate, transmitted by the Hawaiian plume into the upper and lower mantle, is evaluated. With the use of regular patterns of heat transfer in the lower mantle, which is modeled on the horizontal layer, heated from below and cooled from above, the diameter of the plume source, the kinematic viscosity of the melt in the plume conduit, and the velocity of horizontal lower-mantle flows are evaluated and the dependences of the temperature drop, viscosity and Rayleigh number for the lower mantle on the diameter of the plume source are presented.  相似文献   

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

3.
Summary The equations of conservations of momentum and energy scaled with the characteristic values of the mantle indicate the presence of the upper boundary layer to produce the estimated rate of the ocean floor spreading by convection and the importance of the frictional heating. The depth of the upper boundary layer can be estimated from the balance of the viscous force with the horizontal pressure gradient at the sea floor. It is of the orders of 100 km and becomes deeper for the Pacific than for the Atlantic Ocean and also with frictional heating than without it. The frictional heating increases the surface heat flow of the heat conduction by ten to twenty percent for the Pacific Ocean but only by a few percent for the Atlantic Ocean. The similarity solutions are determined for the temperature and horizontal velocity in the upper boundary layer. These solutions are expressed in power series of the variabley x –n , wherex, y, andn are horizontal and vertical coordinates and numerical constant, respectively. Both temperature and horizontal velocity within the boundary layer are higher for the Pacific than for the Atlantic Ocean. When a larger viscosity is applied, it causes the increase of horizontal velocity below the surface because of the surface boundary conditions of the finite velocity and of vanishment of the velocity shear. The higher horizontal velocity generates higher temperature because it advects hotter material from the mid-ocean ridge site. The direct effect of frictional heating on the temperature distribution of the similarity solution is almost negligible, since the shear zone is deep and near the lower boundary of the upper boundary layer. In the similarity solution, the surface heat flow which is increased by the frictional heating is given as the boundary value. The effect of the frictional heating is important below the mid-ocean ridge.  相似文献   

4.
The effect of phenocrysts on convection in magma chambers is investigated experimentally using small heavy particles in convecting fluids. The particles are initially uniformly distributed in a fluid which is either heated from below or cooled from above. The system is allowed to evolve, and temperature and particle concentration profiles are measured as functions of time. When the concentration of particles is sufficiently small, convection is basically unaffected by their presence. When the concentration is above a critical value, however, the convective motion is considerably altered. The effect of particles on the subsequent fluid behaviour is different in the cases of heating from below and cooling from above. In the former case, there are strong convective motions confined to a sedimentary layer of decreasing thickness beneath a clear layer which displays rather weak convective motions. With time, the destabilizing increase of temperature in the lower layer overcomes the stabilizing contribution to the bulk density due to the particles and the layer overturns quite suddenly. In the situation of cooling from above, a critical condition separates a case of continual overturn from a case of no overturn at all, with the sedimentary layer falling unimpeded to the bottom. Theoretical analysis suggests that the critical value is determined primarily by the ratio of the contribution to the bulk density of the suspension due to particles to the change in fluid density due to the thermal effect. The size distribution of the particles can also modify the fluid behaviour. Applying our general results to geological situations, we suggest that the presence of relatively small concentrations of phenocrysts can critically influence the mode of convection in magmas.  相似文献   

5.
The discrete scale spectrum of the convective flows observed on the Sun has not yet received a convincing explanation. Here, an attempt is made to find conditions for the coexistence of convective flows on various scales in a horizontal fluid layer heated from below, where the thermal diffusivity varies with temperature in such a way that the static temperature difference across a thin sublayer near the upper surface of the layer is many times larger than the temperature variation across the remainder of the layer. The equations of two-dimensional thermal convection are solved numerically in an extended Boussinesq approximation, which admits thermal-diffusivity variations. The no-slip conditions are assumed at the lower boundary of the layer; either no-slip or free-slip conditions, at the upper boundary. In the former case, stable large-scale rolls develop, which experience small deformations under the action of small structures concentrated near the horizontal boundaries. In the latter case, the flow structure is highly variable, different flow scales dominate at different heights, the number of large rolls is not constant, and a sort of intermittency occurs: the enhancement of the small-scale flow component is frequently accompanied by the weakening of the large-scale one, and vice versa. The scale-splitting effects revealed here should manifest themselves in one way or another in the structure of solar convection.  相似文献   

6.
Whether in the mantle or in magma chambers, convective flows are characterized by large variations of viscosity. We study the influence of the viscosity structure on the development of convective instabilities in a viscous fluid which is cooled from above. The upper and lower boundaries of the fluid are stress-free. A viscosity dependence with depth of the form ν0 + ν1 exp(?γ.z) is assumed. After the temperature of the top boundary is lowered, velocity and temperature perturbations are followed numerically until convective breakdown occurs. Viscosity contrasts of up to 107 and Rayleigh numbers of up to 108 are studied.For intermediate viscosity contrasts (around 103), convective breakdown is characterized by the almost simultaneous appearance of two modes of instability. One involves the whole fluid layer, has a large horizontal wavelength (several times the layer depth) and exhibits plate-like behaviour. The other mode has a much smaller wavelength and develops below a rigid lid. The “whole layer” mode dominates for small viscosity contrasts but is suppressed by viscous dissipation at large viscosity contrasts.For the “rigid lid” mode, we emphasize that it is the form of the viscosity variation which determines the instability. For steep viscosity profiles, convective flow does not penetrate deeply in the viscous region and only weak convection develops. We propose a simple method to define the rigid lid thickness. We are thus able to compute the true depth extent and the effective driving temperature difference of convective flow. Because viscosity contrasts in the convecting region do not exceed 100, simple scaling arguments are sufficient to describe the instability. The critical wavelength is proportional to the thickness of the thermal boundary layer below the rigid lid. Convection occurs when a Rayleigh number defined locally exceeds a critical value of 160–200. Finally, we show that a local Rayleigh number can be computed at any depth in the fluid and that convection develops below depth zr (the rigid lid thickness) such that this number is maximum.The simple similarity laws are applied to the upper mantle beneath oceans and yield estimates of 5 × 1015?5 × 1016 m2 s?1 for viscosity in the thermal boundary layer below the plate.  相似文献   

7.
A detailed comparison between fully dynamic and kinematic plate formulations has been made in models of mantle convection. Plate velocity is computed self-consistently from fully dynamic plate models with temperature- and stress-dependent viscosity and preexisting mobile faults. In fully dynamic models, the flow is driven solely by internal buoyancy, while in kinematic models the flow is driven by a combination of the prescribed surface velocity and internal buoyancy. Only a temperature-dependent viscosity, close to the effective viscosity determined from the fully dynamic models, is used in the kinematic models. The two types of models give very similar temperature structures and slab evolutionary histories when the effective viscosity and surface velocity are nearly identical. In kinematic plate models, the additional work introduced by the prescribed velocity boundary condition is apparently dissipated within the lithosphere and has little influence on the convection under the lithosphere. In models with periodic lateral boundary conditions, slabs sink into the lower mantle at an oblique angle and this contrasts with the vertical sinking which occurs with reflecting boundary conditions. Models show that we can simulate fully dynamic models with kinematic models under either periodic boundary conditions or reflecting boundary conditions.  相似文献   

8.
Heat and mass transfer processes in the conduit of a thermochemical plume located beneath an oceanic plate far from a mid-ocean ridge (MOR) proceed under conditions of horizontal convective flows penetrating the plume conduit. In the region of a mantle flow approaching the plume conduit (in the frontal part of the conduit), the mantle material heats and melts. The melt moves through the plume conduit at the average velocity of flow v and is crystallized on the opposite side of the conduit. The heat and the chemical dope transferred by the conduit to the mantle flow are carried away by crystallized mantle material at the velocity v. The main equations of heat and mass transfer are obtained for a thermochemical plume interacting with a horizontal convective mantle flow. The joint multiparameter problem of heat and mass transfer is solved for a thermochemical plume located far from an MOR axis. The dope concentration at the base of the plume is found as a function of the Lewis number. The Lewis numbers and, accordingly, the diffusion coefficients of the chemical dope in the plume conduit far from the MOR axis are determined.  相似文献   

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

10.
The lithosphere is interpreted as a thermal boundary layer. Approximate solutions of the boundary layer cooling problem are developed which include mantle radioactivity, partial melt in the asthenosphere, a temperature gradient in the asthenosphere, and a non-zero lithospheric thickness at the ridge crests. The cooling history of oceanic lithosphere is found to be remarkably insensitive to assumptions about the amount of radioactivity in the upper mantle and the extent of melting in the asthenosphere. Determinations of the thickness of oceanic lithosphere and the depths of oceans as a function of age are in excellent agreement with boundary layer predictions which include a heat flux from the asthenosphere. However, the determinations do not resolve how much of the total asthenospheric heat flux might be caused by a temperature gradient in the asthenosphere. Simple thermal arguments indicate that the initial lithospheric thickness, L0, at ridge crests should depend on the local half-spreading rate, V, as L0 = 3 km/V(cm/year).  相似文献   

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

12.
An instantaneous analytical approach is developed to predict sheet flow transport in purely acceleration-skewed oscillatory flow. The approach is derived from exponential approximations of velocity and concentration profiles above a mobile seabed, and it particularly considers factors of phase lead; phase lag (i.e. phase residual and phase shift); acceleration modification; and asymmetries in shear stress, roughness height, and boundary layer development. The approach can predict net boundary layer flow above a mobile seabed, and can revert to the classical bedload model. Instantaneous and net sediment transport rates are studied using the approach. The instantaneous sediment transport rate in an onshore flow stage can be approximated by a power function of velocity in which the exponent is confirmed to range between 1 and 5 with a decrease in the phase residual. The net sediment transport rate predicted using the approach is validated using a considerable amount of measured data, and compared with existing instantaneous and half-period type models that consider the phase lag or acceleration modification. For the net sediment transport rate in purely acceleration-skewed oscillatory flow, the phase residual is less important than the acceleration-skewed boundary layer difference between onshore and offshore acceleration stages.  相似文献   

13.
Sediment movement in the wave boundary layer above a mobile sediment bed is complex.A velocity formula for the boundary layer is proposed for sheet flow induced by asymmetric waves above a mobile sediment bed.The formula consists of a free stream velocity and a defect function which contains a phase-lead,boundary layer thickness and mobile sediment bed.Phase-lag of sediment movement is considered in the formula for the mobile sediment bed.The formula needs six dependent variables about asymmetric wave and sediment characteristics.Asymmetry effects on parameters(orbital amplitude,roughness height,bed shear stress,and boundary layer thickness)are properly considered such that the formula can yield velocity differences among onshore,offshore,acceleration,and deceleration stages.The formula estimates the net boundary layer velocity resulting from the mobile sediment bed and asymmetric boundary layer thickness.In addition,a non-constant phase-lead also contributes to the net boundary layer velocity in asymmetric oscillatory sheet flow.Results of the formula are as good as that of a two-phase numerical model.Sheet flow transport induced by asymmetric waves,and the offshore net sediment transport rate with a large phase-lag under velocity-skewed waves,can be adequately estimated by the formula with a power sediment concentration function.  相似文献   

14.
In this study, we examine the development of topography on a thin dense layer at the base of the lower mantle. The effect of the convecting mantle above is represented as a traction acting on the upper surface of the layer. Topography on the layer boundaries is predicted by a balance of dynamic flow stress and external traction. The nature of boundary topography depends on the magnitude of the driving tractions and the density variation within the layer. If we assume that the layer density is greatest beneath areas of mantle downwelling and decreases to a minimum beneath areas of mantle upwelling (the layer is thermally coupled to the convection in the overlying mantle) then its upper boundary develops a cusp-like peak beneath the upwelling mantle. The height of this peak is potentially much greater than the layer thickness. If, however, the layers are effectively coupled by viscous shear then internal density gradients of the opposite sign may be established. In this case, we observe solutions where the layer is completely swept away beneath areas of mantle downwelling leaving steep-sided ‘islands’ of dense material. This mechanism therefore provides a possible explanation for steep-sided anomalously slow regions at the base of the mantle observed by seismic methods (e.g. beneath south Africa) or for discrete ultralow velocity zones detected at the core-mantle boundary beneath locations of surface hotspots. The magnitude of the upper boundary driving tractions compared to the density gradient within the layer is the key parameter that determines the nature of flow in, and consequently boundary topography of, the layer. The deflection of the core-mantle boundary is small compared with that of the top of the dense layer, but a change in sign of the ratio of these deflections is observed as the magnitude of the driving tractions changes relative to the magnitude of the internal density gradient. We compare seismic measurements of core-mantle boundary topography and D′′ topography with the predictions of this model in an attempt to constrain model parameters, but no clear correlation seems to exist between D′′ thickness and CMB topography.  相似文献   

15.
In order to study the relationship between mantle flow and global tectogenesis, we present a 3-D spherical shell model with incompressible Newtonian fluid medium to simulate mantle flow which fits the global tectogenesis quite well. The governing equations are derived in spherical coordinates. Both the thermal buoyancy force and the self-gravitation are taken into account. The velocity and pressure coupled with temperature are computed, using the finite-element method with a punitive factor. The results show that the lithosphere, as the boundary layer of the earth's thermodynamic system, moves with the entire mantle. Both its horizontal and vertical movements are the results of the earth's thermal motion. The orogenesis occurs not only in the collision zones at the plates' boundaries, but also occurs within the plates. If the core-mantle boundary is impermeable and the viscosity of the lower mantle is considerable, the vertical movement is mostly confined to the upper mantle. The directions of the asthenospheric movements are not fully consistent with those of the lithospheric movements. The depths of spreading movements beneath all ridges are less than 220 km. In some regions, the shear stresses, acting on the base of the lithosphere by the asthenosphere, are the main driving force; but in other regions, the shear stresses are the resisting force.  相似文献   

16.
The wave-induced velocity and pressure fields beneath a large amplitude internal solitary wave of depression propagating over a smooth, flat, horizontal, and rigid boundary in a shallow two-layer fluid are computed numerically. A numerical ocean model is utilised, the set-up of which is designed and tuned to replicate the previously published experimental results of Carr and Davies (Phys Fluids 18(1):016,601–1–016,601–10, 2006). Excellent agreement is found between the two data sets and, in particular, the numerical simulation replicates the finding of a reverse flow along the bed aft of the wave. The numerically computed velocity and pressure gradients confirm that the occurrence of the reverse flow is a consequence of boundary layer separation in the adverse pressure gradient region. In addition, vortices associated with the reverse flow are seen to form near the bed.  相似文献   

17.
Results are presented from both linear stability analysis and numerical simulations of three-dimensional nonlinear convection in a Boussinesq fluid in an annular channel, under experimental boundary conditions, rotating about a vertical axis uniformly heated from below. The focus is placed on the Prandtl number Pr = 7.0, representing liquid water at room temperature. The linear analysis shows that, when the aspect ratio is sufficiently small, there exists only one stationary mode that occupies the whole fluid container. When the aspect ratio is moderate or large, however, there exist three different linear solutions: (i) the outer sidewall-localized traveling wave propagating against the sense of rotation; (ii) the inner sidewall-localized traveling wave propagating in the same sense as rotation; and (iii) both the counter-traveling waves occurring simultaneously. Guided by the result of the linear stability analysis, fully three-dimensional simulations are then performed for a channel with a moderate aspect ratio. It is found that neither the prograde nor the retrograde mode is physically realizable near threshold and beyond. The dynamics of nonlinear convection in a rotating channel are chiefly characterized by the interaction between the sidewall-localized waves and the interior convection cells/rolls, producing an interesting and unusual nonlinear phenomenon. In order to compare with the classical Rayleigh–Bénard problem without vertical sidewalls, we also study linear and nonlinear convection at exactly the same parameters but in an infinitely extended layer with periodic horizontal conditions. This reveals that both the linear instability and nonlinear convection in a rotating channel are characteristically different from those in a rotating layer with periodic horizontal conditions.  相似文献   

18.
裂谷盆地构造-热演化模拟中几个问题的讨论   总被引：2，自引：2，他引：0       下载免费PDF全文
裂谷盆地的构造-热演化模拟是在岩石圈尺度计算裂谷盆地形成演化过程中的热历史和沉降史.拉张模型实现了构造和热的完美结合,在描述裂谷盆地沉降和热流演化方面取得了很大的成功.本文使用二维运动学模型,通过有限元方法,在拉格朗日坐标系下进行拉张背景下的构造热演化模拟,探讨了拉张模型中初始地壳、岩石圈厚度、软流圈对流、模型上边界对构造热演化的影响,以及载水和载沉积物两种情况下盆地侧翼抬升的差异.  相似文献   

19.
A one-dimensional model of flow between a fixed boundary at the bottom and a moving one on top with no net flow through vertical sections is tested for geophysically interesting mantle viscosity-depth functions. Such a model, although simplistic, may help in answering the question to what depth the return flow extends, at least in the case of moving plates measuring many thousand kilometers across, such as the Pacific plate.It the viscosity in the asthenosphere is less than three orders of magnitude smaller than that of the mantle below, the return flow extends to great depth and the asthenosphere is a zone of concentrated shear. If the viscosity contrast is greater, the return flow is concentrated in the asthenosphere. For a wide range of model parameters typical flow velocities below the asthenosphere are about one-tenth of the plate velocity. The pressure gradient required by the mantle flow may be manifest in gravity trends across moving plates, but no excessive gravity anomalies are required by the model if the absolute viscosity values conform to those inferred from post-glacial rebound data. A thinner and lower-viscosity layer is favored over a thicker and more viscous layer if both fit glacial rebound evidence. The present model may not be applicable if down to the core the viscosity is as low as about 1021 N s m–2 with a free-slip bottom boundary.  相似文献   

20.
Introduction Jiashi-Artux area in southwest Xinjiang is one of the most active earthquake provinces at pre-sent in Chinese mainland. In the last century, about 3/4 strong earthquakes in Chinese mainland hit this area, and especially from January 21 to April 16 in 1997, 7 earthquakes with the magnituderanging from 6.0 to 6.9 occurred in a very small area of 9 km×18 km near Jiashi (ZHU et al, 1998). It has never taken place before in Chinese mainland that a series of strong earthquakes shoo…  相似文献   

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

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