pith. sign in

arxiv: 1907.07575 · v1 · pith:6VNIFGPYnew · submitted 2019-07-17 · 🌌 astro-ph.EP

Assessment of the probability of microbial contamination for sample return from Martian moons I: Departure of microbes from Martian surface

Kazuhisa Fujita , Kosuke Kurosawa , Hidenori Genda , Ryuki Hyodo , Shingo Matsuyama , Akihiko Yamagishi , Takashi Mikouchi , Takafumi Niihara This is my paper

Pith reviewed 2026-05-24 20:04 UTC · model grok-4.3

classification 🌌 astro-ph.EP
keywords microbial contaminationMars ejectaPhobosDeimossample returnatmospheric entrymeteoroid impactsmicrobe survival
0
0 comments X

The pith

Mars ejecta smaller than 0.03 m rarely deliver viable microbes to Phobos or Deimos because they decelerate or heat-sterilize in the atmosphere.

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

This paper calculates the microbial content of material blasted off Mars by impacts and whether that material can reach the orbits of Phobos and Deimos. It combines an Earth-analog estimate of surface microbe density, numerical modeling of survival through the impact itself, and trajectory simulations that include drag and heating while the ejecta cross the Martian atmosphere. The central finding is that particles below three centimeters across almost never arrive intact. A sympathetic reader cares because the numbers supply the initial inventory for any later statistical assessment of contamination risk in sample-return missions from the moons.

Core claim

Mars ejecta smaller than 0.03 m in diameter hardly reach the Phobos orbit due to aerodynamic deceleration, or mostly sterilized due to significant aerodynamic heating even though they can reach the Phobos orbit and beyond. A baseline dataset of microbial density in Mars ejecta departing for Martian moons has been presented.

What carries the argument

Numerical analysis of hypervelocity impacts (with and without internal friction and plastic deformation) to estimate microbial survival, combined with computational fluid dynamics of ejecta trajectories through the Martian atmosphere that account for aerodynamic deceleration and heating.

If this is right

  • Only ejecta larger than 0.03 m carry a realistic chance of delivering living microbes to the moons.
  • The baseline microbial-density dataset can be inserted directly into later probability calculations for contamination.
  • Survival rates during launch depend on whether internal friction and plastic deformation are modeled in the impact.
  • Atmospheric drag and heating act as the dominant filter that removes or kills microbes on small particles.

Where Pith is reading between the lines

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

  • Contamination models should concentrate effort on the rare large-ejecta events rather than the far more numerous small ones.
  • The Earth-analog density assumption could be checked against any future Martian surface or returned-sample data.
  • Applying the same methods to larger impacts or a denser early-Martian atmosphere might alter the picture for ancient material transfer.

Load-bearing premise

Microbial density on the Martian surface is estimated by direct analogy to cold arid regions on Earth.

What would settle it

Detection of a substantial fraction of centimeter-scale Mars ejecta that reach Phobos orbit while still carrying viable microbes would falsify the sterilization conclusion.

read the original abstract

Potential microbial contamination of Martian moons, Phobos and Deimos, which can be brought about by transportation of Mars ejecta produced by meteoroid impacts on the Martian surface, has been comprehensively assessed in a statistical approach, based on the most probable history of recent major gigantic meteoroid collisions on the Martian surface. This article is the first part of our study to assess potential microbial density in Mars ejecta departing from the Martian atmosphere, as a source of the second part where statistical analysis of microbial contamination probability is conducted. Potential microbial density on the Martian surface as the source of microorganisms was estimated by analogy to the terrestrial areas having the similar arid and cold environments, from which a probabilistic function was deduced as the asymptotic limit. Microbial survival rate during hypervelocity meteoroid collisions was estimated by numerical analysis of impact phenomena with and without taking internal friction and plastic deformation of the colliding meteoroid and the target ground into consideration. Trajectory calculations of departing ejecta through the Martian atmosphere were conducted with taking account of aerodynamic deceleration and heating by the aid of computational fluid dynamic analysis. It is found that Mars ejecta smaller than 0.03 m in diameter hardly reach the Phobos orbit due to aerodynamic deceleration, or mostly sterilized due to significant aerodynamic heating even though they can reach the Phobos orbit and beyond. Finally, the baseline dataset of microbial density in Mars ejecta departing for Martian moons has been presented for the second part of our study.

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

