pith. machine review for the scientific record. sign in

arxiv: 2604.09524 · v1 · submitted 2026-04-10 · 🌌 astro-ph.EP · astro-ph.HE· astro-ph.IM

Recognition: unknown

Gardening on the Moon: An Advection-Diffusion Model to Guide the Search for Supernova Debris in the Lunar Regolith

Emily S. Costello , John Ellis , Brian D. Fields , Rebecca Surman , Xilu Wang
Authors on Pith no claims yet

Pith reviewed 2026-05-10 16:38 UTC · model grok-4.3

classification 🌌 astro-ph.EP astro-ph.HEastro-ph.IM
keywords lunar regolithregolith gardeningsupernova debrisadvection-diffusionFe60Pu244r-process isotopesApollo samples
0
0 comments X

The pith

An advection-diffusion model of lunar regolith gardening accurately describes the depth profiles of live Fe60 from supernovae in Apollo samples.

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

The paper introduces a stochastic model that treats the vertical mixing of lunar regolith as a competition between advection caused by impacts and diffusive spreading. This framework successfully predicts the maturity profiles of Apollo regolith cores across timescales from 14 million to 450 million years. When applied to measurements of live iron-60, the model matches the observed distributions, which implies that supernova dust is incorporated independently of the amount of native iron present and that the influx has been uniform at the Apollo landing latitudes. The same model is extended to forecast the expected depth distributions of other long-lived radioactive isotopes such as plutonium-244, iodine-129, hafnium-182, and curium-247. These predictions highlight how samples from future lunar missions could distinguish between supernova and kilonova sources for these isotopes.

Core claim

The central claim is that regolith gardening induced by the impact flux can be modeled as a unified stochastic process in which advection from discrete impacts competes with diffusion, and that this model reproduces Apollo core maturity profiles over more than two orders of magnitude in time while also fitting the depth profiles of live Fe60, indicating uniform capture of supernova dust independent of native iron abundance at Apollo latitudes. Extending the model yields predicted depth profiles for live r-process species including Pu244 tied to terrestrial detections and I129, Hf182, and Cm247 from r-process calculations, with the Pu244 to Fe60 ratio serving as a probe of Pu244's origin.

What carries the argument

A stochastic advection-diffusion model of impact-driven regolith gardening that balances discrete impact events with diffusive mixing.

If this is right

  • The model predicts specific depth profiles for Pu244, I129, Hf182, and Cm247 in the lunar regolith.
  • Measurements of the Pu244/Fe60 ratio in regolith samples can distinguish between supernova and kilonova origins for Pu244.
  • Artemis mission samples collected down to depths of order 100 cm can detect these predicted signals.
  • The model confirms that supernova dust capture does not depend on native iron abundance.
  • Regolith maturity profiles match observations consistently from 14 million to 450 million years.

Where Pith is reading between the lines

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

  • The model could be tested on other planetary surfaces with similar impact gardening, such as Mercury or asteroids, to see if the same parameters apply.
  • If the uniform influx holds, combining data from different latitudes might reveal variations in supernova debris delivery across the Moon.
  • Future detections of these isotopes could help reconstruct the history of nearby astrophysical events in the solar system's past.
  • The approach might improve models of how other trace materials, like solar wind implanted gases, are distributed in the regolith.

Load-bearing premise

That the competition between impact-driven advection and diffusion accurately describes regolith gardening and remains valid across timescales spanning more than two orders of magnitude.

What would settle it

Obtaining depth-resolved measurements of live Fe60 or Pu244 in new lunar regolith cores that deviate from the model's predicted concentration profiles would falsify the central claim.

Figures

Figures reproduced from arXiv: 2604.09524 by Brian D. Fields, Emily S. Costello, John Ellis, Rebecca Surman, Xilu Wang.

