Pith. sign in

REVIEW 4 major objections 4 minor 23 references

Investigation of lunar ejecta dynamics: Particles reaching the near-Earth space and their effect on Earth-based observation

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

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

Pith's one-line read Most lunar-ejected dust that impacts Earth arrives quickly—about 70 percent within a year, with the smallest grains arriving in under a week—and carries orbital signatures that separate it from interplanetary dust, giving observers a way…

desk verdict Useful post-processing of an established lunar-ejecta simulation, but the central observational claim ('distinguishable from interplanetary dust') is asserted, not demonstrated. read the letter →

arxiv 2507.14968 v1 pith:MZCG4ZZF submitted 2025-07-20 astro-ph.EP physics.space-ph

Kun Yang , Yu Jiang , Youpeng Liang , Xiaodong Liu This is my paper
classification astro-ph.EPphysics.space-ph
keywords lunarejectadusttorusEarthimpactorsorbitalelementssolarradiationpressurePoynting-RobertsondragzodiacalEarth-basedobservations
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

Most of the lunar dust that eventually hits Earth gets there quickly: about 70 percent of simulated Earth impactors arrive within one year, and the smallest grains (0.2 µm) arrive in under a week 87.2 percent of the time. The paper argues that these incoming grains have distinctive orbital signatures—semi-major axes mostly between 20 and 150 Earth radii, eccentricities peaking near 0.99, and inclinations peaking near 40 and 120 degrees—that separate them from ordinary interplanetary dust. It then projects the simulated particles onto the skies of five Earth-based observatories and finds they form belts and arcs with different orientations, so stars near the ecliptic at low declination are the most likely to be obscured. If the picture holds, it provides a testable map of where and when lunar dust should appear in Earth-based observations and how to tell it apart from other dust.

What carries the argument

The argument is carried by the 216,000-particle dynamical simulation from the paper's preceding study, in which lunar ejecta are launched vertically with power-law size and velocity distributions (exponents 3.7 and 2.2), a bulk density of 3500 kg/$m^{3}$, and a production rate of 0.2 kg/s, then integrated under Earth, Moon, and Sun gravity plus solar radiation pressure and Poynting-Robertson drag until they hit the Moon, hit the Earth, or escape. The central object is the subset of 8,077 Earth impactors extracted from that simulation, which the paper re-analyzes by initial parameter and projects into the topocentric east-north-up frames of five observatories. The selectivity comes from the size- and velocity-dependent dynamics: radiation pressure removes most particles ≤0.5 µm, large grains (≥10 µm) mostly stay in the Earth-Moon torus, and only particles launched at speeds between roughly 1.0 and 1.27 times escape velocity, opposite the Moon's orbital motion, can reach the Earth. That filter, combined with the transfer dynamics, is what produces the observed orbital-element peaks and the rapid arrival times.

What would settle it

A dust detector at the Earth-Moon L1 point or on a near-Earth spacecraft that can measure the orbital elements of individual sub-micron grains would test the predicted bimodal inclination peaks near 40° and 120° and the eccentricity peak near 0.99; observing instead a single broad inclination band or eccentricities far from 0.99 would falsify the distinguishability claim.

Watch

Extended reading notes

Core claim

Using the 8,077 particles from a 216,000-particle simulation of lunar ejecta that are found to strike the Earth, the paper establishes three things. First, the impactors are not a random sample: they favor particle sizes from about 1 to 2 µm, ejection speeds between the lunar escape speed and about 1.27 times that value, and launch directions opposite the Moon's orbital motion. Second, their arrival is fast—most reach Earth within a year, and the sub-micron grains do so within a week—yet their orbital elements at arrival are concentrated: semi-major axes in [20, 150] Earth radii, eccentricities near 0.99, and a bimodal inclination distribution with peaks near 40° (prograde, 54%) and 120° (retrograde, 46%). Third, when these particles are projected onto the sky from five observatories, they take the form of an arc or belt whose position and orientation depend on the observatory's latitude, and the equatorial-coordinate map shows that low-declination stars (Aldebaran, Pollux, Spica, Antares) are most likely to be affected. Together these results support the paper's central claim that lunar-ejected Earth impactors form a distinct, largely time-predictable dust population that can be separated from interplanetary dust.

Load-bearing premise

