{"id":"c2eda163-fa7e-4081-894c-53c5367b0c2e","arxiv_id":"2412.00519","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"3D MHD simulations indicate that terrestrial atmospheric ions reach the Moon mainly inside Earth's magnetotail, that the geodynamo era dominates the non-solar lunar volatile budget, and that Earth's exobase was never below 190 km during implantation.","lead":"This paper argues that Earth's atmosphere, not the solar wind, supplied most of the non-solar nitrogen and noble gases found in lunar soil, and that this delivery happened during Earth's long magnetized (geodynamo-active) history. If correct, buried lunar soils could preserve a multi-billion-year record of Earth's atmosphere and its exobase height.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The MHD atmosphere uses an inflated ~10^3-km scale height with 625-km cells; the Table 2 Earth-wind flux therefore rests on an unresolved, non-Earth-like density gradient, and a realistic-H test could flip the magnetized-case comparison.","rationale":"In good faith, the central mechanism is plausible: terrestrial ions are observed in Earth's magnetotail, the Moon does pass through the tail, and the MHD framework is appropriate at global scales. The paper presents a self-consistent scenario and gives a clear, falsifiable comparison in Table 2. However, the quantitative link from the MHD flux to lunar soil abundances has an unresolved sensitivity. The ~10^3 km scale height is a modeling convenience that can overestimate the neutral density available for pickup, and a 625 km cell cannot resolve an 8 km thermospheric scale height. Because the Table 2 margins are only about 9-25x, this is not a minor numerical detail; a factor-of-ten change in the EW flux would flip the magnetized-case conclusion for N and He. The reader's conditional verdict already reflects this risk, and the proposed resolution/reference-atmosphere test would settle whether the concern actually lands. The paper's own Limitations section does not acknowledge this gap, which reinforces the need for the test.","tokens_in":874,"tokens_out":780,"duration_ms":221730,"concrete_test":"Run the magnetized Case I simulation again with the neutral atmosphere prescribed from NRLMSISE-00 below the exobase and a collisionless exosphere above it, using AMR resolution Delta x <= H/4 near and above the exobase; if full 3D is too costly, use a 1D or 2D axisymmetric surrogate to compute the pickup-ion flux. Recompute the orbit-averaged EW flux at 60 R_E and compare with Table 2: if the EW flux falls below the measured non-solar flux for N or He, the central claim fails, and if it remains above with margin, the resolution concern is retired.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central comparison in Table 2 is the orbit-averaged Earth-wind (EW) flux from the MHD model against the non-solar flux inferred from lunar samples. The paper's own numbers make this comparison load-bearing: for magnetized Case I the EW/non-solar ratios are roughly 12 for N, 9 for He, and 25 for Ne, so an error of about one order of magnitude in the MHD flux would erase the claimed success. Section 4.1.5 states that the finest cell size in the atmosphere is 625 km, and Section 4.2.3 says the scale height used in the MHD atmosphere is ~10^3 km. A real thermosphere has H nearly 8 km below the exobase; with H=10^3 km the implied neutral density at exobase altitudes is orders of magnitude above NRLMSISE-00 values, and the density gradient that controls solar-wind mass loading is not represented. The Supplementary Limitations (S2) discusses rotation and grid artifacts but does not address this. No resolution or convergence test is reported, so the magnitude of the EW flux in Table 2 is not tied to a realistic atmosphere. If a resolved, Earth-like density profile lowers the EW flux by more than ~10x, the magnetized case no longer explains the non-solar N and He fluxes, and the dynamo-history conclusion loses its quantitative support.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper combines 3-D MHD simulations of the solar wind interacting with magnetized (present-day) and unmagnetized (Eoarchean) Earth atmospheres with photoionization and binary-mixing models to interpret light volatile elements in lunar regolith. It argues that terrestrial atmospheric ions are transported to the Moon mainly when the Moon is in Earth's magnetotail, that the non-solar N, He, Ne, and Ar in Apollo samples are best explained by implantation during the long magnetized phase of the geodynamo rather than during a hypothetical unmagnetized epoch, and that the exobase altitude at the time of implantation was never below about 190 km.","tokens_in":34199,"tokens_out":6375,"duration_ms":66283,"significance":"If the central comparison holds, the paper would provide a quantitative framework for reading the lunar regolith as an archive of terrestrial atmospheric escape over billions of years, with direct implications for geodynamo history and paleoatmospheric evolution. The work is constructive: it uses an explicit multi-species MHD model, a photoionization model, and a binary-mixing analysis; it considers multiple parameter cases; and it provides analytic scaling relations and a limitations section in the Supplementary Materials. The main risk is that the quantitative conclusion rests on the orbit-averaged Earth-wind flux from simulations whose atmospheric density gradient is not resolved, and on a one-sided comparison criterion.","major_comments":[{"comment":"The finest atmospheric cell size is 625 km (§4.1.5) and the scale height used in the MHD atmosphere is ~10^3 km (§4.2.3), while the real thermospheric scale height below the exobase is about 8 km. The Earth-wind flux in Table 2 is the load-bearing quantity, and it must exceed the non-solar flux by factors of ~12 (N), ~9 (He), and ~25 (Ne) for the magnetized case to be judged successful. That flux is generated by the solar wind mass loading of an atmosphere whose density gradient is completely unresolved, and no resolution or convergence study is reported. The Supplementary Limitations (S2) discuss rotation and grid artifacts but not this issue. The statement in §4.2.3 that the ion gyroradius (~560 km) is comparable to the MHD scale height justifies the MHD approximation, but it does not address whether the grid resolves the mass-loading gradient. This concern is load-bearing for the central claim and needs to be addressed with either a resolution study or an explicit demonstration that the Earth-wind flux is insensitive to the atmospheric scale height.","section":"§4.1.5, §4.2.3, Table 2"},{"comment":"The criterion used to declare success is that the model Earth-wind flux 'surpasses' the measured non-solar flux. For N, He, and Ne in the magnetized case the excess is roughly one order of magnitude (ratios of ~12, ~9, and ~25, respectively). Without a retention or implantation efficiency, or an uncertainty estimate on both fluxes, an order-of-magnitude excess is not a quantitative match; it is more accurately an upper-limit statement. The paper should either provide a bounded prediction with propagated uncertainties or justify why a large excess is consistent with the observed lunar concentrations.","section":"§2.3, Table 2"},{"comment":"The exobase altitude is effectively a fitted parameter: Figures 6 and ED4 show mixing curves for different exobase heights, and the text identifies ~250 km as the best match for the magnetized case, while the abstract and discussion convert this into a firm statement that the exobase was never below 190 km. This is an inference from the model, not an independent measurement. Furthermore, the Earth-wind fluxes in Table 2 appear to be based on a single exobase value (401 km for the magnetized case and 221 km for the unmagnetized case, as tabulated in Tables ED4 and ED5), which may be inconsistent with the best-fit exobase used in the mixing diagrams. The paper should state explicitly which exobase height underlies each column of Table 2 and recompute the comparison consistently for the best-fit exobase.","section":"Figures 6/ED4, Table 2, Tables ED4/ED5"},{"comment":"The non-solar flux that is compared with the Earth-wind flux is obtained by multiplying the MHD-derived solar-wind flux by a non-solar fraction inferred from the same binary-mixing model that is used to test the Earth-wind hypothesis. This is not fully circular because the Earth-wind flux magnitude is an independent MHD output, but the two sides of the comparison are not independent of the solar-wind flux normalization. The paper should clearly separate quantities inferred from the lunar isotope data from quantities predicted by the model, and should show how the Table 2 conclusion changes when the solar-wind flux is varied within observational uncertainty.","section":"§2.3, Eqs. (29)–(32)"}],"minor_comments":[{"comment":"There is a typo in the definitions of the number densities: the text says 'n_SW denotes the number density of the SW, while n_SW represents that of EW'; the second symbol should be n_EW.","section":"§2.1"},{"comment":"The caption states that the solid mixing curve corresponds to the current exobase at ~401 km, while the main text (§2.3) says the model best matches the data at ~250 km; these statements should be reconciled and the exobase value used for each curve clearly labeled.","section":"Figure 6 caption"},{"comment":"The species-specific Earth-wind fluxes are tabulated for only one exobase height each (401 km and 221 km), even though the mixing diagrams show strong exobase sensitivity. A table or figure giving these fluxes for the full range of exobase heights discussed in Figures 6 and ED4 would make the comparison in Table 2 transparent and reproducible.","section":"Tables ED4/ED5"},{"comment":"The statement that simulation outputs and Python routines are 'available upon request' is weaker than current reproducibility standards; archiving the reduced data and scripts in a permanent repository would strengthen the paper.","section":"Data and Materials Availability"}],"recommendation":"major_revision","confidential_remarks":"The core idea is interesting and the modeling framework is extensive, but the central quantitative conclusion depends on an unresolved atmospheric density gradient in the MHD simulations and on a one-sided comparison criterion. I would be willing to support publication if the authors can supply a resolution or sensitivity study for the atmospheric scale height, recompute Table 2 with a consistent treatment of the exobase, and provide uncertainties for the compared fluxes."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this paper runs 3D MHD simulations with passive tracers to track Earth wind and solar wind separately, for both a magnetized present-day Earth and an unmagnetized Archean Earth, and concludes the non-solar volatile component in lunar soil was probably implanted during the long geodynamo phase, not during a short unmagnetized epoch as Ozima et al. argued in 2005. That is a genuinely new result, and the magnetotail channel they identify is an important idea.\n\nThe paper has real strengths: the setup is explicit, the ionization model and binary mixing framework are standard and well described, and they admit the H data do not fit their model. Using an established AMR code with published tests also helps.\n\nBut the quantitative comparison in Table 2 has a load-bearing weakness. The MHD atmosphere is initialized with a scale height around 10^3 km, about two orders of magnitude larger than the real thermosphere's ~8 km, and the finest cell size is 625 km, so the density gradient that controls mass loading is not resolved. The paper does not report a resolution study, and the S2 Limitations discussion covers rotation and grid artifacts but is silent on this. Since the inferred EW fluxes in Table 2 are within roughly one order of magnitude of the non-solar fluxes for N, He, and Ne, an overestimate of the MHD total flux by even 10x would undo the claim that the magnetized case explains the data. The stress-test concern is therefore solid.\n\nThe other soft spot is the exobase height. The mixing curves are parameterized by exobase altitude, and the data are used to fit it, so the 'never below 190 km' statement is a fitted model output, not an independent constraint. That said, the model self-consistently explains the isotope ratios for the fitted heights, so this is a fair inference but should be labeled as such.\n\nI'd send this to a good referee, but I'd insist on a resolution study or at least a sensitivity test with realistic scale heights, and on redoing Table 2 with that test. The qualitative claim about the magnetotail is worth preserving.\n\nI wouldn't cite the numbers yet; the mechanism might be cited once it's been checked.","headline":"The paper's qualitative claim that a geodynamo, not an unmagnetized early Earth, delivered terrestrial volatiles to the Moon is fresh and worth discussing, but the quantitative flux comparison rests on an inflated MHD atmosphere scale height that needs a resolution study before the numbers can be trusted.","tokens_in":34716,"tokens_out":4007,"would_cite":false,"duration_ms":42502,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper claims that the non-solar light volatiles in nearside lunar soil were implanted as ionized Earth atmosphere carried by the solar wind through Earth's magnetotail during the long magnetized history of the geodynamo, not during a…","keywords":["lunar regolith","Earth wind","solar wind implantation","geodynamo","nitrogen isotopes","noble gases","magnetotail","MHD simulation"],"falsifier":"Recompute the magnetized Case-I Earth-wind flux at higher resolution (finest cell about 312 km and 156 km by adding AMR levels); if the orbit-averaged flux falls by more than about a factor of two below the measured non-solar flux for nitrogen (roughly 2e4 $m^{-2}$ $s^{-1}$), the central match is an artifact of the unresolved atmospheric scale height. Alternatively, measurement of a non-solar nitrogen component in a well-dated lunar sample whose implantation epoch corresponds to an exobase below 190 km would break the stated exobase constraint.","tokens_in":33683,"feed_emoji":"🌙","tokens_out":8091,"duration_ms":75375,"temperature":0.7,"pith_summary":"The paper tries to settle where the non-solar nitrogen and light noble gases in nearside lunar soil came from: the solar wind alone cannot explain their abundances or isotopic ratios, and a long-standing proposal attributes them to ions escaping Earth's atmosphere during an early epoch when Earth had no magnetic field. Using three-dimensional MHD simulations of the solar wind interacting with a present-day magnetized Earth and an Archean unmagnetized Earth, together with a photoionization model and binary isotope-mixing analysis, the authors argue the opposite: the non-solar component was implanted by 'Earth wind' ions during the long magnetized history of the geodynamo, when the Moon passed through Earth's magnetotail. They further conclude that the exobase altitude at the time of implantation was never below 190 km, so lunar soil can serve as a multi-billion-year archive of Earth's atmosphere. If correct, this re-dates the lunar volatile record and changes what the Apollo samples imply about early Earth's magnetic field and atmosphere.","feed_headline":"A magnetized Earth sent its atmosphere to the Moon's soil","feed_subtitle":"3D simulations show the magnetotail carried terrestrial ions into nearside lunar soil, explaining its non-solar nitrogen.","key_machinery":"The load-bearing object is the 'Earth wind' channel: ionized atmospheric species produced above the exobase are picked up by the solar wind and carried antisunward down the magnetotail to the Moon's nearside surface. The argument runs through three linked components: (i) an adaptive-mesh MHD simulation of the solar wind interacting with an isothermal atmosphere matched to a dipole field (exponential density profile below the beta = 1 surface, power-law tail above) and with a fixed-total-mass unmagnetized atmosphere; (ii) a photoionization model using empirical neutral densities and solar EUV spectra to compute species-specific escape fluxes above the exobase, normalized to the MHD bulk flux; and (iii) binary mixing hyperbolae for 15N/14N–D/H, 3He/20Ne, and 3He/40Ar that map the observed Apollo isotope data onto exobase-dependent curves. The decisive comparison is the table of computed Earth-wind flux versus measured non-solar flux for each species in each Earth state.","core_discovery":"Using a single-fluid multi-species MHD model with passive tracers for solar wind and atmospheric material, the paper computes orbit-averaged fluxes of terrestrial ('Earth wind') ions and solar wind ions arriving at the lunar nearside sub-Earth point, for magnetized and unmagnetized Earth configurations. The Earth wind is only appreciable when the Moon is inside the magnetotail; there the planet flux is comparable to or exceeds the measured non-solar fluxes of He, N, Ne, and Ar for the present-day magnetized case, while the Archean unmagnetized case underproduces these species because the stronger early solar wind dilutes the terrestrial contribution. Comparison of Apollo 17 ilmenite isotope data with binary mixing hyperbolae between solar wind and Earth wind end-members places the data on curves corresponding to exobase heights of roughly 190–300 km, yielding the constraint that the exobase never fell below 190 km. When the magnetic field effect is isolated with identical solar wind conditions, the magnetized Earth delivers about an order of magnitude less terrestrial flux to the Moon than the unmagnetized Earth, so the field is protective for this transport channel, yet the much longer magnetized era still dominates the time-integrated implantation.","pith_inferences":["An immediate numerical test is a resolution study: rerunning the magnetized Case-I simulation with one or two additional AMR levels (finest cell 312 km or 156 km instead of 625 km) would show whether the orbit-averaged Earth-wind flux is stable; a large drop would implicate the unresolved density gradient as the driver of the Table 2 match.","The same magnetotail transport geometry should apply to any magnetized planet hosting a close large moon, so terrestrial-ion implantation is a generic mechanism that could seed volatiles on such moons and might be testable with future sample return.","The 190 km exobase floor is a paleoatmospheric thermometer that could be cross-checked against independent reconstructions of Archean atmosphere density from noble-gas escape models."],"forward_implications":["The nearside lunar regolith preserves a record of Earth's atmospheric composition extending over billions of years, not merely a record of solar wind history.","The observed non-solar volatile inventory is consistent with a geodynamo active from at least roughly 3.5–4.2 Ga and with the absence of a long-lived lunar magnetosphere after about 4.36 Ga.","Nitrogen isotope variations in lunar soils are explained as mixing between solar wind and terrestrial end-members, with the data requiring exobase heights below roughly 300 km at implant time.","Isolating the field effect shows the dynamo reduces Earth-wind delivery by about an order of magnitude, clarifying why a magnetic field can be partially protective while still allowing substantial volatile transfer over long times.","The exobase never smaller than 190 km is a testable constraint on the paleo-Earth atmosphere's thermal structure."],"supporting_citations":[{"why":"Proposed that terrestrial nitrogen and noble gases in lunar soil came from Earth's atmospheric ion flow during an early unmagnetized epoch; this is the hypothesis the paper tests and refines.","marker":"[12]"},{"why":"Provided the Apollo ilmenite nitrogen and hydrogen isotope data and the solar-versus-non-solar mixing framework the paper's mixing diagrams compare against.","marker":"[7]"},{"why":"Supplied the helium, neon, and argon isotope measurements in lunar ilmenites and regolith breccias used to build the He-Ne and He-Ar mixing diagrams.","marker":"[45]"},{"why":"Provided the Eoarchean homosphere composition (equal CO2 and N2) and young-Sun XUV flux scaling used in the unmagnetized Earth model.","marker":"[68]"},{"why":"Empirical neutral density model for the present-day thermosphere that feeds the species-specific photoionization escape fluxes.","marker":"[56]"},{"why":"Solar EUV flux model used to compute photoionization rates above the exobase for both present and Archean atmospheres.","marker":"[57]"},{"why":"Paleointensity evidence fixing the earliest geodynamo at roughly 4.2 Ga, defining the time window in which magnetized-Earth implantation could occur.","marker":"[33]"},{"why":"The adaptive mesh refinement MHD code used for the wind-atmosphere interaction simulations; its solvers carry the bulk-flux calculations.","marker":"[49]"}],"fun_headline_variants":["Moon's soil holds Earth's ancient atmosphere, magnetotail delivered it","Earth's magnetic field guided atmosphere into lunar soil","Magnetotail carried Earth's ions to Moon for billions of years","Lunar regolith records Earth's paleoatmosphere via magnetotail","Apollo samples show Earth's magnetic field sent atmosphere to Moon"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The simulations resolve the atmosphere with a finest cell of 625 km while the real pressure scale height is about 8 km, so the numerical atmosphere's density gradient—which controls how much atmospheric mass the solar wind picks up—may not match reality, and the computed Earth-wind flux could be off; the paper does not demonstrate convergence with resolution.","fun_headline_variants_meta":{"raw":{"variants":["Moon's soil holds Earth's ancient atmosphere, magnetotail delivered it","Earth's magnetic field guided atmosphere into lunar soil","Magnetotail carried Earth's ions to Moon for billions of years","Lunar regolith records Earth's paleoatmosphere via magnetotail","Apollo samples show Earth's magnetic field sent atmosphere to Moon"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000582,"raw_usage":{"total_tokens":2786,"prompt_tokens":1036,"completion_tokens":1750,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":652,"completion_tokens_details":{"reasoning_tokens":1661}},"tokens_in":652,"tokens_out":1750,"duration_ms":12090,"temperature":1.0,"reasoning_tokens":1661,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T05:17:05.332159+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the magnetized Case-I Earth-wind flux at higher resolution (finest cell about 312 km and 156 km by adding AMR levels); if the orbit-averaged flux falls by more than about a factor of two below the measured non-solar flux for nitrogen (roughly 2e4 $m^{-2}$ $s^{-1}$), the central match is an artifact of the unresolved atmospheric scale height. Alternatively, measurement of a non-solar nitrogen component in a well-dated lunar sample whose implantation epoch corresponds to an exobase below 190 km would break the stated exobase constraint.","supporting_citations":[{"cited_title":"Ozima, K","cited_arxiv_id":null,"evidence_quote":"Proposed that terrestrial nitrogen and noble gases in lunar soil came from Earth's atmospheric ion flow during an early unmagnetized epoch; this is the hypothesis the paper tests and refines."},{"cited_title":"Hashizume, M","cited_arxiv_id":null,"evidence_quote":"Provided the Apollo ilmenite nitrogen and hydrogen isotope data and the solar-versus-non-solar mixing framework the paper's mixing diagrams compare against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplied the helium, neon, and argon isotope measurements in lunar ilmenites and regolith breccias used to build the He-Ne and He-Ar mixing diagrams."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provided the Eoarchean homosphere composition (equal CO2 and N2) and young-Sun XUV flux scaling used in the unmagnetized Earth model."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Paleointensity evidence fixing the earliest geodynamo at roughly 4.2 Ga, defining the time window in which magnetized-Earth implantation could occur."}],"review_version":1}