Saturn's Evolutionary History and Seismology: Survival of Deep Stably Stratified Regions in Evolutionary Models of Saturn Consistent with Ring Seismology

With recent advances in the modeling of the solar system giant planets, rapid progress has been made in understanding the remaining questions pertaining to their formation and evolution. However, this progress has largely neglected the significant constraints on the interior of Saturn’s structure imposed by the observed oscillation frequencies in its rings. Here, we study initial conditions for Saturn’s evolution that, after $4.56;\mathrm{Gyr}$ of evolution, give rise to planetary structures admitting oscillation frequencies consistent with those observed via Saturn’s ring seismology. Restricting our attention to models without compact rocky cores, we achieve simultaneous good agreement with most observed properties of Saturn at the level of current evolutionary models and with key frequencies in the observed oscillation spectrum. Our preliminary work suggests that Saturn’s interior stably stratified region may be moderately less extended ($\sim 0.4$–$0.5R_{\rm Sat}$) than previously thought, which is important for reconciling the seismic constraints with evolutionary models. We also tentatively find that the deep helium gradients inferred by previous, static structural modelling of Saturn’s ring seismology may not be required to reproduce the observed seismology data.

💡 Research Summary

This paper investigates the interior structure and evolutionary history of Saturn by directly incorporating constraints from Saturn’s ring seismology—specifically the frequencies of low‑degree (ℓ = 2) oscillation modes that are imprinted on the density waves in the rings. The authors aim to identify initial conditions that, after 4.56 Gyr of thermal, compositional, and rotational evolution, produce planetary models whose present‑day bulk properties (equatorial radius, effective temperature, J₂, J₄, and rotation period) match observations while simultaneously reproducing the key ring‑derived mode frequencies (W84.64 and especially W76.44).

To achieve this, they employ the planetary evolution code APPLE, which couples a realistic hydrogen‑helium equation of state (Chabrier & Debras 2021) with an AQUA heavy‑element EOS, includes fourth‑order theory of figures for shape deformation, and conserves angular momentum to evolve Saturn’s rotation rate. Helium immiscibility is modeled using the “scheme B” prescription with a temperature shift of +410 K to reproduce the observed atmospheric helium abundance, and a time‑lag parameter of 5 Myr is introduced to smooth the onset of demixing. The atmospheric boundary condition follows recent work by Chen et al. (2023) and is interpolated as a function of helium and metal fractions at the base of the atmosphere.

The authors explore four families of models: (1) a “gradient” case with an initial helium‑to‑metal fraction (Y′) gradient, (2) a “uniform” case where Y′ is set to the protosolar value throughout the planet, and for each of these two, a variant with a higher initial internal temperature. No compact rocky core is included, allowing the focus to remain on the extended, compositionally stratified region often referred to as a “fuzzy core.”

Oscillation frequencies are computed with a 2‑D pseudo‑spectral solver that takes a smoothed thermodynamic profile and the shape functions from the theory of figures as input. The radial direction is discretized with 120 Gauss‑Lobatto points, and the angular dependence is expanded in spherical harmonics up to ℓ = 30, with a 257‑point Gauss‑Legendre quadrature for latitude integrals. Modes are searched in the frequency interval ω ∈