Outline
Comptes Rendus

Geodynamics / Géodynamique
Simple features of mantle-wide convection and the interpretation of lower-mantle tomograms
Comptes Rendus. Géoscience, Volume 335 (2003) no. 1, pp. 9-22.

Abstracts

Simple fluid dynamic constraints aid in the interpretation of lower mantle tomograms. The geothermal gradient away from slabs and plumes is subadiabatic between the upper (lithosphere) and lower (D″) boundary layers by ∼400 K. Slabs widen proportionally with the square root of the viscosity as they sink through the lower mantle. The time scale for the persistence of slab graveyards at the base of the mantle is comparable to the time for plate motions to reorganize ∼120–200 Myr. At most a few starting plume heads currently exist in the lower mantle. Tomographic inversions may include more sophisticated numerical-modeling versions of these constraints.

Des contraintes simples de dynamique des fluides aident à interpréter les tomogrammes du manteau inférieur. Le gradient de température du manteau, loin des plaques plongeantes et loin des panaches, est sous-adiabatique (d'environ 400 K) entre les couches limites supérieure (la lithosphère) et inférieure (la couche D″). Les plaques plongeantes s'épaississent en fonction de la racine carrée de la viscosité en s'enfonçant dans le manteau. Les cimetières de plaques dans le manteau inférieur ont la même échelle de persistance dans le temps que l'organisation de la dérive des plaques, ∼120–200 Ma. Actuellement, il y a au plus quelques nouvelles têtes de panache dans le manteau inférieur. Les inversions tomographiques peuvent inclure ces contraintes, mais avec des méthodes numériques plus sophistiquées.

Metadata
Received:
Accepted:
Published online:
DOI: 10.1016/S1631-0713(03)00008-7
Keywords: Mantle convection, Slabs, Subduction, Tomography, Geothermal gradient, Plumes, Plume heads
Mot clés : Convection mantélique, Lithosphère plongeante, Subduction, Tomographie, Courbe géothermique, Panaches, Têtes de nouveaux panaches

Norman H. Sleep 1

1 Department of Geophysics, Stanford University, Stanford, CA 94305, USA
@article{CRGEOS_2003__335_1_9_0,
     author = {Norman H. Sleep},
     title = {Simple features of mantle-wide convection and the interpretation of lower-mantle tomograms},
     journal = {Comptes Rendus. G\'eoscience},
     pages = {9--22},
     publisher = {Elsevier},
     volume = {335},
     number = {1},
     year = {2003},
     doi = {10.1016/S1631-0713(03)00008-7},
     language = {en},
}
TY  - JOUR
AU  - Norman H. Sleep
TI  - Simple features of mantle-wide convection and the interpretation of lower-mantle tomograms
JO  - Comptes Rendus. Géoscience
PY  - 2003
SP  - 9
EP  - 22
VL  - 335
IS  - 1
PB  - Elsevier
DO  - 10.1016/S1631-0713(03)00008-7
LA  - en
ID  - CRGEOS_2003__335_1_9_0
ER  - 
%0 Journal Article
%A Norman H. Sleep
%T Simple features of mantle-wide convection and the interpretation of lower-mantle tomograms
%J Comptes Rendus. Géoscience
%D 2003
%P 9-22
%V 335
%N 1
%I Elsevier
%R 10.1016/S1631-0713(03)00008-7
%G en
%F CRGEOS_2003__335_1_9_0
Norman H. Sleep. Simple features of mantle-wide convection and the interpretation of lower-mantle tomograms. Comptes Rendus. Géoscience, Volume 335 (2003) no. 1, pp. 9-22. doi : 10.1016/S1631-0713(03)00008-7. https://comptes-rendus.academie-sciences.fr/geoscience/articles/10.1016/S1631-0713(03)00008-7/

Version originale du texte intégral

1 Introduction

I discuss the thermal structure of the mantle as a guide to interpreting tomograms. That is, what features in the Earth are likely to be detected tomographically. I presume the existence of mantle-wide convection with slabs and plumes. Crudely, the lower mantle contains of descending slabs and graveyards of slabs ponded near its base (Fig. 1). Romanowicz and Gung [30] call the hotter regions between slabs and slabs graveyards superplumes. I retain this term as ‘superplume regions’ for lack of a better one. Hot spots preferentially occur above superplume regions, rather than above slab graveyards. The mantle at 250–500-km depth is also hotter than its surroundings above superplumes. This implies, as expected, that mantle plumes arise from places at the base of the mantle where the basal thermal boundary layer, D″, is thick.

Fig. 1

Schematic cross-section of the Earth, drawn to scale. Slabs sink and pond at the base of the mantle. The subducted material moved upward toward the ridge over 2000 Myr. Plume heads start from cusps in the basal boundary layer.

Schéma en coupe de la Terre, dessiné à l'échelle. Les plaques plongent dans le manteau inférieur et s'accumulent à sa base, où elles stagnent. Le matériel subducté remonte vers la dorsale en 2000 Ma. Les têtes de panache commencent au niveau de rebroussements de la couche limite inférieure.

The question arises as to whether deep superplume regions are hotter than the upper mantle (mid-oceanic ridge basalt, MORB) adiabat, rather than merely the surrounding slab graveyards at great depth. To address this issue, I begin with the average geothermal gradient in the mantle away from slabs and plumes in Section 2. I then discuss subduction with regard to the width of slabs in the lower mantle and the temperature within slab graveyards deep in the lower mantle in Section 3. I then discuss the persistence of lateral temperature anomalies associated with superplume regions in the deep mantle in Section 4. I finally discuss the ascent of starting plume heads as a feature that might be detected by tomography in Section 5 and return to the related problem of the persistence of superplume regions.

I compile and adapt dimensional and analytical methods in this regard. My intent is to make easily understood general conclusions that do not obscure gross observable features of convection with complex mathematics and numerical models. This work serves as a guide for interpreting numerical studies of convection in regard to tomography and vice versa. Tomographers and fluid dynamicists are well aware of the need to act together in a self-consistent manner (e.g., [16,17,22–24,30]).

2 Geothermal gradient at mid-mantle depths

Global tomography images the broad-scale structure of the Earth's mantle. One finds velocities and attenuation of seismic waves and, to some extent, the density and the anisotropy. These measurements constrain the thermal structure as seismic velocity decreases and attenuation increases with temperature. Current tomograms resolve slabs, slab graveyards, and superplume regions. Our lack of knowledge of mineral physics limits what can be directly inferred from these results.

Specifically, common methods of plotting tomograms are potentially misleading. To show these difficulties, I formally represent the shear wave velocity as a function of depth and temperature,

β=β(Z ref ,T ref )+βTZ(T-T ref )+βZS(Z-Z ref )(1)
where T is potential temperature, Tref is a convenient reference potential temperature (like the MORB adiabat), Z is depth, Zref is a reference depth, and S is entropy. Similar formal expressions exist for density, P-wave velocity and the attenuation of both types of waves.

In practice, the depth derivatives are unknown, particularly for the lower mantle. Tomographers thus often explicitly center their plots on the mean value at a given depth. Tomographers are well aware that plots generated with respect to local depth averages need to be carefully interpreted if the actual average thermal gradient in the mantle is subadiabatic at mid-mantle depths as shown later in this section. For example, material hotter than the average at 2400-km depth may be cooler than the shallow MORB adiabat and not able to ascend buoyantly into upper mantle having a cooler color on the plot.

The temperature derivatives may also be unknown, especially if trace phases like partial melt affect seismic velocity and attenuation. To partly adjust for this problem and to show variations at all depths on the same plot, tomographers may explicitly weight color schemes to cover the observed range at each depth. For example, the same color range may represent a 10% variation at 50-km depth but only 1% at 2000-km depth.

It is well known that the thermal gradient of a fluid heated from within is subadiabatic in the middle of the convecting region (e.g., [29]). Whole mantle convection behaves in this way as the major heat sources are radioactive heat generation and internal cooling (e.g., [8,12,13]). The core provides only a modest amount of heating from below [12,13,34]. The function of plate tectonics is to cool the mantle and that of plumes is to cool the core [13].

To estimate the mid-mantle thermal gradient, I begin with the heat flow equation,

ρCTt+ρCV·T=k2T+ρA(2)
where ρ is density, C is heat capacity per mass, t is time, V is material velocity, k is thermal conductivity (constant for simplicity), and A is radioactive heat generation per mass. The conductivity term is important only in the thermal boundary layers, where large thermal gradients exist. These include the lithosphere and slabs formed from it and the basal boundary layer and plumes formed from it (Fig. 1). All mantle material needs to cycle through the lithosphere to lose its heat. This situation contrasts with convection heated from below where only the heated material needs to pass through the upper thermal boundary layer.

I follow a batch of material that was once part of a slab to obtain a quantitative estimate of the mid-mantle geothermal gradient (Fig. 1). The slab material sinks (negatively) buoyantly into the deep mantle. Its temperature anomaly spreads out as the material ponds near the base of the mantle. Thereafter, the conduction term in Eq. (1) is small. The velocity term in Eq. (2) is zero by definition if one sets the coordinate system to follow the material. The heat flow equation (2) then reduces to static heating by radioactive heat generation,

DT Dt =AC(3)
where D/Dt is the time derivative moving with the material. The batch of material ascends on average as new slabs sink below it and pond. That is, shallow material is hotter than deep material because shallow material sank within a slab a long time ago and has had more time to heat up than the deep material that sank recently (Fig. 2).

Fig. 2

Schematic diagram of average mantle properties away from plumes and slabs. Potential temperature has a subadiabatic gradient in the mid-mantle (left). The vertical velocity, normalized as 4πrh2ρVZ, provides the return flow for slabs and plumes (right).

Représentation schématique des propriétés moyennes du manteau à longue distance des panaches et des plaques plongeantes. La courbe de la température potentielle suit une pente sous-adiabatique dans le manteau moyen (à gauche). La vitesse verticale, normalisée à 4πrh2ρVZ, fournit le courant de retour pour les plaques plongeantes et les panaches (à droite).

The interior of the actual Earth is likely to have cooled over geologic time [1,4,18]. In this case, cooling supplies heat to the surface. A subadiabatic gradient exists as with radioactive heating. The material now arriving near the surface has a higher potential temperature than deep material that will arrive say 1000 million years from now. Equivalently, the average temperature of the material that descended at say ∼2000 Myr and is now reaching the surface was higher than the average temperature of material that descended at ∼1000 Myr and will reach the surface in 1000 million years.

Combining the long term cooling rate of the Earth's mantle ∂T/∂t and radioactive heating the departure ΔT of the geotherm across the mantle from adiabatic is

ΔT=-Tt+ACΔt(4)
where the material requires a time Δt to ascend from the deep mantle to the surface. The departure of the thermal gradient from adiabat is:
ΔTM=-Tt+ACΔtM=1VZ-Tt+AC(5)
where M is the thickness of the mantle and VZ is the vertical velocity.

One obtains the subadiabatic gradient (as ΔTt) from Eqs. (4) and (5) from the surface heat flow without knowing the fraction of heat supplied by radioactivity because the sum of these effects is involved. In particular, the magnitude of the departure of the thermal gradient from adiabatic relates to the heat balance in the Earth's mantle. The average surface heat flow coming from the mantle is:

qS=ρCΔTVZ(6)
and the global heat balance is:
4πrE2qS=4πrE2MρCΔT/Δt(7)
where rE is the radius of the Earth and 4πrE2Mρ is the mass of the mantle.

Sclater et al. [33] have estimated the heat flow from the Earth. Following Davies [13], I quickly show the essence of the method. The Earth's surface is 0.6 ocean basins and 0.4 continents. The oceanic heat flow is approximately 500 mW m−2 divided by the square root of age in million years. Currently oceanic crust recycles at a rate of 3 km2 Myr−1, which would replace the ocean basins in 100 Myr. The average heat flow in the ocean basins is thus ∼100 mW m−2. The continental heat flow (not associated with crustal radioactivity) is less well constrained to ∼25 mW m−2. The global average is thus ∼70 mW m−2. About 0.1 of this heat flow comes ultimately from detected plumes [12,34]. This leaves ∼63 mW m−2 from mantle radioactivity and mantle cooling. As the heat capacity of the mantle is ∼1.25 kJ m−3 K−1 [36] and the mass of the mantle is 4×1024 kg, this is equivalent to a heat production of 8×10−12 W kg−1 or to a cooling rate of 0.2 K Myr−1 or a combination there of.

I estimate the return time for the flow Δt from the geometry of plate tectonics. Assuming that the slab and the material entrained with it are 200-km thick, 600 km3 with a mass of 2×1015 kg is subducted each year. The mantle (4×1024 kg) cycles on a time scale Δt of 2000 Myr, implying that ΔT is ∼400 K (from Eq. (7), using qS=63 mW m−2 and C=1.25 kJ m−3 K−1). That is, the thermal gradient in the mantle is likely to be significantly subadiabatic and that tomograms with respect to mean velocity at a given depth should not be construed as diagrams with respect to the MORB adiabat.

To summarize, average geotherm in the mantle (away from slabs and plumes) consists of a superadiabatic conductive geothermal gradient in the lithosphere, a superadiabatic thermal gradient within the D″ thermal boundary layer at the base of the mantle, and a subadiabatic thermal gradient within the intervening mantle (Fig. 2). A chemically denser layer may underlie the convecting part of D″, but is not shown for simplicity. The existence of hot spots associated with mantle plumes indicates that the potential temperature within the convecting part of D″ is hotter than the MORB adiabat.

3 Slab descent

The width of slabs as they descend is detectable by tomography. I quantify the width of slabs sinking into the lower mantle using simple heat and mass transfer methods. I constrain the temperature anomaly associated with slabs in this process.

3.1 Force balance and broadening of the slab with depth

Steinberger [37] dynamically models the descent of slabs into the lower mantle. I use dimensional scaling to represent the essence of this process in the case that viscosity increases strongly with depth. I begin with the velocity that a slab sinks through the lower mantle. For simplicity, my slab sinks vertically (Fig. 3). The slab is able to deform so that its width increases as it sinks into more viscous material (the widening may occur by buckling and folding. I do not consider the mechanical details of widening, as they are not susceptible to simple dimensional analysis and not likely to be resolved in tomograms). The shear traction on the sides of the slab balances the body force producing sinking, dimensionally:

τΔρgW(8)
where W is the width of the slab, Δρ is the density contrast of the slab with its surroundings, and g is the acceleration of gravity. Conduction spreads out this anomaly, but does not affect its magnitude to the first order (I return to conduction in Section 3.2). The shear traction from viscous forces is, dimensionally,
τηVSL(9)
where η is the viscosity of the ordinary mantle surrounding the slab, the slab sinks at a velocity of VS, and L is the length scale for the flow. At steady state, the volume per length flux VSW of slab material is constant with depth. Solving Eqs. (8) and (9) with this constraint (keeping the density contrast constant, as it is likely to change much less than viscosity) yields that the width of the slab increases as
W=W ref ηη ref 1/2(10)
where Wref is the width of the slab at a reference viscosity ηref. The velocity is:
VS=V ref ηη ref -1/2=ΔρgW ref Lη1/2η ref 1/2(11)
where Vref is the velocity of the slab at a reference viscosity ηref. The leading term ΔρgW ref is the (negative) buoyancy per surface area of the slab. This quantity scales with the square root of (subducted) plate age if the reference depth is taken near the top of the lower mantle were the slab is still behaving somewhat rigidly.

Fig. 3

Schematic diagram of subducting slab. The shear traction balances the negative buoyancy of the slab. The slab widens so that both material and forces balance.

Schéma de plaque plongeante. Un équilibre de forces existe entre la contrainte de cisaillement et la poussée d'Archimède (négative) de la plaque plongeante. La plaque plongeante s'élargit pour assurer les équilibres des forces et du bilan de matière.

To obtain closed-formed results, I assume that the viscosity is given by the simple function:

η=η ref exp δTTη exp ZDη(12)
where δT is the potential temperature below a reference adiabat, Tη is the temperature scale for the viscosity, Z is depth below some reference, and Dη is the depth scale for the viscosity. The actual variation of viscosity with depth (as would be determined by a glacial rebound study) includes the effect of a subadiabatic gradient:
δT=ΔTZM(13)
so that the actual variation of viscosity with depth is:
η=η ref exp Δ TZ MT η exp ZDηη ref exp ZD app (14)
where Dapp is the apparent depth scale for viscosity in the mantle. The time for a slab to sink to some depth Z below the reference depth is:
t sink =ZZ0dzV2D app V ref exp Z2D app 2D app η ref ΔρgW ref L exp Z2D app (15)
where Z is a dummy variable for depth and Z0 is the starting depth. The two approximate equalities that apply to the time for a slab to sink into the deep mantle take the lower limit of the integral as −∞.

Eqs. (10) and (15) suffice to illustrate some basic points about the expected properties of slabs within the lower mantle. The increase of viscosity across the lower mantle is unknown but may be significant. For example, Lambeck and Johnson [21] give a constant viscosity of 1022 Pa s in the lower mantle while Forte and Mitrovica [17] give three orders of magnitude change with a maximum around 2000 km depth. Steinberger [37,38] gives an increase of a factor of 40 to 4×1022 Pa s.

To obtain estimates of sinking time and slab width, I assume that viscosity increases with depth as in Eq. (14). I consider increases of 2 and 3 orders of magnitude from 700 to 2700 km depth, to leave 200-km depth for D″. The depth scales Dapp are 434 and 290 km, respectively. The sinking times from Eq. (15) assuming that the vertical plate velocity at ∼700 km, Vref, is 50 km Myr−1, are 173 and 367 Myr, respectively. The material spends 63% of this time in the bottom 2D app of the mantle, so the times to reach this ponding level are 63 and 135 Myr. Slab material sinks through the lower mantle in times that are short compared with the 2000-Myr-cycle time of the mantle. The sinking times for two-orders-of-magnitude contrast are reasonable. The sinking times for three-order-of-magnitude contrast are excessive in that tomograms correlate with the last 120–200 Myr of subduction [7]. The two viscosities imply factors of 10 and 30 broadening of the slab in Eq. (15) or increases in width from ∼200 to ∼2000 and ∼6000 km, respectively. The latter increase is large enough to grossly violate the assumption in the force balance that the slab is thin, compared to the other length scales in the mantle. Rather, the slab would sink through the relatively low-viscosity regions in the uppermost lower mantle and then pond and spreads out slowly in the lowermost lower mantle. The sinking time of 367 Myr is thus an overestimate, indicating that three orders of magnitude of increase are excessive. I thus use the 2-order-of-magnitude increase as a plausible rounded estimate for a significant change in viscosity in subsequent calculations. I return to the gradual spreading of deeply ponded slabs in Section 4.

3.2 Heat conduction and the temperature of deep slabs

Conduction tends to spread out the thermal anomaly associated with the slab. I constrain the average thermal anomaly of dead slabs and hence the effective amount of material that is entrained with slabs. I begin with a slab (that started as a line source) that broadens only by conduction. Mathematically, the temperature is proportional to exp (-x2/4κt), where κk/ρC is thermal diffusivity and x is the distance from the center of the slab (Eq. (4-163) of [43]). The half width of the anomaly is 2κt. The full width of the slab at some time is then:

W=aκt(16)
where the constant a≈4. The change in slab thickness with time relative to its width is
DW Dt =a2κ2W(17)
In the case of interest, the slab broadens mainly mechanically as it penetrates regions of progressively higher viscosity. The relative change in slab width from conduction alone thus determines the average temperature change within the material. That is,
DW W Dt =- DT T Dt =a2κ2W2(18)
where W is the local width of the slab. I express Eq. (18) in terms of depth by noting that:
tZ=1VS(19)
which yields:
WWZ=-TTZ=a2κ2(VW)W=a2κ2(V ref W ref )W(20)
where conservation of material gives the final equality as in Eq. (10). The relative change in temperature by conduction is, using Eqs. (10) and (14):
ln T ref T final =Z0a2κ2V ref W ref 2 exp -z2Z app dza2κZ app V ref W ref 2(21)
where the temperature anomaly is Tref at the reference depth 0, Tfinal is the final temperature anomaly at depth Z beneath the reference depth, and the upper limit of the integral is infinity in the final approximation. Little heat conduction occurs in the lower mantle if viscosity increases with depth. For example, I let the viscosity depth scale Zapp be 434 km, the initial width of the slab Wref be 200 km, the vertical plate velocity at ∼700 km Vref be 50 km Myr−1, and the thermal diffusivity be 10−6 m2 s−1. The relative temperature change is 11%, a modest amount. That is, the average temperature anomaly broadens little from conduction between the time the slab enters the lower mantle until the time where it ponds at near the base of the mantle amongst slightly older slabs.

Viewed in another way, the breadth of the thermal anomaly of a slab and the volume of cool material within it result mostly from conduction while the slab was still oceanic lithosphere rather than conduction after subduction. The same feature exists for mantle plumes [15]. Starting plume heads and plume conduits entrain only minor amounts of material once these leave the boundary layer. This gives some confidence in ignoring entrainment during buckling and folding of slabs because starting plume heads deform greatly during their ascent.

The 400-K, 200-km thick average anomaly in slabs and the cycle time of Δt of 2000 Myr that I obtained in Section 2 are reasonable round estimates since significant entrainment does not occur at midmantle depths. A combination of 300 K, 1500 Myr, and 267 km would be equally good. Geochemical data and our knowledge of past plate tectonics are not now good enough to get a refined estimate.

4 Slow spreading of slabs graveyards

Broad regions of high seismic velocity and low attenuation are prominent features in the deep lower mantle in global tomograms. They occur where many slabs have descended over that last ∼180 Myr for which reliable reconstructions are available [6,26,28,37]. Conversely, superplume regions of low seismic velocity and high attenuation exist elsewhere where slabs have not impinged. Surface hot spots preferentially overlie the slab-free regions [9].

The first order physics of the process are fairly simple. Newly subducted slab material ponds within recently arrived slabs near the base of the mantle. The viscosity of the lower mantle is on average a factor of ∼20 higher than the upper mantle so that the slabs sink significantly slower than surface plate rates [6,44]. The viscosity of the slab graveyard and the surrounding mantle is sufficient to prevent it from quickly spreading laterally along the base of the mantle like dense oceanic bottom water. Rather tomography detects dead slabs in the deep mantle as well as currently sinking slabs.

4.1 Time scale of slab graveyards

The existence of slab graveyards constrains the viscosity of the deep mantle. That is, a graveyard grows as a pile until lateral flow comes into balance with subduction [6] (Fig. 4). I represent the effect by denoting the local elevation of the top of the slab pile relative to its global mean as H. I denote a term of its spherical harmonic expansion of order m as Hm. Assuming a linear fluid where lateral viscosity variations are unimportant, a representative harmonic term evolves as

Hmt=-Hmtm+Fm(22)
where tm is time constant and Fm is the flux (equivalent thickness per time) supplied by slabs for the harmonic. At steady state, the height of the harmonic is
[Hm] ss =Fm/tm(23)
This constrains decay time because the flux Fm is known from reconstructions over the last 200 Myr and the height Hm is known from tomography [28].

Fig. 4

Schematic diagram of slab graveyard. The pile spreads out because it is negatively buoyant and is continually restored by slabs. I represent the relief on the top of the pile with a series of spherical harmonics.

Schéma de cimetière de plaques. Les plaques accumulées à la base du manteau s'étalent, car elles subissent une poussée d'Archimède négative et sont continuellement remplacées par les plaques plongeantes. Je représente le relief de l'empilement par un développement en harmoniques sphériques.

Low-order harmonics, degrees 1 to 3, are prominent on lower mantle tomograms. The harmonic of degree 2 dominates the geoid (the first harmonic of the geoid is by definition zero). The spectrum of hot spot distribution (as points on sphere) peaks mildly at 6 and at low harmonics 1 to 3 [27]. The dominance of lower spherical harmonics in tomograms and the geoid is both a consequence of forcing function Fm and the time scale tm.

Beginning with the forcing function, the spherical harmonic expansion for subducted material expressed as predicted geoidal effect peaks at low harmonics [27]. This occurs because the positions of subduction zones have changed little over 180 Myr compared with the 10000-km length scales of the lowest harmonics so the low-harmonic phase and magnitude of Hm have not changed much [28]. Conversely, plate boundaries have moved significantly relative to higher harmonics over this time although the actual change is not well resolved. Overall, the time scale for the low harmonics cannot be far greater than 180 Myr, because features before that major plate reorganization would be evident on tomograms and not correlate with more recent subduction. In any case, it cannot approach the mantle cycle time of 2000 Myr.

With regard to the time scale, slab piles are within the highest viscosity region of the mantle. The restoration of stable interface of a slab pile is analogous to the growth of an instability on an unstable interface or to the restoration of a stable surface by glacial rebound. The presence of the core causes scale time to be higher at low harmonics (than intermediate ones) even for an isoviscous mantle [2,31,35]. However, a strong depth gradient of viscosity within the lower mantle is not necessary for tomograms to correlate with past subduction. For example, Bunge et al. [7] assume an isoviscous lower mantle, 4×1022 Pa s, and obtain a time scale of 150 Myr and a reasonable fit to observations.

Overall, the persistence of slab piles, the correlation of mid-mantle and deep-mantle slab anomalies with subduction over the last 200 Myr, and the increasing width of slab anomalies with depth in Section 3 all provide evidence that the viscosity of the mantle increases with depth. More numerical calculations like those of Bunge and Richards [6] and Bunge et al. [7], and Steinberger [37] are needed along with better tomographic resolution. Out lack of knowledge of ancient plate motions may remain a difficulty.

4.2 Slab impingement in the deep mantle and its effect of hot spot motion

Flow in the mantle prevents hot spots from being fixed features. Rather, relief on the top of D″, including the cusps beneath plume conduits, and the plume conduits themselves advect with the flow [38,39]. The observed rates of hot spot motion are small compared to surface plate velocities because the viscosity in the mantle above D″ is high compared with that of the upper mantle [25]. Exceptions to the rule are interesting, as they are likely to be detected. I concentrate on an effect that may occur near the base of the mantle.

As noted already, tomograms indicate that the lower-mantle regions between slabs are relatively warm. I consider what happens when subduction begins above such a region. First, the slab sinks rapidly because the surrounding viscosity in Eq. (15) is relatively low. Once to the bottom, the slab material spreads rapidly, again because the surrounding viscosity is relatively low in Eq. (24). The process may cause rapid hot spot wander. The initiation of subduction may somewhat kink a track as the geometry of the return flow immediately adjusts to the new surface boundary condition. The plume cusp may also move as the slab graveyard impinges on it. That is, the cusp advects away from the slab graveyard and toward the adjacent superplume region. Plumes affected in this manner may move significantly relative to the average hot spot frame. However, the mechanism is unlikely to produce a sharp link in a mid-plate track, like the Hawaii-Emperor bend. Once established, a slab graveyard and its buoyancy evolve slowly. There is also no obvious place for a slab graveyard near the Hawaiian hot spot to make the Hawaii-Emperor bend. I do not have a good example of plume deflection by a new slab graveyard.

Conversely, information may be gained by careful tomography. For example, new slabs are potentially distinct from older slabs that impinge on slab piles. New slabs sink into relatively warm, low-viscosity mantle. New slabs should sink rapidly in Eq. (15), widen more slowly in Eq. (10), and have a greater temperature contrast with their surroundings than old slabs. They also lack the underlying negative buoyancy of an established slab pile [37].

4.3 Kinematics of the return flow

The long-term kinematics of flow within the deep mantle are of interest to geochemistry. Crudely, plumes sample subducted material including sediments, oceanic basalt, oceanic gabbro, and depleted residuum with model ages of ∼0.7 to 1.9 Ga (e.g., [5,14]). Components of this age are compatible with the kinematics in Fig. 1 [15]. Old plume components, like 3He, and ‘young’ slab components, indicated by Pb isotopes [10], require more explanation. It is not my purpose to second guess the conclusions of geochemists, rather I show that I expect a mixture of young, intermediate, and old components in plumes supplied by the low-viscosity D″ layer.

To begin with observations, tomograms constrain the gross flow associated with slabs: slabs sink deep into the mantle. As noted in Section 2, the return flow to ridges is on average gradual advection of this material back to the surface. The return flow rate (normalized as 4πr2ρVZ, where r is the local radius from the Earth's center and ρ is the local density) is more-or-less constant between the base of the lithosphere and the depth range where slabs pond. Conversely, plumes tap material from D″ and take it to shallow depths. The return flow into D″ displaces material downward throughout the (sublithospheric) mantle with a uniform normalized flow rate. The normalized flow rate is the sum of the slab and plume effects. It is downward within and immediately above D″ and upward in the rest of the mantle (Fig. 2).

A stagnation region forms a dynamic trap between the upward and downward flowing regions in the deep mantle. Material in this region undergoes a random walk from the vagaries of slab and plume activity. It eventually ends up either in plumes or at ridges. In the former case, it contributes an old component to plumes.

The young component in plumes arises because slabs made from old oceanic lithosphere sink more rapidly and deeply than average slabs [7]. That is, the larger initial width Wref implies a smaller sinking time in Eq. (15). In addition, new slabs made of old lithosphere penetrate relatively low viscosity mantle and sink more deeply. In both cases, slabs impinge on the top of D″ and contribute some of their material to plumes.

5 Ascent of starting plume heads

Starting plume heads are tempting tomographic targets as they should be several hundred kilometers in diameter (Fig. 1). In contrast, plume conduits are narrow and difficult to detect in the lower mantle. Lower mantle ray paths through plume conduits can never be first arrivals and are missed by travel-time tomography. Seismologists cannot easily separate immediately late conduit arrivals from the rest of the coda. They require large arrays to reliably find very late low-amplitude scattered and diffracted arrivals.

5.1 Plume head size and ascent time

I constrain the probable depth range and size range of starting plume heads by modifying a well-known analytical theory reviewed by Schubert et al. [32]. I represent the plume head as a zero-viscosity buoyant sphere with a density contrast Dρ and a radius rh. The Stokes-law viscosity of the sphere is:

Vh=Δρgrh23η(24)
where η is the viscosity of the surrounding material [32, p. 521]. I ignore growth of the plume head by processes, like conduction of heat and entrainment of the surroundings, which leave its thermal and density anomaly Δρrh3 unchanged. The velocity then varies inversely with the radius. For simplicity, the plume head grows at a constant volume rate Q. This implies that:
rh=3tQ4π1/3(25)
where t is the time after the plume head starts filling. The ascent velocity is:
Vh=Δρg3η3tQ4π2/3(26)
In the case that lower mantle viscosity is constant, the time for the plume to ascend to some height H above the base of the mantle is:
tH=5ηHΔρg3/54π3Q2/5(27)
[45, Eq. (21)] and the radius at that time is:
rH=15QHη4πΔρg1/5(28)
Eqs. (27) and (28) suffice to predict what can be expected on tomograms if viscosity does not change much with depth in the lower mantle. Both the ascent time and the radius are weakly dependent on physical parameters. Low-flux plume heads take longer to ascend than high flux ones and should have a higher standing crop in the mantle at any time relative to the number that arrive per time at the surface. Plume heads ascend slowly at first, but the effect is minor. For example, the plume head needs 66% of its ascent time to ascend halfway up. At a random observation time, plume heads are more or less evenly likely over depth in the lower mantle.

I compare the isoviscous situation with a lower mantle where viscosity is strongly dependent on depth. The velocity is then from Eq. (14):

Vh=Δρg3ηB3tQ4π2/3 exp HD app (29)
where the viscosity is ηB near the base of the mantle. I obtain the ascent time by integrating
H0 exp -hD app dh=Δρg3ηB3Q4π2/3tH0t2/3dt(30)
where h is a dummy variable for H. In the case of a large viscosity gradient, the ascent time to the surface is approximately that to H=∞,
t=5ηBD app Δρg3/54π3Q2/5(31)
and the final radius is
r=15QD app ηB4πΔρg1/5(32)
These formulas are homologous to those for an isoviscous mantle in Eqs. (27) and (28), replacing H with Dapp. The weak dependence of plume head radius on physical parameters in the isoviscous case (28) carries through to the case of a strong viscosity gradient.

The likely position of the plume head differs markedly between a strong viscosity gradient and an isoviscous lower mantle. The time to ascend through a viscosity gradient to height H is from Eqs. (30) and (31):

tH=t1- exp -HD app 3/5(33)
The plume head spends 76% of it ascent time in the bottom Dapp of the mantle and 92% of its time in the bottom 2D app . It is quite near its final radius once it ascends above this basal region. That is, a random observation is most likely to find starting plume heads near the base of the mantle and unlikely to find them in shallower regions that they quickly traverse.

The standing crop of plume heads is inversely proportional to the ascent time. I first obtain an estimate for an isoviscous (1022 Pa s) lower mantle. I let the temperature anomaly be 800 K (400 K relative to the MORB adiabat), the density be 5500 kg m−3, and thermal expansion coefficient be 10−5 K−1, the acceleration of gravity be 10.5 m s−2, and the lower mantle ascent distance be 2000 km. I let volume flux be 200 m3 s−1 in the upper range of estimates for individual plumes tabulated by Sleep [34]. This yields an ascent time of 42 Myr and a final radius of 400 km. The ascent time t for Dapp=434 km and ηB=1023 Pa s is 68 Myr and the final radius is 469 km. There is no significant difference between the two cases given the uncertainty in the parameters. My rounded estimates are that the ascent time is ∼50 Myr and that the final head radius is ∼400 km (as shown in Fig. 1).

Current constraints on plume head ascent time are vague because these features are only indirectly detected through their flood basalts. The ascent time needs to be long enough that plume heads have as few 107 years worth of plume flux to produce widespread flood basalts. It cannot be much longer than 108 years without destroying the negative correlation between plume location and slab graveyards.

The expected standing crop of plume heads within our current mantle is small. There have been ∼4 starting plume heads in the last 100 Myr, Yellowstone, East Africa, Réunion Island, and Iceland, or one per 25 Myr. Comparing this with the 50-Myr ascent time, the current standing crop is most likely 2±2. If there is a strong viscosity gradient with depth, all these starting plume heads are most likely near the base of the mantle.

5.2 Ascent of starting plumes from mid-mantle depths

It is conceivable that plumes start from mid-mantle depths, for example, above a thick layer of dense dregs [11,20]. This layer should have significant effects on the thermal and seismic velocity structure of the mantle [42].

I consider its effects on starting plume heads. Plumes should ascend from the peaks of the dreg layer. Such peaks are not resolved in tomograms and a thick layer of dregs may not exist [42]. In particular, there is no evidence that starting plume heads are now forming in such regions or that dregs peaks underlie hot spots.

A potential testable implication is that plumes would ascend from different depths above peaks in a thick dregs layer, but from essentially the same depth if no (or only a thin) dregs layer exists [42]. In the case that the viscosity is not strongly depth dependent, different starting depths weakly affect the ascent time in Eq. (27) and the final radius in Eq. (28) through the ascent distance H. Modest variations in the volume flux Q would mask any real effect in the limited number of real cases. Conversely, the starting height in a mantle with strongly depth-dependent viscosity has strong effects through the viscosity at the starting depth ηB. For example, a change in a factor of 10 of viscosity would change the final volume of the plume head by a factor of 4 in Eq. (32). One would need estimates of both the volume flux Q and the plume head volume to see if the correlation in Eqs. (28) and (32) occurs or if there are effects that might be attributed to variations in the depth to peaks of the dreg layer. At present, these quantities are not precisely constrained even for Iceland and Réunion Island.

5.3 Ascent of superplume regions

The formalism in Section 5.1 is applicable to the ascent of superplume regions detected by tomography. Romanowicz and Gung [30] consider that these regions are buoyant relative to upper mantle MORB adiabat, not just relative to their lateral surroundings. The observed sinking of slabs and ascent of plume heads through the mantle provides some constraint on whether superplume regions are in fact buoyant relative to the MORB adiabat.

The radius and the density contrast differ between superplume regions and plume heads in Stokes-law ascent rate (Eq. (25)). The radius of the superplume regions is a significant part of the lower mantle ∼1000 km compared with ∼400 km for starting plume heads. This implies that the ascent rates would be equal if the density contrast was a factor of 6.25 greater in plume heads than in starting plumes relative to the mantle adiabat. The ascent rate needs to be a factor of few slower than this so that the regions persist on the time scale that slabs form graveyards between them. This implies that the density contrast of superplume regions relative to the mantle adiabat cannot be more than a few percent of that of plume heads. That is, the temperature contrast of superplume regions with the MORB adiabat is less than ∼20 K compared with ∼400 K of plume heads. This situation creates a special circumstances issue. From Section 2, a typical part of the lower mantle heats up 20 K relative to the evolving MORB adiabat in 100 Myr. It is unlikely, though not impossible, that we live in the time that much of the lower mantle is just about to become unstable heating up over ∼2000 Myr. More likely, the superplume regions are cooler than the MORB adiabat, but hotter than their slab-dominated surroundings.

6 Conclusions

Present tomograms resolve the gross features of the lower mantle [3,19]. Resolved features are slabs and warmer regions between them. Unresolved features include starting plume heads and plume tail conduits. Given this situation, how does one proceed with geodynamics and tomography? I have compiled some simple but robust features of mantle-wide convection to aid in this effort.

Plate tectonics cool the mantle (Fig. 1). The heat comes from cooling of the Earth over time and from radioactive decay. The core supplies a small amount of heat that ascends in plumes. Surface plate tectonics does not efficiently stir the mantle because the viscosity of the lower mantle is much greater than that of the upper mantle. The stirring that does occur involves sinking and ponding of slabs in the deep mantle and the eventual return flow of this material to ridge axes on a time scale of 2000 Myr.

A basic feature of this geometry is that the geotherm is subadiabatic between the upper lithospheric boundary layer and the lower D″ boundary layer (Fig. 2). The total subadiabatic temperature change is ∼400 K. The deep mantle above D″ consists of graveyard piles of ponded slabs and somewhat warmer superplume regions between them. The superplume regions are likely to be cooler and negatively buoyant with respect to the upper mantle (MORB) adiabat. Tomograms plotted with respect to the mean velocity at each depth do not represent the temperature structure relative to an adiabat.

The viscosity may change by about two orders of magnitude in the lower mantle, a lot compared with other physical parameters. An observable effect is that slabs should widen proportionally to the square root of viscosity and its descent rate should decrease inversely with the square root of viscosity (Figs. 1 and 3). The temperature anomalies associated with slabs broaden little by conduction while they traverse the lower mantle. The temperature and seismic anomaly should become less distinct as slabs penetrate into (slightly longer ago) subducted material near the base of the mantle.

Graveyards of ponded slabs persist near the base of the mantle on a time scale comparable to that for plate motions to reorganize, ∼120–200 Myr (Fig. 4). Significantly shorter response times preclude slab piles while significantly longer ones would cause tomograms to correlate poorly with relatively recent subduction. The longest wavelengths are most evident because they persist longer and because the source of the material from slabs stays in phase. The viscosity of the lower mantle obtained from these considerations is ∼4×1022 Pa s.

Starting plume heads traverse the lower mantle on a time scale of ∼50 Myr. In the case that viscosity increases with depth, they spend most of their ascent time in the base of the mantle. They traverse an isoviscous mantle at a relatively constant rate. The size of starting plume heads depends weakly on physical parameters including the starting depth. At most a few starting plume heads exist in the present lower mantle.

My inferences provide hints on how to formulate convection calculations. The temperature and depth dependence of viscosity is important to the existence of plate tectonics and plumes in the first place. The widening of slabs with depth is a simple feature that correlates directly with viscosity. Plates need to be included in a ‘damage’ rheology so that weak plate boundaries form and evolve over time [40,41]. This lets plate geometries and the rate of plate turnover be strongly influenced by shallow processes. It lets the subadiabatic thermal gradient in the mid-mantle build-up over time and its stable stratification to affect the flow. Ideally the calculations should resolve D″ and plumes heads and tails. The successful generic models then can be used to formulate models of modern plate motions.

Tomographic modeling may proceed with geodynamic constraints imposed. For example, deep slabs broaden mechanically. They should be relatively isothermal and have thin conductive boundaries. Buoyant regions and negatively buoyant regions not associated with slabs should be rare in the mid-mantle as they already should have ascended or sunk out of this region if they somehow formed [3]. Optimally the density, temperature, and viscosity structure should be consistent with tomography, geoid, long-term fluid dynamics, and short-term glacial rebound. Simple scaling relationships help one get started on this task. They also let one effectively use existing tomograms and dynamic models.

Acknowledgements

This research was in part supported by NSF grants EAR-9902079 and EAR-0000747. I thank Barbara Romanowicz for a helpful review.


References

[1] D.L. Abbott; L. Burgess; J. Longhi; W.H.F. Smith An empirical thermal history of the Earth's upper mantle, J. Geophys. Res., Volume 99 (1994), pp. 13835-13850

[2] D.L. Anderson; R. O'Connell Viscosity of the Earth, Geophys. J. R. Astron. Soc., Volume 14 (1967), pp. 2287-2295

[3] T.W. Becker; L. Boschi A comparison of tomographic and geodynamic mantle models, Geochem. Geophys. Geosyst., Volume 3 (2002) (2001GC000168)

[4] M.J. Bickle Implications of melting for stabilization of the lithosphere and heat loss in the Archean, Earth Planet. Sci. Lett., Volume 80 (1986), pp. 314-324

[5] J. Blichert-Toft; F.A. Frey; F. Albarède Hf isotope evidence for pelagic sediments in the source of Hawaiian basalts, Science, Volume 285 (1999), pp. 879-882

[6] H.-P. Bunge; M.A. Richards The origin of large-scale structure in mantle convection: effects of plate motions and viscosity stratification, Geophys. Res. Lett., Volume 23 (1996), pp. 2987-2990

[7] H.-P. Bunge; M.A. Richards; C. Lithgow-Bertelloni; J.R. Baumgardner; S.P. Grand; B. Romanwicz Time scales and heterogeneous structure in Earth models, Science, Volume 280 (1998), pp. 91-95

[8] H.-P. Bunge; Y. Ricard; J. Matas Non-adiabaticity in mantle convection, Geophys. Res. Lett., Volume 28 (2001), pp. 879-882

[9] C.G. Chase; D.R. Spowl The modern geoid and ancient plate boundaries, Earth Planet. Sci. Lett., Volume 62 (1983), pp. 314-320

[10] G. Class; S.L. Goldstein; S.J.G. Galer; D. Weis Young formation age of a mantle plume source, Nature, Volume 362 (1993), pp. 715-721

[11] A. Davaille Simultanteous generation of hot spots and superswells by convection in a heterogeneous planetary mantle, Nature, Volume 402 (1999), pp. 756-760

[12] G.F. Davies Ocean bathymetry and mantle convection. I. Large-scale flow and hot spots, J. Geophys. Res., Volume 93 (1988), pp. 10467-10480

[13] G.F. Davies Cooling the core and mantle by plume and plate flows, Geophys. J. Int., Volume 115 (1993), pp. 132-146

[14] J. Eisele; M. Sharma; S.J.G. Galer; J. Blichert-Toft; C.W. Devey; A.W. Hofmann The role of sediment recycling in EM-1 inferred from Os, Pb, Hf, Sr isotope and trace element systematics of Pitcarn hot spot, Earth Planet. Sci. Lett., Volume 196 (2002), pp. 197-212

[15] C.G. Farnetani; B. Legras; P.J. Tackley Mixing and deformations in mantle plumes, Earth Planet. Sci. Lett., Volume 196 (2002), pp. 1-15

[16] A.M. Forte Seismic-geodynamic constraints on mantle flow: implications for layered convection, mantle viscosity, and seismic anisotropy in the deep mantle, Geophys. Monogr., Volume 117 (2000), pp. 3-36

[17] A.M. Forte; J.X. Mitrovica Deep-mantle high-viscosity flow and thermochemical structure inferred from seismic and geodynamic data, Nature, Volume 410 (2001), pp. 1049-1056

[18] S.J.G. Galer; K. Mezger Metamorphism denudation and sea level in the Archean and cooling of the Earth, Precambrian Res., Volume 92 (1998), pp. 387-412

[19] E.J. Garnero Heterogeneity of the lowermost mantle, Annu. Rev. Earth Planet. Sci., Volume 28 (2000), pp. 509-537

[20] L.H. Kellogg; B.H. Hager; R.D. van der Hilst Compositional stratification of the deep mantle, Science, Volume 283 (1999), pp. 1881-1884

[21] K. Lambeck; P. Johnson The viscosity of the mantle: Evidence from analyses of glacial-rebound phenomena (I. Jackson, ed.), The Earth's Mantle, Composition, Structure, and Evolution, Cambridge University Press, 1998, pp. 461-502

[22] G. Masters; G. Leske; H. Bolton; A. Dziewonski The relative behavior of shear velocity, bulk sound speed, and compressional velocity in the mantle: implications for chemical and thermal structure, Geophys. Monogr., Volume 117 (2000), pp. 63-87

[23] S.V. Panasyuk; B.H. Hager Inversion for mantle viscosity profiles constrained by dynamic topography and the geoid, and their estimated errors, Geophys. J. Int., Volume 143 (2000), pp. 821-826

[24] G. Pari; W.R. Peltier Global surface heat flux anomalies from seismic tomography-based models of mantle flow: Implications for mantle convection, J. Geophys. Res., Volume 103 (1998), pp. 23743-23780

[25] M.A. Richards Hot spots and the case for a high viscosity lower mantle (R. Sabadini; K. Lambeck; E. Boschi, eds.), Glacial Isostasy, Sea Level and Mantle Rheology, NATO ASI Ser. C, 3344, Kluwer Academic Publishers, Norwell, MA, 1991, pp. 571-587

[26] M.A. Richards; D.C. Engebretson Large scale mantle convection and the history of subduction, Nature, Volume 355 (1992), pp. 437-440

[27] M.A. Richards; B.H. Hager; N.H. Sleep Dynamically supported geoid highs over hot spots: observation and theory, J. Geophys. Res., Volume 93 (1988), pp. 7690-7708

[28] M.A. Richards; Y. Ricard; C. Lithgow-Bertelloni; G. Spada; R. Sabadini An explanation for Earth's long term rotational stability, Science, Volume 275 (1997), pp. 372-375

[29] F.M. Richter; D.P. McKenzie On some consequences and possible causes of layered mantle convection, J. Geophys. Res., Volume 86 (1981), pp. 6133-6142

[30] B. Romanowicz; Y. Gung Superplumes from the core mantle boundary to the lithosphere: Implications for heat flux, Science, Volume 296 (2002), pp. 513-516

[31] S.K. Runcorn Changes in the convection pattern in the Earth's mantle and continental drift: Evidence for a cold origin of the Earth, Phil. Trans. R. Soc. Lond., Volume 258 (1962), p. 228

[32] G. Schubert; D.L. Turcotte; P. Olson Mantle Convection in the Earth and Planets, Cambridge University Press, Cambridge, UK, 2001 (p. 940)

[33] J.G. Sclater; C. Jaupart; D. Galson The heat flow through oceanic and continental crust and the heat loss of the Earth, Rev. Geophys., Volume 18 (1980), pp. 269-311

[34] N.H. Sleep Hot spots and mantle plumes: Some phenomenology, J. Geophys. Res., Volume 95 (1990), pp. 6715-6736

[35] N.H. Sleep Time dependence of mantle plumes: Some simple theory, J. Geophys. Res., Volume 97 (1992), pp. 20007-20019

[36] F.D. Stacey A thermal model of the Earth, Phys. Earth Planet. Inter., Volume 15 (1977), pp. 341-348

[37] B. Steinberger Slabs in the lower mantle: results of dynamic modeling compared with tomographic images and the geoid, Phys. Earth Planet. Inter., Volume 118 (2000), pp. 241-257

[38] B. Steinberger Plumes in a convecting mantle: models and observations for individual hot spots, J. Geophys. Res., Volume 105 (2000), pp. 11127-11152

[39] B. Steinberger; R.J. O'Connell Advection of plumes in mantle flow: implications for hot spot motion, mantle viscosity and plume distribution, Geophys. J. Int., Volume 132 (1998), pp. 412-434

[40] P.J. Tackley Self-consistent generation of tectonic plates in time-dependent, three-dimensional mantle convection simulations. 1. Pseudoplastic yielding, Geochem. Geophys. Geosyst., Volume 1 (2000) (2000GC000036)

[41] P.J. Tackley Self-consistent generation of tectonic plates in time-dependent, three-dimensional mantle convection simulations. 2. Strain weakening and asthenosphere, Geochem. Geophys. Geosyst., Volume 1 (2000) (2000GC000043)

[42] P.J. Tackley Strong heterogeneity caused by deep mantle layering, Geochem. Geophys. Geosyst., Volume 3 (2002) (2001GC000167)

[43] D.L. Turcotte; G. Schubert Mantle Geodynamics Applications of Continuum Physics to Geologogical Problems, John Wiley, Cambridge, UK, 1982 (p. 450)

[44] S. Wang; R. Wang Current plate velocities relative to hot spots: implications for hot spot motion, mantle viscosity and global reference frame, Earth Planet. Sci. Lett., Volume 189 (2001), pp. 133-140

[45] J.A. Whitehead; D.S. Luther Dynamics of laboratory diapir and plume models, J. Geophys. Res., Volume 80 (1975), pp. 705-717


Comments - Politique