Three-dimensional convection of an infinite-Prandtl-number compressible fluid in a basally heated spherical shell

Abstract
A numerical investigation is made of the effects of compressibility on three-dimensional thermal convection in a basally heated, highly viscous fluid spherical shell with an inner to outer radius ratio of approximately 0.55, characteristic of the Earth's whole mantle. Compressibility is implemented with the anelastic approximation and a hydrostatic adiabatic reference state whose bulk modulus is a linear function of pressure. The compressibilities studied range from Boussinesq cases to compressibilities typical of the Earth's whole mantle. Compressibility has little effect on the spatial structure of steady convection when the superadiabatic temperature drop across the shell ΔTsa is comparable to a characteristic adiabatic temperature. When ΔTsa is approximately an order of magnitude smaller than the adiabatic temperature, compressibility is significant. For all the non-Boussinesq cases, the regular polyhedral convective patterns that exist at large ΔTsa break down at small ΔTsa into highly irregular patterns; as ΔTsa decreases convection becomes penetrative in the upper portion of the shell and is strongly time dependent at Rayleigh numbers only ten times the critical Rayleigh number, 〈Ra〉cr. Viscous heating in the compressible solutions is concentrated around the upwelling plumes and is greatest near the top and bottom of the shell. Solutions with regular patterns (and large ΔTsa) remain steady up to fairly high Rayleigh numbers (100〈Ra〉cr), while solutions with irregular convective patterns are time dependent at similar Rayleigh numbers. Compressibility affects the pattern evolution of the irregular solutions, producing fewer upwelling plumes with increasing compressibility.