Start with a mechanical model whose limits can be tested
Pin Beuthe’s 17 October 2014 arXiv version 1410.4735v1, later published in Icarus. Appendix A, equations 105–116, gives the finite homogeneous incompressible shell benchmark. Compare equation 108 with the independently solved tractions in our course. Taking the thin and solid limits tests whether a shell implementation preserves the right physics as its boundary moves.
Appendix C supplies the viscoelastic correspondence and Maxwell law; Appendices E and G connect dissipative work to the full eccentricity tide. Check the optical-libration term before accepting the global power coefficient. Figures 8–9 concern approximation accuracy. Figure 12 introduces thermal transport and rheology assumptions, so its equilibrium thicknesses are not a direct reading of gravity data.
The browser implements one equal-density homogeneous-shell example. It does not reproduce the paper’s depth-dependent ice or convective thermal calculations. That narrow scope makes its energy and solid-sphere limits checkable; the paper shows which additional structure a planetary application needs.
Reading: Beuthe, Tides on Europa: the membrane paradigm, arXiv:1410.4735v1
Read the measured response before the coloured interior
Use the Nature version of record published 17 December 2025, DOI 10.1038/s41586-025-09818-x. Full text, Methods and captions were read from the publisher-deposited open-access record. Figure 1 reports an estimated real gravity response 0.608 with standard error 0.048 and positive loss coefficient 0.135 with standard error 0.035. Its plotted error bars are three-sigma, while these quoted errors are one-sigma.
Figure 2 compares interior-model posteriors with one-sigma data ellipses. Figure 3 depicts an inferred interior. The warm high-pressure ice interpretation uses a stratified structure, material rheology, thermal transport and additional geophysical constraints. It is not the equal-density model, and the schematic does not independently measure liquid absence at every depth.
In Methods, follow the radio tracking analysis separately from the interior inversion. Gravity coefficients and complex tidal response are fitted from tracking; layer structure is inferred with model parameters and priors. Before comparing imaginary responses, translate the positive reported loss coefficient into our negative-imaginary passive convention. Reversing that sign without the convention would reverse the heating conclusion.
Section 2.1 of the 2026 GRL paper compares published real-response estimates near 0.6 with an analysis near 0.4, alongside continuing discussion of processing and formal errors. Keep the 2025 preferred interior attached to its particular gravity solution and structural assumptions. New measurements can test those alternatives as well as reduce statistical errors.
Reading: Petricca et al., Titan’s strong tidal dissipation precludes a subsurface ocean, Nature 648 (2025) · Beuthe, Tides on Europa: the membrane paradigm, arXiv:1410.4735v1 · Petricca & Buccino, A Radio Science Experiment for Dragonfly to Investigate Titan’s Rotation and Interior Structure, GRL 53 (2026)
A 2026 forecast is a proposed experiment, not a new detection
Pin Petricca and Buccino’s published GRL paper of 20 October 2026, DOI 10.1029/2026GL125753. Sections 1–5 and Figures 1–3 were read in full. It simulates a Dragonfly radio experiment to recover rotation and surface tidal displacement. It supplies a formal-covariance forecast; the plotted surface Love response has not thereby been measured.
Read Figure 1 as predicted displacement geometry. In Section 2.2 and Figure 2, compare Doppler-only scenarios with scenarios that add ranging and different ephemeris assumptions. More tracking does not automatically remove parameter correlations. Figure 3 varies structural parameters within the chosen interior models and predicts their response sensitivities.
Section 4 discusses limitations beyond the formal errors, including ephemeris and calibration effects. Its mission duration, tracking schedule and ranging choices are study assumptions. Real observing plans and systematic errors must be reassessed before treating a forecast interval as an achieved uncertainty. The forecast uses a licensed orbit-determination simulation; the public teaching package instead solves the finite-shell response described above.
Ask what would distinguish the competing pictures
Our toy example fixes gravity and displacement Love numbers in one ratio, so adding the latter cannot break its exact degeneracy. A realistic stratified moon need not share that ratio. The relevant research question is how much independent information the new observable carries after correlations, rheological freedom and systematic errors are included.
A useful reading exercise is to trace one arrow at a time: tracking to gravity response, response plus structure to dissipation location, and heat transport to long-term liquid survival. Each arrow needs its own assumptions. A large tidal power does not by itself locate an ocean or establish habitability.
The package reproduces our finite-shell response and conditional inference. It does not reprocess raw Cassini signals, reproduce the Nature posterior, or rerun the Dragonfly forecast. Those are distinct scientific projects. The fixed equations, figure numbers and publication dates here make it possible to read their evidence without silently merging the three tasks.
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) · Petricca & Buccino, A Radio Science Experiment for Dragonfly to Investigate Titan’s Rotation and Interior Structure, GRL 53 (2026)