A benchmark study on mantle convection in a 3‐D spherical shell using CitcomS

As high-performance computing facilities and sophisticated modeling software become available, modeling mantle convection in a three-dimensional (3-D) spherical shell geometry with realistic physical parameters and processes becomes increasingly feasible. However, there is still a lack of comprehensive benchmark studies for 3-D spherical mantle convection. Here we present benchmark and test calculations using a finite element code CitcomS for 3-D spherical convection. Two classes of model calculations are presented: the Stokes' flow and thermal and thermochemical convection. For Stokes' flow, response functions of characteristic flow velocity, topography, and geoid at the surface and core-mantle boundary (CMB) at different spherical harmonic degrees are computed using CitcomS and are compared with those from analytic solutions using a propagator matrix method. For thermal and thermochemical convection, 24 cases are computed with different model parameters including Rayleigh number (7 × 10^3 or 10^5) and viscosity contrast due to temperature dependence (1 to 10^7). For each case, time-averaged quantities at the steady state are computed, including surface and CMB Nussult numbers, RMS velocity, averaged temperature, and maximum and minimum flow velocity, and temperature at the midmantle depth and their standard deviations. For thermochemical convection cases, in addition to outputs for thermal convection, we also quantified entrainment of an initially dense component of the convection and the relative errors in conserving its volume. For nine thermal convection cases that have small viscosity variations and where previously published results were available, we find that the CitcomS results are mostly consistent with these previously published with less than 1% relative differences in globally averaged quantities including Nussult numbers and RMS velocities. For other 15 cases with either strongly temperature-dependent viscosity or thermochemical convection, no previous calculations are available for comparison, but these 15 test calculations from CitcomS are useful for future code developments and comparisons. We also presented results for parallel efficiency for CitcomS, showing that the code achieves 57% efficiency with 3072 cores on Texas Advanced Computing Center's parallel supercomputer Ranger.

[1]  N. Weiss,et al.  Convection in the earth's mantle: towards a numerical simulation , 1974, Journal of Fluid Mechanics.

[2]  M. Gurnis,et al.  Constraints on the lateral strength of slabs from three-dimensional dynamic flow models , 1996 .

[3]  Bradford H. Hager,et al.  The relationship between plate velocity and trench viscosity in Newtonian and power‐law subduction calculations , 1990 .

[4]  Paul J. Tackley,et al.  Effects of strongly temperature‐dependent viscosity on time‐dependent, three‐dimensional models of mantle convection , 1993 .

[5]  Louis Moresi,et al.  Role of temperature‐dependent viscosity and surface plates in spherical shell models of mantle convection , 2000 .

[6]  G. Schubert,et al.  Three-dimensional thermal convection in a spherical shell , 1989, Journal of Fluid Mechanics.

[7]  T. Becker On the effect of temperature and strain-rate dependent viscosity on global mantle flow, net rotation, and plate-driving forces , 2006 .

[8]  Shijie Zhong,et al.  Thermochemical structures within a spherical mantle: Superplumes or piles? , 2004 .

[9]  R. A. Wentzell,et al.  Hydrodynamic and Hydromagnetic Stability. By S. CHANDRASEKHAR. Clarendon Press: Oxford University Press, 1961. 652 pp. £5. 5s. , 1962, Journal of Fluid Mechanics.

[10]  James H. Roberts,et al.  Supercontinent cycles, true polar wander, and very long-wavelength mantle convection , 2007 .

[11]  A. Brooks,et al.  A Petrov-Galerkin Finite Element Formulation for Convection Dominated Flows , 1981 .

[12]  Ulrich R. Christensen,et al.  Convection with pressure- and temperature-dependent non-Newtonian rheology , 1984 .

[13]  Louis Moresi,et al.  Mantle convection with a brittle lithosphere: thoughts on the global tectonic styles of the Earth and Venus , 1998 .

[14]  C. Conrad,et al.  Influence of continental roots and asthenosphere on plate‐mantle coupling , 2006 .

[15]  S. Zhong 7.05 – Numerical Methods for Mantle Convection , 2007 .

[16]  W. Peltier,et al.  The Kinematics and Dynamics of Poloidal–Toroidal Coupling in Mantle Flow: The Importance of Surface Plates and Lateral Viscosity Variations , 1994 .

[17]  James H. Roberts,et al.  Degree-1 convection in the Martian mantle and the origin of the hemispheric dichotomy , 2006 .

[18]  Shijie Zhong,et al.  Mixing in a 3D spherical model of present-day mantle convection , 1999 .

[19]  Ulrich Hansen,et al.  A new method to simulate convection with strongly temperature- and pressure-dependent viscosity in a spherical shell: Applications to the Earth's mantle , 2006 .

[20]  Shijie Zhong,et al.  Effects of plate and slab viscosities on the geoid , 1999 .

[21]  Remo Guidieri Res , 1995, RES: Anthropology and Aesthetics.

[22]  J. T. Ratcliff,et al.  Steady tetrahedral and cubic patterns of spherical shell convection with temperature‐dependent viscosity , 1996 .

[23]  J. T. Ratcliff,et al.  Effects of temperature-dependent viscosity of thermal convection in a spherical shell , 1996 .

[24]  Ulrich R. Christensen,et al.  Some effects of lateral viscosity variations on geoid and surface velocities induced by density anomalies in the mantle , 1993 .

[25]  B. Hager,et al.  Toroidal-Poloidal Partitioning of Lithospheric Plate Motions , 1991 .

