Constraints on the dynamics of mantle plumes from uplift of the Hawaiian Islands

The ∼0.2 mm/yr uplift of Hawaiian islands Lanai and Molokai and Hawaiian swell topography pose important constraints on the structure and dynamics of mantle plumes. We have formulated 3-D models of mantle convection to investigate the effects of plume-plate interactions on surface vertical motions and swell topography. In our models, the controlling parameters are plume radius, excess plume temperature, and upper mantle viscosity. We have found that swell height and swell width constraints limit the radius of the Hawaiian plume to be smaller than 70 km. The additional constraint from the uplift at Lanai requires excess plume temperature to be greater than 400 K. If excess plume temperature is 400 K, models with plume radius between 50 and 70 km and upper mantle viscosity between 10²⁰ and 3×10²⁰ Pa s satisfy all the constraints. Our results indicate that mantle plume in the upper mantle may be significantly hotter than previously suggested. This has important implications for mantle convection and mantle melting. In addition to constraining plume dynamics, our models also provide a mechanism to produce the observed uplift at Lanai and Molokai that has never been satisfactorily explained before.     

Reconciling strong slab pull and weak plate bending: The plate motion constraint on the strength of mantle slabs

Although subducting slabs undergo a bending deformation that resists tectonic plate motions, the magnitude of this resistance is not known because of poor constraints on slab strength. However, because slab bending slows the relatively rapid motions of oceanic plates, observed plate motions constrain the importance of bending. We estimated the slab pull force and the bending resistance globally for 207 subduction zone transects using new measurements of the bending curvature determined from slab seismicity. Predicting plate motions using a global mantle flow model, we constrain the viscosity of the bending slab to be at most ~ 300 times more viscous than the upper mantle; stronger slabs are intolerably slowed by the bending deformation. Weaker slabs, however, cannot transmit a pull force sufficient to explain rapid trenchward plate motions unless slabs stretch faster than seismically observed rates of ~ 10^− 15 s^− 1. The constrained bending viscosity (~ 2 × 10²³ Pa s) is larger than previous estimates that yielded similar or larger bending resistance (here ~ 25% of forces). This apparent discrepancy occurs because slabs bend more gently than previously thought, with an average radius of curvature of 390 km that permits subduction of strong slabs. This gentle bending may ultimately permit plate tectonics on Earth.     

Study on the rheology of subducting slabs

We calculate thermal and phase structures of subducting slabs for different subducting velocities by a modified coupling code of the kinetic phase-transformation equations and the heat-diffusion equation with latent-heat release. Whereafter, we estimate their rheology structures based on the thermal and phase structures from the mineral physical point of view. At shallow depth, the upper layer has a high effective viscosity greater than 10³⁴Pa · s; while the lower layer has a relatively low effective viscosity, which is greater than 10²⁶Pa · s nevertheless. The effective viscosities below the kinetic phase boundary of olivine to wadsleyite decrease obviously, and reach a minimum of 10²²Pa · s. Small areas with higher effective viscosities exist above the depth of about 700 km in subducting slabs, which are produced by lower temperatures that are related with endothermic phase transformation of spinel to perovskite and magnesiowustite. The 1% and 99% isograds of spinel proportion delineate tortuous belts with low effective viscosities, which would affect the geodynamic behavior of subducting slabs.     

Lateral variation in upper mantle viscosity: role of water

Differences in the viscosity of the earth's upper mantle beneath the western US (∼10¹⁸-10¹⁹ Pa s) and global average values based on glacial isostatic adjustment and other data (∼10²⁰-10²¹ Pa s) are generally ascribed to differences in temperature. We compile geochemical data on the water contents of western US lavas and mantle xenoliths, compare these data to water solubility in olivine, and calculate the corresponding effective viscosity of olivine, the major constituent of the upper mantle, using a power law creep rheological model. These data and calculations suggest that the low viscosities of the western US upper mantle reflect the combined effect of high water concentration and elevated temperature. The high water content of the western US upper mantle may reflect the long history of Farallon plate subduction, including flat slab subduction, which effectively advected water as far inland as the Colorado Plateau, hydrating and weakening the upper mantle.     

The mantle convection model with non-Newtonian rheology and phase transitions: The flow structure and stress fields

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.     

