Mantle plume dynamics at the rear of a retreating slab

In this paper, we present 3-D numerical simulations in which a compositional mantle plume rises through a shallow mantle layer at the rear of a retreating slab. The slab–plume buoyancy flux ratio Bs/Bp is varied from 3 to 30 over nine slab–plume simulations. The plume causes an overall decrease of the slab retreat rate from 15 to 7 % in the corresponding range $B_s/B_p=[3,30]$. The retreat rate decrease occurs in two stages: the first decrease occurs remotely when plume and slab are hundreds of kilometres apart; the second decrease is linked to the slab–plume-head impact. Continuous tracking of key positions along the plume head and conduit together with velocity profiling further shows a very close interplay between the conduit and plume head dynamics. In particular, the combination of the slab– and plume(head)–induced flows at the rear of the slab increases the advection and tilt of the conduit and causes its flaring with height in the direction parallel to the trench. As a result, the conduit source slowly drifts away from the slab by hundreds of kilometres and flares by one and a half times its original radius over 30 Myr of plume head spreading. The conduit tilt and flaring result in an increasingly unbalanced azimuthal distribution of the incoming plume flux from the feeding conduit into the head. These changes in the feeding conditions coupled with the mantle flow at the level of the plume head lead to the asymmetric spreading of plume material beneath the plate in the preferential direction that is parallel to the trench. Upon its arrival against the slab, it can be a front of buoyant material set to subduct along the slab width that has widened up to two and a half times its initial dimension. In nature, it is expected to extend from hundreds to thousands of kilometres depending on the slab–plume buoyancy flux ratio. To our knowledge, this study is the first to highlight the inter-relationships between a plume head and its feeding conditions in the plume source region.

[1]  Louis Moresi,et al.  Underworld2: Python Geodynamics Modelling for Desktop, HPC and Cloud , 2020, J. Open Source Softw..

[2]  D. May,et al.  Benchmark of three-dimensional numerical models of subduction against a laboratory experiment , 2018, Physics of the Earth and Planetary Interiors.

[3]  W. Griffin,et al.  Plume-subduction interaction forms large auriferous provinces , 2017, Nature Communications.

[4]  A. Folguera,et al.  Plume overriding triggers shallow subduction and orogeny in the southern Central Andes , 2017 .

[5]  C. Mériaux,et al.  Mantle plumes in the vicinity of subduction zones , 2016 .

[6]  M. Faccenda,et al.  Upper- and mid-mantle interaction between the Samoan plume and the Tonga–Kermadec slabs , 2016, Nature Communications.

[7]  Lijun Liu,et al.  The role of a mantle plume in the formation of Yellowstone volcanism , 2016 .

[8]  C. Mériaux,et al.  A two‐way interaction between the Hainan plume and the Manila subduction zone , 2015 .

[9]  Changyeol Lee,et al.  Short-term and localized plume-slab interaction explains the genesis of Abukuma adakite in Northeastern Japan , 2014 .

[10]  P. Betts,et al.  The influence of a mantle plume head on the dynamics of a retreating subduction zone , 2012 .

[11]  Taras Gerya,et al.  A comparison of numerical surface topography calculations in geodynamic modelling: an evaluation of the ‘sticky air’ method , 2012 .

[12]  R. C. Kerr,et al.  Rayleigh–Taylor instability of an inclined buoyant viscous cylinder , 2011, Journal of Fluid Mechanics.

[13]  P. Betts,et al.  Mesoproterozoic plume‐modified orogenesis in eastern Precambrian Australia , 2009 .

[14]  Gabriele Morra,et al.  A benchmark comparison of spontaneous subduction models – towards a free surface , 2008 .

[15]  Louis Moresi,et al.  Computational approaches to studying non-linear dynamics of the crust and mantle , 2007 .

[16]  R. C. Kerr,et al.  Structure and dynamics of sheared mantle plumes , 2004 .

[17]  D. Wyman,et al.  Assembly of Archean cratonic mantle lithosphere and crust: plume–arc interaction in the Abitibi–Wawa subduction–accretion complex , 2002 .

[18]  G. Brimhall,et al.  Plume-modified orogeny: An example from the western United States , 1998 .

[19]  William Gropp,et al.  Efficient Management of Parallelism in Object-Oriented Numerical Software Libraries , 1997, SciTools.

[20]  N. Sleep.  Lateral flow of hot plume material ponded at sublithospheric depths , 1996 .

[21]  M. Manga.  Mixing of heterogeneities in the mantle: Effect of viscosity differences , 1996 .

[22]  C. Mériaux,et al.  Simple fluid dynamic models of volcanic rift zones , 1995 .

[23]  U. Christensen,et al.  Three-dimensional modeling of plume-lithosphere interaction , 1994 .

[24]  R. W. Griffiths,et al.  Stirring and structure in mantle starting plumes , 1990 .

[25]  R. C. Kerr,et al.  The propagation of two-dimensional and axisymmetric viscous gravity currents at a fluid interface , 1989, Journal of Fluid Mechanics.

[26]  P. Wilmott,et al.  Convective Mixing In the Mantle: the Role of Viscosity Differences.: Convective mixing in the mantle , 1988 .

[27]  M. Richards,et al.  Deflection of plumes by mantle shear flow: experimental results and a simple theory , 1988 .

[28]  H. Huppert.  The propagation of two-dimensional and axisymmetric viscous gravity currents over a rigid horizontal surface , 1982, Journal of Fluid Mechanics.

[29]  Geoffrey Ingram Taylor,et al.  The formation of emulsions in definable fields of flow , 1934 .

[30]  M. Richards,et al.  The fluid dynamics of plume-ridge and plume-plate interactions: An experimental investigation , 1995 .