Assessment of the probability of microbial contamination for sample return from Martian moons II: The fate of microbes on Martian moons
Pith reviewed 2026-05-24 20:01 UTC · model grok-4.3
The pith
The probability of microbial contamination in samples returned from Martian moons is two orders of magnitude below the COSPAR criterion.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Microbes transported to the Martian moons experience severe radiation sterilization, with only 1 ppm surviving on Phobos and 100 ppm on Deimos after accounting for dispersal and shielding in ejecta deposits; this leads to a most likely contamination probability two orders of magnitude lower than the COSPAR limit for sample return missions.
What carries the argument
Impact-transported microbe distribution on moon surfaces combined with depth-dependent radiation sterilization in regolith layers.
If this is right
- Sample return missions collecting 30 g at 10 cm depth qualify as Unrestricted Earth-Return.
- Microbe concentration varies irregularly horizontally but decreases sharply with depth due to radiation.
- The total number of surviving microbes in thick ejecta is 3-4 orders lower than in Mars rock craters.
- Statistical analysis including parameter uncertainties confirms the low contamination probability.
Where Pith is reading between the lines
- Missions could reduce risk further by choosing sampling sites away from crater floors or by going deeper than 10 cm.
- The same transport-plus-sterilization logic could be applied to assess contamination risks for other airless bodies receiving ejecta from habitable worlds.
- Direct laboratory tests of microbial survival under Phobos-like radiation spectra and thermal cycling would tighten the input values used here.
Load-bearing premise
That the microbes originated primarily from the Zunil crater and that laboratory sterilization rates apply directly to the radiation and temperature conditions on the surfaces of Phobos and Deimos.
What would settle it
An experiment exposing relevant microbes to the combined radiation flux, temperature cycles, and vacuum conditions measured or modeled for Phobos and Deimos over multi-year timescales and counting viable survivors.
Figures
read the original abstract
This paper presents a case study of microbe transportation in the Mars-satellites system. We examined the spatial distribution of potential impact-transported microbes on the Martian moons using impact physics by following a companion study (Fujita et al.). We used sterilization data from the precede studies. We considered that the microbes came mainly from the Zunil crater on Mars. We found that 70-80% of the microbes are likely to be dispersed all over the moon surface and are rapidly sterilized due to radiation except for those microbes within a thick ejecta deposit produced by meteoroids. The other 20-30% might be shielded from radiation by thick regolith layers that formed at collapsed layers in craters produced by Mars rock impacts. The total number of potentially surviving microbes at the thick ejecta deposits is estimated to be 3-4 orders of magnitude lower than at the Mars rock craters. The microbe concentration is irregular in the horizontal direction and is largely depth-dependent due to the radiation sterilization. The surviving fraction of transported microbes would be only 1 ppm on Phobos and 100 ppm on Deimos, suggesting that the transport processes and radiation severely affect microbe survival. The microbe sampling probability from the Martian moons was also investigated. We suggest that sample return missions from the Martian moons are classified into Unrestricted Earth-Return missions for 30 g samples and 10 cm depth sampling, even in our conservative scenario. We also conducted a full statistical analysis for sampling the regolith of Phobos to include the effects of uncertainties in input parameters on the sampling probability. The most likely probability of microbial contamination for return samples is estimated to be two orders of magnitude lower than the criterion defined by the planetary protection policy of the Committee on Space Research (COSPAR).
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper models the transport of microbes from Mars (primarily Zunil crater ejecta) to Phobos and Deimos via impacts, their spatial distribution using physics from a companion study, and their survival after applying lab-derived sterilization rates. It partitions 70-80% of microbes into surface layers subject to rapid radiation sterilization and 20-30% into shielded thick deposits, yielding surviving fractions of 1 ppm on Phobos and 100 ppm on Deimos. A statistical analysis of sampling probabilities for regolith leads to the claim that contamination risk for 30 g samples at 10 cm depth is two orders of magnitude below the COSPAR criterion, supporting classification as Unrestricted Earth-Return missions.
Significance. If the survival fractions and resulting probabilities hold, the work supplies a quantitative, uncertainty-aware assessment that could justify relaxed planetary protection requirements for Martian moon sample returns, facilitating mission planning. The use of impact physics from the companion paper and incorporation of parameter uncertainties via full statistical analysis are strengths. The result is directly relevant to COSPAR policy but hinges on the transferability of external sterilization datasets.
major comments (2)
- [Abstract; sterilization and survival modeling sections] Abstract and the section deriving surviving fractions: the 1 ppm (Phobos) and 100 ppm (Deimos) values are obtained by applying sterilization rates from prior laboratory studies directly to moon-surface conditions after partitioning into surface (70-80%) and shielded (20-30%) populations. No adjustment, scaling, or validation is described for differences in radiation spectrum (solar UV + GCR), vacuum, or diurnal temperature cycles versus lab conditions. This assumption directly sets the central survival fractions and is load-bearing for the two-order-of-magnitude margin below the COSPAR criterion; a one-order-of-magnitude change in inactivation rate would eliminate the margin.
- [Sampling probability and statistical analysis] The section on sampling probability and statistical analysis: the conclusion that missions qualify as Unrestricted Earth-Return rests on the above survival fractions. Because those fractions lack demonstrated applicability to Phobos/Deimos environments, the probabilistic claim and policy recommendation require either explicit sensitivity analysis to inactivation-rate uncertainty or additional justification from moon-specific data.
minor comments (2)
- [Abstract] The abstract states 'We used sterilization data from the precede studies' without citing the specific references or tabulating the rates used; add explicit citations and a summary table of input sterilization parameters.
- [Methods] Notation for surviving fractions (ppm) and the partitioning percentages (70-80%, 20-30%) should be defined with equations or a methods subsection to improve traceability.
Simulated Author's Rebuttal
We thank the referee for the constructive comments on our manuscript. We address the major comments below and will revise the paper to include sensitivity analyses on inactivation rates to strengthen the robustness of the survival fractions and sampling probabilities.
read point-by-point responses
-
Referee: [Abstract; sterilization and survival modeling sections] Abstract and the section deriving surviving fractions: the 1 ppm (Phobos) and 100 ppm (Deimos) values are obtained by applying sterilization rates from prior laboratory studies directly to moon-surface conditions after partitioning into surface (70-80%) and shielded (20-30%) populations. No adjustment, scaling, or validation is described for differences in radiation spectrum (solar UV + GCR), vacuum, or diurnal temperature cycles versus lab conditions. This assumption directly sets the central survival fractions and is load-bearing for the two-order-of-magnitude margin below the COSPAR criterion; a one-order-of-magnitude change in inactivation rate would eliminate the margin.
Authors: We agree that the model applies laboratory-derived sterilization rates directly without explicit scaling or validation for Phobos/Deimos-specific conditions such as radiation spectrum details or temperature cycles. These rates represent the best available data from prior studies used in planetary protection contexts, which include vacuum and radiation exposure. To address the load-bearing nature of this assumption, we will add an explicit sensitivity analysis varying the inactivation rates by factors of 0.1–10 and recompute the surviving fractions (1 ppm and 100 ppm) accordingly. revision: yes
-
Referee: [Sampling probability and statistical analysis] The section on sampling probability and statistical analysis: the conclusion that missions qualify as Unrestricted Earth-Return rests on the above survival fractions. Because those fractions lack demonstrated applicability to Phobos/Deimos environments, the probabilistic claim and policy recommendation require either explicit sensitivity analysis to inactivation-rate uncertainty or additional justification from moon-specific data.
Authors: The existing statistical analysis already propagates uncertainties from multiple input parameters into the sampling probabilities. We will extend this framework to include the inactivation-rate sensitivity analysis described above, quantifying the effect on contamination probabilities for 30 g, 10 cm samples. This will show whether the two-order-of-magnitude margin below the COSPAR criterion holds under rate variations, directly supporting the Unrestricted Earth-Return classification. revision: yes
Circularity Check
No significant circularity detected
full rationale
The paper's derivation applies external sterilization rates from prior laboratory studies and transport results from a cited companion study to compute survival fractions (1 ppm on Phobos, 100 ppm on Deimos) and sampling probabilities. These inputs are independent of the outputs; the statistical analysis propagates uncertainties in those external parameters rather than fitting or redefining them internally. No self-definitional steps, fitted inputs renamed as predictions, or load-bearing self-citations that reduce the central claim to its own assumptions by construction appear in the text. The derivation remains self-contained against the stated external benchmarks.
Axiom & Free-Parameter Ledger
free parameters (2)
- surviving fraction on Phobos
- surviving fraction on Deimos
axioms (2)
- domain assumption Microbes transported to the moons originate primarily from the Zunil crater on Mars
- domain assumption Sterilization rates measured in previous laboratory studies apply to conditions on Phobos and Deimos
Reference graph
Works this paper leans on
-
[1]
timate. Even if we employed a = 4.3, the main conclusion of this study does not largely change because it yields a smaller microbial contamination probability than the value in the case with a = 1.8. 2.3. Radiation-induced sterilization Microorganisms transported to the Martian moons are expected to be sterilized over time by solar and galactic cosmic rad...
work page 2017
-
[2]
The results at four different times are shown
Average survival rate due to radiation as a function of depth. The results at four different times are shown. The times and crater names formed at these times are annotated next to the lines. As in Fig. 4, the target depth for sample collection currently used by the JAXA MMX mission are shown by shading. the PDF of the potential microbe density in the Mar...
work page 2007
-
[3]
and was validated via a comparison with an impact experiment (Okamoto et al., in revision). We focused on the Zunil crater (7.7°N, 166°E; 10.1 km in diameter) on Mars as the main source of the microbes potentially living on the current moons. Zunil is the youngest known ray crater with a diameter of >10 km on Mars. Ray systems around a host crater are con...
work page 2004
-
[4]
Distribution of Mars ejecta fragments over the Martian moons In this section we discuss impact bombardment of the microbe-bearing Mars rocks on each of the moons and their distribution after the bombardment. We assumed that all the Mars rocks are spheres of 0.1 m in diameter, because the size–frequency distribution (SFD) of the ejected Mars rocks is highl...
work page 2004
-
[5]
We randomly assigned vimp and qimp to each microbe-bearing rock from the distributions (Fig
was used to produce the random numbers in the Monte Carlo calculations. We randomly assigned vimp and qimp to each microbe-bearing rock from the distributions (Fig. 7). This treatment allowed us to characterize each impact event. The impact outcomes with a given set of impact conditions can be predicted by using the results from previous studies in the fi...
work page 2016
-
[6]
Distributions of the (a) impact velocity and (b) impact angle of Mars rocks colliding with Phobos (red solid lines) and Deimos (blue dashed lines). 3.1. Crater formation Hypervelocity impacts of Mars rocks onto the surfaces of each moon should produce impact craters (e.g., Melosh, 1989). In general, impact-driven gardening process is rather complex becaus...
work page 1989
-
[7]
and mix with the regolith particles (e.g., Daly and Schultz, 2016; Ebert et al., 2014). After transient crater formation, the wall of the transient crater collapses due to gravity, resulting in granular flow towards the crater center. This is often referred to as the “modification stage” in the cratering process (e.g., Melosh, 1989; Barnouin-Jha et al., 2...
work page 2016
-
[8]
The scaling parameters correspond to the values for dry sand (Schmidt and Housen, 1987)
as follows: Df = 1.25Dtr (9) and Dtr = /π6013CD/4π30– β3-ρpρt.13Dp1–βg–β&vimpsinθimp(34, (10) where CD = 1.4, b = 0.17, rp = 2.7 ´ 103 kg/m3, rt = 2 ´ 103 kg/m3, and g = 0.0057 m/s2 for Phobos and g = 0.003 m/s2 for Deimos are a dimensionless scaling constant, dimensionless scaling exponent, projectile density, target density, and gravitational accelerati...
work page 1987
-
[9]
Schematic cross-section of Mars rock impact craters considered in this study. The black dashed line is the profile of the final crater that actually formed in the regolith. The red shaded region corresponds to the collapsed lens where the Mars rock fragments mix with the moon regolith particles. 3.2. Scattered fragments In the previous section, we modeled...
work page 2016
-
[10]
are given by Mdis,p= (1−ψ) Mp, (11) and Mdis,t = CVρtRtr37vesc8gRtr9–3µ (12) where Cv = 0.32, Rtr = 0.5Dtr, and µ = 0.4 are a scaling constant, transient crater radius, and velocity-scaling exponent, respectively, and vesc is the escape velocity from each of the moons (11 m/s for Phobos and 6.9 m/s for Deimos). Equation (12) is derived from the point-sour...
work page 1982
-
[11]
because the moons sweep out their own dust torus. It is difficult to accurately estimate the timescale of re-accumulation without detailed numerical simulations considering the ejection process and orbital evolution of the ejected materials. Nevertheless, the timescale may be much longer than one orbital period of Phobos and Deimos (~10 and ~30 h, respect...
work page 2018
-
[12]
The Dpmax values for Phobos and Deimos in the past 0.1 Myr were estimated to be 2.6 and 1.6 m, respectively. Although there is the possibility that larger impactors have struck the moons, such events are statistically rare. Given that the total number of impacts of bodies smaller than 0.1 m is unknown, we decided to extrapolate Eq. (14) to smaller diamete...
work page 1996
-
[13]
The gray dashed line is a reference power law function with an exponent of –2.5
Cumulative size–frequency impactor distributions (SFDs) pertaining to Mars, Phobos, and Deimos. The gray dashed line is a reference power law function with an exponent of –2.5. The corresponding crater diameters on Mars are shown on the upper x-axis. Given that the cumulative number of impactors N smaller than 0.1 m is unknown, this region is highlighted ...
work page 2017
-
[14]
for both Phobos and Deimos. To investigate the significance of the ejecta deposit covering the thin microbe-bearing layer, knowledge of the rate of decrease of the ejecta deposit thickness with increasing distance along the surface from the host crater is necessary. If the thin microbial layer was covered by an ejecta deposit with a thickness of >3 mm wit...
work page 1983
-
[15]
We used Dp = 0.001, 0.01, 0.1, and 1 m in these calculations
Ejecta thickness as a function of distance from the impact point along the surface of Phobos. We used Dp = 0.001, 0.01, 0.1, and 1 m in these calculations. The grey horizontal line corresponds to an ejecta thickness of 3 mm. probability to the shielded microbial layer Player, where Smoon is the surface area of each moon (1.5 ´ 109 m2 for Phobos and 5.0 ´ ...
work page 2000
-
[16]
performed impact experiments at two different impact angles (40° and 90°) measured from the target surface, and no systematic differences in ξ were observed between the two impact angles. Thus, the average survival fraction ξave was simply obtained by a convolution of the vimp distribution with the survival rate ξ given by Eq. (3). We found that ξave of t...
work page 2018
-
[17]
Thirdly, we consider the covered microbial thin layer. The microbial column density immediately after the formation of the thin global layer σthin0 can be estimated as follows: σthin0 = ξavenMR&1−ψave(Mp,total,aveSmoon. (23) Our model yields σthin0 values on Phobos and Deimos of 2.9 ´ 10–5 and 3.3 ´ 10–5 cells/cm2, respectively, in the fiducial case. The ...
work page 2018
-
[18]
The sum of Ps,crater and Ps,layer is Ps
The sampling probabilities of the microbes from the Mars rock craters and covered microbial thin layer Ps,crater and Ps,layer are respectively given by Ps,crater = Pcraterη(t, H)ncrater0Ms, (28) and Ps,layer = Playerσthin,aveSs, (29) where Ms = ρtSSLs is the sample mass, SS is the sample area, LS is the sample depth, and σthin,ave is the average value of ...
work page 2012
-
[19]
The selected sampling masses are annotated beside each line
Sampling probabilities of microbes from the surfaces of (a) Phobos and (b) Deimos as functions of sampling depth, area, and mass. The selected sampling masses are annotated beside each line. The REQ-10 criterion is shown as a grey dashed horizontal line. The curves are the globally averaged probabilities including the Mars rock craters. The straight lines...
work page 2010
-
[20]
13 except that the calculated results are for an undiscovered crater
Same as Fig. 13 except that the calculated results are for an undiscovered crater. The calculated result in the case of t > 2 kyr (Eq. (31-2)) is shown here. (a) Ms = 60 g and (b) Ms = 100 g. The formation ages are shown beside the lines. 4.4.1. Secondly, we discuss the effects of the sizes of the Martian rocks. Although all the microbe-bearing Mars rocks...
work page 2004
-
[21]
The sampling mass was set to 30 g
Probability density function (PDF) of the microbial contamination probability Ps. The sampling mass was set to 30 g. We calculated Ps for four different sampling depths (i.e., 1, 3, 6, and 10 cm). We obtained the PDF via 100,000 Monte Carlo runs. and Ls = 10 cm. Table 1 lists the expected values of Ps as a function of sampling depth Ls. Although Ps increa...
work page 2004
-
[22]
Conclusions We investigated the surface processes on the moons Phobos and Deimos after the potential transportation and arrival of Martian microbes from a young ray crater (Zunil) on Mars by using the results of Paper 1 (Fujita et al., in submission) as the initial conditions. In this study, we used the orbital parameters of Mars rocks, including impact v...
work page 2010
-
[23]
Planetary and Space Science, 125, 20-26
Survivability of bare, individual Bacillus subtilis spores to high-velocity surface impact: Implications for microbial transfer through space. Planetary and Space Science, 125, 20-26. Barnouin-Jha, O.S., Yamamoto, S., Toriumi, T., Sugita, S., and Matsui, T. (2007) Non-intrusive measurements of crater growth. Icarus 188, 506-521. Bogard, D.D. and Johnson, ...
work page 2007
-
[24]
Survival of bacteria and spores under extreme shock pressures. Mon. Not. R. Astron. Soc., 352, 1273-1278. Chappaz, L., Melosh, H.J., Vaquero, M., and Howell, K.C. (2013) Transfer of impact ejecta material from the surface of Mars to Phobos and Deimos. Astrobiology 13, 963-980. Daly, R.T. and Schultz, P.H. (2016) Delivering a projectile component to the ve...
work page doi:10.1063/1 2013
-
[25]
Genda, H., Kurosawa, K., Okamoto, T., Hydrocode modeling of the spallation process during hypervelocity oblique impacts. in prep. Gillon M., Jehin, E., Lederer, S.M., Delrez, L., de Wit, J., Burdanov, A., Van Grootel, V., Burgasser, A.J., Triaud, A.H.M.J., Opitom, C., Demory, B.-O., Sahu, D.K., Gagliuffi, D.B., Magain, P., and Queloz, D. (2016) Temperate ...
work page 2016
-
[26]
(2005) Martian cratering 8: Isochron refinement and the chronology of Mars
Hartmann, W.K. (2005) Martian cratering 8: Isochron refinement and the chronology of Mars. Icarus 174, 294-320. Hartmann, W.K., Quantin, C., Werner, S.C., and Popova, O. (2010) Do young martian ray craters have ages consistent with the crater count system?. Icarus 208:621-635. Hawke, B.R., Blewett, D.T., Lucey, P.G., Smith, G.A., Bell, J.F., Campbell, B.A...
work page 2005
-
[27]
doi:10.3389/fmicb. 2014.00612. Hazell, P.J., Beveridge, C., Groves, K., Appleby-Thomas, G. (2010) The shock compression of microorganism-loaded broths and emulsions: Experiments and simulations. Int. J. Impact Eng. 37, 433–440. Head, J.N., Melosh, H.J., Ivanov, B.A. (2002) Martian meteorite launch: High-speed ejecta from small craters. Science 298, 1752–1...
-
[28]
On the impact origin of Phobos and Deimos. IV. Volatile depletion. The Astrophysical Journal 860, 150(10pp). Hyodo, R., Kurosawa, K., Genda, H., Fujita, K., and Usui, T. Extensive delivery of Martian ejecta to its moons: the gateway to a time-resolved history of Mars. in prep. Ito, T. and Malhotra, R. (2006) Dynamical transport of asteroid fragments from ...
work page 2006
-
[29]
Physics of the Earth and Planetary Interiors, 129, 131-143
The phase diagram of CaCO3 in relation to shock compression and decomposition. Physics of the Earth and Planetary Interiors, 129, 131-143. Kawakatsu, Y., Kuramoto, K., Ogawa, N., Ikeda, H., Mimasu, Y., One, G., Sawada, H., Yoshikawa, K., Imada, Takane, Otake, Hisashi, Kusano, H., Yamada, K., Otsuki, M., and Baba, M. (2017) Mission concept of Martian Moons...
work page 2017
-
[30]
The rayed crater Zunil and interpretations of small impact craters on Mars. Icarus, 176, 351–381. McKinnon, W.B., Chapman, C.R., and Housen, K.R. (1991) Cratering of the uranian satellites. In: Bergstralh, J.T., Miner, L.D., Matthews, M.S. (Eds.), Uranus. Univ. of Arizona Press, Tucson, pp. 629–692. Melosh, H.J. (1984) Impact ejection, spallation, and the...
-
[31]
Test report on the irradiation inactivation tests results. SterLim-OU-TN15. Patel, M., Pearson, V., Summers, D., Evans, D., Bennet, A., and Truscott, P. (2018) Sterilization Limits for Sample Return Planetary Protection Measures (SterLim),” Presentation to the Committee on the Review of Planetary Protection Requirements for Sample Return from Phobos and D...
work page 2018
-
[32]
Experimental constraints on impact-induced winds. Icarus, 305, 91-104. Ramsley, K. R. and Head, J. W. (2013) Mars impact ejecta in the regolith of Phobos: Bulk concentration and distribution. Planetary and Space Science 87, 115-129. Ramsley, K. R. and Head, J. W. (2017) The Stickney crater ejecta secondary impact crater spike on Phobos: Implications for t...
work page 2013
-
[33]
Journal of Geophysical Research 117, E05004, doi:10.1029/2011JE003966
Database creation, properties, and parameters. Journal of Geophysical Research 117, E05004, doi:10.1029/2011JE003966. Schmedemann, N., Michael, G.G., Ivanov, B.A., Murray, J.B., and Neukum, G. (2014) The age of Phobos and its largest crater, Stickney. Planetary and Space Science 102, 152-163. Schultz, P. H.,
-
[34]
Impact-generated winds on Mars. 292, 86-101. Schmidt, R.M. and Housen, K.R. (1987) Some recent advances in the scaling of impact and explosion cratering. International Journal of Impact Engineering 5:543-560. Shoemaker, E.M. (1962) Interpretation of lunar craters. in Physics and Astronomy of the Moon, edited by Z. Kopal, pp. 283–359, Academic, San Diego,C...
-
[35]
Bugbuster-survivability of living bacteria upon shock compression. Earth Planet. Sci. Lett. 247, 185–196. http://dx.doi.org/10.1016/j.epsl.2006.03.054. Zahnle, K., Schenk, P., Levison, H. and Dones, L. (2003) Cratering rates in the outer solar system. Icarus, 163(2), 263 – 289, doi:http://dx.doi.org/10.1016/S0019-1035(03)00048-4
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.