Tidal Heating

Tidal heating (also termed tidal dissipation or tidal working) is the geophysical process by which orbital and rotational mechanical energy is converted into internal thermal energy via viscoelastic solid-state friction, porous-media fluid pumping, or turbulent ocean dissipation as a planetary body is periodically deformed by a time-varying gravitational tidal field. Whereas terrestrial planets derive their internal heat primarily from primordial accretional energy and the radioactive decay of long-lived isotopes (⁴⁰K, ²³²Th, ²³⁵U, and ²³⁸U), tidal heating is the dominant thermodynamic engine in the outer Solar System and in compact multi-planet systems around low-mass M-dwarf stars (such as TRAPPIST-1).

First quantitatively predicted by Stanton J. Peale, Patrick Cassen, and Ray T. Reynolds in Science in March 1979—just days before Voyager 1 discovered active silicate volcanism on Jupiter's moon Io—tidal heating generates roughly 10¹⁴ W (100,000 GW) of internal power inside Io (a surface heat flux of ~2.5 W m⁻², thirty times Earth's mean geothermal flux). Crucially for astrobiology, tidal dissipation prevents the hydrospheres of icy satellites from freezing solid despite cryogenic surface temperatures (50–110 K), sustaining long-lived subsurface liquid water oceans inside Europa, Enceladus, Ganymede, Titan, Mimas, and Triton, and expanding the canonical circumstellar habitable zone to include ice-covered ocean worlds at any heliocentric distance.

Differential Gravitational Potential and the Tidal Bulge

A satellite of mass m and mean radius R orbiting a primary body of mass M_* at radial distance r experiences a gravitational acceleration GM_*/r² at its center of mass, but a steeper gravitational gradient across its finite diameter. Expanding the primary's gravitational potential to quadrupole (second-degree harmonic) order yields the tide-raising potential U₂ at the satellite's surface:

U₂(r, θ) = −(G M R² / r³) · P₂(cos ψ) = −(G M R² / 2 r³) · (3 cos² ψ − 1)

where ψ is the angle between the sub-planet point and the surface element. This differential pull stretches the satellite into a prolate triaxial ellipsoid whose long axis points toward the primary. If the satellite is in a circular, equatorial, tidally locked (1:1 synchronous spin–orbit) state, the tidal bulge remains stationary in the body frame of the satellite and zero mechanical work is performed.

Eccentricity and Obliquity Tides: The Peale–Cassen–Reynolds Equation

For tidal heating to occur in a synchronously rotating satellite (ω = n, where n = 2π/P is the mean orbital motion), the tidal potential must vary periodically in both magnitude and orientation over each orbital period P:

1. Radial (Eccentricity) Tides: In an eccentric orbit (e > 0), the orbital distance oscillates between periapsis r_min = a(1 − e) and apoapsis r_max = a(1 + e). Because U₂ ∝ r⁻³, the amplitude of the tidal bulge grows and shrinks by a fractional factor of 3*e* each orbit.
2. Librational (Optical) Tides: While the satellite rotates on its axis at a nearly uniform angular rate n, its orbital angular velocity varies according to Kepler's second law (r² dθ/dt = const), causing the sub-planet point to oscillate (librate) east and west in longitude by ±2*e* radians each orbit.
3. Obliquity (Latitudinal) Tides: If the satellite's spin axis is tilted by an obliquity angle θ_obl relative to its orbital normal (as on Triton or early Miranda), the sub-planet point migrates north and south in latitude each orbit.

Because real geological materials are anelastic rather than perfectly elastic, the interior deformation lags behind the forcing potential by a phase angle δ, related to the dimensionless tidal quality factor Q by Q⁻¹ = tan δ ≈ sin δ (representing the fraction of peak stored elastic energy dissipated as heat per radian of cycle). Integrating the viscoelastic stress–strain work over a synchronous orbit to second order in eccentricity e and obliquity θ_obl yields the canonical Peale–Cassen–Reynolds tidal heating rate:

Ė_tidal = −(21 / 2) · (k₂ / Q) · (G M² R⁵ n e² / a⁶) − (3 / 2) · (k₂ / Q) · (G M² R⁵ n sin² θ_obl / a⁶)

Here k₂ is the dimensionless second-degree potential Love number, which measures the structural rigidity and deformability of the body (ranging from k₂ = 0 for an infinitely rigid sphere to k₂ = 3/2 for a homogeneous fluid sphere). Crucially, when a global liquid ocean decouples an outer ice shell from the underlying rocky core, k₂ increases by roughly an order of magnitude (e.g., on Europa, k₂ ≈ 0.02 if frozen solid versus k₂ ≈ 0.25 with a global ocean), creating a strong positive feedback that amplifies tidal dissipation once an ocean forms.

Self-Circularization of Isolated Orbits

Because tidal heating extracts mechanical energy E_orb = −*GM*_*m / (2*a*) from the satellite's orbit while conserving orbital angular momentum L = m √[GM_*a(1 − e²)], dissipation inside an isolated satellite causes both its semi-major axis a and its orbital eccentricity e to decay over time:

de² / dt = (2 a [1 − e²] / G M_* m) · Ė_tidal < 0

