arxiv: 2601.07080 · v2 · submitted 2026-01-11 · 🌌 astro-ph.EP

Magma Ocean Waves and Thermal Variability on Lava Worlds

Mohammad Farhat , Eugene Chiang This is my paper

Pith reviewed 2026-05-16 14:58 UTC · model grok-4.3

classification 🌌 astro-ph.EP

keywords lava worldsmagma oceanstidal heatingthermal variabilitywave interferenceexoplanet interiorsorbital dynamicsplanetary heat transport

0 comments

The pith

Tidal waves in deep magma oceans on lava worlds produce wandering hotspots and aperiodic thermal variability through wave interference.

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

The paper shows that modest orbital eccentricities, sustained by companion planets, can melt deep layers of short-period rocky planets into magma oceans via tidal heating. It models how these oceans respond to the tidal pull with multiple wave modes, bounded by a day-night coastline and slowed by viscous drag. The interfering waves create an uneven dayside temperature pattern that shifts over time, allowing hotspots to drift east or west of the point directly facing the star. Thermal emissions then fluctuate irregularly, with spikes that differ from one orbit to the next and even within a single orbit. This matters because it implies that lava worlds may not show steady, predictable heat signatures but instead complex, time-dependent ones that affect how we interpret observations.

Core claim

Lava tidal waves slosh across deep magma oceans created by tidal heating on short-period rocky planets. The multi-modal ocean response to tidal forcing, with a coastline fixed at the day-night terminator and a parameterized viscous drag, produces wave interference that yields a spatially irregular and highly time-variable dayside heat map. Hotspots wander both east and west of the substellar point, and thermal light curves vary and spike aperiodically from orbit to orbit and within an orbit. Heat deposited by tides reaches steady state through fluid, mushy, and solid-state convection, and for Earth-sized planets with sub-day periods the entire mantle may be tidally liquified.

What carries the argument

The multi-modal response of the magma ocean to tidal forcing under a fixed day-night coastline and parameterized viscous drag, which drives the wave interference responsible for the irregular heat map.

If this is right

Hotspots migrate both east and west of the substellar point on the dayside.
Thermal light curves exhibit aperiodic variations and spikes both within individual orbits and across successive orbits.
For Earth-sized planets with orbital periods shorter than one day, tidal heating can liquify the entire mantle.
Steady-state removal of tidal heat occurs through a combination of fluid, mushy, and solid-state convection in the mantle.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Phase-curve or eclipse observations of lava worlds may require time-resolved modeling rather than equilibrium assumptions to extract surface maps.
Similar wave-driven variability could appear on other bodies with global magma oceans if tidal forcing is present.
The predicted aperiodic spikes offer a potential observational signature to distinguish tidal heating from stellar irradiation alone.
Future missions targeting ultra-short-period planets could test the model by searching for non-repeating thermal patterns over multiple orbits.

Load-bearing premise

The magma ocean response can be captured by a simple parameterized viscous drag together with a coastline fixed exactly at the day-night terminator.

What would settle it

Repeated thermal observations of a known lava world that show no east-west hotspot migration and no orbit-to-orbit or intra-orbit spikes in emission beyond steady-state predictions would falsify the wave-interference mechanism.

read the original abstract

Lava worlds are rocky planets with dayside skins made molten by stellar irradiation. Tidal heating on these shortest-period planets is more than skin deep. We show how orbital eccentricities of just a few percent (within current observed bounds and maintained secularly by exterior companions) can create deep magma oceans. ``Lava tidal waves'' slosh across these oceans; we compute the multi-modal response of the ocean to tidal forcing, subject to a coastline at the day-night terminator and a parameterized viscous drag. Wave interference produces a dayside heat map that is spatially irregular and highly time-variable; hotspots can wander both east and west of the substellar point, and thermal light curves can vary and spike aperiodically, from orbit to orbit and within an orbit. Heat deposited by tides is removed in steady state by a combination of fluid, mushy, and solid-state convection in the mantle. For Earth-sized planets with sub-day periods, the entire mantle may be tidally liquified.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper computes multi-modal tidal wave interference in a magma ocean with a day-night coastline and predicts wandering hotspots plus aperiodic thermal variability, but the results rest on a fixed terminator boundary and constant drag that real rheology may not support.

