Compositional Stratification in the Deep Mantle

See allHide authors and affiliations

Science  19 Mar 1999:
Vol. 283, Issue 5409, pp. 1881-1884
DOI: 10.1126/science.283.5409.1881


A boundary between compositionally distinct regions at a depth of about 1600 kilometers may explain the seismological observations pertaining to Earth's lower mantle, produce the isotopic signatures of mid-ocean ridge basalts and oceanic island basalts, and reconcile the discrepancy between the observed heat flux and the heat production of the mid-ocean ridge basalt source region. Numerical models of thermochemical convection imply that a layer of material that is intrinsically about 4 percent more dense than the overlying mantle is dynamically stable. Because the deep layer is hot, its net density is only slightly greater than adiabatic and its surface develops substantial topography.

Several fundamental constraints must be satisfied by a successful model of the dynamics and thermochemical structure of Earth's mantle. The model must produce sufficient heat, either by radioactive decay or by cooling of the planet, to account for the inferred global heat flux. The model must be capable of producing the rich variety and the observed systematics of geochemical signatures in mantle-derived basalts (1). The model must be consistent with inferences from seismic tomography that some subducted slabs extend to near the base of the mantle (2) and that the lowermost mantle is characterized by long wavelength structure (3, 4) and complex relations between the bulk sound and shear wavespeed (5, 6) [see (7) for an overview]. Finally, the model must be dynamically consistent. Here, we present a model that is dynamically feasible and satisfies the essential geochemical and geophysical observations. It differs from many previous models by placing a boundary between compositionally distinct mantle regions deep in the lower mantle, rather than at a depth of 660 km.

The characteristic isotopic ratios of mid-ocean ridge basalts (MORB) and oceanic island basalts (OIB) provide evidence for a suite of distinct reservoirs in the mantle (1). These reservoirs and signatures are thought to be produced by the formation and recycling of oceanic crust and lithosphere, plus small amounts of recycled continental crust. In addition, 129Xe,3He/4He, and 40Ar contents of the mantle (8–10) indicate that the mantle has not been entirely outgassed.

87Sr/86Sr and143Nd/144Nd isotope ratios of the crust and MORB have been used to estimate the mass of mantle from which the crust was extracted, and hence to infer the mass of the remaining, less depleted component. Estimates for the mass of the depleted mantle range from 25% (11) (coincidentally the mass of the mantle above the 660-km discontinuity) to 90% (1). Similar mass balance arguments are used to determine the amount of mantle that must have been outgassed to produce the 40Ar in Earth's atmosphere (10); these predict a volume of degassed mantle of ∼50%. Uncertainties arise because the K/U ratio of Earth is still under debate (12) and the lower crust or the undegassed parts of the mantle have retained substantial amounts of 40Ar (13), or some Ar may be recycled.

Another fundamental constraint is provided by Earth's heat budget (14, 15). Of the 44 TW (16) of the present-day heat flux out of Earth, 6 TW is generated within the crust by radioactive decay of U, Th, and K, and 38 TW must be provided either by generation of heat within the mantle and core or by cooling of the planet (17). For example, if Earth had the radiogenic heat production of the average chondritic meteorite, the total heat production would be 31 TW; the remaining 13 TW would be provided by cooling of the planet by 65 K per 109 years. Geochemical analyses of basalts, however, show that the source region of MORBs is depleted in heat production by a factor of 5 to 10 relative to a chondritic silicate value (18). Thus, if the MORB source region made up most of the mantle, the mantle heat production would be only 2 to 6 TW, comparable to that of the crust. Matching the observed heat flux would require rapid cooling of the planet by, on average, 175 K per 109 years, which requires excessive internal temperatures during the Archaean (19).

Hence, there must be an extra heat source. The D" region, a layer of anomalous seismic velocities several hundred kilometers thick at the base of the lower mantle, is likely to be compositionally distinct (20), but it can account for only a small fraction of the global heat flux unless there is extreme internal heat production. The latter is unlikely, in particular if theD" layer contains foundered oceanic crust (21), which has a heat production comparable to that of the chondritic silicate Earth (17), or some core material (22), which is likely low in U, Th, and K.

The energy calculations and geochemical mass balances both suggest that the mantle is composed of several reservoirs: a depleted region, which is the source of MORB; a region that, relative to the MORB source, is undepleted in incompatible elements; and smaller volumes containing fragments of subducted lithosphere (including crust). Because mid-ocean ridges sample the shallow mantle over a wide area of the surface and have relatively uniform composition, it is generally assumed that the upper mantle is depleted. As a consequence, a wide range of models have been based on the assumption that the undepleted reservoir corresponds to the lower mantle (1) but that some plumes rise through this region (23) or originate at its lower boundary (1,24). Models involving layering of convective flow with a boundary at 660 km are, however, at odds with the seismological evidence (2).