Possible effects of lateral viscosity variations induced by plate-tectonic mechanism on geoid inferred from numerical models of mantle convection

In a traditional analytical method, the convective features of Earth’s mantle have been inferred from surface signatures obtained by the geodynamic model only with depth-dependent viscosity structure. The moving and subducting plates, however, bring lateral viscosity variations in the mantle. To clarify the effects of lateral viscosity variations caused by the plate-tectonic mechanism, I have first studied systematically instantaneous dynamic flow calculations using new density-viscosity models only with vertical viscosity variations in a three-dimensional spherical shell. I find that the geoid high arises over subduction zones only when the vertical viscosity contrast between the upper mantle and the lower mantle is O(10³) to O(10⁴), which seems to be much larger than the viscosity contrast suggested by other studies. I next show that this discrepancy may be removed when I consider the lateral viscosity variation caused by the plate-tectonic mechanism using two-dimensional numerical models of mantle convection with self-consistently moving and subducting plates, and suggest that the observed geoid anomaly on the Earth’s surface is significantly affected by plate-tectonic mechanism as a first-order effect.     

两种地幔对流模式下俯冲带的热结构

根据准动力学计算方案，通过采用等效热源和等效热传导系数的方法，用有限元法计算了不同俯冲角度，而俯冲速度为８ｃｍ／ａ、年龄为１００Ｍａ的俯冲带在稳定俯冲状态的热结构．计算结果表明俯冲带在接近６７０ｋｍ间断面的最低温度可达到１１００℃．全地幔对流模式热结构的计算结果表明６７０ｋｍ间断面以下可存在最低温度达１０００℃的低温区，相应于有０．７％-３．０％的Ｐ波低速异常存在．双层地幔对流模式表明，在６７０ｋｍ间断面以上可有与周围地幔相差约４００℃的水平舌状低温区存在，相应于０．５％-１．４％的Ｐ波低速异常．     

Stability of the oceanic lithosphere with variable viscosity: an initial-value approach

We have studied the problem concerning the onset of convective instabilities below the oceanic lithosphere. A system of linear partial differential equations, in which the background temperature field is time-dependent, is integrated in time to monitor the evolution of incipient disturbances. Two types of rheologies have been examined. One depends strongly on temperature. The other involves a viscosity which is both temperature- and pressure-dependent. The results from this initial-value approach, in which the viscosity profiles migrate downward with time, reveal the importance of considering temperature- and pressure-dependent rheology in issues regarding the development of local instabilities in upper mantle convection. For temperature-dependent viscosity, viscosities of 0(10²⁰P) are required to produce instabilities with growth-rates of 0(.1/Ma). In contrast, these same growth rates can be attained for a temperature- and pressure-dependent viscosity profile with a mean value close to 0(10²⁰P) in the upper mantle, owing to the presence of a low viscosity zone, 0(10²⁰P), existing right below the lithosphere. Unlike the results of temperature-dependent viscosity, whose growth-rates increase with time, the amplification of disturbances in a fluid medium with temperature- and pressure-dependent rheology reaches a maximum at an early age, < 50 Ma, and decreases thereafter with time. This suggests the potential importance played by initial disturbances in the evolution of the oceanic lithosphere.     

How flat is the lower-mantle temperature gradient?

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⁻¹, uniform internal heating and a viscosity increase of 30 across the phase transition. The global Rayleigh number for basal heating is around 2×10⁶, while an internal heating Rayleigh number as high as 10⁸ 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⁻¹, 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.     

The existence of a thin low-viscosity layer beneath the lithosphere

The horizontal temperature gradient at the base of the lithosphere at an oceanic fracture zone, where plate of different ages is juxtaposed, is expected to drive a local circulation, the characteristics of which can be constrained by the amplitude, wavelength and age-dependence of the geoid. Two-dimensional numerical models of convection in a fluid layer overlain by a solid conducting lid have been used to generate theoretical geoid profiles at right angles to the fracture zone. Only a thin, low-viscosity layer provides a reasonable fit to the data. The best model so far obtained has a fluid layer 150 km thick with viscosity 1.5 × 10¹⁹ Pa s under a 75 km lid. Such a layer, which is incapable of transmitting strong horizontal shear stresses, could provide the decoupling mechanism between plate and deep mantle flow required to balance the forces on the plates.     

