Numerical simulations of mantle convection: Time-dependent,three-dimensional,compressible, spherical shell

Abstract

We describe nonlinear time-dependent numerical simulations of whole mantle convection for a Newtonian, infinite Prandtl number, anelastic fluid in a three-dimensional spherical shell for conditions that approximate the Earth's mantle. Each dependent variable is expanded in a series of 4,096 spherical harmonics to resolve its horizontal structure and in 61 Chebyshev polynomials to resolve its radial structure. A semiimplicit time-integration scheme is used with a spectral transform method. In grid space there are 61 unequally-spaced Chebyshev radial levels, 96 Legendre colatitudinal levels, and 192 Fourier longitudinal levels. For this preliminary study we consider four scenarios, all having the same radially-dependent reference state and no internal heating. They differ by their radially-dependent linear viscous and thermal diffusivities and by the specified temperatures on their isothermal, impermeable, stress-free boundaries. We have found that the structure of convection changes dramatically as the Rayleigh number increases from 10⁵ to 10⁶ to 10⁷. The differences also depend on how the Rayleigh number is increased. That is, increasing the superadiabatic temperature drop, δT, across the mantle produces a greater effect than decreasing the diffusivities. The simulation with a Rayleigh number of 10⁷ is approximately 10,000 times critical, close to estimates of that for the Earth's mantle. However, although the velocity structure for this highest Rayleigh number scenario may be adequately resolved, its thermodynamic structure requires greater horizontal resolution. The velocity and thermodynamic structures of the scenarios at Rayleigh numbers of 10⁵ and 10⁶ appear to be adequately resolved. The 10⁵ Rayleigh number solution has a small number of broad regions of warm upflow embedded in a network of narrow cold downflow regions; whereas, the higher Rayleigh number solutions (with large δT) have a large number of small hot upflow plumes embedded in a broad weak background of downflow. In addition, as would be expected, these higher Rayleigh number solutions have thinner thermal boundary layers and larger convective velocities, temperatures perturbations, and heat fluxes. These differences emphasize the importance of developing even more realistic models at realistic Rayleigh numbers if one wishes to investigate by numerical simulation the type of convection that occurs in the Earth's mantle.