REVIEW 3 major objections 4 minor 35 references
The dynamical surface of Phobos: a morphodynamic atlas
T0 review · 3 major / 4 minor · reviewed 2026-08-02 · deepseek-v4-flash
Pith's one-line read Most of Phobos's surface morphology and its red-blue spectral contrast can be explained as the long-term work of regolith moving along preferred tidal-gravity pathways.
desk verdict A genuinely useful first atlas of Phobos regolith pathways, honestly caveated, but the headline spectral/morphological link depends on a single low-friction threshold that the paper itself cannot justify. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The Regolith Migration Pathway (RMP) — a preferred transport route obtained by integrating surface-constrained trajectories under the local acceleration field with a Coulomb-type friction mobility threshold. The acceleration field includes self-gravity, Martian tides, centrifugal and Coriolis terms; the threshold is explored at 14° (low, revealing past or externally triggered transport) and 30° (high, conservative spontaneous case). The key move is projecting acceleration onto the local tangent plane and using it to route massless tracers, with friction treated as a mobility criterion rather than as a measured material property.
What would settle it
Measure or infer Phobos's effective regolith friction (for example, from MMX lander mechanics or by matching crater-density maps) and test whether predicted RMP endpoints are depleted in small craters. A null result — either an effective friction near 30°, or equal crater densities at endpoints and in quiet areas — would collapse the central claim.
Extended reading notes
Core claim
The central discovery is a network of Regolith Migration Pathways (RMPs) that connects Phobos's morphology to its dynamics. Simulating surface-constrained trajectories under self-gravity, tidal, centrifugal, and Coriolis accelerations with a Coulomb-type friction threshold of 14°, the paper finds that most of the surface is quiet but a sparse set of corridors channels material over kilometre scales. Trajectories from Stickney's eastern outer slope and mid-latitude terrain collect near the sub-Mars point; another dense cluster accumulates in the anti-Mars southern troughs. These termination zones are exactly where images show smooth, mantled, lightly cratered, spectrally neutral surfaces, whe
Load-bearing premise
The load-bearing premise is that a 14° effective friction/mobility threshold reveals the transport corridors that actually shaped Phobos; if the relevant effective resistance were instead close to 30°, the predicted pathways and all their correlations with smooth terrain and spectral units would largely disappear.
Editorial extensions
If this is right
- The sub-Mars MMX sampling site is predicted to contain well-mixed, mature regolith fed from Stickney's outer slope and southern mid-latitudes, while the anti-Mars site is predicted to sample younger, less-weathered material from active troughs.
- Smooth, spectrally neutral, lightly cratered terrain on Phobos marks depositional mantles formed by regolith accumulation along RMP endpoints; rough, blue-spectral terrain marks actively denuded source regions.
- Red-spectral, rough terrain corresponds to dynamically quiet surfaces where negligible regolith motion has occurred, implying they are older and less frequently reworked.
- The 30° result confines present-day spontaneous mobility to steep crater walls and Stickney's outer slope, identifying those as the most plausible sites of ongoing activity.
- A future high-resolution crater-density map can independently test the model: endpoints should be depleted in small craters, and source regions enriched.
Reading between the lines
- Beyond Phobos, the same approach should reveal comparable transport corridors on any tidally locked small body whose mean slope is below its angle of repose; the effective-threshold method sidesteps unknown paleo-topography.
- The correlation between RMPs and grooves along Stickney's outer slope suggests those grooves may be erosional records of long-term regolith migration rather than purely tectonic fractures — a connection the paper raises but leaves open.
- If the true effective friction on Phobos turns out to be close to 30° rather than 14°, the geometric structure of the RMPs would persist but the interpretation of red/blue units as dynamical stratigraphy would need revision.
- A testable extension: run the same trajectory model with a time-varying tidal field, including libration and past eccentricity, to see whether endpoint clustering sharpens or shifts — the paper notes libration is unmodeled.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. Herreros and Charnoz use the RAVEL code to compute the surface acceleration field on a 36-m DTM of Phobos, integrate surface-constrained regolith trajectories under self-gravity, centrifugal, tidal, Coriolis and Coulomb friction, and define Regolith Migration Pathways (RMPs) for two end-member friction angles (30° and 14°). They compare the 14° network with HRSC/HiRISE images and spectral slope maps, reporting qualitative agreement between RMP termini and smooth, spectrally neutral terrains, and between denuded/rough terrains and blue/red units. They interpret the sub-Mars and anti-Mars sites as distinct source/sink regions and use this to frame MMX sample provenance.
Significance. The paper is the first global dynamical transport-atlas calculation for Phobos and provides a falsifiable framework for MMX sample interpretation. Strengths include a clearly documented acceleration model with explicit treatment of tides and eccentricity, sensitivity checks on eccentricity and detachment, and a candid discussion of limitations. The main value is the proposed source-to-sink architecture, which can be tested with MMX data. However, as argued below, the central causal claim currently rests on a single low-friction threshold and on qualitative visual correlations, so the result is promising rather than established.
major comments (3)
- [§4, §5.1.2, App. B.3] The load-bearing threshold choice is not defended against the static/dynamic split. The governing equations (Eq. B.5 and B.6 in App. B) distinguish φ_static and φ_dynamic and explicitly note that dynamic friction is commonly lower than static (App. B.3), yet the 14° case sets φ_static=φ_dynamic=14° (§5.1.2). This means trajectories initiate wherever the present-day dynamical slope exceeds the mean slope (~14°). The paper’s own statement that it does not predict triggering (§1, §4) is in tension with this initiation rule. The skeptic’s scenario—static ≈30° (as in §5.1.1, from [22]) with dynamic ≈14°—would restrict initiation to steep walls and likely erase the RMP correlations. Because the abstract’s causal claim is supported only by the 14° case (§6), please test this split and report whether the correlations survive; if not, the interpretation should be reframed as a diagnostic atlas ra
- [§5.2, Fig. 12] The claimed correlations are qualitative and not quantitatively validated. Sections 5.2.1–5.2.4 rely on visual inspection ('remarkably well', 'good qualitative agreement'), and Fig. 12 compares endpoint clusters with the Basilevsky feature map using different projections without a statistical overlap test. The paper admits (§6) that correlations are reproduced only for φ=14°, but no null or control test is provided (e.g., random endpoints, φ=30° endpoints, or neighboring thresholds). Please add a quantitative metric—endpoint density vs smooth/crater-poor/neutral units, with significance relative to random distributions—and a threshold sweep around 14° to establish that the network is not an artifact of choosing a threshold equal to the mean slope.
- [§6.1, §5.2.3, §5.2.4] The spectral-interpretation claim ('dynamical stratigraphy') exceeds the available evidence. High-resolution spectral data are absent for the anti-Mars region (§5.2.3) and the South Pole correlation is explicitly weak (§5.2.4); moreover, the model cannot distinguish composition from maturity. The manuscript should either soften the claim that 'much of Phobos' surface morphology and spectral heterogeneity can be explained by long-term regolith redistribution' or provide a more explicit mapping of where the correlation is and is not tested.
minor comments (4)
- [Fig. 2] Add axis labels and units; the mean-slope line is useful but the histogram coordinates should be self-contained.
- [App. C.1.4] The weights w1=w2=w3=0.33 are chosen without sensitivity analysis; state that the mobility index is not used in the main RMP calculation, to avoid overinterpretation.
- [App. D] The normal coefficient of restitution is fixed at 0.1; a short sensitivity note would help.
- [§1, App. B.4] The phrase 'the model does not aim to predict triggering' should be reconciled with the θ>φ_static initiation rule; consider renaming φ_static to 'initiation threshold' and clarifying that it is a modeling parameter, not a physical static friction angle.
Circularity Check
No significant circularity: RMPs derive from independent force-field/topography inputs; the 14° threshold is an acknowledged modeling end-member, not a parameter fitted to the observed correlations.
full rationale
The paper's derivation chain is not circular by construction. RMPs are computed from the surface acceleration field (self-gravity, centrifugal, tidal, Coriolis) on a 36 m DTM, using a Coulomb mobility threshold (Appendix B). The observed smooth terrains, crater densities, and red/blue spectral units are external datasets (HiRISE, HRSC, spectral slope maps, Basilevsky et al. 2014) and are not used as inputs to the trajectory calculation. The friction angle is not fitted to these observables: the paper states it explored a suite of friction angles (14°, 20°, 25°, 30°, 35°) and found that trajectory geometry stays consistent while lengths decrease with friction, and it presents 14° as a low end-member close to the mean dynamical slope, not as a measured material property. The central comparison is therefore a posteriori consistency testing, not a renamed fit. The paper explicitly disclaims that the model predicts triggering: 'The model does not aim to predict the triggering of slope failure. Instead, it addresses where material would preferentially move once motion is initiated.' This removes the most obvious circularity. The acknowledged limitation that correlations are reproduced mainly for φ=14°, while the φ=30° case shows no such correspondence, is a model-selection/robustness concern rather than a definitional reduction: the RMP geometry is not logically forced to match the observed units by the choice of threshold. No load-bearing self-citation chain appears; RAVEL is described in the appendices, and external constraints such as Robin et al. (2024) and Ballouz et al. (2019) are used for context. The paper also acknowledges the paleotopography feedback (progressive infill reduces slopes) and frames its low-threshold case as a late-stage pathway-revealing end-member, which is an honest equifinality caveat rather than a circular derivation. Overall, no step reduces to its own inputs by construction.
Assumptions & free parameters
free parameters (3)
- Effective friction angle φ (low/high end-members) =
14°, 30°; 20°, 25°, 35° also explored
- Mobility index weights w1, w2, w3 =
0.33 each
- Normal coefficient of restitution =
0.1
assumptions (6)
- domain assumption Homogeneous interior and uniform mascon gravity
- domain assumption Synchronous rotation; libration neglected
- domain assumption Coulomb friction with static=dynamic for end-members and no cohesion
- domain assumption Surface-constrained motion (no detachment) for main RMPs
- domain assumption Present-day topography used as proxy for past transport
- domain assumption Motion initiation only where dynamical slope exceeds static friction; triggering not modeled
Cite this review
Pith. "Pith review of The dynamical surface of Phobos: a morphodynamic atlas." pith.science (2026). https://pith.science/paper/TNY5J6ZN
@misc{pith2026260710445,
author = {Pith},
title = {Pith review of: The dynamical surface of Phobos: a morphodynamic atlas},
year = {2026},
howpublished = {\url{https://pith.science/paper/TNY5J6ZN}},
note = {Machine review of arXiv:2607.10445}
}
read the original abstract
Phobos evolves in a highly dynamical environment where surface-material motion is controlled by the combined effects of self-gravity, time-dependent Martian tides, and inertial forces. In such a low-gravity regime, the displacement of loose material, cannot be inferred from topographic slope alone, making a dynamical approach essential for interpreting Phobos' surface morphology and for supporting the Martian Moons eXploration (MMX) mission led by JAXA. Here, using our RAVEL code, we apply a dynamical model that combines the surface acceleration field with friction on a digital terrain model of Phobos to compute surface regolith trajectories. The model does not aim to predict the triggering of slope failure. Instead, it addresses where material would preferentially move once motion is initiated. This reveals large scale coherent dynamical regions and a sparse network of preferred regolith transport routes, termed here Regolith Migration Pathways (RMPs). The final positions of the RMPs correlate with smooth, low-relief terrains and spectrally neutral units, consistent with depositional mantles formed by long-term regolith infill, whereas rough, high-standing areas with abundant small craters and blue spectral slopes tend to correspond to dynamically active or denuded source regions. In contrast, spectrally red terrains are generally associated with dynamically quiet, morphologically rough surfaces where our model predicts negligible regolith motion, suggesting older, less frequently reworked units. Taken together, these patterns indicate that much of Phobos' surface morphology and spectral heterogeneity can be explained by long-term regolith redistribution driven by the surface acceleration field along RMPs. We provide a 3D morphodynamic atlas of RMPs across Phobos' surface, which will be useful for constraining the geographical provenance of samples to be collected by the MMX spacecraft.
Figures
Figures from the paper (16 more)
Reference graph
Works this paper leans on
-
[22]
Robin, C.Q., Duchêne, A., Murdoch, N., etal., 2024.Mechanicalpropertiesofrubblepileaster- oids (dimorphos, itokawa, ryugu, and bennu) through surface boulder morphological analysis. Nat. Commun. 15, 6203. https://doi.org/10.1038/s41467-024-50147-w
-
[1]
T., Kawakatsu, Y., Fujimoto, M., 2019
Ballouz, R.-L., Baresi, N., Crites, S. T., Kawakatsu, Y., Fujimoto, M., 2019. Surface refreshing of martian moon phobos by orbital eccentricity-driven grain motion. Nat. Geosci. 12, 229–234. https://doi.org/10.1038/s41561-019-0323-9
-
[2]
Barnouin, O. S., Daly, M. G., Palmer, E. E., et al., 2019. Shape of (101955) bennu indicative of a rubble pile with internal stiffness. Nat. Geosci. 12, 247–252. https://doi.org/10.1038/s41561- 019-0330-x 26 ACCEPTED by Earth and Planetary Science Letters, July 2026
-
[3]
Basilevsky, A. T., Lorenz, C. A., Shingareva, T. V., Head, J. W., Ramsley, K. R., Zubarev, A. E., 2014. The surface geology and geomorphology of phobos. Planet. Space Sci. 102, 95–118. https://doi.org/10.1016/j.pss.2014.04.013
-
[4]
Effect of friction on the angle of repose of elongated particles
Cárdenas-Barrantes, M., Ovalle, C., 2025. Effect of friction on the angle of repose of elongated particles. Powder Technol. 458, 120974. https://doi.org/10.1016/j.powtec.2025.120974
arXiv 2025
-
[5]
D., Olhoeft, G
Carrier, W. D., Olhoeft, G. R., Mendell, W., 1991. Physical properties of the lunar surface. In: Heiken, G. H., Vaniman, D. T., French, B. M. (Eds.), Lunar Sourcebook. Cambridge University Press, pp. 475–594
1991
-
[6]
Ernst, C. M., Daly, R. T., Gaskell, R. W., Barnouin, O. S., Nair, H., Hyatt, B. A., Al Asad, M. M., Hoch, K. K. W., 2023. High-resolution shape models of phobos and deimos from stereophotoclinometry. Earth Planets Space 75, 103. https://doi.org/10.1186/s40623-023- 01814-7
-
[7]
H., Tirsch, D., Matz, K.-D., Doressoundiram, A., Gautier, T., Barucci, M
Fornasier, S., Wargnier, A., Hasselmann, P. H., Tirsch, D., Matz, K.-D., Doressoundiram, A., Gautier, T., Barucci, M. A., 2024. Phobos photometric properties from mars express HRSC observations. Astron. Astrophys. 686 (A203). https://doi.org/10.1051/0004-6361/202449220
Show all 35 references
-
[8]
A., Arvidson, R
Fraeman, A. A., Arvidson, R. E., Murchie, S. L., Rivkin, A. S., 2012. Analysis of disk-resolved OMEGA and CRISM spectral observations of phobos and deimos. J. Geophys. Res. Planets 117, E00J15. https://doi.org/10.1029/2012JE004137
2012 doi
-
[9]
A., Murchie, S
Fraeman, A. A., Murchie, S. L., Arvidson, R. E., et al., 2014. Spectral absorptions on pho- bos and deimos in the visible/near infrared wavelengths and their compositional constraints. ICARUS 229, 196–205. https://doi.org/10.1016/j.icarus.2013.11.021
2014 doi
-
[10]
A., Lainey, V., 2014
Jacobson, R. A., Lainey, V., 2014. Martian satellite orbits and ephemerides. Planet. Space Sci. 102, 35–44. https://doi.org/10.1016/j.pss.2013.06.003
2014 doi
-
[11]
C., Campbell, C
Johnson, B. C., Campbell, C. S., Melosh, H. J., 2016. The reduction of friction in long runout landslides as an emergent phenomenon. J. Geophys. Res. Earth Surf. 121 (5), 881–
2016
-
[12]
P., Oberst, J., Zubarev, A
Karachevtseva, I. P., Oberst, J., Zubarev, A. E., Nadezhdina, I. E., Kokhanov, A. A., Garov, A. S., Uchaev, D. V., Uchaev, Dm. V., Malinnikov, V. A., Klimkin, N. D., 2014. The phobos information system. Planet. Space Sci. 102, 74–85. https://doi.org/10.1016/j.pss.2013.12.015
2014 doi
-
[13]
Signature of Phobos’ interior structure in its gravity field and libration
Le Maistre, S., Rivoldini, A., Rosenblatt, P., 2019. Signature of Phobos’ interior structure in its gravity field and libration. ICARUS 321, 272–290. https://doi.org/10.1016/j.icarus.2018.11.022
2019 doi
-
[14]
C., Dehant, V., Marty, J.-C.,
Le Maistre, S., Rosenblatt, P., Rambaux, N., Castillo-Rogez, J. C., Dehant, V., Marty, J.-C.,
-
[15]
Dynamics of long-runout landslides: a review
Lei, Z., Mao, W., Yu, F., 2025. Dynamics of long-runout landslides: a review. Appl. Sci. 15 (21). https://doi.org/10.3390/app152111300
2025 doi
-
[16]
P., 2014
Lucas, A., Mangeney, A., Ampuero, J. P., 2014. Frictional velocity-weakening in landslides on Earth and on other planetary bodies. Nat. Commun. 5, 3417. https://doi.org/10.1038/ncomms4417 27 ACCEPTED by Earth and Planetary Science Letters, July 2026
2014 doi
-
[17]
K., Scott, R
Mitchell, J. K., Scott, R. F., Houston, W. N., et al., 1972. Mechanical properties of lunar soil: density, porosity, cohesion, and angle of internal friction. In: Proceedings of the Lunar Science Conference. Vol. 3, pp. 3235–3253
1972
-
[18]
Surface environment of phobos and phobos simulant UTPS
Miyamoto, H., Niihara, T., Wada, K., Ogawa, K., Senshu, H., Michel, P., et al., 2021. Surface environment of phobos and phobos simulant UTPS. Earth Planets Space 73, 214. https://doi.org/10.1186/s40623-021-01406-3
2021 doi
-
[19]
Rotational motion of phobos
Rambaux, N., Castillo-Rogez, J.-C., Le Maistre, S., Rosenblatt, P., 2012. Rotational motion of phobos. Astron. Astrophys. 548, A14. https://doi.org/10.1051/0004-6361/201219710
2012 doi
-
[20]
R., Head, J
Ramsley, K. R., Head, J. W., 2019. Origin of phobos grooves: testing the stickney crater ejecta model. Planet. Space Sci. 165, 137–147. https://doi.org/10.1016/j.pss.2018.11.004
2019 doi
-
[21]
R., Head, J
Ramsley, K. R., Head, J. W., 2021. The origins and geological histories of deimos and phobos: hypotheses and open questions. Space Sci. Rev. 217 (86). https://doi.org/10.1007/s11214-021- 00864-1
2021 doi
-
[23]
Inte- grated method for crater detection from topography and optical images and the new PH9224GT catalogue of Phobos impact craters
Salamunićcar, G., Lončarić, S., Pina, P., Bandeira, L., Saraiva, J., 2014. Inte- grated method for crater detection from topography and optical images and the new PH9224GT catalogue of Phobos impact craters. Adv. Space Res. 53 (12), 1798–1809. https://doi.org/10.1016/j.asr.2013.11.006
2014 doi
-
[24]
J., Hartzell, C
Scheeres, D. J., Hartzell, C. M., Sánchez, P., Swift, M., 2010. Scaling forces to asteroid surfaces: the role of cohesion. ICARUS 210 (2), 968–984. https://doi.org/10.1016/j.icarus.2010.07.009
2010 doi
-
[25]
Evolution of phobos’ orbit, tidal forces, dynamical topography, and related surface modification processes
Shi, X., Willner, K., Oberst, J., 2013. Evolution of phobos’ orbit, tidal forces, dynamical topography, and related surface modification processes. In: 44th Lunar and Planetary Science Conference. Abstract 1889
2013
-
[26]
Mapping the thermal inertia of phobos using MGS-TES observations and thermophysical modeling
Smith, N., Edwards, C., Mommert, M., Trilling, D., Glotch, T., 2018. Mapping the thermal inertia of phobos using MGS-TES observations and thermophysical modeling. In: Lunar and Planetary Science Conference. Vol. 49. LPI Contribution No. 2083
2018
-
[27]
T., Herkenhoff, K
Thomas, N., Stelter, R., Ivanov, A., Bridges, N. T., Herkenhoff, K. E., McEwen, A. S., 2011. Spectral heterogeneity on phobos and deimos: HiRISE observations and comparisons to mars pathfinder results. Planet. Space Sci. 59 (13), 1281–1292. https://doi.org/10.1016/j.pss.2010.04.018
2011 doi
-
[28]
The surfaces of phobos and deimos
Veverka, J., 1978. The surfaces of phobos and deimos. Vistas Astron. 22 (2), 163–192. https://doi.org/10.1016/0083-6656(78)90014-4
1978 doi
-
[29]
A new topographic image atlas of phobos
Wählisch, M., Willner, K., Oberst, J., Matz, K.-D., Scholten, F., Roatsch, T., Hoffmann, H., Semm, S., Neukum, G., 2010. A new topographic image atlas of phobos. Earth Planet. Sci. Lett. 294, 547–553. https://doi.org/10.1016/j.epsl.2009.11.003 28 ACCEPTED by Earth and Planetar...
2010 doi
-
[32]
For each starting pointi, we consider the endpoints of the pericenter and apocenter trajecto- ries
-
[33]
If two distinct endpoints exist, we define their mutual separation in 3D, dee =∥e 1 −e 2∥.(C.10)
-
[34]
We compute the start–to–endpoint distances, d1 =∥e 1 −p 0∥, d 2 =∥e 2 −p 0∥,(C.11) and their mean dse = 1 2 (d1 +d 2).(C.12)
Letp 0 be the initial position. We compute the start–to–endpoint distances, d1 =∥e 1 −p 0∥, d 2 =∥e 2 −p 0∥,(C.11) and their mean dse = 1 2 (d1 +d 2).(C.12)
-
[35]
When evaluated with the paired pericenter/apocenter trajectories, the eccentricity–effect index exhibits a broad, strongly right-skewed distribution
The eccentricity–effect index is then defined as the dimensionless ratio Ei = dee dse ≤2,(C.13) with an additional safeguard: starting points withdse below a small tolerance are treated as immobile and are omitted from this analysis. When evaluated with the paired pericenter/a...
2026
-
[889]
https://doi.org/10.1002/2015JF003751
-
[2013]
Phobos interior from librations determination using Doppler and star tracker measure- ments. Planet. Space Sci. 85, 106–122. https://doi.org/10.1016/j.pss.2013.06.015
2013 doi
Reviewed August 2, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.