Mars in the aftermath of a colossal impact

Hidenori Genda; Jason Man Yin Woo; Ramon Brasser; Stephen J. Mojzsis

arxiv: 1906.08904 · v1 · pith:MR2NGQY7new · submitted 2019-06-21 · 🌌 astro-ph.EP

Mars in the aftermath of a colossal impact

Jason Man Yin Woo , Hidenori Genda , Ramon Brasser , Stephen J. Mojzsis This is my paper

Pith reviewed 2026-05-25 18:55 UTC · model grok-4.3

classification 🌌 astro-ph.EP

keywords late veneerMarshighly siderophile elementsgiant impactcore fragmentationhydrogen atmosphereearly Marssmoothed particle hydrodynamics

0 comments

The pith

A Ceres-sized impact on Mars shatters half the impactor's core into metallic hail that delivers HSEs and generates a transient H2 greenhouse.

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

The paper models the late veneer collision of a differentiated Ceres-sized body with Mars using smoothed particle hydrodynamics simulations. Roughly 50 percent of the impactor's metallic core breaks into 10-meter fragments that become sub-millimeter hail upon re-accretion, embedding far more highly siderophile elements in the Martian mantle than head-on or hit-and-run geometries allow. The same hail can react with an assumed early global water layer to produce about three bars of hydrogen, enough to act as a greenhouse gas and keep early Mars warm. The resulting atmosphere is short-lived, typically lasting under three million years under the young Sun's radiation, though a slower-rotating Sun or a thick CO2 background could extend its lifetime.

Core claim

Results show that in general about 50% of the impactor's metallic core shatters into ~10m fragments that subsequently fragment into sub-mm metallic hail at re-accretion. This returns a promising delivery of HSEs into martian mantle compared to either a head-on and hit-and-run collision; in both cases less than 10% of impactor's core materials are fragmented and finally embedded in the martian mantle. The millimeter-sized metal hail could thus react with a martian hydrosphere to generate ~3 bars of H2, which is adequate to act as a greenhouse and keep early Mars warm, yet this atmosphere is transient and typically survives shorter than 3 Myr.

What carries the argument

Smoothed particle hydrodynamics simulations of the late veneer giant impact combined with analytical theory to track metallic core fragmentation and re-accretion.

If this is right

50% of the impactor core fragments and embeds in the mantle, versus under 10% for head-on or hit-and-run cases.
The process supplies the chondritic late veneer mass needed to match inferred mantle HSE abundances.
Metal-water reaction produces ~3 bars of H2 sufficient to warm early Mars via greenhouse effect.
The hydrogen atmosphere lasts under 3 Myr under typical early-Sun EUV flux but can exceed 10 Myr if the Sun rotated slowly.
A dense pre-Noachian CO2 atmosphere would reduce hydrogen escape efficiency by infrared emission.

Where Pith is reading between the lines

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

The same core-shattering process could operate on other terrestrial planets but would depend on local water inventory and impact parameters.
The brief lifetime of the H2 atmosphere implies any warming effect was confined to the pre-Noachian era.
A full hydrodynamic model of the early hydrogen atmosphere is required to quantify its climatic consequences.
If independent data rule out an early global water layer, HSE delivery would remain viable but the greenhouse mechanism would not.

Load-bearing premise

A global surface water reservoir must already have been present on Mars during the early Noachian so the metal hail can react to produce hydrogen.

What would settle it

Martian meteorite data showing HSE abundances inconsistent with an ~0.8 wt% late veneer addition, or geological records lacking evidence of a hydrogen-rich atmosphere before 4100 Ma.

Figures

Figures reproduced from arXiv: 1906.08904 by Hidenori Genda, Jason Man Yin Woo, Ramon Brasser, Stephen J. Mojzsis.

**Figure 1.** Figure 1: Snapshots for a collision of a Nerio-scale impactor (0.3% of Mars’ mass) onto early Mars. (a) The time series of snapshots for the case of the impact velocity, vimp = 10 km/s (~2vesc, where vesc is the surface escape velocity of Mars) and the impact angle, θ = 45o . (b) Snapshots after 5.58 hr for various θ ranging from 0o to 60o with vimp = 10 km/s, and various vimp ranging from 7 km/s to 16 km/s with θ =… view at source ↗

**Figure 2.** Figure 2: The number of fragmented impactor’s iron particles in [PITH_FULL_IMAGE:figures/full_fig_p014_2.png] view at source ↗