For an isolated satellite without external perturbations, the eccentricity damping timescale τ_e = e / |*de*/*dt*| is geologically brief: ~10⁶ years for Io, ~10⁷ years for Enceladus, and ~10⁸ years for Europa. Left to itself, tidal dissipation circularizes the orbit (e → 0) and extinguishes itself.

Mean-Motion Resonance as a Gravitational Pump

Sustained tidal heating over billions of years requires a mean-motion orbital resonance with at least one neighboring satellite. In the Jovian system, Io, Europa, and Ganymede are locked in the 4:2:1 Laplace resonance (n_Io − 3*n*_Eur + 2*n*_Gan = 0), while at Saturn, Enceladus and Dione occupy a 2:1 mean-motion resonance.

In this coupled dynamical engine, the rapid rotation of the giant planet (Jupiter's 9.9-hour spin or Saturn's 10.7-hour spin) raises a tidal bulge on the planet that leads the inner satellite in longitude, exerting a positive gravitational torque that transfers rotational kinetic energy and angular momentum from the planet to the satellite's orbit, driving outward migration (da/dt > 0). As the inner moon migrates outward into the resonant commensurability with the outer moon, periodic gravitational conjunctions at the same orbital longitude force the orbital eccentricity (e_forced) upward, establishing a quasi-steady equilibrium where resonant eccentricity pumping balances viscoelastic tidal circularization.

Viscoelastic Rheology of Warm Ice (Maxwell and Andrade Models)

Within an icy ocean world, solid-state tidal dissipation occurs predominantly in the warm, ductile convective basal layer of the Ice I_h shell (temperatures T ≈ 240–273 K), where the Maxwell relaxation time τ_M = η / μ (the ratio of dynamic shear viscosity η to elastic shear modulus μ ≈ 3.5 GPa) approaches the orbital forcing period P (~1.4 to 7 days). Because dynamic viscosity obeys an Arrhenius exponential law (η(T) = η₀ exp[E_a(T_m/T − 1)]), tidal dissipation peaks sharply when ωτ_M ≈ 1. Modern geodynamic models employ the Andrade viscoelastic rheology, which incorporates transient anelastic grain-boundary sliding (ε ∝ t^α, with α ≈ 0.2–0.33) and predicts broader dissipation peaks across lower temperatures than classical Maxwell rheology.

Porous Core Pumping and Turbulent Ocean Tides

On small, low-gravity moons such as Enceladus (R = 252.1 km, central pressure ~20–30 MPa), ice-shell dissipation alone (~1–5 GW) cannot match the 15 to 40 GW of conductive and advective heat radiated from the south polar region. Choblet et al. (2017) demonstrated that if Enceladus possesses an unconsolidated, porous rocky core (porosity ϕ ≈ 20–30%), periodic tidal deformation squeezes the core matrix and drives cyclic Darcy fluid flow alongside grain-to-grain frictional sliding. This porous-core tidal dissipation generates 10 to 30 GW of heat, heating ocean water to >90 °C and driving upwelling hydrothermal plumes that fuel serpentinization and thin the polar ice shell above the Tiger Stripes. A third mechanism—resonant Rossby and gravity wave dissipation directly within the liquid ocean layer—becomes significant when the ocean is thin (0.1 to 2 km) or when obliquity tides force large-amplitude barotropic currents against rough seafloor topography.

How does tidal heating keep subsurface oceans liquid far beyond the habitable zone?

In the outer Solar System (5 to 30 AU), solar irradiance is 25 to 900 times weaker than at Earth, freezing exposed surface water ice to 50–110 K. However, when a moon in an eccentric or inclined orbit flexes within its planet's gravitational gradient, internal viscoelastic friction generates continuous heat according to the Peale–Cassen–Reynolds equation (Ė_tidal ∝ k₂ e² / Q a⁶). Beneath an insulating ice shell, this internal power maintains global liquid water oceans on Europa, Enceladus, and Titan independent of stellar sunlight.

Why would tidal heating stop without an orbital resonance?

Tidal dissipation inside a satellite draws energy directly from its orbital mechanical energy while conserving angular momentum, which causes the orbit to circularize (de/dt < 0) on timescales of 10⁶ to 10⁸ years. Once the eccentricity drops to zero (e = 0), a synchronously locked moon experiences a static tidal bulge and heating ceases. A mean-motion orbital resonance (such as the 4:2:1 Laplace resonance of Io–Europa–Ganymede or the 2:1 Enceladus–Dione resonance) continuously pumps eccentricity using rotational energy transferred from the host giant planet.

Why does forming a subsurface ocean increase tidal heating by tenfold?

The tidal dissipation rate is directly proportional to the potential Love number k₂, which measures how much a body deforms under a tidal force. When an icy moon's ice shell is frozen directly to its rigid silicate core, the entire body resists flexing and k₂ is small (~0.01–0.02). Once a global liquid water layer melts between the core and the ice shell, the outer ice shell becomes mechanically decoupled from the rocky interior and flexes freely by meters (e.g., ~30 m radial tide on Europa), raising k₂ to ~0.20–0.26 and amplifying tidal heating by an order of magnitude.

1. Peale, S. J., Cassen, P., & Reynolds, R. T. (1979). "Melting of Io by Tidal Dissipation." Science, 203(4383), 892–894. https://doi.org/10.1126/science.203.4383.892
2. Tobie, G., Mocquet, A., & Sotin, C. (2005). "Tidal dissipation within large icy satellites: Applications to Europa and Titan." Icarus, 177(2), 534–549. https://doi.org/10.1016/j.icarus.2005.04.006
3. Choblet, G., Tobie, G., Sotin, C., et al. (2017). "Powering prolonged hydrothermal activity inside Enceladus." Nature Astronomy, 1(12), 841–847. https://doi.org/10.1038/s41550-017-0289-8
4. Nimmo, F., & Pappalardo, R. T. (2016). "Ocean worlds in the outer solar system." Journal of Geophysical Research: Planets, 121(8), 1378–1399. https://doi.org/10.1002/2016JE005081