We present a set of four global Eulerian/semi-Lagrangian fluid solver (EULAG) hydrodynamical (HD) and magnetohydrodynamical (MHD) simulations of solar convection, two of which are restricted to the nominal convection zone, and the other two include an underlying stably stratified fluid layer. While all four simulations generate reasonably solar-like latitudinal differential rotation profiles where the equatorial region rotates faster than the polar regions, the rotational isocontours vary significantly among them. In particular, the purely HD simulation with a stable layer alone can break the Taylor–Proudman theorem and produce approximately radially oriented rotational isocontours at medium to high latitudes. We trace this effect to the buildup of a significant latitudinal temperature gradient in the stable fluid immediately beneath the convection zone, which imprints itself on the lower convection zone. It develops naturally in our simulations as a consequence of convective overshoot and rotational influence of rotation on convective energy fluxes. This favors the establishment of a thermal wind balance that allows evading the Taylor–Proudman constraint. A much smaller latitudinal temperature gradient develops in the companion MHD simulation that includes a stable fluid layer, reflecting the tapering of deep convective overshoot that occurs at medium to high latitudes, which is caused by the strong magnetic fields that accumulate across the base of the convection zone. The stable fluid layer also has a profound impact on the large-scale magnetic cycles developing in the two MHD simulations. Even though both simulations operate in the same convective parameter regime, the simulation that includes a stable layer eventually loses cyclicity and transits to a non-solar, steady quadrupolar state.
展开▼
