Pith. sign in

REVIEW 3 major objections 4 minor 1 cited by

Thermal-orbital evolution of Eris

T0 review · 3 major / 4 minor · reviewed 2026-08-05 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Eris's synchronous spin is best explained by a subsurface ocean that decoupled its ice shell from its rocky interior.

desk verdict A careful thermal-orbital model that makes a plausible case for a subsurface ocean on Eris, but the inference is conditional on the disputed synchronous-rotation assumption. read the letter →

arxiv 2508.16532 v1 pith:OBA4UQTJ submitted 2025-08-22 astro-ph.EP

Ryunosuke Akiba , Francis Nimmo This is my paper · ORCID
classification astro-ph.EP
keywords KuiperbeltErisDysnomiaorbitalevolutionsubsurfaceoceanthermalhistorytidaldissipationAndraderheology
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper argues that the only practical way to spin Eris down from a fast post-impact rotation to its current 378.862-hour synchronous lock with Dysnomia within 4.5 billion years is for Eris to have had a subsurface ocean. A warm, convecting ice shell can do it only with a thick insulating layer, a thick hydrosphere, and very dissipative anelastic ice; once the ice's dissipation is capped at experimentally supported levels, nearly all successful simulations contain an ocean. The reason is mechanical: an ocean decouples the ice shell from the stiff rocky core, lowering the tidal response by roughly two orders of magnitude and letting Dysnomia's tides drain Eris's spin efficiently. The result matters because it would add a Kuiper belt object to the family of ocean worlds that are kept warm not by giant-planet tides but by radioactive decay and insulation. The paper also finds these oceans tend to freeze over by the present day unless porosity, gas clathrates, or antifreeze preserves them.

What carries the argument

The load-bearing mechanism is the decoupling of the ice shell from the rocky interior by a subsurface ocean: a few tens of kilometres of liquid water changes the tidal response from that of a stiff, cold body to a viscoelastic shell sliding over a fluid layer, lowering Q/k2 by roughly a factor of 100. The supporting machinery is a coupled thermal-orbital model: a 1-D finite-difference heat equation with mixing-length-theory convection, porosity evolution, and ocean growth/refreezing, feeding temperature- and rheology-dependent tidal Love numbers into the standard tidal spin-orbital equations. The Andrade beta parameter — the anelastic term in the ice compliance that sets how dissipative cold

What would settle it

A definitive determination of Eris's rotation period — for example, high-cadence space or ground-based photometry that confirms or dismisses the 18.85-hour signal as an artifact or a close-in satellite — would settle whether the synchronous assumption holds. Separately, measuring the ice's Andrade beta parameter above 3e-11 Pa^-1 s^-0.25 would reopen the no-ocean branch, and a spectroscopic detection of ammonia or methane clathrates at Eris's surface would corroborate a present-day ocean.

Watch

Extended reading notes

Core claim

Coupled 1-D thermal evolution of a differentiated Eris — rocky core plus ice shell, with convection, porosity, ocean formation, clathrate/antifreeze insulation — is integrated with the tidal spin-orbital evolution of the Eris-Dysnomia system, using the tidal response computed from the evolving internal structure. The central discovery is that the observed doubly synchronous state requires Eris to be dissipative, and the only structure that reliably delivers that dissipation is a subsurface ocean. When an ocean forms, it decouples the ice shell from the rigid interior and drops the tidal response parameter Q/k2 by about a factor of 100, letting Eris reach the present-day state by about 1.3 Gy

Load-bearing premise

Eris is actually synchronously rotating with Dysnomia at the observed 378.862-hour period; the paper sets aside the 18.85-hour periodicity reported in Gaia photometry, which could be Eris's true rotation or an undiscovered close-in satellite. If Eris is not synchronous, the 4.5-Gyr despin constraint that drives the ocean conclusion disappears.

Editorial extensions

