Recent precise measurements at Jupiter's and Saturn's gravity fields constrain the properties of the zonal flows in the outer envelopes of these planets. A simplified dynamic equation, sometimes called the thermal wind or thermo-gravitational wind equation, establishes a link between zonal flows and the related buoyancy perturbation, which in turn can be exploited to yield the dynamic gravity perturbation. Whether or not the action of the dynamic gravity perturbation needs to be explicitly included in this equation, an effect we call the Dynamic Self Gravity (DSG), has been a matter of intense debate. We show that, under reasonable assumptions, the equation can be solved (semi) analytically. This allows us to quantify the impact of the DSG on each gravity harmonic, practically independent of the zonal flow or the details of the planetary interior model. The impact decreases with growing spherical harmonic degree l. For degrees l=2 to about l=4, the DSG is a first order effect and should be taken into account in any attempt of inverting gravity measurements for zonal flow properties. For degrees of about l=5 to roughly l=10, the relative impact of DSG is about 10% and thus seems worthwhile to include, in particular since this comes at little extra costs with the method presented here. For yet higher degrees, is seems questionable whether gravity measurements or interior models will ever reach the required precision equivalent of the DSG impact of only a few percent of less.

The new data delivered by NASA's Juno spacecraft significantly increase our understanding of Jupiter's internal dynamics. The gravity data constrain the depth of the zonal flows observed at cloud level and suggest that they slow down considerably at a depth of about $0.96\,r_J$, where $r_J$ is the mean radius at the one bar level. Juno's magnetometer reveals the planet's internal magnetic field. We combine the new zonal flow and magnetic field models with an updated electrical conductivity profile to assess the zonal wind induced dynamo action, concentrating on the outer part of Jupiter's molecular hydrogen region where the conductivity increases very rapidly with depth. Dynamo action remains quasi-stationary and can thus reasonably be estimated where the magnetic Reynolds number remains smaller than one, which is roughly the region above $0.96\,r_J$. We calculate that the locally induced radial magnetic field reaches rms values of about $10^{-6}\,$T in this region and may just be detectable by the Juno mission. Very localized dynamo action and a distinct pattern that reflects the zonal wind system increases the chance to disentangle this locally induced field from the background field. The estimates of the locally induced currents also allow calculating the zonal flow related Ohmic heating and associated entropy production. The respective quantities remain below new revised predictions for the total dissipative heating and total entropy production in Jupiter for any of the explored model combinations. Thus neither Ohmic heating nor entropy production offer additional constraints on the depth of the zonal winds.

Celestial objects host interfaces between convective and stable stratified interior regions. The interaction between both, e.g., the transfer of heat, mass, or angular momentum depends on whether and how flows penetrate into the stable layer. Powered from the unstable, convective regions, radial flows can pierce into the stable region depending on their inertia (overshooting). In rapidly rotating systems, the dynamics are strongly influenced by the Coriolis force and radial flows penetrate in stratified regions due to the geostrophic invariance of columnar convection even in the limit of vanishing inertia. Within this study, we numerically investigate both mechanisms and hence explore the nature of penetrative convection in rapidly rotating spherical shells. The study covers a broad range of system parameters, such as the strength of the stratification relative to the Coriolis force or the inertia. Guided by the application to Saturn, we model a sandwiched stable stratified layer (SSL) surrounded by two convective zones. A comprehensive analysis of the damping behavior of convective flows at the edges of the SSL showed that the mean penetration depth is controlled by the ratio of stratified and unstratified buoyancy gradients and is hence independent of rotation. A scaling law is derived and suggests that the penetration depth decreases with the square root of the ratio of unstabilizing and stabilizing entropy gradients. The influence of the Coriolis force, however, is evident by a modulation of the penetration depth along latitude, since convective columns are elongated vertically and hence pierce predominantly into the SSL around mid-latitudes and outside the tangent cylinder. Our result also show that the penetration depth decreases linearly with the flow length scale (low pass filter), confirming predictions from the linear theory of rotating partially stratified convection.