3 major / 2 minor

Summary. The manuscript is the first part of a two-part study on microbial contamination risk for Phobos and Deimos via Mars ejecta. It estimates surface microbial density via analogy to terrestrial arid/cold sites and derives an asymptotic probabilistic function; models microbial survival during hypervelocity impacts via numerical analysis (with and without internal friction/plastic deformation); and computes atmospheric trajectories of ejecta using CFD to include aerodynamic deceleration and heating. The central result is that ejecta <0.03 m diameter either fail to reach Phobos orbit or are sterilized, and a baseline microbial density dataset for departing ejecta is supplied for the contamination-probability analysis in part II.

Significance. If the numerical results hold, the work supplies a quantitative source-term dataset that is directly usable for planetary-protection calculations in future Phobos/Deimos sample-return missions. The combination of impact hydrodynamics and CFD-based trajectory modeling constitutes a methodical treatment of the departure phase. The terrestrial-analogy source term, however, remains the dominant uncertainty and limits the robustness of any downstream probability estimates.

major comments (3)
  1. [microbial density estimation section] The probabilistic microbial-density function (abstract and the section describing estimation of potential microbial density on the Martian surface) is obtained solely by mapping terrestrial arid/cold analog sites onto Mars. Because this function is the linear source term for the entire baseline dataset, the absence of any Martian-specific bounding constraint (e.g., perchlorate abundance, cumulative radiation dose, or desiccation history) makes the 0.03 m threshold result and the final dataset inherit the full, unquantified transferability error.
  2. [trajectory calculations section] The trajectory calculations that establish the 0.03 m diameter cutoff (abstract and the section on trajectory calculations of departing ejecta) are presented without reported sensitivity tests to atmospheric-density profile variations, ejecta shape factors, or lift/drag coefficient uncertainties. These parameters directly control whether the aerodynamic-deceleration conclusion is robust.
  3. [impact phenomena numerical analysis section] The numerical analysis of impact survival rates (abstract and the section on microbial survival rate during hypervelocity meteoroid collisions) does not state the sampled ranges or distributions of impact velocity, angle, and target/impactor material properties used to generate the survival fractions that enter the baseline dataset. Without this information the sterilization claim for the 0.03 m population cannot be reproduced or stress-tested.
minor comments (2)
  1. [abstract] The abstract states the 0.03 m threshold but supplies no accompanying numerical values, error estimates, or figure references; adding one sentence that points to the relevant table or figure would improve traceability.
  2. [microbial density estimation section] Notation for the probabilistic density function should be introduced once with a clear symbol and units when it is first defined, rather than only in the final dataset description.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We appreciate the referee's detailed review and constructive suggestions. We address each of the major comments below. Revisions will be made to enhance the transparency of our methods and to discuss uncertainties more explicitly.

read point-by-point responses

  1. Referee: [microbial density estimation section] The probabilistic microbial-density function (abstract and the section describing estimation of potential microbial density on the Martian surface) is obtained solely by mapping terrestrial arid/cold analog sites onto Mars. Because this function is the linear source term for the entire baseline dataset, the absence of any Martian-specific bounding constraint (e.g., perchlorate abundance, cumulative radiation dose, or desiccation history) makes the 0.03 m threshold result and the final dataset inherit the full, unquantified transferability error.

    Authors: We acknowledge the limitations of relying on terrestrial analogs for estimating microbial densities on Mars. Direct measurements from Mars are not available, and terrestrial arid and cold environments provide the closest proxies as used in prior planetary protection assessments. The 0.03 m threshold is derived from the trajectory calculations and is independent of the microbial density estimate. The density function serves as the source term for part II. In the revised version, we will include an expanded discussion on potential Martian-specific factors such as perchlorates, radiation, and desiccation, noting that these would likely result in lower densities, rendering our baseline conservative. We will also quantify the transferability uncertainty where possible based on available literature. revision: yes

  2. Referee: [trajectory calculations section] The trajectory calculations that establish the 0.03 m diameter cutoff (abstract and the section on trajectory calculations of departing ejecta) are presented without reported sensitivity tests to atmospheric-density profile variations, ejecta shape factors, or lift/drag coefficient uncertainties. These parameters directly control whether the aerodynamic-deceleration conclusion is robust.

    Authors: The original manuscript used the standard Mars atmospheric model and assumed spherical particles with drag and lift coefficients drawn from established literature values for high-speed entry. Explicit sensitivity tests were not included in the submitted version. To address this, we will add a new subsection in the revised manuscript presenting sensitivity analyses for atmospheric density variations, non-spherical shape factors, and ranges of drag/lift coefficients. This will confirm the robustness of the 0.03 m cutoff conclusion. revision: yes

  3. Referee: [impact phenomena numerical analysis section] The numerical analysis of impact survival rates (abstract and the section on microbial survival rate during hypervelocity meteoroid collisions) does not state the sampled ranges or distributions of impact velocity, angle, and target/impactor material properties used to generate the survival fractions that enter the baseline dataset. Without this information the sterilization claim for the 0.03 m population cannot be reproduced or stress-tested.

    Authors: The impact survival calculations were based on a Monte Carlo sampling of parameters drawn from literature values for Martian surface conditions and meteoroid properties. Specific ranges included impact velocities of 2-12 km/s, impact angles of 15-75 degrees from horizontal, and material strengths consistent with basaltic regolith and chondritic impactors. We will revise the manuscript to explicitly document these parameter ranges, distributions, and the number of simulations performed, allowing for reproducibility and further stress-testing. revision: yes

