Lagrangian structures and stirring in the Earth’s mantle

Abstract In this paper we investigate three Lagrangian methods that have been recently proposed to quantify mixing in chaotic and time-aperiodic geophysical flows. The analytical method proposed by Haller [e.g., G. Haller, Chaos 10 (2000) 99–108] has a strong mathematical foundation and seeks to determine the location of stable (i.e., most attracting) and unstable (i.e., most repelling) material lines. The main hyperbolic lines describe the spatial organization of chaotic mixing. The finite-size Lyapunov exponents estimate the local mixing properties of the flow using the finite dispersion of particles, while the finite-time Lyapunov exponents have been often used to locate dynamically distinguished regions in geophysical flows. Mantle stirring is induced by the repeated action of stretching and folding, thus requiring to follow the strain history of a fluid element along a trajectory. We calculate the trajectories of more than half a million passive tracers forward and backward in time. The Eulerian velocity field is computed using a finite element code for solid state convection. We focus on a thermochemical model with a chemically denser layer at the base of the Earth’s mantle. This case allows us to test the ability of the Lagrangian techniques to detect the location of a dynamical barrier that inhibits mass exchanges and delimits domains characterized by different efficiency of stirring. We find that the Lagrangian techniques provide a satisfactory description of the main structures governing stirring, and enlighten different and complementary aspects: the methods based on the Lyapunov exponents provide a clear picture of mantle domains characterized by different strength of stirring, while the method proposed by Haller identifies the skeleton of the main structures around which stirring is organized. Our paper builds toward a more rigorous analysis of the stirring processes in the Earth’s mantle, which is required to understand the existence of geochemical reservoirs under a dynamical prospective.

[1]  A. Ōkubo.  Horizontal dispersion of floatable particles in the vicinity of velocity singularities such as convergences , 1970 .

[2]  Ulrich R. Christensen,et al.  Mixing by time-dependent convection , 1989 .

[3]  Comment on "Finding finite-time invariant manifolds in two-dimensional velocity fields" [Chaos 10, 99 (2000)]. , 2001, Chaos.

[4]  Sylvaine Ferrachat,et al.  Regular vs. chaotic mantle mixing , 1998 .

[5]  George Haller,et al.  Geometry of Cross-Stream Mixing in a Double-Gyre Ocean Model , 1999 .

[6]  S. Hart,et al.  Mantle Plumes and Entrainment: Isotopic Evidence , 1992, Science.

[7]  Bernard Legras,et al.  Hyperbolic lines and the stratospheric polar vortex. , 2002, Chaos.

[8]  B. Legras,et al.  Relation between kinematic boundaries, stirring, and barriers for the Antarctic polar vortex , 2002 .

[9]  R. Boehler,et al.  Thermal expansivity in the lower mantle , 1992 .

[10]  J. M. Ottino,et al.  Morphological structures produced by mixing in chaotic flows , 1988, Nature.

[11]  N. Hoffman,et al.  The destruction of geochemical heterogeneities by differential fluid motions during mantle convection , 1985 .

[12]  B. Legras,et al.  Mixing and deformations in mantle plumes , 2002 .

[13]  G. Haller.  Distinguished material surfaces and coherent structures in three-dimensional fluid flows , 2001 .

[14]  M. Gurnis,et al.  Mixing in numerical models of mantle convection incorporating plate kinematics , 1986 .

[15]  M. Gurnis,et al.  The effect of depth‐dependent viscosity on convective mixing in the mantle and the possible survival of primitive mantle , 1986 .

[16]  J. Gollub,et al.  Experimental measurements of stretching fields in fluid mixing. , 2002, Physical review letters.

[17]  Miss A.O. Penney.  (b) , 1974, The New Yale Book of Quotations.

[18]  G. Haller,et al.  Lagrangian coherent structures and mixing in two-dimensional turbulence , 2000 .

[19]  Guido Boffetta,et al.  Detecting barriers to transport: a review of different techniques , 2001, nlin/0102022.

[20]  George Haller,et al.  Lagrangian structures and the rate of strain in a partition of two-dimensional turbulence , 2001 .

[21]  A. Hofmann,et al.  Segregation of subducted oceanic crust in the convecting mantle , 1994 .

[22]  David A. Yuen,et al.  Comparison of mixing properties in convection with the Particle‐Line Method , 1998 .

[23]  J. Schmalzl,et al.  Mixing in vigorous, time-dependent three-dimensional convection and application to Earth's mantle , 1996 .

[24]  J. Ottino.  The Kinematics of Mixing: Stretching, Chaos, and Transport , 1989 .

[25]  A. Mcnamara,et al.  Cooling of the Earth: A parameterized convection study of whole versus layered models , 2000 .

[26]  Bradford H. Hager,et al.  Conman: vectorizing a finite element code for incompressible two-dimensional convection in the Earth's mantle , 1990 .

[27]  H. Samuel,et al.  Thermochemical convection and helium concentrations in mantle plumes , 2003 .

[28]  Raymond T. Pierrehumbert,et al.  Large-scale horizontal mixing in planetary atmospheres , 1991 .

[29]  Louise H. Kellogg,et al.  Mixing in the Mantle , 1992 .

[30]  G. Batchelor,et al.  The effect of homogeneous turbulence on material lines and surfaces , 1952, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[31]  A. Provost,et al.  Topology in isotopic multispace and origin of mantle chemical heterogeneities , 1987 .

[32]  HIROSHI AKIMA,et al.  Algorithm 760: Rectangular-grid-data surface fitting that has the accuracy of a bicubic polynomial , 1996, TOMS.

[33]  G. Haller.  Finding finite-time invariant manifolds in two-dimensional velocity fields. , 2000, Chaos.