The whole analysis rests on the assumed initial conditions for lunar ejecta—straight-up launch, steep power-law size and speed distributions, a grain density of 3500 kg/$m^{3}$, and a production rate of 0.2 kg/s—and the paper reports no sensitivity tests, so if any of these misrepresents real lunar ejecta, every reported impactor fraction, arrival time, and sky map inherits the error.

Editorial extensions

If this is right

  • A dust detector capable of determining grain orbits would see a population with semi-major axes under about 150 Earth radii and eccentricities near 0.99, distinct from the typical interplanetary dust background.
  • Transient ejection events on the Moon—a fresh cratering impact—should produce a measurable spike of sub-micron grains at Earth within about a week, giving a direct test of the arrival-time statistics.
  • Surveys that observe low-declination stars such as Aldebaran, Pollux, Spica, and Antares are the most likely to have their images contaminated by lunar-dust scattering, so photometric pipelines for those fields should account for a structured background.
  • The roughly 30 percent of impactors that take more than a year are the main builders of the steady-state Earth-Moon dust torus, so the fast-arrival and slow-torus components are two observable faces of the same lunar ejecta process.
  • The mass input of about 2.3e-4 kg/s from lunar ejecta is a non-negligible part of the near-Earth dust environment, comparable to earlier estimates and worth including in micrometeoroid flux models.

Reading between the lines

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

  • If the vertical-ejection assumption is too simple, the launch-angle preference found here suggests that a more realistic inclined ejecta distribution would likely reduce the absolute Earth-impact fraction, but the orbital-element fingerprint at arrival—set by dynamics after leaving the Moon's Hill sphere—might remain broadly similar, so the distinguishability claim could survive.
  • The sky-map method could be extended to predict how the dust belt shifts with the lunar phase and the observing season, since the torus plane does not coincide with the Earth's equator; such a time-dependent map would let observers schedule around the dust.
  • The same simulation infrastructure could be used to estimate how much lunar-ejected dust accumulates on Earth's upper atmosphere or on space-based detectors in low Earth orbit, connecting the orbital distributions to deposition rates.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

4 major / 4 minor

Summary. This manuscript is a post-processing study of a 216,000-particle simulation of lunar impact ejecta from the authors' previous work (Yang et al. 2022). The authors select the 8,077 simulated particles that hit the Earth and analyze their impactor fractions as functions of size, launch angle, and ejection velocity; their orbital-element distributions at encounter; travel-time distributions; and sky projections toward five Earth-based observatories. The headline results are that about 70% of lunar-ejected Earth impactors arrive within one year, that small particles (0.2 and 0.5 micrometer) mostly arrive within one week, and that a 'large proportion' of these impactors can be distinguished from interplanetary dust on the basis of orbital element differences. The paper also presents optical-depth and number-density maps of the lunar dust torus.

Significance. The strength of the paper is that it extracts, from an already published dynamical simulation, a concrete set of observationally relevant products: particle fractions, orbital-element histograms, and observatory-specific sky maps. The authors are transparent about the underlying model and its parameters (Table 1), and the previous simulation is a reasonable basis for this analysis. However, the headline observational claim of distinguishability from interplanetary dust is not established, because no interplanetary dust population is introduced and no quantitative separation metric is computed. In addition, the orbital-element histograms used for that claim deliberately exclude particles that leave the Earth's Hill sphere, without stating how many are excluded or how the distributions change. If these gaps are filled, the paper would provide a useful reference for assessing lunar dust contamination of ground-based observations; at present it is mostly a descriptive extension of Yang et al. (2022).