If this is right

  • A subsurface ocean is the most probable explanation for Eris's current synchronous state; without one, successful despinning requires ice dissipation values at or beyond the upper end of experimental measurements.
  • Present-day oceans are not guaranteed: in pure-ice models every ocean refreezes by today; only porosity, clathrate lids, or antifreeze like ammonia keeps an ocean alive.
  • Eris's spin constraint restricts composition: no successful simulations with rock densities below about 3050 kg m^-3 (hydrosphere thinner than roughly 110 km) occur except with thick surface clathrates, limiting how much low-density organic-rich material Eris can hold.
  • A convecting ice shell or past/ongoing ocean implies relaxed topography and possible cryovolcanism, consistent with Eris's bright surface and the D/H ratio of its methane ice.
  • If the rotational state is as observed, the ice's Andrade beta value becomes a decisive parameter for distinguishing ocean versus no-ocean histories.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the 18.85-hour periodic signal recently reported in Gaia photometry of Eris is actually Eris's true rotation or an undiscovered close-in satellite, the synchronous assumption collapses and the ocean inference does not follow; the paper's own neglect of this signal is the main observational risk.
  • The same coupled framework could be applied to other binary Kuiper belt objects with measured spin-orbital states, such as Orcus–Vanth or Salacia–Actaea, to test whether ocean-favored despinning is a general feature of large, differentiated trans-Neptunian objects.
  • A testable prediction is that if Eris's ocean is still present, volatiles such as ammonia or methane clathrates should be detectable at the surface, and future geophysical observations of shape or moment of inertia could distinguish a frozen from a liquid hydrosphere.
  • The paper's neglect of tidal heating is quantitatively safe for Eris, but for a more massive satellite the thermal-orbital coupling would need to be two-way, with tidal heating included as an interior heat source.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper couples a 1-D thermal evolution model of a differentiated Eris (conductive/convective ice shell, possible ocean, porous or clathrate insulation, radiogenic heating) with a two-body tidal spin-orbit model to test whether Eris can be despun into its reported doubly synchronous state with Dysnomia within 4.5 Gyr. The authors find that successful spin-down is possible in a minority of their tested grid, that a subsurface ocean is present in the large majority of successful cases (77–100% depending on the insulation/heating case, rising to >98% when the Andrade β parameter is restricted to experimentally plausible values), and that oceans usually freeze by the present day unless porosity, clathrates, or antifreeze are present. The key physical mechanism is that an ocean decouples the ice shell from the rocky interior, greatly reducing Q/k2, whereas a convecting ice shell alone is generally too weakly dissipative to despin Eris unless the ice is anomalously anelastic (high β).

Significance. If the result holds, it is significant: it identifies a plausible subsurface ocean on the most massive known dwarf planet, extends the population of candidate ocean worlds to the distant Kuiper belt, and makes testable predictions about surface relaxation, shape, and D/H-based internal activity. The study is a genuine forward model: the ocean is not imposed a priori, no parameter is fitted to the target spin state, and the tidal response is computed with an open-source code (California Planetary Geophysics Code). The sensitivity coverage is unusually broad (rheology, heating rate, porosity, clathrate geometry, antifreeze, initial temperature), which strengthens the robustness of the central mechanism. The main weaknesses are that the headline statistics are grid fractions over an ad hoc parameter space rather than posterior probabilities, and that the entire inference is contingent on the disputed synchronous-rotation interpretation of Eris's photometry.

major comments (3)
  1. [§1 (para. 2), §3.1, §5] The central inference is conditional on Eris being doubly synchronous at 378.862 h. The authors explicitly set aside the 18.85 h periodicity reported by Ortiz et al. (2025), which could be Eris's true rotation or a close-in satellite. If the former, the 4.5 Gyr despin constraint disappears and the ocean preference has no observational anchor; if the latter, the two-body orbital model, mass/density estimate, and tidal history all change. Because this assumption is load-bearing, the abstract and conclusions should state the result as conditional on the Szakáts/Bernstein synchronous interpretation, or the authors should provide a quantitative robustness test. As written, the headline 'subsurface ocean is preferred' overstates the support.
  2. [§3.2–3.3, Table 3] The percentages '77%' and '>98%' are counts over a finite grid in (rho_rock, eta_ref, beta) plus selected insulation modes. The parameter bounds are ad hoc (e.g., beta up to 1e-10 Pa^-1 s^-0.25 is outside the experimental range), and no priors are defined. These numbers are therefore grid fractions, not posterior probabilities. The text usually says 'of successful simulations,' but the abstract's 'Oceans make up 77–100% of successful models' and the conclusion 'oceans are preferred' invite a probabilistic reading. Recommend explicitly labeling these as fractions of the tested parameter grid and, if a probabilistic claim is intended, integrating over stated priors.
  3. [Table 1; §2.3, Eq. (9)] Equation (9) depends on the secondary mass M_j through M_j^2, but Table 1 lists only an upper bound f=0.0084 for the mass ratio, and the text says 'we use the central values' without defining a central Dysnomia mass. If the models adopt f as the nominal value, the tidal torque is maximized, making the success rates optimistic. Please state the adopted Dysnomia mass explicitly and, ideally, test a lower mass (e.g., f=0.004 or 0.002) to show how the ocean fraction depends on this assumption.
