REVIEW 2 major objections 5 minor 48 references
IceCube neutrinos map Earth's density shells and recover its mass and moment of inertia via weak interactions alone.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · grok-4.5
2026-07-12 08:01 UTC pith:4RH5OXHB
load-bearing objection Solid order-of-magnitude upgrade of IceCube Earth tomography; mass and MOI posteriors are new, consistent with PREM/gravity, and the analysis is carefully done. the 2 major comments →
High-Energy Neutrino Tomography of the Earth's Interior with IceCube
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
From the posterior densities of five (and eight) concentric shells fitted to the zenith- and energy-dependent suppression of upgoing atmospheric and astrophysical neutrinos, the Earth's mass and polar moment of inertia are recovered with the highest precision yet achieved through weak interactions; both quantities are statistically consistent with PREM and gravitational reference values.
What carries the argument
Energy- and zenith-dependent neutrino attenuation through a concentric uniform-density shell model, jointly fitted with 37 nuisance parameters via MCMC; the resulting density samples are integrated to give mass and moment of inertia under spherical symmetry.
Load-bearing premise
The Standard-Model neutrino-nucleon cross section is assumed known to a few percent, so residual energy-dependent biases are absorbed into the shell densities rather than treated as free physics.
What would settle it
A future higher-statistics measurement (or an independent cross-section determination) that drives the neutrino-derived mass or moment of inertia outside the present 95 percent highest-posterior-density intervals while PREM remains fixed would falsify the claim of consistency.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports a measurement of the Earth’s radial density profile from the energy- and zenith-dependent attenuation of upgoing muon neutrinos in 10.7 years of IceCube data (500 GeV–100 TeV). The Earth is parametrized as five (or eight) concentric uniform-density shells whose densities are free parameters with flat nonnegative priors; 37 nuisance parameters cover flux, CSMS NLO cross-section, detector, and ice systematics. MCMC sampling of the joint posterior yields shell-density constraints that are statistically consistent with PREM-averaged values. From the same posterior samples the authors derive the Earth’s mass (mode 7.25 imes10^24 kg, 68 % HPD [6.31, 8.31] imes10^24 kg) and polar moment of inertia (mode 1.05 imes10^38 kg m^{2}, 68 % HPD [9.04 imes10^37, 1.24 imes10^38] kg m^{2}), both of which contain the gravitational/PREM reference values inside the 95 % HPD. Blind staged unblinding, a posterior-predictive p-value of 16 %, and an explicit comparison of fits with and without near-horizontal events are used to support the claim that these are the most precise weak-interaction determinations of the two global quantities to date.
Significance. If the result holds, the work supplies an independent, weak-interaction probe of planetary bulk properties that is complementary to gravity and seismology. The order-of-magnitude increase in statistics relative to Donini et al. (2019), the blind analysis protocol, the documented MCMC convergence criteria, and the posterior-predictive goodness-of-fit test constitute a clear methodological advance. The derived mass and moment-of-inertia posteriors are falsifiable against future gravitational or seismic determinations and against next-generation neutrino-telescope data. The measurement therefore strengthens the case that high-energy neutrino attenuation can become a precision geophysical tool once detector and cross-section systematics improve.
major comments (2)
- Methods and §I: the claim that nuclear shadowing and non-isoscalarity are “adequately described” by the two attenuation nuisance parameters (νAtt, ν̄Att) rests on a single EPPS21 study whose quantitative impact on the final mass and MOI posteriors is not shown. A short table or figure quantifying the shift in the mass and MOI modes when the nuclear-PDF correction is turned on (or when a non-isoscalar target is used) would make the residual systematic fully transparent and would strengthen the central consistency claim.
- Fig. 3 and Table I: the posterior modes for mass and MOI lie ~20–30 % above the gravitational/PREM references, with only 6–9 % of the posterior mass below those values. While the 95 % HPD still covers the references, the paper should state more explicitly whether this mild upward bias is expected from the violin ensembles of PREM-based pseudo-data or whether it could indicate a residual energy-dependent cross-section or flux mismodelling that is partially absorbed into the outer-shell densities.
minor comments (5)
- Fig. 1 caption: the phrase “exaggerated by a factor of 10 for illustrative purposes” should specify whether the exaggeration is applied to the rate difference or to the absolute rate; the present wording is ambiguous.
- Table I: the PREM densities used for normalization are given only in parentheses in the row labels; a separate column would improve readability.
- Supplementary Fig. 9: the caption states that mass or MOI is fixed to the reference values, but the numerical values used are not restated; adding them would make the figure self-contained.
- Eq. (1): the structure-function definitions in Eq. (2) are given only for neutrinos at leading order; a brief note that the NLO CSMS calculation used for the actual attenuation includes the corresponding antineutrino and higher-order terms would avoid confusion.
- Methods: the statement that tau regeneration is “negligibly small compared to uncertainties on the inferred Earth densities” would be more convincing if a one-sentence quantitative bound (e.g., “<2 % on column depth for E < 100 TeV”) were supplied.
Circularity Check
No significant circularity: shell densities are free parameters with flat priors; mass and MOI are derived posteriors, not inputs; PREM is used only for comparison and pseudo-data ensembles.
full rationale
The central claim is an inference of Earth's radial density profile from the energy- and zenith-dependent attenuation of IceCube upgoing muon tracks, followed by derivation of mass and polar moment of inertia from those density posteriors. The paper parametrizes the Earth as concentric uniform-density shells with nonnegative uniform priors, jointly samples them with 37 nuisance parameters (flux, CSMS cross-section scalings, ice/DOM systematics) via MCMC, and reports posterior modes and HPD intervals. PREM densities appear only as a comparison curve, as the mean of the statistical-ensemble violins generated from PREM-averaged shells, and as a normalization reference in Table I; they are not imposed as constraints on the primary fit (the constrained-mass/MOI exercise of Fig. 9 is presented separately as a diagnostic). Mass and MOI are computed after the fact from the free density samples under spherical symmetry and are compared to independent gravitational/PREM values; the gravitational mass and PREM MOI lie inside the 95 % HPDs. The CSMS NLO cross section is an external Standard-Model input whose residual nuclear effects are shown (via EPPS21 tests) to be absorbed by the two attenuation nuisance parameters; this is ordinary model dependence, not circularity. Self-citations to earlier IceCube flux, cross-section, and reconstruction papers supply the data sample and nuisance model but do not force the density result. Blind staging, posterior-predictive p = 16 %, and the order-of-magnitude statistics gain over Donini et al. further confirm that the result is data-driven rather than definitional. No step reduces a claimed prediction to a fitted input or to a self-citation uniqueness theorem by construction.
Axiom & Free-Parameter Ledger
free parameters (3)
- five (or eight) shell densities
- 37 nuisance parameters (Norm, DOMeff, Ice A0-4, Ice Phs1-4, Hole Ice, DaemonFlux hadronic and GSF parameters, astrophysi
- astrophysical flux normalization and spectral indices
axioms (4)
- domain assumption CSMS NLO neutrino-nucleon cross section (HERA PDFs) is accurate to a few percent in the 0.5-100 TeV range
- ad hoc to paper Earth density can be approximated by concentric uniform-density spherical shells
- domain assumption Atmospheric and astrophysical neutrino fluxes are adequately described by DaemonFlux + broken power law once nuisance parameters are floated
- domain assumption Fixed nuclear composition (isoscalar nucleons) for converting mass density to nucleon number density
read the original abstract
The Earth's interior reflects its geological evolution, from accretion to present-day dynamics. Its structure drives the geodynamo in the outer core, generating the magnetic field that shields the surface from charged cosmic radiation. The primary observables of the Earth's interior are its radial density distribution and derived quantities such as its mass and moment of inertia. These have traditionally been inferred from gravity and seismic wave propagation, which probe the macroscopic response of matter to gravitational and elastic forces. Here we instead constrain the Earth's density profile using high-energy neutrinos observed by the IceCube Neutrino Observatory at the South Pole. We analyze 10.7 years of predominantly muon-neutrino data spanning 500 GeV--100 TeV, including atmospheric neutrinos produced by cosmic-ray interactions in the Earth's atmosphere and the diffuse astrophysical neutrino flux. Neutrino attenuation depends on both the traversed column density and neutrino energy. By measuring the zenith- and energy-dependent flux suppression, we infer the Earth's radial density profile by fitting a concentric uniform-density shell model that incorporates neutrino fluxes, interaction cross sections, detector response, and glacial-ice systematic uncertainties. From the resulting density posteriors, we derive the Earth's mass and polar moment of inertia as measured by neutrinos. These are the most precise weak-interaction measurements of these quantities to date and are consistent with the Preliminary Reference Earth Model and independent gravitational determinations. Our results demonstrate that neutrinos provide a novel probe of planetary interiors via a distinct physical interaction, complementing gravity and seismology. With improved detectors and precision, neutrinos will further contribute to a multifaceted understanding of the Earth's structure.
Figures
Reference graph
Works this paper leans on
-
[1]
E. D. Williamson and L. H. Adams, Journal of the Washington Academy of Sciences13, 413 (1923)
1923
-
[2]
S. P. Clark Jr and A. E. Ringwood, Reviews of Geophysics2, 35 (1964)
1964
-
[3]
A. M. Dziewonski and D. L. Anderson, Physics of the Earth and Planetary Interiors25, 297 (1981)
1981
-
[4]
Gilbert and A
F. Gilbert and A. M. Dziewonski, Philosophical Transactions of the Royal Society of London, Series A: Mathematical and Physical Sciences278, 187 (1975)
1975
-
[5]
B. L. N. Kennett, E. R. Engdahl, and R. Buland, Geophysical Journal International122, 108 (1995)
1995
-
[6]
A. Garcia, R. Gauld, A. Heijboer, and J. Rojo, JCAP09, 025, arXiv:2004.04756 [hep-ph]
Pith/arXiv arXiv 2004
-
[7]
M. H. Reno, Ann. Rev. Nucl. Part. Sci.73, 181 (2023)
2023
-
[8]
M. C. Gonzalez-Garcia, F. Halzen, M. Maltoni, and H. K. M. Tanaka, Phys. Rev. Lett.100, 061802 (2008), arXiv:0711.0745 [hep-ph]
Pith/arXiv arXiv 2008
-
[9]
A. C. Vincent, C. A. Argüelles, and A. Kheirandish, JCAP11, 012, arXiv:1706.09895 [hep-ph]
-
[10]
R. Abbasiet al.(IceCube), Phys. Rev. D104, 022001 (2021), arXiv:2011.03560 [hep-ex]
arXiv 2021
-
[11]
M. G. Aartsenet al.(IceCube), JINST12(03), P03012, [Erratum: JINST 19, E05001 (2024)], arXiv:1612.05093 [astro- ph.IM]
Pith/arXiv arXiv 2024
-
[12]
Placci and E
A. Placci and E. Zavattini, CERN Internal Report (1973)
1973
-
[13]
L. V. Volkova and G. T. Zatsepin, Izv. Akad. Nauk Ser. Fiz.38N5, 1060 (1974)
1974
-
[14]
M. M. Reynoso and O. A. Sampayo, Astropart. Phys.21, 315 (2004), arXiv:hep-ph/0401102
Pith/arXiv arXiv 2004
-
[15]
Hoshina and H
K. Hoshina and H. K. M. Tanaka, inEGU General Assembly Conference Abstracts, EGU General Assembly Conference Abstracts (2012) p. 3246
2012
-
[16]
A. Donini, S. Palomares-Ruiz, and J. Salvado, Nature Phys.15, 37 (2019), arXiv:1803.05901 [hep-ph]
Pith/arXiv arXiv 2019
-
[17]
R. Abbasiet al.(IceCube), Nucl. Instrum. Meth. A618, 139 (2010), arXiv:1002.2442 [astro-ph.IM]
Pith/arXiv arXiv 2010
-
[18]
Navaset al.(Particle Data Group), Phys
S. Navaset al.(Particle Data Group), Phys. Rev. D110, 030001 (2024)
2024
-
[19]
Devenish and A
R. Devenish and A. Cooper-Sarkar,Deep inelastic scattering(Oxford University Press, 2004)
2004
-
[20]
Y. L. Dokshitzer, Sov. Phys. JETP46, 641 (1977)
1977
-
[21]
V. N. Gribov and L. N. Lipatov, Sov. J. Nucl. Phys.15, 438 (1972)
1972
-
[22]
V. N. Gribov and L. N. Lipatov, Sov. J. Nucl. Phys.15, 675 (1972)
1972
-
[23]
Altarelli and G
G. Altarelli and G. Parisi, Nucl. Phys. B126, 298 (1977)
1977
-
[24]
F. D. Aaronet al.(H1, ZEUS), JHEP01, 109, arXiv:0911.0884 [hep-ex]
-
[25]
Abramowiczet al.(H1, ZEUS), Eur
H. Abramowiczet al.(H1, ZEUS), Eur. Phys. J. C75, 580 (2015), arXiv:1506.06042 [hep-ex]
Pith/arXiv arXiv 2015
-
[26]
A. Cooper-Sarkar, P. Mertsch, and S. Sarkar, JHEP08, 042, arXiv:1106.3723 [hep-ph]
-
[27]
R. Abbasiet al.(IceCube), Phys. Rev. D110, 092009 (2024), arXiv:2405.08077 [hep-ex]
arXiv 2024
-
[28]
R. Abbasiet al.(IceCube), Phys. Rev. Lett.133, 201804 (2024), arXiv:2405.08070 [hep-ex]
arXiv 2024
-
[29]
Goodman and J
J. Goodman and J. Weare, Commun. Appl. Math. Comput. Sc.5, 65 (2010)
2010
-
[30]
D. Foreman-Mackey, D. W. Hogg, D. Lang, and J. Goodman, PASP125, 306 (2013), arXiv:1202.3665 [astro-ph.IM]
Pith/arXiv arXiv 2013
-
[31]
Luzum, N
B. Luzum, N. Capitaine, A. Fienga, W. Folkner, T. Fukushima, J. Hilton, C. Hohenkerk, G. Krasinsky, G. Petit, E. Pitjeva, M. Soffel, and P. Wallace, Celestial Mechanics and Dynamical Astronomy110, 293 (2011)
2011
-
[32]
Ishihara (IceCube), PoSICRC2019, 1031 (2021), arXiv:1908.09441 [astro-ph.HE]
A. Ishihara (IceCube), PoSICRC2019, 1031 (2021), arXiv:1908.09441 [astro-ph.HE]
Pith/arXiv arXiv 2021
-
[33]
M. G. Aartsenet al.(IceCube-Gen2), J. Phys. G48, 060501 (2021), arXiv:2008.04323 [astro-ph.HE]
arXiv 2021
-
[34]
Adrian-Martinezet al.(KM3NeT), J
S. Adrian-Martinezet al.(KM3NeT), J. Phys. G43, 084001 (2016), arXiv:1601.07459 [astro-ph.IM]
Pith/arXiv arXiv 2016
-
[35]
F. Halzen and D. Saltzberg, Phys. Rev. Lett.81, 4305 (1998), arXiv:hep-ph/9804354
Pith/arXiv arXiv 1998
-
[36]
J. P. Yañez and A. Fedynitch, Phys. Rev. D107, 123037 (2023), arXiv:2303.00022 [hep-ph]
Pith/arXiv arXiv 2023
-
[37]
J. M. Picone, A. E. Hedin, D. P. Drob, and A. C. Aikin, Journal of Geophysical Research (Space Physics)107, 1468 (2002). 12
2002
-
[38]
F. Riehn, R. Engel, A. Fedynitch, T. K. Gaisser, and T. Stanev, Phys. Rev. D102, 063002 (2020), arXiv:1912.03300 [hep-ph]
Pith/arXiv arXiv 2020
-
[39]
S. R. Klein, S. A. Robertson, and R. Vogt, Phys. Rev. C102, 015808 (2020), arXiv:2001.03677 [hep-ph]
Pith/arXiv arXiv 2020
-
[40]
P. L. R. Weigel, J. M. Conrad, and A. Garcia-Soto, Phys. Rev. D111, 043044 (2025), arXiv:2408.05866 [hep-ph]
Pith/arXiv arXiv 2025
-
[41]
K. J. Eskola, P. Paakkinen, H. Paukkunen, and C. A. Salgado, Eur. Phys. J. C82, 413 (2022), arXiv:2112.12462 [hep-ph]
Pith/arXiv arXiv 2022
-
[42]
M. G. Aartsenet al.(IceCube), JCAP10, 048, arXiv:1909.01530 [hep-ex]
arXiv 1909
-
[43]
Abbasiet al.(IceCube), (2025), arXiv:2506.04491 [hep-ex]
R. Abbasiet al.(IceCube), (2025), arXiv:2506.04491 [hep-ex]
arXiv 2025
-
[44]
Gelman, X.-L
A. Gelman, X.-L. Meng, and H. Stern (Institute of Statistical Science, Academia Sinica, 1996) pp. 733–760
1996
-
[45]
M. A. Aceroet al.(NOvA), Phys. Rev. D110, 012005 (2024), arXiv:2311.07835 [hep-ex]
Pith/arXiv arXiv 2024
-
[46]
M. G. Aartsenet al.(IceCube), Phys. Rev. Lett.125, 121104 (2020), arXiv:2001.09520 [astro-ph.HE]
arXiv 2020
-
[47]
R. Abbasiet al.(IceCube), Phys. Rev. D104, 022002 (2021), arXiv:2011.03545 [astro-ph.HE]
arXiv 2021
-
[48]
R. Abbasiet al.(IceCube), Astrophys. J.928, 50 (2022), arXiv:2111.10299 [astro-ph.HE]. 13 IV. METHODS The data selection and reconstruction procedures are described in detail in Refs. [27, 28]; we summarize the key elements relevant for this analysis here. We analyze muon tracks produced in muon-neutrino charged-current interactions. These events are clas...
arXiv 2022
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.