The Juno mission is delivering spectacular data of Jupiter's magnetic field, while the gravity measurements finally allow constraining the depth of the winds observed at cloud level. However, to which degree the zonal winds contribute to the planet's dynamo action remains an open question. Here we explore numerical dynamo simulations that include an Jupiter-like electrical conductivity profile and successfully model the planet's large scale field. We concentrate on analyzing the dynamo action in the Steeply Decaying Conductivity Region (SDCR) where the high conductivity in the metallic Hydrogen region drops to the much lower values caused by ionization effects in the very outer envelope of the planet. Our simulations show that the dynamo action in the SDCR is strongly ruled by diffusive effects and therefore quasi stationary. The locally induced magnetic field is dominated by the horizontal toroidal field, while the locally induced currents flow mainly in the latitudinal direction. The simple dynamics can be exploited to yield estimates of surprisingly high quality for both the induced field and the electric currents in the SDCR. These could be potentially be exploited to predict the dynamo action of the zonal winds in Jupiter's SDCR but also in other planets.

We present thermodynamic material and transport properties for the extreme conditions prevalent in the interiors of massive giant planets and brown dwarfs. They are obtained from extensive \textit{ab initio} simulations of hydrogen-helium mixtures along the isentropes of three representative objects. In particular, we determine the heat capacities, the thermal expansion coefficient, the isothermal compressibility, and the sound velocity. Important transport properties such as the electrical and thermal conductivity, opacity, and shear viscosity are also calculated. Further results for associated quantities including magnetic and thermal diffusivity, kinematic shear viscosity, as well as the static Love number $k_2$ and the equidistance are presented. In comparison to Jupiter-mass planets, the behavior inside massive giant planets and brown dwarfs is stronger dominated by degenerate matter. We discuss the implications on possible dynamics and magnetic fields of those massive objects. The consistent data set compiled here may serve as starting point to obtain material and transport properties for other substellar H-He objects with masses above one Jovian mass and finally may be used as input for dynamo simulations.

The Juno mission will measure Jupiter's magnetic field with unprecedented precision and provide a wealth of additional data that will allow to constrain the planet's interior structure and dynamics. Here we analyse 66 numerical simulations in order to explore the sensitivity of the dynamo-generated magnetic field to the planets interior properties. The degree l=4 field model VIP4 and up-to-date interior models based on ab initio simulations serve as benchmarks. Our results suggest that VIP4-like magnetic fields can be found for a number of different models. We find that whether we assume an ideal gas or use the more realistic interior model based on ab initio simulations makes no difference. However, two other factors are important. Low Rayleigh number leads to strong axial dipole contribution while the axial dipole dominance is lost when the convective driving is too strong. The required intermediate range that yields Jupiter-like magnetic fields depends on the other system properties. The second factor is the convective magnetic Reynolds number profile Rmc(r), a product of the non-axisymmetric flow velocity and electrical conductivity. We find that the depth where Rmc exceeds about 50 is a good proxy for the top of the dynamo region. When the dynamo region sits too deep, the axial dipole is once more too dominant due to geometric reasons. We extrapolated our results to Jupiter and the result suggests that the Jovian dynamo extends to 95% of the planetary radius. The zonal flows in our simulations are dominated by an equatorial jet largely confined to the molecular layer. Where the jet reaches down to higher electrical conductivities, it is gives rise to a secondary alpha-Omega dynamo that modifies the dipole-dominated field produced deeper in the planet. This secondary dynamo can lead to strong magnetic field patches at low latitudes that seem compatible with the VIP4 model.

MESSENGER magnetometer data show that Mercury's magnetic field is not only exceptionally weak but also has a unique geometry. The internal field resembles an axial dipole that is offset to the North by 20% of the planetary radius. This implies that the axial quadrupol is particularly strong while the dipole tilt is likely below 0.8 degree. The close proximity to the sun in combination with the weak internal field results in a very small and highly dynamic Hermean magnetosphere. We review the current understanding of Mercury's internal and external magnetic field and discuss possible explanations. Classical convection driven core dynamos have a hard time to reproduce the observations. Strong quadrupol contributions can be promoted by different measures, but they always go along with a large dipole tilt and generally rather small scale fields. A stably stratified outer core region seems required to explain not only the particular geometry but also the weakness of the Hermean magnetic field. New interior models suggest that Mercury's core likely hosts an iron snow zone underneath the core-mantle boundary. The positive radial sulfur gradient likely to develop in such a zone would indeed promote stable stratification. However, even dynamo models that include the stable layer show Mercury-like magnetic fields only for a fraction of the total simulation time. Large scale variations in the core-mantle boundary heat flux promise to yield more persistent results but are not compatible with the current understanding of Mercury's lower mantle.