read the letter

The main thing here is the explicit calculation of how tidal waves slosh in a magma ocean bounded by the day-night line, producing interfering modes that make dayside heating irregular and time-variable. Hotspots wander east and west of the substellar point, and the model yields light curves that spike and change from orbit to orbit. They also tie modest eccentricities to full mantle melting for the shortest-period cases. This is a straightforward extension of standard tidal theory into a new setting, and the concrete predictions for observable variability are the useful output. The work is aimed squarely at people who build interior and atmosphere models for lava worlds, and it gives them a mechanism they can test against data. The soft spot is exactly the one flagged in the stress-test note. The coastline is held fixed at the terminator and drag is a single free parameter. Magma viscosity changes by orders of magnitude with temperature and melt fraction, and the molten region itself should adjust with local energy balance rather than stay pinned to a geometric boundary. Replacing either assumption with something more self-consistent could shift the modal amplitudes and phases enough to remove the reported wandering and aperiodicity. The abstract supplies no equations or validation runs, so it is hard to judge how sensitive the numerics actually are. This is worth sending to referees. The central claim is a plausible mechanism under the stated assumptions, the math is standard, and the implications are direct for the subfield. Reviewers can check the derivations and test the boundary conditions without needing to overhaul the whole approach.

Referee Report

2 major / 2 minor

Summary. The paper claims that orbital eccentricities of a few percent can drive deep magma oceans on short-period lava worlds via tidal heating. The multi-modal response of these oceans to tidal forcing, subject to a fixed coastline at the day-night terminator and a parameterized viscous drag, produces wave interference that yields spatially irregular, highly time-variable dayside heat maps. Hotspots wander east and west of the substellar point, generating aperiodic spikes in thermal light curves both within and across orbits. For Earth-sized planets with sub-day periods, the entire mantle may be tidally liquified, with heat removed by a combination of fluid, mushy, and solid-state convection.

Significance. If the central variability result is robust, the work identifies an intrinsic, tide-driven mechanism for aperiodic thermal fluctuations on lava worlds that could affect interpretation of their observed light curves and phase curves. It also supplies a pathway to complete mantle melting on the shortest-period rocky planets, with implications for their long-term thermal evolution and interior structure.

major comments (2)

[Modeling approach (abstract and §3)] The central claim of wandering hotspots and aperiodic light-curve variability is load-bearing on the choice of a rigid, fixed coastline exactly at the day-night terminator together with a spatially uniform, constant viscous drag coefficient. The manuscript provides no sensitivity tests replacing either assumption with a temperature- or strain-rate-dependent rheology; if the drag term or boundary condition is altered, the modal amplitudes and phase relations that produce the reported interference pattern can change qualitatively (see skeptic note on rheology).
[Numerical results (§4)] The shallow-water or weakly nonlinear formulation is used to compute the ocean response, yet no validation against known analytic limits (e.g., equilibrium tide or single-mode solutions) or comparison to existing magma-ocean models is shown. Without such checks, it is unclear whether the reported east-west wandering and orbit-to-orbit spikes are physical or artifacts of the linear drag parameterization.

minor comments (2)

[Abstract] The abstract states the modeling approach but supplies no governing equations, numerical resolution, or validation metrics; adding a concise equation block or table of parameters would improve accessibility.
[Figures] Figure captions for the heat maps and light curves should explicitly state the orbital period, eccentricity, and drag coefficient values used, as well as the time interval over which variability is shown.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive comments, which have helped clarify the robustness of our results. We address each major point below and have revised the manuscript accordingly to incorporate additional tests and validations.

