A numerical model of convective heat transfer in Titan’s subsurface ocean

• Published in: Icarus (New York, N.Y. 1962) Vol. 376; p. 114853
• Main Authors: Kvorka, Jakub, Čadek, Ondřej
• Format: Journal Article
• Language: English
• Published: Elsevier Inc 01-04-2022
• Subjects: Heat flux variations; Ocean dynamics; Thermal convection; Titan; Thermal convection; Ocean dynamics; Titan; Heat flux variations
• ISSN: 0019-1035; 1090-2643
• DOI: 10.1016/j.icarus.2021.114853

It has been suggested that the long-wavelength topography of Titan is related to lateral variations in the heat flux from the ocean. Recent studies of the heat transfer in Titan’s ocean agree that the time-averaged heat flux can vary in latitude by tens of percent, but they predict different distributions of heat flux anomalies at the upper boundary of the ocean. In order to clarify this issue, we perform 115 numerical simulations of thermal convection in a rotating spherical shell, varying the mechanical boundary conditions and dimensionless input parameters (Rayleigh, Ekman and Prandtl numbers) by at least one order of magnitude. The results of the simulations are examined in terms of the modified transitional number, RG∗=RaEk12/7Pr−1. Depending on the relative importance of rotation, the heat flux maximum is located either at the equator (equatorial cooling) or at the poles (polar cooling). We demonstrate that equatorial cooling occurs when RG∗<1 or RG∗>10 while polar cooling occurs when RG∗∈〈1,10〉. Based on this result, we predict that Titan’s ocean is in the polar cooling mode and the heat flux distribution is controlled by zonal degree 2 and 4 harmonics. The predicted heat flux shows a high degree of similarity with the axisymmetric part of Titan’s long-wavelength topography, indicating a strong relationship between ocean dynamics and the processes in the ice shell.
•We investigate lateral variations in heat flux coming from Titan’s ocean.•The heat flux peaks at the poles and is dominated by zonal degree 2 and 4 terms.•The predicted heat flux explains a significant portion of Titan’s topography.•The numerical model is sensitive to boundary conditions and the Prandtl number.•A new scaling law to estimate the heat flux distribution is presented.