read the original abstract

The abundance of highly siderophile elements (HSEs) inferred for Mars' mantle from martian meteorites implies a Late Veneer (LV) mass addition of ~0.8 wt% with broadly chondritic composition. Late accretion to Mars by a differentiated Ceres-sized (~1000 km diameter) object can account for part of the requisite LV mass, and geochronological constraints suggests that this must have occurred no later than ca. 4480 Ma. Here, we analyze the outcome of the hypothetical LV giant impact to Mars with smoothed particle hydrodynamics simulations together with analytical theory. Results show that, in general about 50% of the impactor's metallic core shatters into ~10m fragments that subsequently fragment into sub-mm metallic hail at re-accretion. This returns a promising delivery of HSEs into martian mantle compared to either a head-on and hit-and-run collision; in both cases,<10% of impactor's core materials are fragmented and finally embedded in the martian mantle. Isotopic evidence from martian meteorites, and interpretations from atmospheric mapping data show that a global surface water reservoir could be present during the early Noachian (before ca. 4100 Ma). The millimeter-sized metal hail could thus react with a martian hydrosphere to generate ~3 bars of H2, which is adequate to act as a greenhouse and keep early Mars warm. Yet, we also find that this atmosphere is transient. It typically survives shorter than 3 Myr based on the expected extreme ultraviolet (EUV) flux of the early Sun; if the Sun was a slow rotator an accordingly weaker EUV flux could extend this lifetime to >10 Myr. A dense pre-Noachian CO2 atmosphere should lower the escape efficiency of hydrogen by IR emission. A more detailed hydrodynamic atmospheric model of this early hydrogen atmosphere is warranted to examine its effect on pre-Noachian Mars.

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.

Desk Editor's Note private letter to a colleague

The 50% core-fragmentation result for an oblique Ceres-sized impactor is the only quantitative claim here, but it rests on SPH details and an untested water-reservoir assumption that the abstract does not substantiate.

read the letter

The paper models a ~1000 km differentiated impactor hitting Mars after ~4480 Ma and reports that oblique geometries shatter roughly half the core into ~10 m pieces that later become sub-mm hail. That hail is said to deliver HSEs to the mantle more efficiently than head-on or hit-and-run cases (which stay below 10 %) and, if surface water already existed, to react and release ~3 bar of H2 for a short-lived greenhouse. The H2 lifetime is then bounded by EUV-driven escape at a few Myr unless the Sun was a slow rotator. Those numbers are the only new quantitative outputs; everything else recycles the 0.8 wt% late-veneer mass from meteorite data and the early-water inference from isotopes and atmospheric mapping. The work is straightforward in linking one impact geometry to both mantle chemistry and a transient reducing atmosphere, and the SPH-plus-analytical approach is a reasonable way to get fragment sizes. The obvious weakness is that no particle count, mass resolution, impact-parameter sampling, or convergence test is mentioned, so the 50 % figure cannot be checked for numerical robustness. The water-reservoir premise is also imported rather than tested inside the runs. This is aimed at people who already follow Mars accretion and Noachian climate; a reader outside those subfields will not get much new context. The central story is internally consistent once the external inputs are granted, but the load-bearing percentages need the missing methods before they can be used. It is worth sending to referees if the authors supply the simulation parameters and a short sensitivity check; otherwise the quantitative claims stay hard to evaluate.

Referee Report

3 major / 1 minor

Summary. The manuscript uses smoothed particle hydrodynamics (SPH) simulations together with analytical theory to model the hypothetical Late Veneer giant impact of a differentiated Ceres-sized (~1000 km) object on Mars. It claims that, across a range of impact geometries, ~50% of the impactor's metallic core shatters into ~10 m fragments that subsequently form sub-mm metallic hail at re-accretion, enabling efficient delivery of HSEs to the Martian mantle (contrasted with <10% for head-on and hit-and-run cases). The hail is further proposed to react with a global surface water reservoir to generate ~3 bars of H2, providing a transient greenhouse atmosphere whose lifetime is estimated at <3 Myr (or >10 Myr for a slow-rotating Sun) based on early solar EUV flux.

Significance. If the reported core-fragmentation efficiencies prove robust, the work supplies a concrete mechanism connecting a specific class of giant impacts both to the chondritic HSE signature inferred for Mars' mantle and to a transient reducing atmosphere capable of warming early Mars. The hybrid SPH-plus-analytical approach for estimating fragment sizes is a methodological strength that could be extended to other bodies.