Rayleigh-B\'enard convection in rotating spherical shells can be considered as a simplified analogue of many astrophysical and geophysical fluid flows. Here, we use three-dimensional direct numerical simulations to study this physical process. We construct a dataset of more than 200 numerical models that cover a broad parameter range with Ekman numbers spanning $3\times 10^{-7} \leq E \leq 10^{-1}$, Rayleigh numbers within the range $10^3 < Ra < 2\times 10^{10}$ and a Prandtl number unity. We investigate the scaling behaviours of both local (length scales, boundary layers) and global (Nusselt and Reynolds numbers) properties across various physical regimes from onset of rotating convection to weakly-rotating convection. Close to critical, the convective flow is dominated by a triple force balance between viscosity, Coriolis force and buoyancy. For larger supercriticalities, a subset of our numerical data approaches the asymptotic diffusivity-free scaling of rotating convection $Nu\sim Ra^{3/2}E^{2}$ in a narrow fraction of the parameter space delimited by $6\,Ra_c \leq Ra \leq 0.4\,E^{-8/5}$. Using a decomposition of the viscous dissipation rate into bulk and boundary layer contributions, we establish a theoretical scaling of the flow velocity that accurately describes the numerical data. In rapidly-rotating turbulent convection, the fluid bulk is controlled by a triple force balance between Coriolis, inertia and buoyancy, while the remaining fraction of the dissipation can be attributed to the viscous friction in the Ekman layers. Beyond $Ra \simeq E^{-8/5}$, the rotational constraint on the convective flow is gradually lost and the flow properties vary to match the regime changes between rotation-dominated and non-rotating convection. The quantity $Ra E^{12/7}$ provides an accurate transition parameter to separate rotating and non-rotating convection.

Zonal flows in rapidly-rotating celestial objects such as the Sun, gas or ice giants form in a variety of surface patterns and amplitudes. Whereas the differential rotation on the Sun, Jupiter and Saturn features a super-rotating equatorial region, the ice giants, Neptune and Uranus harbour an equatorial jet slower than the planetary rotation. Global numerical models covering the optically thick, deep-reaching and rapidly rotating convective envelopes of gas giants reproduce successfully the prograde jet at the equator. In such models, convective columns shaped by the dominant Coriolis force typically exhibit a consistent prograde tilt. Hence angular momentum is pumped away from the rotation axis via Reynolds stresses. Those models are found to be strongly geostrophic, hence a modulation of the zonal flow structure along the axis of rotation, e.g. introduced by persistent latitudinal temperature gradients, seems of minor importance. Within our study we stimulate these thermal gradients and the resulting ageostrophic flows by applying an axisymmetric and equatorially symmetric outer boundary heat flux anomaly ($Y_{20}$) with variable amplitude and sign. Such a forcing pattern mimics the thermal effect of intense solar or stellar irradiation. Our results suggest that the ageostrophic flows are linearly amplified with the forcing amplitude $q^\star$ leading to a more pronounced dimple of the equatorial jet (alike Jupiter). The geostrophic flow contributions, however, are suppressed for weak $q^\star$, but inverted and re-amplified once $q^\star$ exceeds a critical value. The inverse geostrophic differential rotation is consistently maintained by now also inversely tilted columns and reminiscent of zonal flow profiles observed for the ice giants. Analysis of the main force balance and parameter studies further foster these results.

The amount and spatial pattern of heat extracted from cores of terrestrial planets is ultimately controlled by the thermal structure of the lower rocky mantle. Using the most common model to tackle this problem, a rapidly rotating and differentially cooled spherical shell containing an incompressible and viscous liquid is numerically investigated. To gain the physical basics, we consider a simple, equatorial symmetric perturbation of the CMB heat flux shaped as a spherical harmonic $Y_{11}$. The thermodynamic properties of the induced flows mainly depend on the degree of nonlinearity parametrised by a horizontal Rayleigh number $Ra_h=q^\ast Ra$, where $q^\ast$ is the relative CMB heat flux anomaly amplitude and $Ra$ is the Rayleigh number which controls radial buoyancy-driven convection. Depending on $Ra_h$ we characterise three flow regimes through their spatial patterns, heat transport and flow speed scalings: in the linear conductive regime the radial inward flow is found to be phase shifted $90^\circ$ eastwards from the maximal heat flux as predicted by a linear quasi-geostrophic model for rapidly rotating spherical systems. The advective regime is characterised by an increased $Ra_h$ where nonlinearities become significant, but is still subcritical to radial convection. There the upwelling is dispersed and the downwelling is compressed by the thermal advection into a spiralling jet-like structure. As $Ra_h$ becomes large enough for the radial convection to set in, the jet remains identifiable on time-average and significantly alters the global heat budget in the convective regime. Our results suggest, that the boundary forcing not only introduces a net horizontal heat transport but also suppresses the convection locally to such an extent, that the net Nusselt number is reduced by up to 50\%, even though the mean CMB heat flux is conserved.