Heat transport by variable viscosity convection and implications for the Earth's thermal evolution

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 10⁴?10⁷ 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.     

Rheology of the lower mantle

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 10² to 10³ 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.     

On the thermal evolution of the earth

The two principal contributions to the surface heat flow of the earth are the cooling of the earth and the heat production of radioactive isotopes. As the rate of heat production decreases with time the temperature of the interior of the earth also decreases. The rate of decrease is determined by the ability of solid-state mantle convection to transport the heat to the surface. The dominant effect is the exponential temperature dependence of the mantle viscosity. The non-dimensional mantle temperature can be parameterised in terms of the Rayleigh number for mantle convection. It is found that the mantle is currently cooling at a rate of 36°K/10⁹ years and that three billion years ago the mean temperature was 150°K higher than it is today; 83% of the present surface heat flow is attributed to the decay of radioactive isotopes and 17% to the cooling of the earth. The corresponding mean concentration of uranium in the mantle is 32 ppb.     

地幔对流的数值模拟及其与表面观测的关系

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

Mantle rheology: radial and lateral viscosity variations inferred from microphysical creep laws

Inferences on the rheology of the mantle based on theoretical and experimental rate equations for steady state creep are discussed and compared with results from geophysical models. The radial increase of viscosity by one to three orders of magnitude across the mantle, required by inversion of postglacial rebound and geodynamic data, is confirmed by microphysical models based on the estimation of continuous and discontinuous changes of creep parameters with depth. The upper mantle (viscosity 10²⁰–10²¹ Pa s) is likely to show non-Newtonian rheology (power-law creep) for average grain sizes larger than 0.1 mm as an order of magnitude. Given the variability of both grain size and stress conditions, local regions of linear rheology can be present. The rheology of transition zone and lower mantle (viscosity 10²²–10²⁴ Pa s) cannot be definitely resolved at present. Estimation of creep parameters leads to possible nonlinear or mixed rheology, if grain sizes are not lower than 0.1 mm and flow conditions can be approximated by a constant strain rate of about 10⁻¹⁵ s⁻¹. This conclusion can be modified by different flow conditions (e.g. a decrease in strain rate or constant viscous dissipation). Furthermore, experiments on fine-grained garnetites and perovskite analogues have shown that diffusion creep is predominant at laboratory conditions. However, the pressure dependence of creep in these phases is unknown, and therefore direct extrapolation to lower mantle conditions is necessarily speculative. Lateral variations of viscosity, largest in the upper and lowermost mantle (up to 2–4 orders of magnitude) are predicted by models based on lateral temperature anomalies derived from seismic tomographic models.     

Convective instabilities in a variable viscosity fluid cooled from above

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 ν₀ + ν₁ 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 10⁷ and Rayleigh numbers of up to 10⁸ are studied.For intermediate viscosity contrasts (around 10³), 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 z_r (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 × 10¹⁵?5 × 10¹⁶ m² s^?1 for viscosity in the thermal boundary layer below the plate.     

How valid are dynamic models of subduction and convection when plate motions are prescribed?

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.     

Tomographic structure of East Asia:Ⅱ.Stagnant slab above 660 km discontinuity and its geodynamic implications

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.     

Rheological properties of deep subducted oceanic lithosphere and their geodynamic implications

According to the experimental studies on the rheology of two important mantle rocks (eclogite and harzburgite), the rheological properties of the deep subducted oceanic lithosphere are investigated by assuming a simplified harzburgite type slab model with moderate thickness of basaltic layer. When the mantle convergence rate is small or the subducting slab has been trapped in the mantle for an enough long time, the strength profile of the slab is characterized by a strong subducting crustal component lying on a weak subducting upper mantle. However, if the convergence rate is large enough, the subducting slab will be featured only by a rigid cold center. Our study suggests that the detachment of the subducting crust component from the underlying upper mantle is only likely to happen in hot slow subducting slabs, but not the cold fast subducting lithosphere. Rheological properties of the harzburgitic and the eclogitic upper mantle vary with depths. The eclogitic upper mantle is stronger than the peridotitic upper mantle across the upper mantle. Transition zone is the high strength and high viscosity layer in the upper mantle except the lithosphere.