minor comments (4)
  1. [§2.3] The notation Q/k2 is nonstandard and easy to confuse with the usual dissipation factor k2/Q. Please define it explicitly at first use and state that lower Q/k2 means more tidal dissipation.
  2. [Abstract, §4.4, Table 3] The unqualified '77–100%' in the abstract conflicts with the reduced-heating rows in Table 3 (67.7% and 30.0% ocean fractions). Add 'excluding the reduced-heating cases' or otherwise qualify the range in the abstract.
  3. [§2.3] The assumption of zero eccentricity is made despite the nonzero eccentricity reported by Holler et al. (2021). A one-sentence justification of why eccentric tides cannot qualitatively change the conclusion would help.
  4. [Data availability] The statement 'Codes ... available upon request' is weaker than a permanent repository. Please deposit the thermal-orbital code and the tidal code version used, with version identifiers, to improve reproducibility.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the ocean preference emerges from a forward thermal-orbital model matched to an independent observed spin state.

full rationale

The paper builds a forward chain: choose interior parameters (rock density, ice viscosity, Andrade beta, insulation mode); integrate a 1-D thermal evolution model; compute tidal Love numbers from the resulting viscoelastic structure with a benchmarked open-source code; then integrate the spin-orbit equations and accept a model only if it reproduces the observed synchronous spin-orbit state within 4.5 Gyr. The conclusion that oceans are preferred is not imposed or fitted to the target observable: it arises because an ocean decouples the ice shell from the rocky interior and lowers Q/k2 in the forward tidal calculation. The spin-down constraint itself is taken from Nimmo & Brown (2023), a peer-reviewed, observation-based inference that is external to the present model and not re-fit here; this self-citation is independent support rather than a circular premise. The Andrade beta restriction is based on experimental/computational compilations (Bierson 2024), not on the success criterion. The neglect of the ~18.85 h periodicity (Ortiz et al. 2025) is an observational-assumption risk affecting robustness, but it is not a circular reduction because the model equations do not assume the ocean conclusion. No equation or fitted parameter is definitionally equivalent to the claimed result.

Assumptions & free parameters 6 free parameters · 7 assumptions · 0 invented entities

The model introduces no new physical entities. Its conclusions rest on standard interior assumptions, a suite of rheological and thermal parameters drawn from prior literature, and a key observational assumption (synchronous rotation). The main free parameters are explored as grids rather than fitted, but the success fractions are conditional on the ranges chosen.

free parameters (6)
  • Ice reference viscosity eta_ref = 1e13 - 1e15 Pa s (grid)
    Controls convective heat transport and tidal dissipation; varied over a grid, with outcome strongly depending on its value.
  • Andrade beta parameter = 1e-12 to 1e-10 Pa^-1 s^-0.25; main conclusion uses beta <= 3e-11
    Sets anelastic dissipation in ice; the no-ocean successful cases mostly require the highest values, so the experimental upper bound is load-bearing.
  • Rock density rho_rock = 2800-3500 kg m^-3
    Determines ice fraction and hydrosphere thickness (70-180 km), which strongly affects ocean formation and spin-down success.
  • Initial temperature = 100 K (150 K tested)
    Cold start assumption; warmer start slightly expands ocean/convection regions but has small effect on conclusions.
  • Clathrate insulation layer parameters = 5 or 10 km thickness; conductivity 0.6 or 1.0 W m^-1 K^-1
    Assigned values for insulating layers, which strongly influence ocean survival and success rates.
  • Antifreeze melting point depression = -20 K (10 wt% ammonia)
    Assumed shift of the ice-water phase curve; antifreeze cases produce 100% ocean successes.
