---
title: "An icy world can warm as it bends"
description: "A changing tide can turn orbital work into heat. Follow the lag, then see why a strong gravity response alone cannot identify an underground ocean."
canonical_url: "https://sandboxphysics.com/stories/ocean-worlds/how-tides-warm-an-icy-world/"
language: "en"
page_type: "story"
alternate_languages:
  zh-CN: "https://sandboxphysics.com/zh/stories/ocean-worlds/how-tides-warm-an-icy-world/"
last_modified: "2026-10-28"
---

[Ocean worlds from material memory to tidal heat](https://sandboxphysics.com/stories/ocean-worlds/)

AC12 · Start

# An icy world can warm as it bends

Sandbox Physics · 2026-10-28

An eccentric orbit makes the tide change even when a moon keeps one face toward its planet. If the interior responds late, the planet does net work on the deforming moon. Some of that work becomes heat. Whether the heat sustains liquid water depends on where it is deposited and how quickly the interior loses it.

## Keeping one face toward the planet does not freeze the tide

Imagine standing on a small frozen moon. The planet stays near one part of the sky, but along an eccentric orbit its distance changes. Its apparent direction also rocks slightly because orbital speed changes while uniform synchronous spin does not. Both effects alter the stretching force. A perfectly circular, exactly synchronous orbit removes this particular changing tide.

The force that deforms the moon is a difference: the planet pulls on its near side, centre and far side by different amounts. Much of the common pull merely carries the entire moon around its orbit. Removing that common acceleration leaves the tide. In the Lab the arrows show the changing part of this differential pull, while the steady tidal bulge has been removed.

The shell window makes the response visible: an ice layer surrounds liquid, and both boundaries move according to the same displacement solution. The world is deliberately synthetic, with a radius of five hundred kilometres and equal density in both layers. Real ice and water differ in density, and real moons can contain rock, several ice phases and oceans with their own dynamics.

Reading: Beuthe, Tides on Europa: the membrane paradigm, arXiv:1410.4735v1

## A late response leaves work behind

An ideal spring stores energy when stretched and gives it back when released. Over a closed cycle it keeps none. A material with viscosity also resists motion. Its deformation does not exactly follow the changing load, so the force and displacement trace a loop. That loop represents a net energy transfer per cycle.

The first data plot divides force and displacement by their separate peaks so that their timing can be compared. Their similar heights therefore do not mean equal physical units or equal energy. The second plot follows one radial surface component. Total heating requires the whole surface, all the changing tidal components and their physical scales; the power readout performs that calculation.

Try ideal elasticity. The body still bends, but the phase lag and heating vanish. Then choose the circular preset: the displayed changing deformation vanishes because its driver has gone. These are two different routes to zero power. One removes dissipation from the material; the other removes this time-dependent forcing.

Reading: Beuthe, Tides on Europa: the membrane paradigm, arXiv:1410.4735v1

## A large tide is not a photograph of an ocean

A subsurface liquid can make an outer shell easier to deform. But a warm, soft solid can also bend strongly. The gravity response measures a mass redistribution, not the state of every layer. Within this Lab’s equal-density model, changing the shell thickness while rescaling stiffness and viscosity can preserve the entire external response. The solid preset demonstrates that exact ambiguity.

This matters when reading actual spacecraft results. A 2025 reanalysis of Cassini tracking estimated Titan’s complex gravity response and interpreted it with layered interior models. Its preferred warm high-pressure ice structure follows from those observations together with material, thermal and structural assumptions. The coloured interior diagram in that paper is an inference, not a radar image through the moon.

The next useful measurement is one that responds differently to competing interiors. Surface displacement, rotation and gravity can complement each other in a realistic stratified moon. A 2026 Dragonfly study explores such radio measurements through simulations. Its simulations estimate the precision of a proposed surface-response measurement.

Reading: Petricca et al., Titan’s strong tidal dissipation precludes a subsurface ocean, Nature 648 (2025) · Petricca & Buccino, A Radio Science Experiment for Dragonfly to Investigate Titan’s Rotation and Interior Structure, GRL 53 (2026)

## Heat has to reach the right place

The default synthetic world produces about one hundred gigawatts in its viscoelastic ice. That number is a model prediction, not a measured moon luminosity. The thickness plot places it beside a separate conductive reference. That reference assumes fixed temperatures at the two boundaries and a constant thermal conductivity; it is not a coupled thermal evolution calculation.

Where the curves cross, the assumed heating balances the assumed conductive loss. Warming can change viscosity, convection can carry more heat, melting can move the boundary, and eccentricity itself can evolve. Those feedbacks decide whether liquid survives. The small experiment helps identify the terms that a more complete calculation must connect.

The apparent shape changes are magnified, and the screen states the gain. For each thickness and material the gain is referenced to the same small eccentricity, so increasing eccentricity visibly increases the changing bulge. Read the metre-valued displacement when comparing absolute sizes across different materials. No colour in the cutaway is a calculated temperature map.

Reading: Beuthe, Tides on Europa: the membrane paradigm, arXiv:1410.4735v1 · Petricca et al., Titan’s strong tidal dissipation precludes a subsurface ocean, Nature 648 (2025)

## [Bend the same world in three different ways](https://sandboxphysics.com/labs/tidal-locking/)

1. Watch the delayed-response preset, then try ideal elasticity without changing eccentricity. Compare the loop and total power.
2. Compare the shell and rescaled solid presets. Their interiors differ while the external Love response agrees.
3. Drag the liquid boundary, then switch to Orbit view and change eccentricity with the highlighted periapsis grip.

[Next: Material memory turns a tidal cycle into heat](https://sandboxphysics.com/stories/ocean-worlds/from-a-maxwell-material-to-tidal-lag/)

## Sources and calculation scope

1. [Beuthe, Tides on Europa: the membrane paradigm, arXiv:1410.4735v1](https://arxiv.org/html/1410.4735v1)

   Fixed 17 October 2014 preprint; Icarus publication DOI 10.1016/j.icarus.2014.10.027. Appendix A equations 105–116 give the finite-thickness homogeneous incompressible shell; Appendix C gives Maxwell correspondence; Appendix E derives work and power; Appendix G fixes the eccentricity/libration tide. Figures 8–9 and 12–13 expose approximation and thermal-model dependence. Full relevant sections and appendices were read. Our uniform-density Maxwell example is much simpler than a real moon.

2. [Petricca et al., Titan’s strong tidal dissipation precludes a subsurface ocean, Nature 648 (2025)](https://www.nature.com/articles/s41586-025-09818-x)

   Version of record, 17 December 2025, DOI 10.1038/s41586-025-09818-x, CC BY 4.0; full publisher-deposited text and Methods, PMC12711566, read. Figure 1 reports gravity-response estimates, Figure 2 compares conditional interior posteriors and Figure 3 is an inferred-interior schematic. Figure 1 error bars are three-sigma; Figure 2 ellipses are one-sigma. Positive reported imaginary response denotes their loss convention; our positive-time-exponential convention uses negative imaginary Love numbers for passive loss.

3. [Petricca & Buccino, A Radio Science Experiment for Dragonfly to Investigate Titan’s Rotation and Interior Structure, GRL 53 (2026)](https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2026GL125753?af=R)

   Published 20 October 2026, DOI 10.1029/2026GL125753. Full sections 1–5 and Figures 1–3 read. This is a proposed experiment and simulated formal-covariance forecast, not a measured displacement Love number. Figure 2 compares Doppler/range and mission assumptions; Figure 3 contains model-dependent structural sensitivities. Ephemeris and calibration systematics are not all captured by formal uncertainties.