We discuss here a purely hydrodynamical mechanism to invert the sign of the kinetic helicity, which plays a key role in determining the direction of propagation of cyclical magnetism in most models of dynamo action by rotating convection. Such propagation provides a prominent, and puzzling constraint on dynamo models. In the Sun, active regions emerge first at mid-latitudes, then appear nearer the equator over the course of a cycle, but most previous global-scale dynamo simulations have exhibited poleward propagation (if they were cyclical at all). Here, we highlight some simulations in which the direction of propagation of dynamo waves is altered primarily by an inversion of the kinetic helicity throughout much of the interior, rather than by changes in the differential rotation. This tends to occur in cases with a low Prandtl number and internal heating, in regions where the local density gradient is relatively small. We analyse how this inversion arises, and contrast it to the case of convection that is either highly columnar (i.e., rapidly rotating) or locally very stratified; in both of those situations, the typical profile of kinetic helicity (negative throughout most of the northern hemisphere) instead prevails.

Within the fluid iron cores of terrestrial planets, convection and the resulting generation of global magnetic fields are controlled by the overlying rocky mantle. The thermal structure of the lower mantle determines how much heat is allowed to escape the core. Hot lower mantle features, such as the thermal footprint of a giant impact or hot mantle plumes, will locally reduce the heat flux through the core mantle boundary (CMB), thereby weakening core convection and affecting the magnetic field generation process. In this study, we numerically investigate how parametrised hot spots at the CMB with arbitrary sizes, amplitudes, and positions affect core convection and hence the dynamo. The effect of the heat flux anomaly is quantified by changes in global flow symmetry properties, such as the emergence of equatorial antisymmetric, axisymmetric (EAA) zonal flows. For purely hydrodynamic models, the EAA symmetry scales almost linearly with the CMB amplitude and size, whereas self-consistent dynamo simulations typically reveal either suppressed or drastically enhanced EAA symmetry depending mainly on the horizontal extent of the heat flux anomaly. Our results suggest that the length scale of the anomaly should be on the same order as the outer core radius to significantly affect flow and field symmetries. As an implication to Mars and in the range of our model, the study concludes that an ancient core field modified by a CMB heat flux anomaly is not able to heterogeneously magnetise the crust to the present-day level of north--south asymmetry on Mars. The resulting magnetic fields obtained using our model either are not asymmetric enough or, when they are asymmetric enough, show rapid polarity inversions, which are incompatible with thick unidirectional magnetisation.

We simulate numerically Boussinesq convection in non-rotating spherical shells for a fluid with a unity Prandtl number and Rayleigh numbers up to $10^9$. In this geometry, curvature and radial variations of the gravitationnal acceleration yield asymmetric boundary layers. A systematic parameter study for various radius ratios (from $\eta=r_i/r_o=0.2$ to $\eta=0.95$) and gravity profiles allows us to explore the dependence of the asymmetry on these parameters. We find that the average plume spacing is comparable between the spherical inner and outer bounding surfaces. An estimate of the average plume separation allows us to accurately predict the boundary layer asymmetry for the various spherical shell configurations explored here. The mean temperature and horizontal velocity profiles are in good agreement with classical Prandtl-Blasius laminar boundary layer profiles, provided the boundary layers are analysed in a dynamical frame, that fluctuates with the local and instantaneous boundary layer thicknesses. The scaling properties of the Nusselt and Reynolds numbers are investigated by separating the bulk and boundary layer contributions to the thermal and viscous dissipation rates using numerical models with $\eta=0.6$ and a gravity proportional to $1/r^2$. We show that our spherical models are consistent with the predictions of Grossmann \& Lohse's (2000) theory and that $Nu(Ra)$ and $Re(Ra)$ scalings are in good agreement with plane layer results.

