Fast asthenosphere motion in high‐resolution global mantle flow models

A variety of geologic observations point to fast upper mantle flow that may exceed plate tectonic velocities by an order of magnitude. At the same time there is mounting evidence from seismology for flow-like structures in the upper 100–200 km of the mantle. Here we present a set of geodynamic simulations to link these observations. In a synthetic setting, we include asthenospheric channels of varying thickness, with an extreme case of 100 km, and a significant viscosity contrast of up to 4 orders of magnitude relative to the deeper mantle. Using our new global high-resolution code TERRA-NEO, we obtain an increase in velocity by a factor of 10 between a 1000 km thick and the very thin channel, translating into velocities of ∼ 20 cm/a within the narrow asthenosphere. We further present and verify a simple Poiseuille flow model, predicting that the upper mantle velocity scales with the inverse of the asthenosphere thickness.

[1]  Martin Kronbichler,et al.  High accuracy mantle convection simulation through modern numerical methods , 2012 .

[2]  FAST TRACK PAPER: A simple model of high Prandtl and high Rayleigh number convection bounded by thin low-viscosity layers , 2006 .

[3]  A. Jones.  OBSERVATIONS OF THE ELECTRICAL ASTHENOSPHERE BENEATH SCANDINAVIA , 1982 .

[4]  S. Stein,et al.  Implications of oceanic intraplate seismicity for plate stresses, driving forces and rheology , 1985 .

[5]  W. Ryan,et al.  Constraints on asthenospheric flow from the depths of oceanic spreading centers: The East Pacific Rise and the Australian‐Antarctic Discordance , 2009 .

[6]  Georg Stadler,et al.  Large-scale adaptive mantle convection simulation , 2013 .

[7]  Mark A. Richards,et al.  A sensitivity study of three-dimensional spherical mantle convection at 108 Rayleigh number: Effects of depth-dependent viscosity, heating mode, and an endothermic phase change , 1997 .

[8]  Adrian Lenardic,et al.  Three‐dimensional mantle convection simulations with a low‐viscosity asthenosphere and the relationship between heat flow and the horizontal length scale of convection , 2008 .

[9]  J. Mitrovica,et al.  Haskell [1935] revisited , 1996 .

[10]  N. White,et al.  Depth, age and dynamic topography of oceanic lithosphere beneath heavily sedimented Atlantic margins , 2009 .

[11]  Richard G. Gordon,et al.  Geologically current plate motions , 2010 .

[12]  Bramley J. Murton,et al.  A continuous 55-million-year record of transient mantle plume activity beneath Iceland , 2014 .

[13]  Ulrich Rüde,et al.  A Massively Parallel Multigrid Method for Finite Elements , 2006, Computing in Science & Engineering.

[14]  P. Cobbold,et al.  Episodic burial and exhumation in NE Brazil after opening of the South Atlantic , 2012 .

[15]  D. L. Anderson,et al.  Partial melting in the upper mantle , 1970 .

[16]  H. Bunge,et al.  Rapid plate{motion variations through geological time: observations serving geodynamic interpretation , 2015 .

[17]  S. Karato,et al.  Water, partial melting and the origin of the seismic low velocity and high attenuation zone in the upper mantle , 1998 .

[18]  B. Schuberth,et al.  Stability of the rotation axis in high‐resolution mantle circulation models: Weak polar wander despite strong core heating , 2009 .

[19]  N. A. Haskell.  The Motion of a Viscous Fluid Under a Surface Load , 1935 .

[20]  Sri Widiyantoro,et al.  Global seismic tomography: A snapshot of convection in the Earth: GSA Today , 1997 .

[21]  Andreas Fichtner,et al.  Full waveform tomography of the upper mantle in the South Atlantic region: Imaging a westward fluxing shallow asthenosphere? , 2013 .

[22]  M. Richards,et al.  Depth-dependent viscosity and mantle stress amplification: implications for the role of the asthenosphere in maintaining plate tectonics , 2012 .

[23]  V. Solomatov,et al.  Scaling of temperature‐ and stress‐dependent viscosity convection , 1995 .

[24]  B. Hager,et al.  Long-wavelength variations in Earth’s geoid: physical models and dynamical implications , 1989, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

[25]  Ulrich Rüde,et al.  Peta-Scale Hierarchical Hybrid Multigrid Using Hybrid Parallelization , 2013, LSSC.

[26]  Lars Stixrude,et al.  Thermodynamics of mantle minerals – I. Physical properties , 2005 .

[27]  A. Lenardic,et al.  Long wavelength convection, Poiseuille–Couette flow in the low-viscosity asthenosphere and the strength of plate margins , 2010 .

[28]  Nicholas J. White,et al.  Transient convective uplift of an ancient buried landscape , 2010 .

[29]  Irina M. Artemieva,et al.  Global 1°×1° thermal model TC1 for the continental lithosphere: Implications for lithosphere secular evolution , 2006 .

[30]  Georg Stadler,et al.  Parallel scalable adjoint-based adaptive solution of variable-viscosity Stokes flow problems , 2009 .

[31]  Patrick Wu,et al.  Rheology of the Upper Mantle: A Synthesis , 1993, Science.

[32]  B. Kennett,et al.  Global azimuthal seismic anisotropy and the unique plate-motion deformation of Australia , 2005, Nature.

[33]  A. Jellinek,et al.  Viscous coupling at the lithosphere‐asthenosphere boundary , 2011 .

[34]  Barbara I. Wohlmuth,et al.  Performance and Scalability of Hierarchical Hybrid Multigrid Solvers for Stokes Systems , 2015, SIAM J. Sci. Comput..

[35]  John H. Woodhouse,et al.  S40RTS: A degree-40 shear-velocity model for the mantle from new Rayleigh wave dispersion, teleseismic traveltime and normal-mode splitting function measurements , 2011 .

[36]  Paul J. Tackley,et al.  Effects of strongly variable viscosity on three‐dimensional compressible convection in planetary mantles , 1996 .

[37]  H. Bunge,et al.  Rapid South Atlantic spreading changes and coeval vertical motion in surrounding continents: Evidence for temporal changes of pressure‐driven upper mantle flow , 2014 .

[38]  P. Japsen,et al.  Elevated, passive continental margins: Not rift shoulders, but expressions of episodic, post-rift burial and exhumation , 2012 .

[39]  Barbara Romanowicz,et al.  Waveform Tomography Reveals Channeled Flow at the Base of the Oceanic Asthenosphere , 2013, Science.

[40]  Andreas Fichtner,et al.  Full seismic waveform tomography for upper-mantle structure in the Australasian region using adjoint methods , 2009 .

[41]  A. Paulson,et al.  On the resolution of radial viscosity structure in modelling long‐wavelength postglacial rebound data , 2009 .

[42]  P. R. Vogt,et al.  Asthenosphere motion recorded by the ocean floor south of Iceland , 1971 .

[43]  R. Müller,et al.  Geological, tomographic, kinematic and geodynamic constraints on the dynamics of sinking slabs , 2014 .

[44]  C. G. Chase.  Asthenospheric counterflow: a kinematic model , 1979 .