We propose instead that the transition in seismological heterogeneity observed at ∼1600 km depth (7) is the manifestation of a boundary between the depleted MORB source region and enriched, intrinsically dense material in the deeper mantle (Fig. 1). A compositionally distinct layer in the deep mantle would provide an explanation for long-lived geochemical reservoirs; with about 20 to 30% of the mantle volume, it is sufficiently massive to contain the 40Ar and 3He that have not been lost to the atmosphere. In the transition zone above the dense layer, recycled subducted material, mixed with depleted and enriched mantle, would form the complex family of geochemical components seen in oceanic basalts. In this model, the deep, dense layer provides the missing heat production needed to explain the observed heat flux. Such a layer could develop during the early differentiation of Earth, by processes associated with a deep magma ocean, or by formation and recycling of a mafic crust in the Archaean (25).

Figure 1

Diagram illustrating the possible dynamics of an intrinsically dense layer in the lower mantle. Depth to the top of the layer ranges from ∼1600 km to near the CMB, where it is deflected by downwelling slabs. Internal circulation within the layer is driven by internal heating and by heat flow across the CMB. A thermal boundary layer develops at the interface, and plumes arise from local high spots, carrying recycled slab and some primordial material.

To determine whether this model is dynamically plausible, we used a finite-element model of mantle convection (26) modified for thermochemical convection (27) to examine the interaction of a subducting slab with an intrinsically dense layer in the deep mantle. Viscosity is temperature- and pressure-dependent (28), and plate-like behavior is induced in the cold surface boundary layer by providing weak zones at the upper corners of the box (29). An isothermal boundary condition at the core-mantle boundary (CMB) heats the mantle from below at a Rayleigh number of 2 × 107 (30). The mantle is also heated from within; in the upper layer, the internal heating Rayleigh number is Ra h = 1.16 × 108 [corresponding to a U concentration of 7 ppb with Th/U = 2.5 and K/U = 104 (18, 31)]. The deep layer is enriched in heat-producing elements withRa h = 5.12 × 108, corresponding to a U concentration of 25.6 ppm with Th/U = 4.

Figure 2 shows a snapshot of the resulting temperature and composition (Fig. 2A) and viscosity (Fig. 2B) for the mantle after a subducting slab has penetrated the intrinsically dense layer (32). In this simulation—for which we used a buoyancy number B = ΔρC0αTΔT = 0.9 (corresponding to an intrinsic density contrast of ∼4%), where ΔρC is density change due to composition alone, ρ0 is the reference density, αT is the coefficient of thermal expansion, and ΔT is the temperature drop across the mantle—the dense layer is stable for the duration of the calculation. Because B < 1, this model is in the regime of penetrative convection (33, 34). The high-viscosity downwelling slab penetrates to the CMB, and steep, large-amplitude topography develops on the interface (35). Because the deep layer contains more heat-producing elements, the heat flow across the CMB is lower than it would need to be if the entire mantle had a heat production corresponding to the MORB source region. In this model, the vertical thermal gradient across the thermal boundary layer at the interface between the two domains is larger than that across the thermal boundary layer at the CMB. The interface contains a mix of depleted and subducted material. Plumes arise from high spots on the boundary and carry a mix of depleted material, plus a small amount of entrained material from the less depleted lower layer (36).

Figure 2

Numerical simulation of a dense layer, enriched in heat-producing elements, at the base of the mantle. (A) Temperature (color) is elevated in the dense layer; a thermal boundary layer is visible on the interface between the two materials. Contours of composition (lines) illustrate large-scale topography on the interface; in this model, with B = 0.9, a subducting slab on the right side (blue) nearly reaches the CMB. A small amount of dense material is entrained in the hot plume (white) rising from the peak of the layer. (B) The cold slab retains a high viscosity ν (28), whereas the higher temperature in the lower layer tends to counteract the effect of pressure on viscosity. (C) Density Δρ, relative to the mean density at a depth of 1100 km (above the compositional boundary). On this nondimensional scale, ±0.5 corresponds to ±2 to 4%. (D) Laterally averaged densities for the model in (A) to (C) and for a reference model containing no compositional variations (dashed line) are plotted with respect to their values at a depth of 1100 km. Net density varies by less than 1% compared to the reference model. ΔρTis density change due to temperature alone, ΔρC is density change due to composition alone, and ΔρTotal is the net density for the chemically stratified model.

Because of numerical diffusion, our convection model overestimates the rate at which material in the undepleted layer is entrained by plumes; thus, the rate of entrainment and mixing between layers that we calculate is an upper bound. Even with this small artificial diffusivity, the layers in our calculations remain compositionally distinct for at least 10 overturn times (billions of years). Scaling of mixing rates and dynamic topography with viscosities suggests that mixing was less efficient and the layer more stable in the past (37), so a compositionally distinct layer with the intrinsic density contrast that we propose should be stable over the age of Earth.