Field reversals are the most spectacular changes in the geomagnetic field but remain little understood. Paleomagnetic data primarily constrain the reversal rate and provide few additional clues. Reversals and excursions are characterized by a low in dipole moment that can last for some 10kyr. Some paleomagnetic records also suggest that the field decreases much slower before an reversals than it recovers afterwards and that the recovery phase may show an overshoot in field intensity. Here we study the dipole moment variations in several extremely long dynamo simulation to statistically explored the reversal and excursion properties. The numerical reversals are characterized by a switch from a high axial dipole moment state to a low axial dipole moment state. When analysing the respective transitions we find that decay and growth have very similar time scales and that there is no overshoot. Other properties are generally similar to paleomagnetic findings. The dipole moment has to decrease to about 30% of its mean to allow for reversals. Grand excursions during which the field intensity drops by a comparable margin are very similar to reversals and likely have the same internal origin. The simulations suggest that both are simply triggered by particularly large axial dipole fluctuations while other field components remain largely unaffected. A model at a particularly large Ekman number shows a second but little Earth-like type of reversals where the total field decays and recovers after some time.

Spacecraft data reveal a very Earth-like Jovian magnetic field. This is surprising since numerical simulations have shown that the vastly different interiors of terrestrial and gas planets can strongly affect the internal dynamo process. Here we present the first numerical dynamo that manages to match the structure and strength of the observed magnetic field by embracing the newest models for Jupiter's interior. Simulated dynamo action primarily occurs in the deep high electrical conductivity region while zonal flows are dynamically constrained to a strong equatorial jet in the outer envelope of low conductivity. Our model reproduces the structure and strength of the observed global magnetic field and predicts that secondary dynamo action associated to the equatorial jet produces banded magnetic features likely observable by the Juno mission. Secular variation in our model scales to about 2000 nT per year and should also be observable during the one year nominal mission duration.

The surface winds of Jupiter and Saturn are primarily zonal. Each planet exhibits strong prograde equatorial flow flanked by multiple alternating zonal winds at higher latitudes. The depth to which these flows penetrate has long been debated and is still an unsolved problem. Previous rotating convection models that obtained multiple high latitude zonal jets comparable to those on the giant planets assumed an incompressible (Boussinesq) fluid, which is unrealistic for gas giant planets. Later models of compressible rotating convection obtained only few high latitude jets which were not amenable to scaling analysis. Here we present 3-D numerical simulations of compressible convection in rapidly-rotating spherical shells. To explore the formation and scaling of high-latitude zonal jets, we consider models with a strong radial density variation and a range of Ekman numbers, while maintaining a zonal flow Rossby number characteristic of Saturn. All of our simulations show a strong prograde equatorial jet outside the tangent cylinder. At low Ekman numbers several alternating jets form in each hemisphere inside the tangent cylinder. To analyse jet scaling of our numerical models and of Jupiter and Saturn, we extend Rhines scaling based on a topographic $\beta$-parameter, which was previously applied to an incompressible fluid in a spherical shell, to compressible fluids. The jet-widths predicted by this modified Rhines length are found to be in relatively good agreement with our numerical model results and with cloud tracking observations of Jupiter and Saturn.

Mars Global Surveyor measurements revealed that the Martian crust is strongly magnetized in the southern hemisphere while the northern hemisphere is virtually void of magnetization. Two possible reasons have been suggested for this dichotomy: A once more or less homogeneously magnetization may have been destroyed in the northern hemisphere by, for example, resurfacing or impacts. The alternative theory we further explore here assumes that the dynamo itself produced a hemispherical field. We use numerical dynamo simulations to study under which conditions a spatial variation of the heat flux through the core-mantle boundary (CMB) may yield a strongly hemispherical surface field. We assume that the early Martian dynamo was exclusively driven by secular cooling and we mostly concentrate on a cosine CMB heat flux pattern with a minimum at the north pole, possibly caused by the impacts responsible for the northern lowlands. This pattern consistently triggers a convective mode which is dominated by equatorially anti-symmetric and axisymmetric (EAA) thermal winds. Convective up- and down-wellings and thus radial magnetic field production then tend to concentrate in the southern hemisphere which is still cooled efficiently while the northern hemisphere remains hot. The dynamo changes from an alpha^2- for a homogeneous CMB heat flux to an alpha-Omega-type in the hemispherical configuration. These dynamos reverse on time scales of about 10 kyrs. This too fast to allow for the more or less unidirectional magnetization of thick crustal layer required to explain the strong magnetization in the southern hemisphere.