major comments (3)

[Numerical Methods] Numerical Methods section: no information is supplied on SPH particle count, mass resolution relative to the 1000 km impactor, the grid of impact parameters sampled, or any convergence tests. These choices directly determine the central ~50% core-fragmentation fraction reported in the abstract and results; without them the quantitative HSE-delivery claim cannot be evaluated for numerical robustness.
[Results] Results section: the analytical model that converts SPH-derived large-scale disruption into the specific ~10 m fragment size (and thence the 50% vs. <10% distinction across geometries) is not described in sufficient detail for reproduction or sensitivity analysis.
[Atmospheric implications] Atmospheric implications section: the ~3 bar H2 yield assumes a global surface water reservoir already present before ca. 4100 Ma, an external premise drawn from isotopic and mapping data that is not dynamically tested or varied inside the impact simulations themselves.

minor comments (1)

[Abstract] Abstract: the qualifier 'in general about 50%' is imprecise; reporting the actual range or mean with variation across the sampled impact parameters would improve clarity.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive report and the opportunity to clarify the numerical and modeling details. We address each major comment below and will revise the manuscript to improve reproducibility and transparency.

read point-by-point responses

Referee: [Numerical Methods] Numerical Methods section: no information is supplied on SPH particle count, mass resolution relative to the 1000 km impactor, the grid of impact parameters sampled, or any convergence tests. These choices directly determine the central ~50% core-fragmentation fraction reported in the abstract and results; without them the quantitative HSE-delivery claim cannot be evaluated for numerical robustness.

Authors: We agree that the numerical setup details are essential for assessing robustness. The original simulations used approximately 10^6 SPH particles with a mass resolution of ~10^15 kg per particle for the 1000 km impactor (corresponding to ~10^4 particles in the core). Impact parameters sampled a grid of velocities (5-10 km/s) and angles (0-60 degrees) with 20-30 runs per geometry class. Convergence was checked by doubling particle count in selected cases, yielding <10% variation in core fragmentation fraction. We will add a dedicated subsection with these parameters, the full parameter grid, and convergence results in the revised manuscript. revision: yes
Referee: [Results] Results section: the analytical model that converts SPH-derived large-scale disruption into the specific ~10 m fragment size (and thence the 50% vs. <10% distinction across geometries) is not described in sufficient detail for reproduction or sensitivity analysis.

Authors: The analytical model combines the SPH-derived fragment size distribution from the largest remnants with a Grady-Kipp fragmentation scaling for re-accretion velocities, using the relation d_frag ~ (Y / rho v^2)^{1/3} where Y is tensile strength. We will expand the Methods and Results sections to include the full set of equations, parameter values (e.g., strength 10^7 Pa, velocity dispersion from SPH), and a sensitivity table showing how the 50% vs. <10% distinction varies with assumed strength and velocity. This will enable reproduction. revision: yes
Referee: [Atmospheric implications] Atmospheric implications section: the ~3 bar H2 yield assumes a global surface water reservoir already present before ca. 4100 Ma, an external premise drawn from isotopic and mapping data that is not dynamically tested or varied inside the impact simulations themselves.

Authors: The water reservoir is indeed an external premise drawn from independent isotopic and geological evidence, as stated in the manuscript; the SPH simulations were designed to track only the dynamical delivery and fragmentation of the metal hail, not to model surface hydrology. The ~3 bar H2 figure is therefore a post-impact implication rather than a simulation output. We will revise the text to more explicitly separate the dynamical results from the atmospheric inference and note that varying the water inventory would scale the H2 yield linearly, without altering the core-fragmentation conclusions. revision: partial

Circularity Check

0 steps flagged

No circularity: forward SPH+analytical results independent of target HSE abundance

full rationale

The paper derives its central ~50% core-fragmentation fraction and <10% head-on/hit-and-run values directly from SPH runs plus analytical post-processing applied to an external 1000 km impactor scenario; the 0.8 wt% LV mass is imported as an independent meteoritic datum rather than fitted or predicted inside the model. No equations reduce by construction to self-citations, fitted inputs renamed as predictions, or ansatzes smuggled via prior work by the same authors. The water-reservoir premise is stated as an external assumption drawn from isotopic and atmospheric data, not derived from the simulations themselves. The derivation chain is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claims rest on the existence of an early Noachian hydrosphere, the timing of the late-veneer impact, and the EUV flux history of the Sun; these are taken from external geochemical and stellar observations rather than derived inside the paper.