assumptions (7)
  • domain assumption Eris is fully differentiated into a rocky core and an ice shell, with no porosity affecting bulk density.
    Motivates the interior model; supported by high bulk density and prior work, but is an assumption about formation.
  • domain assumption Ice viscosity follows a Newtonian Arrhenius law with reference viscosity at 270 K.
    Standard in ice shell modeling; Section 2.1, Eq. 3.
  • domain assumption Solid-state convection is parameterized with mixing length theory (Kamata 2018).
    Section 2.1: the MLT parameterization is used as a stand-in for full 3-D convection, with a fixed mixing length parameter x=1.35.
  • domain assumption Tidal response is computed for an incompressible, spherically symmetric body; ocean dissipation is neglected.
    Section 2.3: the tidal Love number calculation assumes incompressibility and ignores any dissipation in the ocean itself.
  • domain assumption Radiogenic heating follows CI chondrite abundances; no 26Al.
    Section 2.1: abundances of radioisotopes in Eris are uncertain; reduced heating cases are tested as sensitivity.
  • domain assumption Eris is currently synchronously rotating with Dysnomia, neglecting the 18.85-hour periodicity reported by Ortiz et al. (2025).
    Section 1: the whole spin-down constraint depends on this observational choice, which is flagged but not modeled.
  • domain assumption The Andrade beta experimental upper bound is 3e-11 Pa^-1 s^-0.25.
    Section 4.3: used to compute the >98% ocean fractions; if higher beta values were physically allowed, more no-ocean cases would succeed.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Thermal-orbital evolution of Eris." pith.science (2026). https://pith.science/paper/OBA4UQTJ

@misc{pith2026250816532,
  author       = {Pith},
  title        = {Pith review of: Thermal-orbital evolution of Eris},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/OBA4UQTJ}},
  note         = {Machine review of arXiv:2508.16532}
}
abstract

The large Kuiper Belt object (KBO) Eris is nearly as big as Pluto and has a small moon, Dysnomia. Constraints on the system's spin and orbit characteristics were recently used to argue for a dissipative Eris, requiring a differentiated structure but not necessarily a subsurface ocean. Here, we model the thermal history of Eris coupled to its spin-orbital evolution, finding a subsurface ocean is preferred in order for Eris to be sufficiently dissipative. Spinning down Eris without an ocean is difficult, requiring a warm convecting ice shell protected by a thick insulating layer and very dissipative anelastic behavior in ice. Oceans make up 77-100% of successful thermal-orbital evolution models, depending on the parameters assumed, which increases to >98% when the Andrade $\beta$ parameter for ice is restricted to $\beta\leq3\times10^{-11}$ Pa$^{-1}$ s$^{-0.25}$. Oceans freeze over by the present day unless insulation (porosity, gas clathrates) or antifreeze are present.

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Synchronous Rotation in the (120347) Salacia-Actaea System

    astro-ph.EP 2025-09 conditional novelty 6.0 of 10

    Observations show that Salacia and Actaea are likely in fully synchronous rotation, with Salacia's albedo-variation lightcurve matching the 5.49389-day mutual orbital period.

Reference graph

Works this paper leans on