To quantify the effects of intrinsic density contrasts on the net properties, we compare the average density as a function of depth for these models with that of a reference model with the same viscosity law but with uniform composition. The reference model with uniform composition exhibits a subadiabatic density increase of ∼0.1ρ0αTΔT that results from the accumulation of cold material in the lower mantle (Fig. 2D, dashed line). For the compositionally stratified model, assuming a volumetric coefficient of thermal expansion for the lower mantle of αT = 2 × 10–5 K–1 and ΔT = 1800°C, the increase in density due to composition alone (green, Fig. 2D) is ∼4%. This density increase is not abrupt but is distributed over the bottom half of the domain in a way that mimics the topography of the boundary between the compositionally distinct layers. There is a nearly compensating decrease in density due to increasing temperature (red, Fig. 2D). For the compositionally layered case (blue, Fig. 2D), the total density excess over that for the uniform composition reference model is never more than ∼1%.

When mineral physics estimates of density and seismic velocity as a function of depth for a particular model composition are compared with estimates from analyses of the normal modes of Earth's free oscillation, the results suggest that the lower mantle is ≤2% more dense than would be expected if it were made of the composition estimated for the upper mantle (38, 39). This range encompasses the net density increase in our models (Fig. 2C). There are, however, a number of complicating factors. First, although the density averaged over depths on the order of 1000 km can be determined from normal mode data to an accuracy of better than 0.3% (40), there is little resolving power to determine the variations over shorter distances. Earth models typically start from the hypothesis that density variations are adiabatic, but variations of at least 0.5% over a depth range of ∼100 km are possible (41). Second, mineral physics estimates of density in the lower mantle require knowledge of both the temperature and the thermal expansion coefficient, neither of which is well determined. Third, our convection calculations are for an incompressible model with constant thermal expansivity, which complicates comparisons with Earth. Within these uncertainties, however, our model appears to be consistent with the constraints from seismology and mineral physics.

If a stable, dense layer exists, the amplitude and pattern ofS- and P-velocity anomalies deep in the mantle would depend on whether the increase in density results from enrichment in iron or silicon. To first order, Fe enrichment (Fig. 3A) does not affect the elastic modulus of perovskite, so the main effect is a wavespeed reduction due to increased density (42). Changing the Si content alters the relative proportion of perovskite and magnesiowüstite, leading to a lower layer that is seismically fast even though it is hot (Fig. 3B). Thus, adding some Si to an Fe-enriched region would reduce the large velocity variations of Fig. 3A to reasonable values. Both components may be important in Earth.

Figure 3

Model seismic velocities (lateral variations in ∂V s/V s) assuming that the intrinsic density increase in the lower layer is due to enrichment in (A) Fe or (B) Si (39,42). Contour intervals are 1%.

We do not have precise constraints on the volume of the depleted mantle. The change in seismic signature characterized by the radial correlation function (7) and an investigation of possible nonadiabatic density jumps (41) suggest that this transition should begin at a depth in excess of 1500 km. The high amplitude of topography on the interface may lead to incoherent scattering and complicate the interpretation of radial correlation functions, which average out signals over a wide lateral region. Because of the dynamic topography, the interface would more likely manifest itself as regional discontinuities of variable depths, locally perhaps as shallow as 1200 km, and the interface might therefore lead to incoherent scattering. Through analysis of scattered waves, seismologists may be able to identify the nature and morphology of this feature. Simulations of the substantially thinner D" layer in two dimensions (43) and in three dimensions (34) show that compositionally dense material piles up under upwellings [possibly explaining the superplume structures seen under Africa, for example (34)] and that some dense material is entrained in upwellings (34, 36). These simulations exhibit varying degrees of topography onD", depending on the density of the layer. Further studies will be required to determine how convection in the relatively thick, dense layer that we propose relates to the D" region.

Because our calculations are for an incompressible medium and a Cartesian geometry, we cannot assign great significance to the exact depth of the interface. The important concept is that the anomalous mantle region should be voluminous enough to generate substantial heat and is not easily swept out of the way by the convection in the overlying layer associated with plate motions. Also, the dynamics of convection seem to favor massive, long-wavelength deflections of the compositional boundary, rather than highly localized penetrative convection (44).

The seismic signature of the downwelling slab is apparent through the upper two-thirds of the mantle. At greater depths, where the slab has encountered the dense layer, the signature of the slab becomes weaker, broader, and more complex. Although it is clear from the temperature and composition that the slab reaches to near the CMB, its seismic signal is diffuse. In Earth, we would expect to see variability along the strike of the slab as well. Models of downwellings in a spherical shell indicate that slab-like structures may be expected to break up into drips in the deep mantle (45). Lateral variations in the age of slabs originating from different tectonic settings will also affect the ability of a slab to deflect the boundary. We therefore envision slab penetration to near the CMB as being a relatively rare and probably short-lived event.

Fossil subducted slabs would collect at the interface between the two materials. A thermal boundary layer at the interface would generate plumes that rise to the surface, carrying components of ancient slab as well as sometimes carrying components of the dense layer itself. Plumes would tend to arise from local high spots on the interface; these high spots act to stabilize the location of plumes. Because the surface of the interface varies with depth, plumes could originate at depths ranging from 1600 km to perhaps the CMB.


View Abstract

Navigate This Article