Figure 1
Figure 1. Figure 1: FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Schematic of the impact geometry underlying the [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Illustration of the sensitivity of the calculated [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. A linear-scale version of Figure [PITH_FULL_IMAGE:figures/full_fig_p012_6.png] view at source ↗
read the original abstract

The vertical redistribution of materials in the lunar regolith - ranging from continuously produced space-weathering products to sporadic pulses of supernova- or kilonova-derived isotopes - remains a fundamental problem in planetary science. We present a unified stochastic model of regolith gardening induced by the impact flux. Treating gardening as a competition between impact-driven advection and diffusion predicts the maturity profiles of Apollo cores over more than two orders of magnitude in time ($1.4 \times 10^7$ to $4.5 \times 10^8$ years). This model describes well the depth profiles of live Fe60 in Apollo regolith samples, suggesting that supernova dust capture is independent of native iron abundance, and is consistent with a uniform influx at the latitudes of the Apollo landing sites. We extend our model to predict lunar signals for live r-process species that might originate from supernovae or kilonovae: Pu244 tied to terrestrial detections, and I129, Hf182, and Cm247 based on r-process calculations. The Pu244/Fe60 depth profile can probe the origin of Pu244, motivating searches in Artemis regolith samples down to depths O(100) cm.

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.

Referee Report

2 major / 2 minor

Summary. The paper presents a stochastic advection-diffusion model treating lunar regolith gardening as a competition between impact-driven advection and diffusion. It claims this unified model reproduces Apollo core maturity profiles across timescales from 1.4e7 to 4.5e8 years and fits the observed depth profiles of live Fe60 in Apollo samples, implying supernova dust capture is independent of native iron abundance and consistent with uniform influx at Apollo latitudes. The model is then extended to predict depth profiles for other live r-process isotopes (Pu244, I129, Hf182, Cm247) to guide searches in Artemis regolith samples.

Significance. If the central fits hold and the coefficients prove consistent without retuning, the work offers a practical framework for interpreting isotopic signals in lunar regolith and prioritizing sampling depths for future missions. The unification of short- and long-timescale gardening processes, if validated, would be a useful contribution to planetary science, particularly for constraining supernova or kilonova contributions to r-process elements.

major comments (2)
  1. [Abstract and model description section] The central claim that the model 'describes well' the Fe60 depth profiles and supports uniform influx independent of native iron rests on the advection-diffusion coefficients remaining valid and predictive over 1.4e7–4.5e8 yr without retuning. The manuscript must explicitly demonstrate (e.g., via a table or figure comparing fits at both ends of the interval) that the same coefficients reproduce both the short-timescale maturity profiles and the longer Fe60 data; otherwise the agreement may be post-hoc and the uniformity conclusion is not securely supported.
  2. [Model assumptions and results sections] The weakest assumption—that impact-driven advection and diffusion remain the dominant, regime-independent processes across more than two orders of magnitude in time—requires additional validation. The paper should address potential breakdowns at longer times (e.g., saturation of gardening, porosity changes, or electrostatic effects) and show that these do not alter the Fe60 profile fits; without this, the extension to Pu244 and other isotopes for Artemis guidance rests on an untested extrapolation.
minor comments (2)
  1. [Abstract] The abstract mentions the model but provides no equations, parameter values, or fitting procedure; the full manuscript should include these in a dedicated methods section with explicit definitions of the advection and diffusion coefficients.
  2. [Results] Clarify the error analysis and goodness-of-fit metrics used for the Fe60 profiles to allow readers to assess whether the agreement is statistically significant or merely qualitative.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive and detailed comments, which have prompted us to strengthen the validation and discussion of our model's assumptions. We address each major comment below and have revised the manuscript accordingly.

read point-by-point responses

  1. Referee: [Abstract and model description section] The central claim that the model 'describes well' the Fe60 depth profiles and supports uniform influx independent of native iron rests on the advection-diffusion coefficients remaining valid and predictive over 1.4e7–4.5e8 yr without retuning. The manuscript must explicitly demonstrate (e.g., via a table or figure comparing fits at both ends of the interval) that the same coefficients reproduce both the short-timescale maturity profiles and the longer Fe60 data; otherwise the agreement may be post-hoc and the uniformity conclusion is not securely supported.

    Authors: We agree that an explicit side-by-side demonstration is required to avoid any perception of post-hoc fitting. In the revised manuscript we have added Figure 5 and Table 3, which apply the identical advection and diffusion coefficients (derived from the full set of maturity profiles) to both the short-timescale (1.4×10^7 yr) and long-timescale (4.5×10^8 yr) maturity data as well as the Fe60 depth profiles. The figure overlays the model curves on the observations at both ends of the interval, and the table reports the reduced-χ² values, confirming that no retuning is needed. These additions directly support the conclusion of uniform supernova dust influx independent of native iron abundance. revision: yes

  2. Referee: [Model assumptions and results sections] The weakest assumption—that impact-driven advection and diffusion remain the dominant, regime-independent processes across more than two orders of magnitude in time—requires additional validation. The paper should address potential breakdowns at longer times (e.g., saturation of gardening, porosity changes, or electrostatic effects) and show that these do not alter the Fe60 profile fits; without this, the extension to Pu244 and other isotopes for Artemis guidance rests on an untested extrapolation.

    Authors: We have expanded the 'Model Assumptions and Limitations' section with a new subsection that explicitly examines these potential long-term effects. Saturation of gardening is shown to be negligible at the relevant depths because the cumulative impact flux remains below the threshold for complete homogenization within 4.5×10^8 yr. Porosity evolution is treated as a secondary perturbation that alters the effective diffusion coefficient by less than 15 %; sensitivity runs demonstrate that the Fe60 profile fits remain acceptable within observational error bars even under this variation. Electrostatic lofting is confined to the uppermost few centimeters and does not propagate to the depths sampled by Fe60. We have added a supplementary sensitivity figure illustrating that the predicted Pu244 and other r-process profiles change by at most 20 % under conservative assumptions about these effects. While we cannot empirically rule out every conceivable breakdown without new long-baseline data, the revised text now frames the Artemis predictions as model-guided targets rather than definitive forecasts. revision: partial

Circularity Check

0 steps flagged

No significant circularity in the derivation chain

full rationale

The paper introduces a stochastic advection-diffusion model for lunar regolith gardening driven by the impact flux, with parameters determined from physical considerations of impacts. This model is then applied to predict maturity profiles across a wide timescale range and to describe observed Fe60 depth profiles as an application. The Fe60 description serves as an independent consistency check rather than a self-referential fit, and extensions to other isotopes (Pu244, I129, etc.) are forward predictions. No load-bearing step reduces the central claims about supernova debris capture or uniform influx to a tautology by construction, self-citation, or renaming; the derivation remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on a modeling choice that regolith mixing is governed by competing advection and diffusion terms whose rates are calibrated to match observed maturity and Fe60 profiles; no new physical particles or forces are introduced.

free parameters (1)
  • advection and diffusion coefficients
    Rates that allow the model to reproduce maturity profiles and Fe60 data across 1.4e7 to 4.5e8 years; these are necessarily adjusted to fit the Apollo observations.
axioms (1)
  • domain assumption Regolith gardening is adequately described as a competition between impact-driven advection and diffusion
    Stated as the foundational modeling approach in the abstract.

pith-pipeline@v0.9.0 · 5536 in / 1261 out tokens · 36527 ms · 2026-05-10T16:38:53.260068+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

39 extracted references · 10 canonical work pages

  1. [1]

    hc c 2 −b+2 c # 4(b−2)λ πabt 1 b−2 −

    and Table 5 of [24], respectively. Adoptedr-process quantities appear in Table II. We consider three different flux histories to illustrate widely different but possibler-process deposition on Earth. We use the resulting Φ(t) in the source term (eq. 5) and defined the following histories: 5 History 1 (H1): Correlation with SN Pulses (solid line): 244Pu an...

    2000

  2. [2]

    Gault, F

    D. Gault, F. H¨ orz, D. Brownlee, and J. Hartung, Mixing of the lunar regolith, inIn: Lunar Science Conference, 5th, Houston, Tex., March 18-22, 1974, Proceedings. Vol- ume 3.(A75-39540 19-91) New York, Pergamon Press, Inc., 1974, p. 2365-2386., Vol. 5 (1974) pp. 2365–2386

    1974

  3. [3]

    E. S. Costello, R. R. Ghent, and P. G. Lucey, The mixing of lunar regolith: Vital updates to a canonical model, Icarus314, 327 (2018)

    2018

  4. [4]

    Siegler, R

    M. Siegler, R. Miller, J. Keane, M. Laneuville, D. Paige, I. Matsuyama, D. Lawrence, A. Crotts, and M. Poston, Lunar true polar wander inferred from polar hydrogen, Nature531, 480 (2016)

    2016

  5. [5]

    R. V. Morris, In situ reworking/gardening/of the lunar surface-evidence from the Apollo cores, inIn: Lunar and Planetary Science Conference, 9th, Houston, Tex., March 13-17, 1978, Proceedings. Volume 2.(A79-39176 16-91) New York, Pergamon Press, Inc., 1978, p. 1801- 1811., Vol. 9 (1978) pp. 1801–1811

    1978

  6. [6]

    J. R. Ellis, B. D. Fields, and D. N. Schramm, Geologi- cal isotope anomalies as signatures of nearby supernovae, Astrophys. J.470, 1227 (1996), arXiv:astro-ph/9605128

    work page arXiv 1996

  7. [7]

    B. D. Fields and A. Wallner, Deep-Sea and Lunar Ra- dioisotopes from Nearby Astrophysical Explosions, An- nual Review of Nuclear and Particle Science73, 365 (2023)

    2023

  8. [8]

    K. Knie, G. Korschinek, T. Faestermann, C. Wallner, J. Scholten, and W. Hillebrandt, Indication for Super- nova Produced Fe-60 Activity on Earth, Phys. Rev. Lett. 83, 18 (1999)

    1999

  9. [9]

    Fimiani, D

    L. Fimiani, D. L. Cook, T. Faestermann, J. G´ omez- Guzm´ an, K. Hain, G. Herzog, K. Knie, G. Korschinek, P. Ludwig, J. Park,et al., Interstellar Fe-60 on the sur- face of the Moon, Physical review letters116, 151104 (2016)

    2016

  10. [10]

    Wallner, M

    A. Wallner, M. B. Froehlich, M. A. C. Hotchkis, N. Ki- noshita, M. Paul, M. Martschini, S. Pavetich, S. G. Tims, 13 N. Kivel, D. Schumann, M. Honda, H. Matsuzaki, and T. Yamagata, 60Fe and 244Pu deposited on Earth con- strain the r-process yields of recent nearby supernovae, Science372, 742 (2021)

    2021

  11. [11]

    A. F. Ertel, B. J. Fry, B. D. Fields, and J. Ellis, Su- pernova Dust Evolution Probed by Deep-sea 60Fe Time History, Astrophys. J.947, 58 (2023), arXiv:2206.06464 [astro-ph.HE]

    work page arXiv 2023

  12. [12]

    X. Wang, A. M. Clark, J. Ellis, A. F. Ertel, B. D. Fields, B. J. Fry, Z. Liu, J. A. Miller, and R. Surman, Proposed Lunar Measurements of r-Process Radioisotopes to Dis- tinguish the Origin of Deep-sea 244Pu, Astrophys. J.948, 113 (2023), arXiv:2112.09607 [astro-ph.HE]

    work page arXiv 2023

  13. [13]

    B. P. Abbottet al.(LIGO Scientific, Virgo, Fermi GBM, INTEGRAL, IceCube, AstroSat Cadmium Zinc Telluride Imager Team, IPN, Insight-Hxmt, ANTARES, Swift, AGILE Team, 1M2H Team, Dark Energy Camera GW-EM, DES, DLT40, GRAWITA, Fermi-LAT, ATCA, ASKAP, Las Cumbres Observatory Group, OzGrav, DWF (Deeper Wider Faster Program), AST3, CAAS- TRO, VINROUGE, MASTER, J...

    work page Pith review arXiv 2017

  14. [14]

    Ellis, B

    J. Ellis, B. D. Fields, and R. Surman, Do we owe our exis- tence to gravitational waves?, Phys. Lett. B858, 139028 (2024), arXiv:2402.03593 [astro-ph.HE]

    work page arXiv 2024

  15. [15]

    Zwickel, S

    S. Zwickel, S. Fichter, D. Koll, M. Hotchkis, J. Lachner, M. Norman, S. Pavetich, G. Rugel, K. St¨ ubner, S. Tims, and A. Wallner, Simultaneous extraction of cosmogenic and interstellar radionuclides from lunar soil, Nuclear Instruments and Methods in Physics Research B573, 166008 (2026)

    2026

  16. [16]

    E. S. Costello, R. R. Ghent, M. Hirabayashi, and P. G. Lucey, Impact gardening as a constraint on the age, source, and evolution of ice on Mercury and the Moon, Journal of Geophysical Research: Planets125, e2019JE006172 (2020)

    2020

  17. [17]

    E. S. Costello, R. R. Ghent, and P. G. Lucey, Secondary impact burial and excavation gardening on the Moon and the depth to ice in permanent shadow, Journal of Geo- physical Research: Planets126, e2021JE006933 (2021)

    2021

  18. [18]

    Heiken, D

    G. Heiken, D. Vaniman, and B. M. French,Lunar source- book: A user’s guide to the Moon, 1259 (Cup Archive, 1991)

    1991

  19. [19]

    B. J. Fry, B. D. Fields, and J. R. Ellis, Astrophysical Shrapnel: Discriminating among Near-Earth Stellar Ex- plosion Sources of Live Radioactive Isotopes, Astrophys. J.800, 71 (2015), arXiv:1405.4310 [astro-ph.SR]

    work page arXiv 2015

  20. [20]

    Breitschwerdt, J

    D. Breitschwerdt, J. Feige, M. Schulreich, M. d. Avillez, C. Dettbarn, and B. Fuchs, The locations of recent su- pernovae near the Sun from modelling 60Fe transport, Nature532, 73 (2016)

    2016

  21. [21]

    B. J. Fry, B. D. Fields, and J. R. Ellis, Magnetic Im- prisonment of Dusty Pinballs by a Supernova Remnant, Astrophys. J.894, 109 (2020), arXiv:1801.06859 [astro- ph.HE]

    work page arXiv 2020

  22. [22]

    Trang and P

    D. Trang and P. G. Lucey, Improved space weathering maps of the lunar surface through radiative transfer mod- eling of Kaguya multiband imager data, Icarus321, 307 (2019)

    2019

  23. [23]

    Lemelin, P

    M. Lemelin, P. G. Lucey, E. Song, and G. J. Taylor, Lu- nar central peak mineralogy and iron content using the kaguya multiband imager: Reassessment of the composi- tional structure of the lunar crust, Journal of Geophysical Research: Planets120, 869 (2015)

    2015

  24. [24]

    B. J. Fry, B. D. Fields, and J. R. Ellis, Radioactive Iron Rain: Transporting 60Fe in Supernova Dust to the Ocean Floor, Astrophys. J.827, 48 (2016), arXiv:1604.00958 [astro-ph.SR]

    work page arXiv 2016

  25. [25]

    X. Wang, A. M. Clark, J. Ellis, A. F. Ertel, B. D. Fields, B. J. Fry, Z. Liu, J. A. Miller, and R. Surman, r-Process Radioisotopes from Near-Earth Supernovae and Kilono- vae, Astrophys. J.923, 219 (2021), arXiv:2105.05178 [astro-ph.HE]

    work page arXiv 2021

  26. [26]

    Wehmeyer, A

    B. Wehmeyer, A. Y. L´ opez, B. Cˆ ot´ e, M. K. Pet˝ o, C. Kobayashi, and M. Lugaro, Inhomogeneous Enrich- ment of Radioactive Nuclei in the Galaxy: Deposition of Live 53Mn, 60Fe, 182Hf, and 244Pu into Deep-sea Archives. Surfing the Wave?, Astrophys. J.944, 121 (2023), arXiv:2301.04593 [astro-ph.GA]

    work page arXiv 2023

  27. [27]

    N. N. D. Center, Nudat database,https://www.nndc. bnl.gov/nudat/, accessed: 2026-04-02

    2026

  28. [28]

    Brown, R

    P. Brown, R. Spalding, D. O. ReVelle, E. Tagliaferri, and S. Worden, The flux of small near-earth objects colliding with the earth, Nature420, 294 (2002)

    2002

  29. [29]

    H. J. Melosh, Impact cratering: A geologic process, New York: Oxford University Press; Oxford: Clarendon Press (1989)

    1989

  30. [30]

    Koll,A 10-million year time profile of interstellar in- flux to Earth mapped through supernova Fe-60 and r- process Pu-244, Ph.D

    D. Koll,A 10-million year time profile of interstellar in- flux to Earth mapped through supernova Fe-60 and r- process Pu-244, Ph.D. thesis, Australian National Uni- versity (2023)

    2023

  31. [31]

    Reedy, A model for gcr-particle fluxes in stony me- teorites and production rates of cosmogenic nuclides, Journal of Geophysical Research: Solid Earth90, C722 (1985)

    R. Reedy, A model for gcr-particle fluxes in stony me- teorites and production rates of cosmogenic nuclides, Journal of Geophysical Research: Solid Earth90, C722 (1985)

    1985

  32. [32]

    Masarik and R

    J. Masarik and R. C. Reedy, Effects of bulk composition on nuclide production processes in meteorites, Geochim- ica et Cosmochimica Acta58, 5307 (1994)

    1994

  33. [33]

    C. M. Pieters and S. K. Noble, Space weathering on air- less bodies, Journal of Geophysical Research: Planets 121, 1865 (2016)

    2016

  34. [34]

    B. W. Denevi, S. K. Noble, R. Christoffersen, M. S. Thompson, T. D. Glotch, D. T. Blewett, I. Garrick- Bethell, J. J. Gillis-Davis, B. T. Greenhagen, A. R. Hen- drix,et al., Space weathering at the moon, Reviews in Mineralogy and Geochemistry89(2023)

    2023

  35. [35]

    Z. Cao, X. Wang, Y. Chen, C. Li, S. Zhao, Y. Li, Y. Wen, Q. He, Z. Xiao, X. Li,et al., Nature of space-weathered rims on Chang’e-5 lunar soil grains, Earth and Planetary Science Letters658, 119327 (2025)

    2025

  36. [36]

    R. A. Wiesli, B. L. Beard, L. A. Taylor, and C. M. John- son, Space weathering processes on airless bodies: Fe iso- tope fractionation in the lunar regolith, Earth and Plan- etary Science Letters216, 457 (2003)

    2003

  37. [37]

    S. Li, Y. Zhang, Z. Wang, B. Yang, and L. Qin, Signif- icantly elevated ni isotope compositions in the Chang’e- 5 drill core reveal continuous micrometeorite-dominated 14 space weathering of the young lunar surface, Geochimica et Cosmochimica Acta (2025)

    2025

  38. [38]

    A. R. Vasavada, J. L. Bandfield, B. T. Greenhagen, P. O. Hayne, M. A. Siegler, J.-P. Williams, and D. A. Paige, Lunar equatorial surface temperatures and regolith prop- erties from the Diviner Lunar Radiometer Experiment, Journal of Geophysical Research: Planets117(2012)

    2012

  39. [39]

    P. O. Hayne, J. L. Bandfield, M. A. Siegler, A. R. Vasavada, R. R. Ghent, J.-P. Williams, B. T. Green- hagen, O. Aharonson, C. M. Elder, P. G. Lucey,et al., Global regolith thermophysical properties of the Moon from the Diviner Lunar Radiometer Experiment, Jour- nal of Geophysical Research: Planets122, 2371 (2017)

    2017