Dynamics and internal structure of a lower mantle plume conduit
Introduction
Geochemical studies of Hawaiian basalts have yielded a wealth of data relevant to the understanding of the chemical and isotopic internal structure of mantle plumes. However, the interpretation of these data has remained highly controversial, partly because there are significant uncertainties as to how the volcanic plumbing system samples the underlying plume, but also because there is little agreement as to how the existing geochemical heterogeneities might be spatially arranged within the plume conduit. The pioneering model by Hauri et al. (1994), suggested a concentrically zoned plume conduit, whereby geochemical variations occur only in the radial direction, due to entrainment of surrounding mantle. Entrainment is governed by thermal diffusion and is expected to decrease with increasing plume flux (Hauri et al., 1994). In spite of the high buoyancy flux of the Hawaiian plume (Sleep, 1990) the concentric model has been extensively invoked to explain the time and space variability of isotope ratios e.g., (DePaolo and Stolper, 1996, Hauri et al., 1996, Lassiter et al., 1996). However, recent data from the Hawaiian Scientific Drilling Project HSDR-2, have shown spatial and temporal variations for Pb, Nd, Hf, and Sr isotope ratios (see DePaolo et al., 2001, Blichert-Toft et al., 2003, Eisele et al., 2003, Abouchami et al., 2005, Bryce et al., 2005) and references therein) that are inconsistent with concentric zoning and have prompted the suggestion of new conceptual models.
Presenting a contrasting view of the internal plume structure, Blichert-Toft et al. (2003) propose a plume structure that minimizes radial variations across the conduit by assuming a ‘plug-flow’, whereby the vertical velocity remains nearly constant across the main part of the conduit and decreases sharply only at its periphery. The ‘plug-flow’ velocity profile contrasts with the exponentially varying velocity profile derived analytically by Olson et al. (1993). According to Blichert-Toft et al. (2003) the geochemical cross section of the plume does not need to be time invariant, so that the conduit structure may be similar to a stack of nearly horizontal thin ‘layers’ with distinct isotope ratios.
In yet another interpretation (Eisele et al., 2003, Abouchami et al., 2005) the conduit internal structure is similar to a bundle of vertically elongated filaments, successively sampled by different volcanoes as the Pacific plate moves over the Hawaiian plume. The underlying physical process is stretching of deep heterogeneities due to velocity gradients within the plume conduit.
The above interpretations are purely based on geochemical observations and there is a clear need to investigate the internal structure of a plume conduit from a fluid dynamics perspective. In this paper we use numerical simulations of a thermal plume to address some fluid dynamics questions expected to be relevant for the geochemical interpretation of internal plume structure. First, what is the relation between length-scales of heterogeneities across the source region (assumed to be the D″ layer overlying the core–mantle boundary) and length- and time-scales of geochemical variations once heterogeneities are upwelling in the plume conduit? Second, deep heterogeneities may be thought of as isolated finite-size volumes, or they may be distributed in space so as to form an internal architecture in the D″ region. In this case, what would be the internal structure of a plume conduit fed by such a complex source region? Third, are there fundamental differences between the central and the peripheral zones of the plume conduit, for example due to different deformations undergone during upwelling and/or different parts of the source region being sampled? Finally, which aspects of the three different (and to some extent mutually exclusive) models of conduit structure presented above are consistent with fluid dynamics constraints?
To our knowledge, previous fluid dynamics studies have not explored the relation between finite-size heterogeneities in the thermal boundary layer (TBL) and the conduit internal structure. Although it is beyond the scope of our paper to review the important contributions of laboratory experiments and numerical models for the understanding of mantle plumes, we note that early laboratory experiments, for example by Griffiths and Campbell (1990), focused on the plume head rather than on the long lived conduit. Moreover, laboratory experiments have rigid bottom boundary conditions and the hot plume is often injected in the tank, thereby precluding any investigation of the dynamics and internal deformations in the basal TBL. More recent experiments by Kerr and Mériaux (2004) on plumes upwelling in a mantle wind generated by surface plate motion, concluded that the Hawaiian plume conduit should be zoned azimuthally, rather than concentrically. However, Kerr and Mériaux (2004) did not consider the deformations of finite-size heterogeneities, nor did they quantify strain rates across the conduit and the TBL.
Our work is based on a purely thermal flow within an idealized TBL and the associated plume conduit. Certainly, the D″ region is more complex than a thermal boundary layer (see for example Jellinek and Manga, 2004, Garnero and McNamara, 2008, and references therein), and our model does not include most of the recent findings, namely the post-perovskite phase transition (Murakami et al., 2004) and the existence of chemical heterogeneities in the deep mantle. While the role of the post-perovskite on mantle plumes is still a matter of debate, it is now clear that chemically denser material plays an important role on plumes shape and dynamics (Tackley, 1998, Davaille, 1999, Farnetani and Samuel, 2005), composition (Christensen and Hofmann, 1994, Samuel and Farnetani, 2003) and excess temperature (Farnetani, 1997). We are aware that considering the D″ region as a purely thermal boundary layer is a gross simplification; nevertheless, from a modeling stand point, we consider it important and instructive to start with the most basic fluid flow.
It is also difficult to characterize the geochemical fingerprint of the D″ region. The base of the Earth's mantle is likely to be a slab ‘graveyard’. In such a case, denser subducted crust may segregate and evolve geochemically (Christensen and Hofmann, 1994) before being recycled in mantle plumes (Hofmann and White, 1982). But the D″ region, or parts of it, could also host a ‘hidden reservoir’, with a distinct geochemical composition, remnant from the early differentiation of the Earth's mantle (Boyet and Carlson, 2005, Tolstikhin and Hofmann, 2005). Deliberately, we do not attempt to reproduce this complex geochemistry and geochemical evolution over time. Instead, we conduct a ‘thought experiment’ to investigate the relation between a hypothetically heterogeneous D″ and the internal structure of a plume conduit.
Finally, we had to restrict our modeling to a few examples of kilometer scale geochemical heterogeneities with simple shapes and basic spatial orientations, which are certainly poor representations of the great variety created by mantle stirring. As reviewed by Hofmann (2003) the spatial scale of heterogeneities may vary from 102–104 km, to the decimeter scale (Allègre and Turcotte, 1986) and less. On the other hand, the kilometer-scale explored here is likely to be particularly relevant to mantle plumes, particularly in view of the rapid fluctuations in isotopic composition observed in the stratigraphic record of the Hawaiian Scientific Drilling Project (Eisele et al., 2003).
Our high resolution numerical model allows us to calculate deformations in a plume conduit with a detail that is not achieved in global mantle mixing models. Although such models have improved our understanding of convective stirring (see for example Kellogg, 1992, Davies, 2002, van Keken et al., 2003, Tackley, 2007, and references therein), it is important to focus on deformations inside hot, low viscosity and fast upwelling narrow plumes that are likely to represent unique places in the lower mantle where strain rates are particularly high and the associated stretching is more rapid and intense than elsewhere in the lower mantle.
This paper focuses on the basal thermal boundary layer feeding the plume and on the lower mantle part of the conduit, while, in a companion paper (in prep.) we model the dynamics of an upper mantle conduit sheared by a fast moving plate, allowing us to investigate how volcanoes sample a spatially heterogeneous plume conduit.
Lastly, we point out that although our results are used to elucidate the internal structure of the Hawaiian conduit, they may be of interest to other plumes, for example Galápagos, where isotope and incompatible element ratios have a complex spatial distribution, namely a horseshoe pattern with enriched signatures around a depleted center (White et al., 1993, Harpp and White, 2001). Interestingly, the observed asymmetrical spatial zonation of the Galápagos hotspot has persisted for several millions of years (Hoernle et al., 2000), suggesting the involvement of large scale heterogeneities in the source region feeding the plume.
Section snippets
Numerical model
We use the two dimensional code ConmanCYL (Farnetani, 1997, King et al., 1990) in cylindrical geometry to solve the equations of conservation of mass, conservation of momentum and conservation of energy for an incompressible viscous fluid at infinite Prandtl number. The size of the model domain is 1500 × 2890 km in the radial and vertical directions, respectively. The element size is 4 × 4 km in the regions of interest (i.e., the TBL and the plume conduit), elsewhere in the lower mantle the element
Results
Starting from the initial condition described above we run the calculation until the transient plume head reaches the base of the lithosphere and the velocity field approaches steady state. Fig. 1a (solid line) shows the vertical velocity component vz as a function of radial distance r across the plume conduit at 2000 km depth: at the axis vaxisz = 50 cm/yr, at r = 38 km vz = vaxisz/2 and at r = 70 km vz = vaxisz / 10. The radial dependence of vz is quite invariant throughout the height of the plume and it
Comparison with previous models of conduit structure and conclusion
The first conceptual model of conduit structure was suggested nearly fifteen years ago and is known as the ‘concentric zoning’ (Hauri et al., 1994, DePaolo and Stolper, 1996) whereby the plume center samples the deepest part of the TBL, while the conduit periphery either samples the upper part of the TBL (Farnetani et al., 2002, Bryce et al., 2005) or, in case of entrainment, it samples also the surrounding mantle (Hauri et al., 1996). Our simulations confirm one basic idea of the concentric
Acknowledgements
We thank Bill White, Ross Kerr and an anonymous reviewer for their constructive and useful reviews. Erik Hauri provided thorough comments on an earlier version of the manuscript. CGF thanks Bernard Legras for discussions of the Supplementary Material and Claude Jaupart for support. We also wish to thank Rick Carlson for his editorial handling. IPGP Contribution No. 2489; LDEO Contribution No. 7247.
References (42)
Stirring geochemistry in mantle convection models with stiff plates and slabs
Geochim. Cosmochim. Acta
(2002)- et al.
Mixing and deformations in mantle plumes
Earth Planet. Sci. Lett.
(2002) - et al.
Stirring and structure in mantle starting plumes
Earth Planet. Sci. Lett.
(1990) - et al.
Mantle plumes from ancient oceanic crust
Earth Planet. Sci. Lett.
(1982) - et al.
ConMan: vectorizing a finite element code for incompressible two-dimensional convection in the Earth's mantle
Phys. Earth Planet. Inter.
(1990) - et al.
Ultrafast mantle plumes and implications for flood basalt volcanism in the Northern Atlantic Region
Tectonophysics
(1999) - et al.
The dynamical and thermal structure of deep mantle plumes
Phys. Earth Planet. Inter.
(1983) - et al.
Thermochemical convection and helium concentrations in mantle plumes
Earth Planet. Sci. Lett.
(2003) Mantle Geochemical Geodynamics
- et al.
Hawaiian double volcanic chain triggered by an episodic involvement of recycled material: constraints from temporal Sr–Nd–Hf–Pb isotopic trend of the Loa-type volcanoes
Earth Planet. Sci. Lett.
(2008)
Early crust on top of the Earth's core
Phys. Earth Planet. Inter.
Pb isotopes reveal bilateral asymmetry and vertical continuity in the Hawaiian plume
Nature
Algorithm 760: rectangular-grid-data surface fitting that has the accuracy of a bicubic polynomial
ACM Trans. Math. Softw.
Implications of a two component marble-cake mantle
Nature
Hawaiian hot spot dynamics as inferred from the Hf and Pb isotope evolution of Mauna Kea volcano
Geochem. Geophys. Geosyst.
142Nd evidence for early (> 4.53 Ga) global differentiation of the silicate Earth
Science
Geochemical structure of the Hawaiian plume: Sr, Nd, and Os isotopes in the 2.8 km HSDP-2 section of Mauna Kea volcano
Geochem. Geophys. Geosyst.
Segregation of subducted oceanic crust in the convecting mantle
J. Geophys. Res.
Simultaneous generation of hotspots and superswells by convection in a heterogeneous planetary mantle
Nature
Models of Hawaiian volcano growth and plume structure: implications of results from the Hawaii Scientific Drilling Project
J. Geophys. Res.
Isotopic evolution of Mauna Loa and the chemical structure of the Hawaiian plume
Geochem. Geophys. Geosyst.
Cited by (77)
Mid-mantle seismic scatterers beneath the Samoan hotspot
2023, Physics of the Earth and Planetary InteriorsDynamics of a chemically pulsing mantle plume
2020, Earth and Planetary Science LettersCitation Excerpt :An effect of a vertically elongate interface would be to prolong the transition to an incoming heterogeneity, or in the case of a fixed transition time (0.3 Myr), it would require a plume ascent rate proportionally faster than 167 km Myr−1. If the velocity at 0.7r from the axis is ∼10 times slower than at the centreline (Farnetani and Hofmann, 2009) a rise from 700 km depth would result in >500 km of relief on the boundary, requiring either high-velocity ascent or a protracted transition period. An exponentially declining radial velocity profile is likely to result in different strain rates across the plume stem.
Intraplate Volcanism
2020, Encyclopedia of Geology: Volume 1-6, Second EditionOsmium isotope evidence for a heterogeneous <sup>3</sup>He/<sup>4</sup>He mantle plume beneath the Juan Fernandez Islands
2019, Geochimica et Cosmochimica ActaCitation Excerpt :Weis et al. (2011) proposed that the geochemical differences between the Kea and Loa trends reflect preferential sampling of two distinct sources of deep mantle material: the Loa trend sampling LLSVP (Large Low-Shear Velocity Province) material, and the Kea trend sampling lower mantle material. Alternatively, Farnetani & Hofmann (2009, 2010) have proposed that a continuous deep layer overlained by a bilaterally varying layer best reproduces the “bilateral asymmetry” observed for Pb isotopic differences between the Kea and Loa-trends of Hawaiian volcanoes and for the concentric distribution of 3He/4He (DePaolo et al., 2001). It has also been suggested that high-3He/4He materials could be distributed as isolated primitive blobs in the lower mantle (Becker et al., 1999), or in an isolated primitive reservoir in the lower mantle (Allègre et al., 1983; Tolstikhin and Hofmann, 2005) or even in the core (Bouhifd et al., 2013), to explain the isotopic variability observed within individual hotspot tracks.
A review of large low shear velocity provinces and ultra low velocity zones
2019, Tectonophysics