Coupling of thermal evolution and despinning of early Iapetus

The Cassini mission revealed two spectacular characteristics of Iapetus: (1) a geologically old and high equatorial ridge, which is unique in the Solar System and (2) a large flattening of 35 km consistent with the equilibrium figure for a hydrostatic body rotating with a period of 16 h, whereas the current spin period is 79.33 days. This study describes three-dimensional simulations of solid-state convection within an undifferentiated Iapetus. It investigates the implications for the evolution of the interior thermal structure and its spin rate and global shape using radially layered viscoelastic models. The role of the concentration in the short-lived radiogenic element [ 26 Al], just after accretion is completed, is specifically addressed. The first result is to show that whatever the [ 26 Al] value, convection occurs. As suggested by Castillo-Rogez et al. [Castillo-Rogez, J., Matson, D., Sotin, C., Johnson, T., Lunine, J., Thomas, P. [2007] Icarus, 190, 179–202], convection reduces the warming of the interior compared to the conductive evolution and therefore limits the conditions for despinning. In our calculations, two conceptual linear viscoelastic models are used. When considering a Maxwell rheology, the interior temperature (viscosity) never reaches a value high (low) enough to induce despinning. In order to promote dissipation at low temperature, a Burgers rheology, which includes an additional dissipation peak, is introduced. For favorable parameter values, this latter rheology leads to despinning. However, only models associated with large amounts of short-lived radiogenic elements ( [ 26 Al ] ⩾ 25 ppb ) lead to the observed flattening. This suggests that the accretion process needs to be completed shortly after the formation of CAIs (Calcium–Aluminum-rich Inclusions) (⩽4 Myr). For [ 26 Al] varying between 72 and 46 ppb, the observed flattening is obtained only for a limited range of initial spin period, between 9.5 and 10.2 h. For [ 26 Al] ranging between 30 and 15 ppb, initial spin rates smaller than 8.5 h are required. For smaller values of [ 26 Al], the body is too cold and viscous to acquire a significant flattening even if a rotation period close to the body disruption limit is considered. Even with a thin lithosphere during the early stage, our simulations show that Iapetus never reaches the equilibrium figure for a hydrostatic body due to the non-zero rigidity of the lithosphere. The 35 km value of the flattening is the result of the partial relaxation of an ancient larger flattening ranging between 45 and 80 km, depending on the evolution of the lithosphere thickness mainly controlled by the radiogenic content. A thin lithosphere is consistent with an early building of the equatorial ridge. The lithosphere thickening due to interior cooling can explain the preservation of the ridge throughout the remaining evolution of Iapetus.