read point-by-point responses

Referee: [Modeling approach (abstract and §3)] The central claim of wandering hotspots and aperiodic light-curve variability is load-bearing on the choice of a rigid, fixed coastline exactly at the day-night terminator together with a spatially uniform, constant viscous drag coefficient. The manuscript provides no sensitivity tests replacing either assumption with a temperature- or strain-rate-dependent rheology; if the drag term or boundary condition is altered, the modal amplitudes and phase relations that produce the reported interference pattern can change qualitatively (see skeptic note on rheology).

Authors: We agree that the fixed coastline and constant drag are simplifying assumptions central to the model. The coastline location is physically motivated by the sharp day-night temperature contrast that solidifies the nightside, consistent with prior lava-world studies. To address sensitivity, we have added new simulations in revised §4.2 varying the dimensionless drag coefficient over two orders of magnitude (10^{-3} to 10^{-1}). The east-west hotspot wandering and aperiodic light-curve spikes remain qualitatively robust, although quantitative amplitudes shift modestly. We have included these results as a new Figure 6 and expanded the caveats on rheology in §5. A fully temperature- or strain-rate-dependent rheology would require a different numerical framework beyond the present shallow-water approach and is noted as a limitation for future work. revision: yes
Referee: [Numerical results (§4)] The shallow-water or weakly nonlinear formulation is used to compute the ocean response, yet no validation against known analytic limits (e.g., equilibrium tide or single-mode solutions) or comparison to existing magma-ocean models is shown. Without such checks, it is unclear whether the reported east-west wandering and orbit-to-orbit spikes are physical or artifacts of the linear drag parameterization.

Authors: We have added the requested validations in the revised manuscript. Appendix A now demonstrates that the numerical scheme recovers the equilibrium tide solution in the zero-drag, small-eccentricity limit and matches analytic single-mode solutions for a closed basin. We have also included a direct comparison in §4.1 to time-averaged heat fluxes from existing magma-ocean models (e.g., those using parameterized convection), finding order-of-magnitude consistency. These checks confirm that the multi-modal interference and resulting variability are physical features of the tidal response rather than artifacts of the linear drag term. revision: yes

Circularity Check

0 steps flagged

No circularity: forward computation from standard tidal shallow-water model with explicit assumptions

full rationale

The paper derives the dayside heat map and light-curve variability by solving the multi-modal response of a shallow-water ocean to tidal forcing, subject to an imposed coastline at the terminator and a constant viscous drag coefficient. This is a direct numerical or analytic integration of the forced wave equation under those boundary conditions; the resulting interference pattern, east-west hotspot wandering, and orbit-to-orbit spikes are outputs of the integration rather than re-statements of the inputs. No self-citation chain, uniqueness theorem, fitted parameter renamed as prediction, or ansatz smuggled via prior work is present in the derivation chain. The central claim therefore remains independent of its own assumptions and is self-contained against external tidal theory benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claim rests on a parameterized viscous drag whose value is not derived from first principles and on two domain assumptions (eccentricity maintained by companions, coastline at terminator) whose validity is not independently demonstrated in the abstract.

free parameters (1)

viscous drag parameter
Used to close the ocean wave response model; value not specified or derived in abstract.

axioms (2)

domain assumption Orbital eccentricities of a few percent are maintained secularly by exterior companions
Invoked to keep the tidal forcing active within observed bounds.
domain assumption Magma ocean has a coastline at the day-night terminator
Defines the boundary condition for wave reflection.

pith-pipeline@v0.9.0 · 5461 in / 1386 out tokens · 52023 ms · 2026-05-16T14:58:07.105810+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We utilize the classical Laplace tidal equations (LTEs) in the strongly viscous, creep flow limit... subject to a coastline at the day-night terminator and a parameterized viscous drag.
IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

Wave interference produces a dayside heat map that is spatially irregular and highly time-variable

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.