major comments (4)
  1. [Abstract and Sect. 5] The central claim that a large proportion of lunar-ejected Earth impactors can be distinguished from interplanetary dust is not supported by any comparison population. Neither the abstract nor Sect. 5 defines an interplanetary dust population, provides its orbital-element distribution, or computes a quantitative overlap, confusion rate, or separation metric. To substantiate the claim, the authors should introduce a concrete IDP model (e.g., based on the Grun et al. 1985 or Dikarev et al. 2005 populations cited in the introduction) and quantify what fraction of lunar particles fall outside the IDP envelope, or provide a classification/confusion analysis.
  2. [Sect. 3.2 and Fig. 3] The orbital-element distributions in Fig. 3, which underlie the distinguishability claim, include only particles that never leave the Earth's Hill sphere, but the manuscript does not state how many of the 8,077 impactors are in that subset or how the excluded particles' orbital elements differ. Since Sect. 5 notes that about 30% of impactors take longer than one year and contribute greatly to the dust torus, the excluded subset is not negligible. The authors should report the number and fraction of excluded particles and either include them or demonstrate that their orbital distribution does not alter the quoted peaks and ranges.
  3. [Sect. 5 and Figs. 6-12] The conclusion concedes that the uncertainties in the solid-angle fraction and number density are 'at least one order of magnitude', yet no figure or quantitative value in Sects. 3 and 4 carries an error bar, a range, or a propagated uncertainty. A reader cannot judge which of the reported fractions, arrival times, or sky-projection features are robust to the stated uncertainty in mass production rate and initial parameter distributions. The authors should propagate this uncertainty into the figures or perform a sensitivity analysis and state which conclusions survive.
  4. [Sect. 2.1 and Table 1] The principal quantitative results rest on assumptions that are not tested: vertical ejection, power-law size and velocity exponents of 3.7 and 2.2, bulk density 3500 kg/m3, and mass production rate 0.2 kg/s. No sensitivity analysis is provided for any of these inputs. Since the paper itself acknowledges order-of-magnitude uncertainties, the authors should show, at minimum, how the '70% within one year' and the 'large proportion distinguishable from interplanetary dust' statements change under plausible variations of the size exponent and ejection angular distribution.
minor comments (4)
  1. [Title and abstract] There are typographical errors in the displayed title/abstract: 'reachin g' and 'observati on' should be 'reaching' and 'observation'; also 'Aollo 15 and 17' in the introduction should be 'Apollo 15 and 17'.
  2. [Fig. 2] The color scale of Fig. 2 has an ambiguous annotation '10^-3' on the right; the figure should explicitly label the colorbar with the quantity and units (fraction of impactors per bin).
  3. [Fig. 3] The caption should state the number of particles included in the subset of particles that never leave the Earth's Hill sphere, and give the bin widths used for the histograms of a, e, and i.
  4. [Fig. 4] The caption notes that only particles with travel times less than one year are shown, but it does not state that these constitute about 70% of all impactors; adding this information would help readers interpret the omitted 30% tail.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity; the paper is a readout of a forward simulation anchored to external measurements, with no fitted parameter renamed as a prediction.

full rationale

The paper analyzes outputs of a prior forward simulation (Yang et al. 2022) that is anchored to LDEX measurements and literature values (Table 1; Sects. 2.1, 3.1). No parameter is fitted to the Earth-impact statistics, travel times, or sky projections that are later reported, and no reported output is fed back into the model. The authors' self-citation for the simulation is load-bearing, but the cited simulation is an independent numerical experiment with stated assumptions that do not contain the target results. The main caveats—Fig. 3 excludes particles that leave Earth's Hill sphere, and Sect. 5 reports ~30% long-term impactors not shown there—are scope and completeness limitations, not circular reductions. The abstract and Sect. 5 claim that a large proportion of impactors can be distinguished from interplanetary dust, but no interplanetary-dust population or separation metric is introduced; this is an unsupported claim, not a circular derivation. The concluding 'at least one order of magnitude' uncertainty caveat also bears on the strength of the conclusions but does not make the derivation circular. No equation or fitted parameter reduces to an input by construction.

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

The central results depend entirely on the initial condition model inherited from Yang et al. 2022: power-law size and velocity distributions, vertical ejection, fixed density, and a production rate from the literature. These are adopted, not fitted here, but they are free parameters of the simulation. The analysis introduces no new physical entities.

free parameters (6)
  • Size distribution power-law exponent = 3.7
    Controls the fraction of impactors per size bin (Fig. 1); adopted from Yang et al. 2022, not fitted or varied in this paper.
  • Velocity distribution power-law exponent = 2.2
    Shapes the fraction of impactors vs. ejection velocity (Fig. 2); adopted from prior work.
  • Mass production rate of lunar ejecta = 0.2 kg/s
    Converted to the Earth impact rate 2.3e-4 kg/s; taken from Pokorny et al. 2019.
  • Bulk density = 3500 kg/m^3
    Affects solar radiation pressure and mass; taken from Solomon 1974.
  • Ejection velocity grid = 0.95 to 2.00 vesc (10 discrete values)
    The claim that only [vesc, 1.27 vesc] particles impact Earth is based on this discrete sampling; intermediate velocities are not simulated.
  • Particle radius grid = 0.2 to 100 microns (9 discrete values)
    Fraction per size bin depends on this binning; ranges from Table 1.