11 extracted references · 11 canonical work pages · cited by 1 Pith paper

  1. [3]

    Geophysical Research Letters 46, 1963–1972

    Conditions for the long-term preservation of a deep brine reservoir in Ceres. Geophysical Research Letters 46, 1963–1972. Castillo-Rogez, J.C., Kalousová, K.,

    work page 1963

  2. [11]

    Cambridge university press

    Geodynamics. Cambridge university press. Vance,S.,Harnmeijer,J.,Kimura,J.,Hussmann,H.,DeMartin,B.,Brown,J.M.,2007. Hydrothermalsystemsinsmalloceanplanets. Astrobiology 7, 987–1005. : Preprint submitted to Elsevier Page 15 of 15

    work page 2007

  3. [1972]

    Methods in computational physics 11, 217–295

    Seismic surface waves. Methods in computational physics 11, 217–295. Truong,N.,Glein,C.R.,Lunine,J.I.,2024. AbroadsetofsolarandcosmochemicaldataindicateshighCNOabundancesforthesolarsystem. The Astrophysical Journal 976,

    work page 2024

  4. [1985]

    Geophysical & Astrophysical Fluid Dynamics 33, 63–96

    A mathematical model of magma transport in the asthenosphere. Geophysical & Astrophysical Fluid Dynamics 33, 63–96. Glein, C.R., Grundy, W.M., Lunine, J.I., Wong, I., Protopapa, S., Pinilla-Alonso, N., Stansberry, J.A., Holler, B.J., Cook, J.C., Souza-Feliciano, A.C.,2024. ModerateD/HratiosinmethaneiceonErisandmakemakeasevidenceofhydrothermalormetamorphic...

    work page 2024

  5. [2000]

    Science 287, 107–109

    Evidence for crystalline water and ammonia ices on Pluto’s satellite Charon. Science 287, 107–109. Castillo-Rogez,J.,Weiss,B.,Beddingfield,C.,Biersteker,J.,Cartwright,R.,Goode,A.,MelwaniDaswani,M.,Neveu,M.,2023. Compositionsand interior structures of the large moons of Uranus and implications for future spacecraft observations. Journal of Geophysical Rese...

    work page 2023

  6. [2002]

    Icarus 157, 219–227

    Monte carlo modeling of the thermal conductivity of porous cometary ice. Icarus 157, 219–227. Sicardy,B.,Ortiz,J.,Assafin,M.,Jehin,E.,Maury,A.,Lellouch,E.,Hutton,R.G.,Braga-Ribas,F.,Colas,F.,Hestroffer,D.,etal.,2011. APluto-like radius and a high albedo for the dwarf planet Eris from an occultation. Nature 478, 493–496. Szakáts,R.,Kiss,C.,2023. Rotational...

    work page 2011

  7. [2006]

    Icarus 185, 258–273

    Subsurface oceans and deep interiors of medium-sized outer planet satellites and large trans-neptunian objects. Icarus 185, 258–273. Journaux,B.,Brown,J.M.,Pakhomova,A.,Collings,I.E.,Petitgirard,S.,Espinoza,P.,BoffaBallaran,T.,Vance,S.D.,Ott,J.,Cova,F.,etal.,2020. Holistic approach for studying planetary hydrospheres: Gibbs representation of ices thermody...

    work page 2020

  8. [2017]

    Nature Astronomy 1, 841–847

    Powering prolonged hydrothermal activity inside enceladus. Nature Astronomy 1, 841–847. Choukroun,M.,Grasset,O.,2010. Thermodynamicdataandmodelingofthewaterandammonia-waterphasediagramsupto2.2GPaforplanetary geophysics. The Journal of chemical physics

    work page 2010

Show all 11 references
  1. [2019]

    Nature Astronomy 3, 802–807

    Early formation of moons around large trans-neptunian objects via giant impacts. Nature Astronomy 3, 802–807. Arakawa,S.,Kamata,S.,Genda,H.,2025. Lowrockmassfractionwithintrans-neptunianobjectsinferredfromthespin–orbitevolutionofOrcus– Vanth and Salacia–Actaea. Journal of Geop...

    2025

  2. [2024]

    URL:https://doi.org/10.5281/zenodo.14029354

    California planetary geophysics code. URL:https://doi.org/10.5281/zenodo.14029354. Fowler, A.,

    doi

  3. [2025]

    Monthly Notices of the Royal Astronomical Society: Letters 538, L37–L42

    Analysis of the Gaia DR3 photometry of the dwarf planet Eris. Monthly Notices of the Royal Astronomical Society: Letters 538, L37–L42. Park,R.,Jacobson,R.,GomezCasajus,L.,Nimmo,F.,Ermakov,A.,Keane,J.,McKinnon,W.,Stevenson,D.,Akiba,R.,Idini,B.,etal.,2025. Io’s tidal response pr...

    2025

Pith tools

Reviewed August 5, 2026 · model on record in the stance chip above.