Planetary and stellar magnetic fields are thought to be sustained by helical motions ($\alpha$-effect) and, if present, differential rotation ($\Omega$-effect). In the Sun, the strong differential rotation in the tachocline is responsible for an efficient $\Omega$-effect creating a strong axisymmetric azimuthal magnetic field. This is a prerequisite for Parker dynamo waves that may be responsible for the solar cycle. In the liquid iron cores of terrestrial planets, the Coriolis force organizes convection into columns with a strong helical flow component. These likely dominate magnetic field generation while the $\Omega$-effect is of secondary importance. Here we use numerical simulations to show that the planetary dynamo scenario may change when the heat flux through the outer boundary is higher in one hemisphere than in the other. A hemispherical dynamo is promoted that is dominated by fierce thermal wind responsible for a strong $\Omega$-effect. As a consequence Parker dynamo waves are excited equivalent to those predicted for the Sun. They obey the same dispersion relation and propagation characteristics. We suggest that Parker waves may therefore also play a role in planetary dynamos for all scenarios where zonal flows become an important part of convective motions.

Stellar differential rotation can be separated into two main regimes: solar-like when the equator rotates faster than the poles and anti-solar when the polar regions rotate faster than the equator. We investigate the transition between these two regimes with 3-D numerical simulations of rotating spherical shells. We conduct a systematic parameter study which also includes models from different research groups. We find that the direction of the differential rotation is governed by the contribution of the Coriolis force in the force balance, independently of the model setup (presence of a magnetic field, thickness of the convective layer, density stratification). Rapidly-rotating cases with a small Rossby number yield solar-like differential rotation, while weakly-rotating models sustain anti-solar differential rotation. Close to the transition, the two kinds of differential rotation are two possible bistable states. This study provides theoretical support for the existence of anti-solar differential rotation in cool stars with large Rossby numbers.

Observations of active M dwarfs show a broad variety of large-scale magnetic fields encompassing dipole-dominated and multipolar geometries. We detail the analogy between some anelastic dynamo simulations and spectropolarimetric observations of 23 M stars. In numerical models, the relative contribution of inertia and Coriolis force in the global force balance -estimated by the so-called local Rossby number- is known to have a strong impact on the magnetic field geometry. We discuss the relevance of this parameter in setting the large-scale magnetic field of M dwarfs.

The observed surface dynamics of Jupiter and Saturn is dominated by a banded system of zonal winds. Their depth remains unclear but they are thought to be confined to the very outer envelopes where hydrogen remains molecular and the electrical conductivity is small. The dynamo maintaining the dipole-dominated magnetic fields of both gas giants likely operates in the deeper interior where hydrogen assumes a metallic state. Here, we present numerical simulations that attempt to model both the zonal winds and the interior dynamo action in an integrated approach. Using the anelastic version of the MHD code MagIC, we explore the effects of density stratification and radial electrical conductivity variation. The electrical conductivity is mostly assumed to remain constant in the thicker inner metallic region and it decays exponentially towards the outer boundary throughout the molecular envelope. Our results show that the combination of stronger density stratification and weaker conducting outer layer is essential for reconciling dipole dominated dynamo action and a fierce equatorial zonal jet. Previous simulations with homogeneous electrical conductivity show that both are merely exclusive, with solutions either having strong zonal winds and multipolar magnetic fields or weak zonal winds and dipole-dominated magnetic fields. All jets tend to be geostrophic and therefore reach right through the convective shell in our simulations. The particular setup explored here allows a strong equatorial jet to remain confined to the weaker conducting outer region where it does not interfere with the deeper seated dynamo action. The flanking mid to high latitude jets, on the other hand, have to remain faint to yield a strongly dipolar magnetic field. The fiercer jets on Jupiter and Saturn only seem compatible with the observed dipolar fields when they remain confined to a weaker conducting outer layer.