assumptions (4)
  • domain assumption The gravitational and non-gravitational force model (Earth and Moon point masses, solar gravity, solar radiation pressure, Poynting-Robertson drag) is sufficient for the simulated orbital evolution.
    Invoked in Sect. 2.1; omits Earth oblateness and lunar mascons, and no validation against an independent ephemeris is shown.
  • domain assumption Dust particles are ejected vertically from the lunar surface.
    Stated in Sect. 2.1; real ejecta have a distribution of ejection angles, so launch-angle trends in Fig. 2 are not directly sampled.
  • ad hoc to paper The finite grid of starting regions, times, sizes, and velocities is representative of the real continuous distribution.
    Used to compute fractions; no convergence test with denser grids is reported.
  • ad hoc to paper Impactors that never leave the Earth's Hill sphere are representative of all Earth impactors for orbital element statistics.
    Fig. 3 and Sect. 3.2 are restricted to this subset, but the excluded population is not characterized.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Investigation of lunar ejecta dynamics: Particles reaching the near-Earth space and their effect on Earth-based observation." pith.science (2026). https://pith.science/paper/MZCG4ZZF

@misc{pith2026250714968,
  author       = {Pith},
  title        = {Pith review of: Investigation of lunar ejecta dynamics: Particles reaching the near-Earth space and their effect on Earth-based observation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/MZCG4ZZF}},
  note         = {Machine review of arXiv:2507.14968}
}
abstract

Aims. Particles ejected from the lunar surface via hypervelocity impacts form a torus between the Earth and the Moon. According to our previous study (Yang et al., A\&A, 659, A120), among them about $2.3\times10^{-4}\,\mathrm{kg/s}$ particles impact the Earth after long-term orbital evolution. We mainly focus on these Earth impactors, analyze their orbital element distribution, and estimate their influence on Earth-based observations. Methods. In previous work we simulated the long-term orbital evolution of particles ejected from the lunar surface, and obtained their steady-state spatial distribution in the Earth-Moon system. In this work, we analyze the simulation results about the Earth impactors, including the fraction of impactors with different initial parameters among all impactors, the orbital element distribution, and the projection of particles onto several Earth-based observatories. Results. Particles ejected from the lunar surface are more likely to impact the Earth within a certain range of initial parameters. Most of these lunar-ejected impactors ($\sim70\%$) reach the Earth within one year, while most of the small ones ($87.2\%$ of $0.2\,\mathrm{\mu m}$ particles and $64.6\%$ of $0.5\,\mathrm{\mu m}$ particles) reach the Earth within one week. A large proportion of lunar-ejected Earth impactors can be distinguished from interplanetary dust particles according to the differences in their orbital distributions. Besides, lunar-ejected particles may exhibit distinct configurations and orientations from the perspectives of different Earth-based observatories.

Figures

Figures reproduced from arXiv: 2507.14968 by the authors.