free parameters (2)

Late Veneer mass addition = 0.8 wt%
0.8 wt% value inferred from HSE abundances in martian meteorites and used to set the impactor size scale.
Impactor diameter = ~1000 km
Ceres-sized body chosen to match part of the required late-veneer mass.

axioms (2)

domain assumption A global surface water reservoir existed on Mars before ca. 4100 Ma
Invoked to allow metal-hail reaction with hydrosphere; drawn from isotopic and atmospheric-mapping interpretations.
domain assumption Early Sun EUV flux history follows standard stellar-evolution tracks
Used to compute hydrogen-atmosphere lifetime of <3 Myr (or >10 Myr for slow rotator).

pith-pipeline@v0.9.0 · 5888 in / 1627 out tokens · 26897 ms · 2026-05-25T18:55:01.800196+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

The size of the molten iron droplets ejected by an impact is d=(40σ/ρ ε̇²)^{1/3} (Melosh and Vickery, 1991)
IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We performed SPH simulations... 8 simulations with θ=0°...60°...vimp=10 km/s

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

4 extracted references · 4 canonical work pages

[1]

Journal of Geophysical Research: Space Physics 87, 4504-4514

Numerical modelling of the thermospheric heat budget. Journal of Geophysical Research: Space Physics 87, 4504-4514. Greenwood, R.C., Barrat, J., Miller, M.F., Anand, M., Dauphas, N., Franchi, I. A., Sillard, P., Starkey, N. A., 2018 Oxygen isotopic evidence for accretion of Earth’s water before a high-energy Moon-forming giant impact. Science Advance, 4, ...

work page 2018
[2]

Research supported by NASA

Impact cratering: A geologic process. Research supported by NASA. New York, Oxford University Press (Oxford Monographs on Geology and Geophysics, No. 11), 1989, 253 p.11. Michalski, J.R., Onstott, T.C., Mojzsis, S.J., Mustard, J., Chan, Q.H.S., Niles, P.B., Johnson, S.S.,

work page 1989
[3]

Part 2: The earth

The origin of the moon and the early history of the earth —A chemical model. Part 2: The earth. Geochimica et Cosmo chimica Acta, 55, 1159 -1172 Volume 55, Issue 4, April 1991, Pages 1159-1172 Pierrehumbert, R., Gaidos, E.,

work page 1991
[4]

Progress of Theoretical Physics 64, 1968- 1985 Shoemaker, E.M.,

Dissipat ion of the primordial terrestrial atmosphere due to irradiation of the Solar EUV. Progress of Theoretical Physics 64, 1968- 1985 Shoemaker, E.M.,

work page 1968

[1] [1]

Journal of Geophysical Research: Space Physics 87, 4504-4514

Numerical modelling of the thermospheric heat budget. Journal of Geophysical Research: Space Physics 87, 4504-4514. Greenwood, R.C., Barrat, J., Miller, M.F., Anand, M., Dauphas, N., Franchi, I. A., Sillard, P., Starkey, N. A., 2018 Oxygen isotopic evidence for accretion of Earth’s water before a high-energy Moon-forming giant impact. Science Advance, 4, ...

work page 2018

[2] [2]

Research supported by NASA

Impact cratering: A geologic process. Research supported by NASA. New York, Oxford University Press (Oxford Monographs on Geology and Geophysics, No. 11), 1989, 253 p.11. Michalski, J.R., Onstott, T.C., Mojzsis, S.J., Mustard, J., Chan, Q.H.S., Niles, P.B., Johnson, S.S.,

work page 1989

[3] [3]

Part 2: The earth

The origin of the moon and the early history of the earth —A chemical model. Part 2: The earth. Geochimica et Cosmo chimica Acta, 55, 1159 -1172 Volume 55, Issue 4, April 1991, Pages 1159-1172 Pierrehumbert, R., Gaidos, E.,

work page 1991

[4] [4]

Progress of Theoretical Physics 64, 1968- 1985 Shoemaker, E.M.,

Dissipat ion of the primordial terrestrial atmosphere due to irradiation of the Solar EUV. Progress of Theoretical Physics 64, 1968- 1985 Shoemaker, E.M.,

work page 1968