A model of force balance in Saturns magnetodisc

We present calculations of magnetic potential associated with the perturbation of Saturn’s magnetic field by a rotating, equatorially-situated disc of plasma. Such structures are central to the dynamics of the rapidly rotating magnetospheres of Saturn and Jupiter. They are fed' internally by sources of plasma from moons such as Enceladus (Saturn) and Io (Jupiter). We use a scaled form of Euler potentials for the Jovian magnetodisc field (Caudal, 1986). In this formalism, the magnetic field is assumed to be azimuthally symmetric about the planet's axis of rotation, and plasma temperature is constant along a field line. We perturb the dipole potential by using simplified distributions of plasma pressure and angular velocity for both planets, based on observations by Cassini (Saturn) and Voyager (Jupiter). Our results quantify the degree of radial stretching’ exerted on the dipolar field lines through the plasma’s rotational motion and pressure. A simplified version of the field model, the `homogeneous disc’, can be used to easily estimate the distance of transition in the outer magnetosphere between pressure-dominated and centrifugally-dominated disc structure. We comment on the degree of equatorial confinement as represented by the scale height associated with disc ions of varying mass and temperature. For Saturn, we identify the principal forces which contribute to the magnetodisc current and make comparisons between the field structure predicted by the model and magnetic field measurements from Cassini. For Jupiter, we reproduce Caudal’s original calculation in order to validate our model implementation. We also show that compared to Saturn, where plasma pressure gradient is, on average, weaker than centrifugal force, the outer plasmadisc of Jupiter is clearly a pressure-dominated structure.

💡 Research Summary

The paper presents a quantitative model of the force balance that shapes the magnetodisc of Saturn, and by extension, that of Jupiter. The authors adopt the Euler‑potential formalism originally developed by Caudal (1986) for the Jovian magnetodisc, which assumes azimuthal symmetry about the planetary rotation axis and a constant plasma temperature along each magnetic field line. Within this framework the magnetic field is expressed as the gradient of a scalar potential; the dipolar planetary field is then perturbed by adding contributions from a rotating, equatorially confined plasma disc.

To make the model applicable to Saturn, the authors incorporate observationally derived radial profiles of plasma pressure and angular velocity obtained from the Cassini spacecraft. The pressure profile reflects the mixture of cold water‑group ions sourced from Enceladus and hotter, more tenuous plasma from the inner magnetosphere, decreasing steeply with radial distance. The angular velocity profile is close to rigid corotation (Ω≈Ωₚ) in the inner region but shows a modest sub‑corotational lag farther out. By inserting these profiles into the Euler‑potential equations, the dipole potential is “stretched” outward, producing a more disc‑like magnetic geometry.

The degree of stretching is governed by the competition between two forces: the centrifugal force (ρ Ω² r) that tends to fling plasma outward, and the pressure gradient force (∇P) that pushes plasma inward. For Saturn, the authors find that the centrifugal term dominates on average, making the outer disc a centrifugal‑dominated structure. This dominance is quantified by a transition radius rₜ at which the two forces become comparable. To facilitate rapid estimates, the authors introduce a “homogeneous disc” approximation in which pressure and density are assumed constant with radius. Within this simplification, rₜ can be expressed analytically as a function of plasma temperature, ion mass, and rotation rate.

The model also predicts the vertical scale height H of the disc, given approximately by H≈(k T / m Ω²)¹ᐟ², where k is Boltzmann’s constant, T the plasma temperature, m the mean ion mass, and Ω the angular velocity. Because Saturn’s ion population is relatively light (dominantly H⁺ and water‑group ions) and its temperature modest, H is of order a few thousand kilometres, yielding a comparatively thin disc.

To validate the implementation, the authors reproduce Caudal’s original Jovian calculation using Voyager plasma data. In the Jovian case the pressure gradient exceeds the centrifugal term, so the outer magnetodisc is pressure‑dominated. The higher temperatures and heavier ion species (S⁺, O⁺) supplied by Io lead to a much larger scale height (tens of thousands of kilometres) and a thicker disc.

A detailed comparison with Cassini magnetometer (MAG) measurements shows good agreement in regions where the disc is strongly developed. The model captures the observed enhancement of the north‑south magnetic component (Bₙ) and the rapid change in the radial and azimuthal components (Bᵣ, B_φ) near the transition radius. By decomposing the disc current into contributions from the pressure‑gradient current and the centrifugal current, the authors identify the principal forces that sustain the magnetodisc current system at Saturn.

The paper concludes that the Euler‑potential approach, combined with realistic pressure and angular‑velocity profiles, provides a robust and computationally inexpensive tool for estimating the structure of planetary magnetodiscs. It clarifies why Saturn’s disc is primarily shaped by centrifugal forces, whereas Jupiter’s is governed by plasma pressure. The homogeneous‑disc model offers a quick method to locate the pressure‑to‑centrifugal transition and to estimate the vertical confinement of ions of different mass and temperature. These results have direct relevance for interpreting in‑situ magnetic field data, planning future missions to the outer planets, and informing global magnetohydrodynamic simulations of rapidly rotating magnetospheres.