Icy Worlds

Rotating Convection With A Melting Boundary: An Application To The Icy Moons

By Keith Cowing
Status Report
astro-ph.EP
December 17, 2024
Filed under , , , , , , , , ,
Rotating Convection With A Melting Boundary: An Application To The Icy Moons
Three-dimensional renderings of snapshots of four selected simulations. For each simulation, the inner sphere shows the temperature at π‘Ÿ = π‘Ÿπ‘– + 0.02 atop the inner thermal boundary layer, while the outer surface corresponds to the melting radius π‘ŸM. Equatorial and meridional slices show the temperature in the liquid and solid phase with two different separated colormaps. — astro-ph.EP

A better understanding of the ice-ocean couplings is required to better characterise the hydrosphere of the icy moons. Using global numerical simulations in spherical geometry, we have investigated here the interplay between rotating convection and a melting boundary. To do so, we have implemented and validated a phase field formulation in the open-source code MagIC.

We have conducted a parameter study varying the influence of rotation, the vigour of the convective forcing and the melting temperature. We have evidenced different regimes akin to those already found in previous monophasic models in which the mean axisymmetric ice crust transits from pole-ward thinning to equator-ward thinning with the increase of the rotational constraint on the flow. The derivation of a perturbative model of heat conduction in the ice layer enabled us to relate those mean topographic changes to the underlying latitudinal heat flux variations at the top of the ocean.

The phase change has also been found to yield the formation of sizeable non-axisymmetric topography at the solid-liquid interface with a typical size close to that of the convective columns. We have shown that the typical evolution timescale of the interface increases linearly with the crest-to-trough amplitude and quadratically with the mean melt radius. In the case of the largest topographic changes, the convective flows become quasi locked in the topography due to the constructive coupling between convection and ice melting.

The tentative extrapolation to the planetary regimes yields O(102 βˆ’ 103 ) meters for the amplitude of non-axisymmetric topography at the base of the ice layer of Enceladus and O(103 βˆ’ 104 ) meters for Titan.

T. Gastine, B. Favier

Comments: 24 pages, 17 figures, 2 tables, accepted for publication in Icarus
Subjects: Earth and Planetary Astrophysics (astro-ph.EP); Fluid Dynamics (physics.flu-dyn); Geophysics (physics.geo-ph)
Cite as: arXiv:2412.09700 [astro-ph.EP] (or arXiv:2412.09700v1 [astro-ph.EP] for this version)
https://doi.org/10.48550/arXiv.2412.09700
Focus to learn more
Submission history
From: Thomas Gastine
[v1] Thu, 12 Dec 2024 19:42:10 UTC (8,064 KB)
https://arxiv.org/abs/2412.09700
Astrobiology

Explorers Club Fellow, ex-NASA Space Station Payload manager/space biologist, Away Teams, Journalist, Lapsed climber, Synaesthete, Na’Vi-Jedi-Freman-Buddhist-mix, ASL, Devon Island and Everest Base Camp veteran, (he/him) πŸ––πŸ»