The ice shell on Enceladus, an icy moon of Saturn, exhibits strong asymmetry between the northern and southern hemispheres, with all known geysers concentrated over the south pole, even though the expected pattern of tidal-rotational deformation should be symmetric between the north and south poles. Using an idealized ice evolution model, we demonstrate that this asymmetry may form spontaneously, without any noticeable a priori asymmetry (such as a giant impact or a monopole structure of geological activity), in contrast to previous studies. Infinitesimal asymmetry in the ice shell thickness due to random perturbations are found to be able to grow indefinitely, ending up significantly thinning the ice shell at one of the poles, thereby allowing fracture formation there. Necessary conditions to trigger this hemispheric symmetry breaking mechanism are found analytically. A rule of thumb we find is that, for Galilean and Saturnian icy moons, the ice shell can undergo hemispheric symmetry breaking only if the mean shell thickness is around 10-30~km.

Ionization occurs in the upper atmospheres of hot Jupiters and in the interiors of Gas Giant Planets, leading to magnetohydrodynamic (MHD) effects which couple the momentum and the magnetic field, thereby significantly altering the dynamics. In regions of moderate temperatures the gas is only partially ionized, which also leads to interactions with neutral molecules. To explore the turbulent dynamics of these regions we utilize Partially-Ionized MHD (PIMHD), a two-fluid model -- one neutral and one ionized -- coupled by a collision term proportional to the difference in velocities. Motivated by planetary settings where rotation constrains the large-scale motions to be mostly two-dimensional, we perform a suite of simulations to examine the parameter space of 2D PIMHD turbulence and pay particular attention to collisions and their role in the dynamics, dissipation, and energy exchange between the two species. We arrive at, and numerically confirm, an expression for the energy loss due to collisions in both the weakly and strongly collisional limits, and show that, in the latter limit, the neutral fluid couples to the ions and behaves as an MHD fluid. Finally, we discuss some implications of our findings to current understanding of gas giant planet atmospheres.

We report on direct numerical simulations of two-dimensional, horizontally periodic Rayleigh-B\'enard convection, focusing on its ability to drive large-scale horizontal flow that is vertically sheared. For the Prandtl numbers ($Pr$) between 1 and 10 simulated here, this large-scale shear can be induced by raising the Rayleigh number ($Ra$) sufficiently, and we explore the resulting convection for $Ra$ up to $10^{10}$. When present in our simulations, the sheared mean flow accounts for a large fraction of the total kinetic energy, and this fraction tends towards unity as $Ra\to\infty$. The shear helps disperse convective structures, and it reduces vertical heat flux; in parameter regimes where one state with large-scale shear and one without are both stable, the Nusselt number of the state with shear is smaller and grows more slowly with $Ra$. When the large-scale shear is present with $Pr\lesssim2$, the convection undergoes strong global oscillations on long timescales, and heat transport occurs in bursts. Nusselt numbers, time-averaged over these bursts, vary non-monotonically with $Ra$ for $Pr=1$. When the shear is present with $Pr\gtrsim3$, the flow does not burst, and convective heat transport is sustained at all times. Nusselt numbers then grow roughly as powers of $Ra$, but the growth rates are slower than any previously reported for Rayleigh-B\'enard convection without large-scale shear. We find the Nusselt numbers grow proportionally to $Ra^{0.077}$ when $Pr=3$ and to $Ra^{0.19}$ when $Pr=10$. Analogies with tokamak plasmas are described.

Three-dimensional studies of convection in deep spherical shells have been used to test the hypothesis that the strong jet streams on Jupiter, Saturn, Uranus, and Neptune result from convection throughout the molecular envelopes. Due to computational limitations, these simulations must adopt viscosities and heat fluxes many orders of magnitude larger than the planetary values. Several numerical investigations have identified trends for how the mean jet speed varies with heat flux and viscosity, but no previous theories have been advanced to explain these trends. Here, we show using simple arguments that if convective release of potential energy pumps the jets and viscosity damps them, the mean jet speeds split into two regimes. When the convection is weakly nonlinear, the equilibrated jet speeds should scale approximately with F/nu, where F is the convective heat flux and nu is the viscosity. When the convection is strongly nonlinear, the jet speeds are faster and should scale approximately as (F/nu)^{1/2}. We demonstrate how this regime shift can naturally result from a shift in the behavior of the jet-pumping efficiency with heat flux and viscosity. Moreover, the simulations hint at a third regime where, at sufficiently small viscosities, the jet speed becomes independent of the viscosity. We show based on mixing-length estimates that if such a regime exists, mean jet speeds should scale as heat flux to the 1/4 power. Our scalings provide a good match to the mean jet speeds obtained in previous Boussinesq and anelastic, three-dimensional simulations of convection within giant planets over a broad range of parameters. When extrapolated to the real heat fluxes, these scalings suggest that the mass-weighted jet speeds in the molecular envelopes of the giant planets are much weaker--by an order of magnitude or more--than the speeds measured at cloud level.