A Spectral Method to Compute the Tides of Laterally-Heterogeneous Bodies

Marc Rovira-Navarro,Isamu Matsuyama,Alexander Berne

DOI: https://doi.org/10.3847/PSJ/ad381f

2024-06-04

Abstract:Body tides reveal information about planetary interiors and affect their evolution. Most models to compute body tides rely on the assumption of a spherically-symmetric interior. However, several processes can lead to lateral variations of interior properties. We present a new spectral method to compute the tidal response of laterally-heterogeneous bodies. Compared to previous spectral methods, our approach is not limited to small-amplitude lateral variations; compared to finite element codes, the approach is more computationally-efficient. While the tidal response of a spherically-symmetric body has the same wave-length as the tidal force; lateral heterogeneities produce an additional tidal response with an spectra that depends on the spatial pattern of such variations. For Mercury, the Moon and Io the amplitude of this signal is as high as $1\%-10\%$ the main tidal response for long-wavelength shear modulus variations higher than $\sim 10\%$ the mean shear modulus. For Europa, Ganymede and Enceladus, shell-thickness variations of $50\%$ the mean shell thickness can cause an additional signal of $\sim 1\%$ and $\sim 10\%$ for the Jovian moons and Encelaudus, respectively. Future missions, such as $ \textit{BepiColombo}$ and $\textit{JUICE}$, might measure these signals. Lateral variations of viscosity affect the distribution of tidal heating. This can drive the thermal evolution of tidally-active bodies and affect the distribution of active regions.

Earth and Planetary Astrophysics,Geophysics

What problem does this paper attempt to address?

This paper proposes a new spectral method to calculate the tidal response of celestial bodies with lateral heterogeneity. Traditional tidal models usually assume spherical symmetry within the celestial bodies, but there are lateral structural changes in the Earth, Moon, and other planets. Compared with previous spectral methods, this approach is not limited by small-amplitude lateral variations and has higher computational efficiency compared with finite element codes. The paper shows that for celestial bodies such as Mercury, the Moon, and Io, an additional tidal response signal can reach 1% to 10% of the main response when there is a lateral shear modulus variation of more than 10%. For Europa, Ganymede, and Enceladus, a 50% variation in shell thickness may result in similar signals. The paper also discusses the impact of viscosity lateral variations on the tidal heat distribution, thereby affecting the internal evolution of celestial bodies and the distribution of surface active areas. The authors developed a software package called LOV3D to calculate the tidal response of celestial bodies with lateral heterogeneity. The research method involves transforming the viscoelastic-gravitational problem into mass, momentum, and Poisson equations and solving them using spectral methods while considering boundary conditions. The paper concludes with numerical comparisons and practical application examples.