The Influence of Heat Flux Boundary Heterogeneity on Heat Transport in Earth's Core

Rotating convection in planetary systems can be subjected to large lateral variations in heat flux from above; for example, due to the interaction between the metallic cores of terrestrial planets and their overlying silicate mantles. The boundary anomalies can significantly reorganise the pattern of convection and influence global diagnostics such as the Nusselt number. We have conducted a suite of numerical simulations of rotating convection in a spherical shell geometry comparing convection with homogeneous boundary conditions to that with two patterns of heat flux variation at the outer boundary: one hemispheric pattern, and one derived from seismic tomographic imaging of Earth's lower mantle. We consider Ekman numbers down to 10-6 and flux-based Rayleigh numbers up to 800 times critical. The heterogeneous boundary conditions tend to increase the Nusselt number relative to the equivalent homogeneous case by altering both the flow and temperature fields, particularly near the top of the convecting region. The enhancement in Nusselt number tends to increase as the amplitude and wavelength of the boundary heterogeneity is increased and as the system becomes more supercritical. In our suite of models, the increase in Nusselt number can be as large as 25%. The slope of the Nusselt-Rayleigh scaling also changes when boundary heterogeneity is included, which has implications when extrapolating to planetary conditions. Additionally, regions of effective thermal stratification can develop when strongly heterogeneous heat flux conditions are applied at the outer boundary.