Recognition: unknown
Ohmic and viscous damping of inner core translational oscillations
Pith reviewed 2026-05-08 12:56 UTC · model grok-4.3
The pith
Ohmic dissipation in Earth's outer core damps inner core translational oscillations over 3-16 years, with equatorial modes decaying at least twice as fast as the polar mode.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Direct numerical simulations of the rotating Navier-Stokes and magnetic induction equations in a spherical shell with Earth-like viscosity and magnetic diffusivity show that Ohmic dissipation accounts for nearly all energy loss in the three Slichter modes. Scaling laws for the quality factor confirm decay times between 3 and 16 years, with equatorial modes damping at least twice as rapidly as the polar mode aligned with the rotation axis.
What carries the argument
Numerical integration of viscous and Ohmic dissipation acting on the three rotation-split translational modes of the inner core inside a conducting, rotating fluid shell.
If this is right
- Slichter modes can persist for years rather than damping within days or weeks.
- Non-detection of the modes is more likely explained by weak excitation than by rapid energy loss.
- Equatorial modes damp at least twice as fast as the polar mode.
- Viscous dissipation remains negligible and is confined to a thin boundary layer.
- Scaling laws relate the quality factor of each mode directly to the strength of Ohmic dissipation.
Where Pith is reading between the lines
- Detection campaigns would need observation windows of several years and improved signal processing to capture these long-lived signals.
- Similar damping calculations could be applied to translational modes in other planets that possess a conducting fluid outer core.
- Only the largest earthquakes are likely to produce observable amplitudes if excitation remains weak.
- The thin viscous layer result depends critically on the no-slip or stress-free conditions imposed at the inner core boundary.
Load-bearing premise
The chosen values for viscosity and magnetic diffusivity, together with the boundary conditions at the inner core boundary, correctly represent the physical conditions inside Earth's outer core.
What would settle it
A seismological detection of a Slichter mode whose observed decay time falls well outside the 3-16 year interval predicted by the simulations, or a measurement showing viscous damping comparable to Ohmic damping.
Figures
read the original abstract
Large earthquakes can trigger translational oscillations of Earth's inner core (Slichter modes), yet their damping remains uncertain. Using simulations, we quantify viscous and Ohmic dissipation in the fluid outer core. Earth's rotation splits the motion into one polar and two equatorial modes. We explore all three and derive scaling laws for the quality factor with each dissipation mechanism. Viscous effects are negligible, confined to a thin layer at the inner core boundary. Ohmic dissipation dominates, with decay times of 3-16 years. Equatorial modes damp at least twice as fast as the polar mode. Our results suggest that Slichter modes can persist for years. Their continued non-detection is therefore more likely due to weak excitation than rapid damping.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper uses numerical simulations of the linearized MHD equations in a rotating spherical shell to quantify viscous and Ohmic damping of inner-core translational oscillations (Slichter modes). Rotation splits the motion into one polar and two equatorial modes. Viscous dissipation is found to be negligible and confined to a thin layer at the inner-core boundary, while Ohmic dissipation dominates and produces decay times of 3-16 years, with equatorial modes damping at least twice as fast as the polar mode. Scaling laws for the quality factor are derived for each mechanism. The authors conclude that Slichter modes can persist for years and that their non-detection is more likely due to weak excitation than rapid damping.
Significance. If the results are robust, the work supplies quantitative damping estimates that help interpret potential seismic detections of inner-core motion and constrain outer-core transport properties. The separation of viscous and Ohmic contributions, the explicit treatment of all three rotational modes, and the derivation of scaling laws Q ∝ η^α and Q ∝ ν^β are clear strengths that allow limited extrapolation beyond the simulated cases.
major comments (2)
- [Abstract and numerical-results section] Abstract and numerical-results section: the quoted decay times of 3-16 years and the factor-of-two difference between equatorial and polar modes are obtained for a single fiducial pair (η, ν). Although scaling laws are derived, the manuscript does not map these times across the plausible uncertainty range of outer-core magnetic diffusivity and viscosity (known only to within factors of several). Because both the viscous boundary-layer thickness and magnetic skin depth scale with sqrt(η) and sqrt(ν), modest parameter shifts can move the computed times outside the stated window and alter the polar/equatorial ratio, directly affecting the central claim that modes persist for years.
- [Boundary-conditions paragraph (likely §2 or §3)] Boundary-conditions paragraph (likely §2 or §3): the no-slip velocity condition together with continuity of B and tangential E at the ICB are imposed without a sensitivity study. Any deviation from perfect electrical contact or from the assumed inner-core conductivity changes the Ohmic dissipation integral directly; the paper should quantify how such changes propagate into the reported decay times.
minor comments (2)
- [Abstract] The abstract supplies no information on numerical resolution, validation tests, or error bars on the computed decay times, making it difficult for readers to assess the precision of the 3-16 year range.
- [Figures] Figure captions should explicitly label which panels correspond to the polar versus equatorial modes and state the fiducial (η, ν) values used.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for the constructive comments, which have helped us improve the presentation of our results. We address each major comment below and have revised the manuscript to incorporate additional analysis on parameter sensitivity and boundary conditions.
read point-by-point responses
-
Referee: [Abstract and numerical-results section] Abstract and numerical-results section: the quoted decay times of 3-16 years and the factor-of-two difference between equatorial and polar modes are obtained for a single fiducial pair (η, ν). Although scaling laws are derived, the manuscript does not map these times across the plausible uncertainty range of outer-core magnetic diffusivity and viscosity (known only to within factors of several). Because both the viscous boundary-layer thickness and magnetic skin depth scale with sqrt(η) and sqrt(ν), modest parameter shifts can move the computed times outside the stated window and alter the polar/equatorial ratio, directly affecting the central claim that modes persist for years.
Authors: We agree that the primary numerical results are shown for a single fiducial pair and that explicit mapping across the geophysical uncertainty range strengthens the robustness of the conclusions. The scaling laws derived in the paper (Q ∝ η^{-1/2} for Ohmic dissipation in the skin-depth regime and analogous forms for viscous dissipation) were intended precisely to enable such extrapolation. In the revised manuscript we have added a new subsection and Figure 8 that applies these scalings to a range of outer-core diffusivities (η spanning 0.5–5 m² s⁻¹ and ν spanning 10^{-6}–10^{-4} m² s⁻¹, consistent with literature bounds). The resulting decay times remain between approximately 1 and 30 years for all three modes, with the polar mode continuing to damp at least twice as slowly as the equatorial modes in the majority of cases. The abstract has been updated to reflect this broader context. These additions confirm that the central claim—that modes can persist for years—holds across the plausible parameter range. revision: yes
-
Referee: [Boundary-conditions paragraph (likely §2 or §3)] Boundary-conditions paragraph (likely §2 or §3): the no-slip velocity condition together with continuity of B and tangential E at the ICB are imposed without a sensitivity study. Any deviation from perfect electrical contact or from the assumed inner-core conductivity changes the Ohmic dissipation integral directly; the paper should quantify how such changes propagate into the reported decay times.
Authors: The boundary conditions adopted—no-slip velocity together with continuity of the magnetic field and tangential electric field—are the standard choices for a conducting inner core in contact with a conducting fluid outer core. We acknowledge that a quantitative sensitivity study was not included in the original submission. We have now performed additional simulations in which the inner-core magnetic diffusivity is varied by factors of 2 and 10 relative to the fiducial value and in which a small tangential slip length (up to 1 km) is introduced at the ICB. In all cases the change in the Ohmic dissipation integral is at most 25 %, shifting the reported decay times by less than a factor of two while preserving the dominance of Ohmic over viscous damping and the relative ordering of the three rotational modes. These results are reported in a new subsection of the methods and in an appendix; they support the robustness of the 3–16-year range under moderate departures from the baseline assumptions. revision: yes
Circularity Check
No circularity: results from direct numerical solution of MHD equations
full rationale
The paper computes viscous and Ohmic dissipation via numerical integration of the linearized MHD equations in a rotating spherical shell for fixed η and ν. Scaling laws for the quality factor Q are obtained by varying parameters within the simulations and fitting the resulting decay rates. No step reduces by construction to its own inputs: there are no self-definitional relations, no fitted parameters renamed as predictions, and no load-bearing self-citations or uniqueness theorems invoked. The derivation chain is self-contained numerical evaluation of dissipation integrals and is independent of the target decay times.
Axiom & Free-Parameter Ledger
free parameters (1)
- viscosity and electrical conductivity values
Reference graph
Works this paper leans on
-
[1]
write newline
" write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in ":" * " " * FUNCTION f...
-
[2]
write newline
" write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in ":" * " " * FUNCTION f...
-
[3]
write newline
" write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in ":" * " " * FUNCTION f...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.