The surface zonal winds observed in the giant planets form a complex jet pattern with alternating prograde and retrograde direction. While the main equatorial band is prograde on the gas giants, both ice giants have a pronounced retrograde equatorial jet. We use three-dimensional numerical models of compressible convection in rotating spherical shells to explore the properties of zonal flows in different regimes where either rotation or buoyancy dominates the force balance. We conduct a systematic parameter study to quantify the dependence of zonal flows on the background density stratification and the driving of convection. We find that the direction of the equatorial zonal wind is controlled by the ratio of buoyancy and Coriolis force. The prograde equatorial band maintained by Reynolds stresses is found in the rotation-dominated regime. In cases where buoyancy dominates Coriolis force, the angular momentum per unit mass is homogenised and the equatorial band is retrograde, reminiscent to those observed in the ice giants. In this regime, the amplitude of the zonal jets depends on the background density contrast with strongly stratified models producing stronger jets than comparable weakly stratified cases. Furthermore, our results can help to explain the transition between solar-like and "anti-solar" differential rotations found in anelastic models of stellar convection zones. In the strongly stratified cases, we find that the leading order force balance can significantly vary with depth (rotation-dominated inside and buoyancy-dominated in a thin surface layer). This so-called "transitional regime" has a visible signature in the main equatorial jet which shows a pronounced dimple where flow amplitudes notably decay towards the equator. A similar dimple is observed on Jupiter, which suggests that convection in the planet interior could possibly operate in this regime.

We solve the equations of motion of a one-dimensional planar Heisenberg (or Vaks-Larkin) model consisting of a system of interacting macro-spins aligned along a ring. Each spin has unit length and is described by its angle with respect to the rotational axis. The orientation of the spins can vary in time due to random forcing and spin-spin interaction. We statistically describe the behaviour of the sum of all spins for different parameters. The term "domino model" in the title refers to the interaction among the spins. We compare the model results with geomagnetic field reversals and find strikingly similar behaviour. The aggregate of all spins keeps the same direction for a long time and, once in a while, begins flipping to change the orientation by almost 180 degrees (mimicking a geomagnetic reversal) or to move back to the original direction (mimicking an excursion). Most of the time the spins are aligned or anti-aligned and deviate only slightly with respect to the rotational axis (mimicking the secular variation of the geomagnetic pole with respect to the geographic pole). Reversals are fast compared to the times in between and they occur at random times, both in the model and in the case of the Earth's magnetic field.

Context: observations of rapidly rotating M dwarfs show a broad variety of large-scale magnetic fields encompassing dipole-dominated and multipolar geometries. In dynamo models, the relative importance of inertia in the force balance -- quantified by the local Rossby number -- is known to have a strong impact on the magnetic field geometry. Aims: we aim to assess the relevance of the local Rossby number in controlling the large-scale magnetic field geometry of M dwarfs. Methods: we explore the similarities between anelastic dynamo models in spherical shells and observations of active M-dwarfs, focusing on field geometries derived from spectropolarimetric studies. To do so, we construct observation-based quantities aimed to reflect the diagnostic parameters employed in numerical models. Results: the transition between dipole-dominated and multipolar large-scale fields in early to mid M dwarfs is tentatively attributed to a Rossby number threshold. We interpret late M dwarfs magnetism to result from a dynamo bistability occurring at low Rossby number. By analogy with numerical models, we expect different amplitudes of differential rotation on the two dynamo branches.

Context: dynamo action in giant planets and rapidly rotating stars leads to a broad variety of magnetic field geometries including small scale multipolar and large scale dipole-dominated topologies. Previous dynamo models suggest that solutions become multipolar once inertia becomes influential. Being tailored for terrestrial planets, most of these models neglected the background density stratification. Aims: we investigate the influence of the density stratification on convection-driven dynamo models. Methods: three-dimensional nonlinear simulations of rapidly rotating spherical shells are employed using the anelastic approximation to incorporate density stratification. A systematic parametric study for various density stratifications and Rayleigh numbers allows to explore the dependence of the magnetic field topology on these parameters. Results: anelastic dynamo models tend to produce a broad range of magnetic field geometries that fall on two distinct branches with either strong dipole-dominated or weak multipolar fields. As long as inertia is weak, both branches can coexist but the dipolar branch vanishes once inertia becomes influential. The dipolar branch also vanishes for stronger density stratifications. The reason is the concentration of the convective columns in a narrow region close to the outer boundary equator, a configuration that favors non-axisymmetric solutions. In multipolar solutions, zonal flows can become significant and participate in the toroidal field generation. Parker dynamo waves may then play an important role close to onset of dynamo action leading to a cyclic magnetic field behavior. Conclusion: Our simulations also suggest that the fact that late M dwarfs have dipolar or multipolar magnetic fields can be explained in two ways. They may differ either by the relative influence of inertia or fall into the regime where both types of solutions coexist.