[26]  David J. Stevenson,et al.  Effects of an endothermic phase transition at 670 km depth in a spherical model of convection in the Earth's mantle , 1993, Nature.

[27]  D. Yuen,et al.  The interaction of a subducting lithospheric slab with a chemical or phase boundary , 1984 .

[28]  D. Silvester,et al.  ITERATIVE METHODS FOR STABILIZED MIXED VELOCITY PRESSURE FINITE-ELEMENTS , 1992 .

[29]  David A. Yuen,et al.  The influences of lower mantle viscosity stratification on 3D spherical-shell mantle convection , 1995 .

[30]  S. Zhong Role of ocean-continent contrast and continental keels on plate motion, net rotation of lithosphere, and the geoid , 2001 .

[31]  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.

[32]  T. Becker,et al.  Azimuthal seismic anisotropy constrains net rotation of the lithosphere , 2008 .

[33]  M. Richards,et al.  Role of a low‐viscosity zone in stabilizing plate tectonics: Implications for comparative terrestrial planetology , 2001 .

[34]  Paul J. Tackley,et al.  Self‐consistent generation of tectonic plates in time‐dependent, three‐dimensional mantle convection simulations , 2000 .

[35]  Shijie Zhong,et al.  Interaction of weak faults and non-newtonian rheology produces plate tectonics in a 3D model of mantle flow , 1996, Nature.

[36]  Akira Kageyama,et al.  Application of the Yin-Yang grid to a thermal convection of a Boussinesq fluid with infinite Prandtl number in a three-dimensional spherical shell , 2004 .

[37]  Carl W. Gable,et al.  3D convection at infinite Prandtl number in Cartesian geometry — a benchmark comparison , 1994 .

[38]  D. Yuen,et al.  Numerical Methods for Mantle Convection , 2007 .

[39]  C. Dumoulin,et al.  ŒDIPUS: a new tool to study the dynamics of planetary interiors , 2007 .

[40]  S. Honda,et al.  Generation of plumes under a localized high viscosity lid in 3‐D spherical shell convection , 1999 .

[41]  Manfred Koch,et al.  A benchmark comparison for mantle convection codes , 1989 .

[42]  G. Glatzmaier Numerical simulations of mantle convection: time-dependent, three-dimensional, compressible, spherical shell , 1988 .

[43]  D. Helmberger,et al.  Dynamics of a phase change at the base of the mantle consistent with seismological observations , 1999 .

[44]  Ulrich Hansen,et al.  A finite-volume solution method for thermal convection and dynamo problems in spherical shells , 2005 .

[45]  Gregory M. Hulbert,et al.  Accurate determination of surface normal stress in viscous flow from a consistent boundary flux method , 1993 .

[46]  David Bercovici,et al.  The generation of plate tectonics from mantle convection , 2003 .

[47]  C. Sotin,et al.  Turbulent 3‐D thermal convection in an infinite Prandtl number, volumetrically heated fluid: implications for mantle dynamics , 1994 .

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

[49]  A. Wathen,et al.  Iterative solution techniques for the stokes and Navier‐Stokes equations , 1994 .

[50]  Paul J. Tackley,et al.  Testing the tracer ratio method for modeling active compositional fields in mantle convection simulations , 2003 .

[51]  Marie-Pierre Doin,et al.  A comparison of methods for the modeling of thermochemical convection , 1997 .

[52]  Ulrich R. Christensen,et al.  A one-plume model of martian mantle convection , 1996, Nature.

[53]  Louis Moresi,et al.  Numerical investigation of 2D convection with extremely large viscosity variations , 1995 .

[54]  G. Hirth,et al.  Rheologic controls on slab dynamics , 2007 .

[55]  R. Young Finite-amplitude thermal convection in a spherical shell , 1974, Journal of Fluid Mechanics.

[56]  J. Tromp,et al.  Theoretical Global Seismology , 1998 .

[57]  G. Schubert,et al.  A benchmark comparison of numerical methods for infinite Prandtl number thermal convection in two-dimensional Cartesian geometry , 1990 .

[58]  Mark A. Richards,et al.  Effect of depth-dependent viscosity on the planform of mantle convection , 1996, Nature.

[59]  B. Hager,et al.  Entrainment of a dense layer by thermal plumes , 2003 .

[60]  M. Gurnis A reassessment of the heat transport by variable viscosity convection with plates and lids , 1989 .

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

[62]  J. Davies,et al.  Numerical investigation of layered convection in a three‐dimensional shell with application to planetary mantles , 2004 .

[63]  J. Hunen,et al.  New evidence for dislocation creep from 3-D geodynamic modeling of the Pacific upper mantle structure , 2005 .

[64]  Louis Moresi,et al.  Role of faults, nonlinear rheology, and viscosity structure in generating plates from instantaneous mantle flow models , 1998 .

[65]  Shijie Zhong,et al.  Dynamic feedback between a continentlike raft and thermal convection , 1993 .

[66]  Bradford H. Hager,et al.  A simple global model of plate dynamics and mantle convection , 1981 .

[67]  Louis Moresi,et al.  The accuracy of finite element solutions of Stokes's flow with strongly varying viscosity , 1996 .

[68]  B. Parsons,et al.  Convective instabilities in a variable viscosity fluid cooled from above , 1985 .

[69]  John R. Baumgardner,et al.  Three-dimensional treatment of convective flow in the earth's mantle , 1985 .

[70]  J. Z. Zhu,et al.  The finite element method , 1977 .