Figure 1
Figure 1. Fraction of impactors vs. grain size [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Fraction of impactors with different launch angles and ejection velocities. 3.2. Orbital distribution and evolution The orbital elements of lunar-ejected particles at the moment of reaching the Earth are shown in [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Distribution of a, e and i for lunar-ejected particles at the moment of reaching the Earth. The symbols of a, e and i denote the semi-major axis, the eccentricity and the inclination of parti￾cles, which are measured in the ECI frame. Panel a: Distribution of a. Panel b: Distribution of e. Panel c: Distribution of i. continuously, while the apogee remains nearly constant around one Earth radius after a rapid descent… view at source ↗
Figures from the paper (7 more)
Figure 4
Figure 4. Figure 4: Distribution of travel time of Earth impactors. [PITH_FULL_IMAGE:figures/full_fig_p004_4.png]
Figure 5
Figure 5. Figure 5: Evolution of e, i, q and Q for a 0.2 µm impactor from ini￾tial retrograde orbit, where e, i, q and Q denote the eccentricity, the inclination, the perigee and the apogee of the particle, which are measured in the geocentric lunar orbital frame. The dashed lines denote …
Figure 6
Figure 6. Figure 6: Fraction of solid angle subtended by particles per an [PITH_FULL_IMAGE:figures/full_fig_p005_6.png]
Figure 9
Figure 9. Figure 9: Fraction of solid angle subtended by particles per an [PITH_FULL_IMAGE:figures/full_fig_p005_9.png]
Figure 10
Figure 10. Figure 10: Fraction of solid angle subtended by particles per a [PITH_FULL_IMAGE:figures/full_fig_p006_10.png]
Figure 11
Figure 11. Figure 11: Geometric Optical depth and number density of lunar [PITH_FULL_IMAGE:figures/full_fig_p006_11.png]
Figure 12
Figure 12. Figure 12: Fraction of solid angle subtended by particles in th [PITH_FULL_IMAGE:figures/full_fig_p007_12.png]

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

23 extracted references · 21 canonical work pages

  1. [1]

    , " * write output.state after.block = add.period write newline

    ENTRY address archiveprefix author booktitle chapter edition editor howpublished institution eprint journal key month note number organization pages publisher school series title type volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 ...

  2. [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 format....

  3. [3]

    1984, Advances in Space Research, 4, 27

    Alexander, M., Anz, P., Hyde, T., et al. 1984, Advances in Space Research, 4, 27

    work page 1984

  4. [4]

    2005, Advances in Space Research, 35, 1282

    Dikarev, V., Gr \"u n, E., Baggaley, J., et al. 2005, Advances in Space Research, 35, 1282

    work page 2005

  5. [5]

    D., Glenar, D

    Feldman, P. D., Glenar, D. A., Stubbs, T. J., et al. 2014, Icarus, 233, 106

    work page 2014

  6. [6]

    2022, Investigations of the dynamical evolution of lunar ejecta

    Fladeland, L. 2022, Investigations of the dynamical evolution of lunar ejecta

    work page 2022

  7. [7]

    2012, Icarus, 221, 1183

    Fritz, J. 2012, Icarus, 221, 1183

    work page 2012

  8. [8]

    1983, in LUNAR AND PLANETARY SCIENCE XIV, P

    Gault, D. 1983, in LUNAR AND PLANETARY SCIENCE XIV, P. 243-244. Abstract., Vol. 14, 243--244

    work page 1983

Show all 23 references
  1. [9]

    J., Burns, J

    Gladman, B. J., Burns, J. A., Duncan, M. J., & Levison, H. F. 1995, Icarus, 118, 302

    1995

  2. [10]

    A., Fechtig, H., & Giese, R

    Gr \"u n, E., Zook, H. A., Fechtig, H., & Giese, R. 1985, Icarus, 62, 244

    1985

  3. [11]

    Hartzell, C. M. & Scheeres, D. J. 2013, Journal of Geophysical Research: Planets, 118, 116

    2013

  4. [12]

    2015, Nature, 522, 324

    Hor \'a nyi, M., Szalay, J., Kempf, S., et al. 2015, Nature, 522, 324

    2015

  5. [13]

    1996, Advances in Space Research, 17, 177

    Iglseder, H., Uesugi, K., & Svedhem, H. 1996, Advances in Space Research, 17, 177

    1996

  6. [14]

    & Brown, P

    Jones, J. & Brown, P. 1993, Monthly Notices of the Royal Astronomical Society, 265, 524

    1993

  7. [15]

    Meeus, J. H. 1991, Astronomical algorithms (Willmann-Bell, Incorporated)

    1991

  8. [16]

    O'Brien, B. J. 2011, planetary and Space Science, 59, 1708

    2011

  9. [17]

    2019, Journal of Geophysical Research: Planets, 124, 752

    Pokorn \`y , P., Janches, D., Sarantos, M., et al. 2019, Journal of Geophysical Research: Planets, 124, 752

    2019

  10. [18]

    & Criswell, D

    Rennilson, J. & Criswell, D. R. 1974, The Moon, 10, 121

    1974

  11. [19]

    Solomon, S. C. 1974, The Moon, 9, 147

    1974

  12. [20]

    Szalay, J. R. & Horányi, M. 2016, Geophysical Research Letters, 43, 4893

    2016

  13. [21]

    R., Pokorn \`y , P., Sternovsky, Z., et al

    Szalay, J. R., Pokorn \`y , P., Sternovsky, Z., et al. 2019, Journal of Geophysical Research: Planets, 124, 143

    2019

  14. [22]

    & Mukai, T

    Yamamoto, S. & Mukai, T. 1996, in International Astronomical Union Colloquium, Vol. 150, Cambridge University Press, 47--50

    1996

  15. [23]

    2022, Astronomy & Astrophysics, 659, A120

    Yang, K., Schmidt, J., Feng, W., & Liu, X. 2022, Astronomy & Astrophysics, 659, A120

    2022

Pith tools

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