The banded structures observed on the surfaces of the gas giants are associated with strong zonal winds alternating in direction with latitude. We use three-dimensional numerical simulations of compressible convection in the anelastic approximation to explore the properties of zonal winds in rapidly rotating spherical shells. Since the model is restricted to the electrically insulating outer envelope, we therefore neglect magnetic effects. A systematic parametric study for various density scaleheights and Rayleigh numbers allows to explore the dependence of convection and zonal jets on these parameters and to derive scaling laws. While the density stratification affects the local flow amplitude and the convective scales, global quantities and zonal jets properties remain fairly independent of the density stratification. The zonal jets are maintained by Reynolds stresses, which rely on the correlation between zonal and cylindrically radial flow components. The gradual loss of this correlation with increasing supercriticality hampers all our simulations and explains why the additional compressional source of vorticity hardly affects zonal flows. All these common features may explain why previous Boussinesq models were already successful in reproducing the morphology of zonal jets in gas giants.

Geomagnetism is characterized by intermittent polarity reversals and rapid fluctuations. We have recently proposed a coupled macro-spin model to describe these dynamics based on the idea that the whole dynamo mechanism is described by the coherent interactions of many small dynamo elements. In this paper, we further develop this idea and construct a minimal model for magnetic variations. This simple model naturally yields many of the observed features of geomagnetism: its time evolution, the power spectrum, the frequency distribution of stable polarity periods, etc. This model has coexistent two phases; i.e. the cluster phase which determines the global dipole magnetic moment and the expanded phase which gives random perpetual perturbations that yield intermittent polarity flip of the dipole moment. This model can also describe the synchronization of the spin oscillation. This corresponds to the case of sun and the model well describes the quasi-regular cycles of the solar magnetism. Furthermore, by analyzing the relevant terms of MHD equation based on our model, we have obtained a scaling relation for the magnetism for planets, satellites, sun, and stars. Comparing it with various observations, we can estimate the scale of the macro-spins.

Various possibilities are currently under discussion to explain the observed weakness of the intrinsic magnetic field of planet Mercury. One of the possible dynamo scenarios is a dynamo with feedback from the magnetosphere. Due to its weak magnetic field Mercury exhibits a small magnetosphere whose subsolar magnetopause distance is only about 1.7 Hermean radii. We consider the magnetic field due to magnetopause currents in the dynamo region. Since the external field of magnetospheric origin is antiparallel to the dipole component of the dynamo field, a negative feedback results. For an alpha-omega-dynamo two stationary solutions of such a feedback dynamo emerge, one with a weak and the other with a strong magnetic field. The question, however, is how these solutions can be realized. To address this problem, we discuss various scenarios for a simple dynamo model and the conditions under which a steady weak magnetic field can be reached. We find that the feedback mechanism quenches the overall field to a low value of about 100 to 150 nT if the dynamo is not driven too strongly.

The transition from liquid metal to silicate rock in the cores of the terrestrial planets is likely to be accompanied by a gradient in the composition of the outer core liquid. The electrical conductivity of a volatile enriched liquid alloy can be substantially lower than a light-element-depleted fluid found close to the inner core boundary. In this paper, we investigate the effect of radially variable electrical conductivity on planetary dynamo action using an electrical conductivity that decreases exponentially as a function of radius. We find that numerical solutions with continuous, radially outward decreasing electrical conductivity profiles result in strongly modified flow and magnetic field dynamics, compared to solutions with homogeneous electrical conductivity. The force balances at the top of the simulated fluid determine the overall character of the flow. The relationship between Coriolis and Lorentz forces near the outer boundary controls the flow and magnetic field intensity and morphology of the system. Our results imply that a low conductivity layer near the top of Mercury's liquid outer core is consistent with its weak magnetic field.