Circularity Check

0 steps flagged

No circularity; derivation uses external terrestrial analogy and independent numerical models

full rationale

The microbial source density is obtained from terrestrial arid/cold site analogy (an external input, not fitted to Martian ejecta outcomes). Survival rates and trajectories are computed via separate numerical impact and CFD analyses. No equation or step reduces the final baseline dataset to the source term by construction, and no self-citation chains or ansatzes are invoked to force results. The paper is self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The assessment rests on terrestrial environmental analogies for microbial density and on numerical models for impact survival and atmospheric trajectories whose internal parameters and validation are not shown in the abstract.

free parameters (1)
  • microbial density probabilistic function
    Deduced as the asymptotic limit from analogy to terrestrial arid and cold environments
axioms (2)
  • domain assumption Microbial density on Mars can be estimated by analogy to terrestrial arid and cold areas
    Used as source for microorganisms in the statistical approach
  • domain assumption Numerical analysis of impact phenomena with and without internal friction and plastic deformation accurately estimates microbial survival rate
    Basis for survival rate during hypervelocity collisions

pith-pipeline@v0.9.0 · 5830 in / 1441 out tokens · 26176 ms · 2026-05-24T20:04:11.070350+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

27 extracted references · 27 canonical work pages

  1. [1]

    Introduction The planet Mars and its moons, Phobos and Deimos, have attracted continual scientific interests from viewpoints of the process of chemical evolution and the origins of life. Motivated by such scientific interests, Japan Aerospace Exploration Agency (JAXA) is currently entertaining a sample return mission from Martian moons for launch in 2024 ...

    work page 2024

  2. [2]

    As a whole, all the processes are quantitatively examined to obtain an estimate of the potential microbial density on the surface of the Martian moons, from which microbial contamination probability in samples collected on the surface of the Martian moons is assessed. Fig. 1 Potential scenario of microbial contamination of Martian moons. 5 2.2. Comparison...

    work page 2011

  3. [3]

    Potential Microbial Density on Martian Surface Estimation of microbial density on the Martian surface is of primary importance since the final result of microbial contamination probability linearly changes with the estimated microbial density. Since it is reasonable to get insight from microbial density in the terrestrial areas having the similar arid and...

    work page 2003

  4. [4]

    Maier et al

    However, it should be noted that microbial densities measured in the Yungay area, which is the most arid zone in the Atacama Desert, are not greater than 107 CFU/kg. Maier et al. (2004) reported microbial densities ranging from 106 to 108 CFU/kg in the Yungay area as well. Their results show good agreement with those of Navarro-González et al. (2003). Sin...

    work page 2004

  5. [5]

    Sterilization during Mars Ejecta Formation Mars ejecta transported to the Martian moons are produced by hypervelocity impacts of meteoroids against the Martian surface. In such the impact events that can transport the ejected rocklets to the Martian moons, the incident meteoroids must have considerably high impact velocities, since the ejected rocklets sh...

    work page 2018

  6. [6]

    Strength parameters are taken from the work for basalt (Ivanov et al., 2010)

    with taking internal friction and plastic deformation into consideration. Strength parameters are taken from the work for basalt (Ivanov et al., 2010). The numerical suite constructed by Kurosawa et al. (2018) was employed in this simulation, but we newly introduced a constitutive model for basaltic materials. The general setup parameters for the simulati...

    work page 2010

  7. [7]

    when internal friction and plastic deformation are not present (Fig. 5a). In contrast to this, the ejecta having velocities greater than 3.8 km/s are Table 2 General setup parameters for the 2-D iSALE calculation. The detailed descriptions for values listed here can be found in Collins et al. (2016). Computational geometry Cylindrical coordinates Number o...

    work page 2016

  8. [8]

    (2005) a) b) Fig

    2 Modified by Collins and Melosh (2014) after Pierazzo et al. (2005) a) b) Fig. 5 Temperature of ejected particles at ejection in a hypervelocity impact in relation to ejection velocity; a) without internal friction and plastic deformation, and b) with internal friction and plastic deformation. The two grey dashed lines indicate the temperatures required ...

    work page 2014

  9. [9]

    At the same time, however, Nyquist et al

    This fact supports the above argument that the microbial survival rate is sufficiently small in Mars ejecta formation. At the same time, however, Nyquist et al. (2001) reported that a considerable portion of Martian meteorites do not have a sign of shock-heating. This might come from the fact that the geochemical barometers, such as shock melted glasses (...

    work page 2001

  10. [10]

    Nevertheless, the results by Nyquist et al

    and atomic diffusion (e.g., Takenouchi et al., 2017), for the estimation of the temperature rise have relatively low sensitivities at a moderate temperature (several hundred K). Nevertheless, the results by Nyquist et al. (2001) suggest that the Mars ejecta can be sometimes accelerated to the escape velocity and above without sufficient shock-heating. Fur...

    work page 2017

  11. [11]

    From an orbital mechanics viewpoint, Mars ejecta need to have at least an initial ejection velocity higher than 3.8 km/s to reach the orbits of the Martian moons

    Sterilization by Aerodynamic Heating Mars ejecta produced by hypervelocity impacts on the Martian surface must undergo aerodynamic deceleration and consequent aerodynamic heating along their flight trajectory in Martian atmosphere. From an orbital mechanics viewpoint, Mars ejecta need to have at least an initial ejection velocity higher than 3.8 km/s to r...

    work page 2005

  12. [12]

    In this analysis, Mars ejecta are assumed to be in a spherical form, having diameters ranging from 0.01 to 10 m

    At first, a representative aerodynamic drag coefficient and a distribution of aerodynamic heat transfer rate around Mars ejecta are estimated by CFD calculations. In this analysis, Mars ejecta are assumed to be in a spherical form, having diameters ranging from 0.01 to 10 m. Atmospheric temperature and pressure are assumed as 210 K and 0.8 kPa, respective...

    work page 2006

  13. [13]

    Time history of aerodynamic heat transfer rate at the forebody stagnation point is determined by Tauber’s formula (Tauber et al., 1990; Tauber and Sutton,

    is used. Time history of aerodynamic heat transfer rate at the forebody stagnation point is determined by Tauber’s formula (Tauber et al., 1990; Tauber and Sutton,

    work page 1990

  14. [14]

    Distributions of aerodynamic heat transfer rate around Mars ejecta are then estimated by using the representative distribution obtained by the CFD analysis

    along the flight trajectory. Distributions of aerodynamic heat transfer rate around Mars ejecta are then estimated by using the representative distribution obtained by the CFD analysis. According to time history of aerodynamic heat transfer rate distributions around Mars ejecta, unsteady heat conduction analysis about Mars ejecta is conducted to obtain ev...

    work page 2005

  15. [15]

    The results of trajectory analysis are summarized in Figs. and 10, in terms of the threshold velocity to reach the Phobos orbit and the corresponding total heat transferred to the stagnation point for Mars ejecta reaching the Phobos orbit, as a function of the ejection angle and the diameter (it should be noted that the ejection angle is denoted as the fl...

    work page 1991

  16. [16]

    Nature, 555, 643–646

    Timing of oceans on Mars from shoreline deformation. Nature, 555, 643–646. https://doi.org/10.1038/nature26144. Collins, G. S., Elbeshausen, D., Davison, T. M., Wünnemann, K., Ivanov, B., and Melosh, H. J.,

    work page doi:10.1038/nature26144

  17. [17]

    https://doi.org/ 10.6084/m9.figshare.3473690.v2

    iSALE-Dellen manual, Figshare. https://doi.org/ 10.6084/m9.figshare.3473690.v2. Collins, G. S. and Melosh, H. J.,

    work page doi:10.6084/m9.figshare.3473690.v2

  18. [18]

    AIAA Paper 2005-1239, 43rd AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada, USA

    Global Reference Atmospheric Model (GRAM) Series for Aeroassist Applications. AIAA Paper 2005-1239, 43rd AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada, USA. Fritz, J., Greshake, A., and Fernandes, V. A.,

    work page 2005

  19. [19]

    AIAA Paper 2012-3001, 43rd Thermophysics Conference, New Orleans, Louisiana, USA

    Prediction of Forebody and Aftbody Heat Transfer Rate for Mars Aerocapture Demonstrator. AIAA Paper 2012-3001, 43rd Thermophysics Conference, New Orleans, Louisiana, USA. Genda, H., Kurosawa, K., and Okamoto, T., Hydrocode modeling of the spallation process during hypervelocity oblique impacts. in prep. Gilichinsky, D. A., Wilson, G. S., Friedmann, E. I.,...

    work page 2012

  20. [20]

    Earth and Planetary Science Letters, 185, 1–5

    Detecting pyrolysis products from bacteria on Mars. Earth and Planetary Science Letters, 185, 1–5. https://doi.org/10.1016/S0012-821X(00)00370-8. Glavin, D. P., Cleaves, J., Schubert, M., Aubrey, A., and Bada, J. L.,

    work page doi:10.1016/s0012-821x(00)00370-8

  21. [21]

    AIAA Paper 2005-4106, 41st AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit, Tucson, Arizona, USA

    Atmospheric Models for Mars Aerocapture. AIAA Paper 2005-4106, 41st AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit, Tucson, Arizona, USA. Klein, H. P.,

    work page 2005

  22. [22]

    Icarus, 34, 666–674

    The Viking biological experiments on Mars. Icarus, 34, 666–674. https://doi.org/10.1016/0019-1035(78)90053-2. Kurosawa, K. and Genda, H.,

    work page doi:10.1016/0019-1035(78)90053-2

  23. [23]

    Icarus, 301, 219–234

    Hydrocode modeling of the spallation process during hypervelocity impacts: Implications for the ejection of Martian meteorites. Icarus, 301, 219–234. https://doi.org/10.1016/j.icarus.2017.09.015. 23 Lester, E. D., Satomi, M., and Ponce, A.,

    work page doi:10.1016/j.icarus.2017.09.015 2017

  24. [24]

    Meteoritics and Planetary Science, 42(12), 2079–2098

    A hydrocode equation of state for SiO2. Meteoritics and Planetary Science, 42(12), 2079–2098. https://doi.org/10.1111/ j.1945-5100.2007.tb01009.x. Melosh, H. J.,

    work page arXiv 2079

  25. [25]

    Analytical Chemistry, 79, 7249-7256

    The Viking GC/MS and the Search for Organics on Mars. Analytical Chemistry, 79, 7249-7256. https://doi.org/10.1021/ac071972t. Navarro-González, R., Rainey, F. A., Molina, P., Bagaley, D. R., Hollen, B. J., de la Rosa, J., Small, A. M., Quinn, R. C., Grunthaner, F. J., Cáceres, L., Gomez-Silva, B., and McKay, C. P.,

    work page doi:10.1021/ac071972t

  26. [26]

    Astrobiology, 19, 655–668

    A maximum subsurface biomass on Mars from untapped free energy: CO and H2 as potential antibiosignatures. Astrobiology, 19, 655–668. https://doi.org/10.1089/ast.2018.1835. Summers, D.,

    work page doi:10.1089/ast.2018.1835 2018

  27. [27]

    Meteoritics & Planetary Sciences, 52, 2491‒2504

    Mineralogical study of brown olivine in Northwest Africa 1950 shergottite and implications for the formation mechanism of iron nanoparticles. Meteoritics & Planetary Sciences, 52, 2491‒2504. Tauber, M. E., Bowles, J. V., and Yang, L.,

    work page 1950