---
title: Stochastic Model of Regolith Gardening
url: https://www.emergentmind.com/topics/stochastic-model-of-regolith-gardening
type: topic
---

# Stochastic Model of Regolith Gardening

Searching arXiv for relevant papers on stochastic models of regolith gardening and adjacent observational/constitutive constraints.
A stochastic model of regolith gardening treats regolith evolution on airless bodies as the cumulative outcome of discrete, spatially random impact events that excavate, bury, and overlap one another through time. In the formulation developed for small, simple lunar craters by Hirabayashi et al., the central quantity is not a single bulk regolith thickness but the probability distribution of thickness at a point, expressed as the fraction of area whose regolith thickness exceeds a specified depth \(h\) [1802.03664]. Within this framework, “gardening” is modeled primarily as regolith production and spatial accumulation by impact cratering, with ejecta blanketing outside crater rims and breccia-lens infilling inside crater cavities as the two explicit regolith-generating processes [1802.03664]. The resulting picture is intrinsically heterogeneous: regolith thickness varies laterally because impacts occur randomly in space, and continued cratering can thicken regolith even after crater equilibrium is reached [1802.03664].

## 1. Conceptual definition and scope

In the stochastic formulation of regolith gardening, repeated impacts into an initially intact, flat surface generate a thickness field that is probabilistic rather than spatially uniform. The modeled processes are excavation and deposition of ejecta blankets, regolith infilling of the transient crater cavity represented as a breccia lens beneath the final crater floor, overlap and superposition of multiple impacts, lateral variability in regolith thickness, and temporal evolution of regolith thickness as crater number increases through time [1802.03664].

This scope is deliberately narrower than many broader uses of the term “gardening.” The model does not explicitly include topographic diffusion of crater forms, regolith mixing or turnover rates in the sense of vertical overturn statistics, particle size evolution, compaction, melting, shock-damage halos outside the adopted crater-regolith units, or spatially heterogeneous distal rays [1802.03664]. A common misconception is that a regolith-gardening model necessarily constitutes a complete maturity model. In this case it does not: the formulation is directed at impact-generated regolith production and accumulation by small, simple craters rather than grain comminution, agglutinate production, or optical maturity [1802.03664].

The stochastic aspect is central because deterministic average-thickness models cannot represent the fact that thickness at a given point depends on the random history of crater overlaps of different sizes and positions [1802.03664]. This is closely related, but not identical, to later continuum treatments in which discrete impact statistics are coarse-grained into effective vertical transport. A 2026 lunar model recasts gardening as a competition between impact-driven advection and stochastic diffusion in depth, arguing that pure diffusion cannot reproduce sharp gradients and displaced centroids in maturity and radionuclide profiles [2604.09524]. This suggests that “stochastic regolith gardening” has at least two complementary meanings in current literature: lateral stochastic coverage by crater-generated regolith units [1802.03664], and stochastic-impact-informed continuum transport in depth [2604.09524].

## 2. Crater geometry, assumptions, and single-event regolith production

The analytical model of Hirabayashi et al. is restricted to populations of small, simple craters [1802.03664]. Each crater contributes two idealized regolith units: an ejecta blanket outside the final rim and a breccia lens inside the crater cavity. For the Apollo 15 application, craters with radius \(r \le 381\ \mathrm{m}\) contribute via ejecta blanketing, whereas craters with radius \(r \le 100\ \mathrm{m}\) contribute via breccia-lens infilling [1802.03664]. The asymmetry reflects the distinction between local cavity infilling and ejecta contributions that can come from larger or more distant craters.

The model assumes that impacts occur randomly on an initially intact, flat surface; newly formed craters are circular; each crater produces regolith only in the ejecta blanket and breccia lens; no distinction is made between regolith particle sizes in ejecta and breccia lens; and no melting, compaction, topographic evolution, explicit mixing law, thermal fatigue, or volcanism is included in the analytical formulation [1802.03664]. Crater production follows a cumulative size-frequency distribution of power-law form, crater overlap is treated statistically rather than geometrically, the minimum crater radius is taken effectively to zero for analytical simplicity, and maximum contributing crater radii are imposed separately for infilling and ejecta blanketing [1802.03664]. The model also does not assume that crater equilibrium halts regolith growth [1802.03664].

Outside the crater, ejecta thickness is prescribed by a power law,
\[
h_{\mathrm{out}} \sim \sigma r_c \left(\frac{r}{r_c}\right)^{-\kappa}, \qquad r>r_c,
\]
where \(r_c\) is the final crater radius, \(\sigma\) is the ejecta thickness ratio at the rim, and \(\kappa\) is the ejecta thickness power-law slope [1802.03664]. For most of the paper, the adopted values are \(\sigma=0.014\) and \(\kappa=3\), following Sharpton (2014) [1802.03664].

Inside the crater, breccia-lens thickness is approximated by
\[
h_{\mathrm{in}}=
\begin{cases}
\delta r_c\left(1-\dfrac{r^2}{r_c^2}\right), & r\le r_c,\\[6pt]
0, & r>r_c,
\end{cases}
\]
where \(\delta\) is the ratio of maximum breccia-lens thickness at crater center to crater radius [1802.03664]. For lunar fresh small simple craters, the inferred range is \(\delta \approx 0.30\text{--}0.34\), and for the Apollo 15 application the adopted value is \(\delta=0.32\pm0.02\) [1802.03664]. The physical significance of this geometry is that the breccia lens is much thicker than the ejecta blanket near crater center [1802.03664].

## 3. Probabilistic formulation of regolith thickness

The stochastic model asks what fraction of area has regolith thickness greater than a threshold \(h\) [1802.03664]. For craters of a single size \(i\), the area fraction at depth \(h\) affected by one crater is
\[
P_i(h,r_{c,i})=\frac{S_i(h,r_{c,i})}{A},
\]
where \(S_i\) is the area at depth \(h\) affected by a crater of radius \(r_{c,i}\), and \(A\) is the total target area [1802.03664]. If \(n_i\) craters of that size are randomly emplaced, the empty fraction approaches the Poisson limit
\[
\left(1-\frac{S_i}{A}\right)^{n_i}\simeq \exp\left(-\frac{S_i n_i}{A}\right),
\]
and extension to all crater sizes yields
\[
P(h)=1-\exp\left(-\frac{1}{A}\sum_i S_i n_i\right),
\]
where \(P(h)\) is the fraction of surface area where regolith thickness exceeds \(h\) [1802.03664].

This quantity has several equivalent interpretations. It is an areal thickness exceedance function, the survival function of regolith thickness at a randomly chosen point, and a probabilistic description that naturally incorporates crater overlap [1802.03664]. Because the formulation returns a full thickness distribution rather than only a mean, it predicts that some locations are buried beneath thick breccia lenses and multiple ejecta blankets while others are affected only by small thin deposits [1802.03664].

Crater production is represented through a cumulative size-frequency distribution
\[
C=A\,\xi\,X\,r_c^{-\eta},
\]
where \(\xi\) is the normalization coefficient, \(X\) is a dimensionless normalized crater-production number proportional to exposure time for constant impact flux, and \(\eta\) is the cumulative power-law slope [1802.03664]. In the Apollo 15 application, the inferred parameters are \(\eta=3.2\), \(\xi=2.21\ \mathrm{m}^{1.2}\), and \(X=1\) for the present produced crater CSFD [1802.03664].

For a given threshold \(h\), the area affected by a single crater is computed separately for ejecta and breccia-lens infill. For ejecta,
\[
S_{\mathrm{out}}=\pi r_c^2\left(\frac{\sigma r_c}{h}\right)^{2/\kappa}-\pi r_c^2,
\]
and for the breccia lens,
\[
S_{\mathrm{in}}=\pi r_c^2\left(1-\frac{h}{\delta r_c}\right),
\]
for \(h\le \delta r_c\) [1802.03664]. These relations are then integrated over crater sizes to obtain the continuous stochastic coverage law [1802.03664].

Because the model yields a thickness distribution over area, the expected regolith thickness for a randomly selected point is defined by
\[
\bar{h}=-\int \frac{dP}{dh}\,h\,dh.
\]
This expected value is the quantity compared with empirical constraints such as Apollo seismic regolith thickness [1802.03664].

## 4. Breccia lenses as the dominant regolith source

The central conclusion of the analytical model is that breccia-lens infilling dominates impact-generated regolith production [1802.03664]. In the adopted geometry, the breccia lens is a paraboloid whose maximum thickness at crater center is \(\delta r_c\), tapering quadratically to zero at the rim [1802.03664]. The parameter \(\delta\) is inferred from crater geometry using a transient crater depth/diameter ratio of \(1/3\), a final crater diameter \(=1.25\) times transient diameter, and observed final crater depth-to-diameter ratios of \(\sim 0.13\) to \(0.15\), leading to \(\delta \approx 0.34\) for final \(d/D=0.13\) and \(\delta \approx 0.30\) for final \(d/D=0.15\) [1802.03664].

At Apollo 15, the quantitative contrast between ejecta blanketing and breccia-lens production is strong. Using local crater counts, the ejecta-only expected thickness is \(0.41\ \mathrm{m}\), whereas the ejecta-plus-breccia-lens expected thickness is \(3.85\pm0.25\ \mathrm{m}\) locally [1802.03664]. Sensitivity tests reinforce this conclusion: varying the ejecta thickness parameter \(\sigma\) has little effect, whereas sensitivity to \(\delta\) is strong [1802.03664]. The reported \(\delta\)-dependence is \(\bar{h}\approx 1.2\ \mathrm{m}\) for \(\delta=0.1\), \(\bar{h}\approx 2.4\ \mathrm{m}\) for \(\delta=0.2\), and \(\bar{h}\approx 4.1\ \mathrm{m}\) for \(\delta=0.34\) [1802.03664].

Volume arguments point in the same direction. For the nominal Apollo 15 parameters, total regolith volume generated by one crater is partitioned approximately as 15% ejecta and 85% breccia lens [1802.03664]. Even under an approximate volume-conserving adjustment of the rim and ejecta assumptions, the partition becomes 42% ejecta and 58% breccia lens, while predicted total regolith thickness remains nearly unchanged: \(3.91\ \mathrm{m}\) for the volume-conserving case versus \(3.85\ \mathrm{m}\) for the original case [1802.03664].

A frequent misconception in earlier crater-gardening discussions is that ejecta blanketing alone explains observed lunar regolith thicknesses. The Hirabayashi model explicitly argues the opposite: ejecta blanketing alone is far too weak, and breccia-lens formation controls total regolith production [1802.03664]. This conclusion is specific to the adopted crater class and regolith units, but within that domain it is the decisive result.

## 5. Apollo 15 calibration, regolith growth, and lunar thermophysical context

The Apollo 15 application computes regolith thickness using the observed crater size-frequency distribution of small, simple lunar craters, with \(r\le381\ \mathrm{m}\) for ejecta blanketing and \(r\le100\ \mathrm{m}\) for regolith infilling [1802.03664]. Allowing for some amount of regolith coming from outside the area, the modeled result is consistent with the empirical result from the Apollo 15 seismic experiment [1802.03664]. The paper also concludes that the timescale of regolith growth is longer than that of crater equilibrium, implying that crater equilibrium on a surface does not imply a fixed regolith thickness [1802.03664].

This distinction matters for interpretation of lunar surface states. Crater equilibrium is a statement about crater populations, whereas the stochastic regolith-thickness distribution continues to evolve through additional impacts [1802.03664]. A plausible implication is that surfaces classified as saturated in crater counts can still be evolving in subsurface physical structure.

Independent lunar thermophysical observations provide a complementary near-surface constraint. Diviner mapping shows that the upper regolith fines layer is globally remarkably uniform, with globally averaged thermal inertia at \(273\ \mathrm{K}\) of \(55\pm2\ \mathrm{J\,m^{-2}\,K^{-1}\,s^{-1/2}}\), global mean \(H\approx6.8\ \mathrm{cm}\), and global-scale variations over \(10^3\ \mathrm{km}\) scales of \(<10\%\) [1711.00977]. Hayne et al. interpret this uniformity as evidence that the upper \(\sim 10\ \mathrm{cm}\) has been rapidly homogenized by impact gardening and lateral mixing on timescales \(<1\ \mathrm{Gyr}\) [1711.00977]. The diurnally active near-surface layer is \(\sim 4\text{--}7\ \mathrm{cm}\), and density increases from \(1100\ \mathrm{kg\,m^{-3}}\) at the surface to \(1800\ \mathrm{kg\,m^{-3}}\) at about \(1\ \mathrm{m}\) depth [1711.00977].

The stochastic thickness model and the Diviner thermophysical results refer to different aspects of regolith structure. The former addresses cumulative production and lateral variability at meter scales and greater [1802.03664], whereas the latter constrains the thermally active fines layer in the upper few centimeters [1711.00977]. Taken together, they indicate that laterally heterogeneous regolith production at depth is compatible with large-scale homogenization of the very shallow fines layer.

## 6. Extensions, analogues, and limitations

A later lunar development generalizes stochastic gardening into an effective one-dimensional transport model in which impact statistics produce a net gardening depth \(D_z(t)\), an effective advection velocity
\[
v_z(t)=\frac{1}{b-2}\frac{D_z(t)}{t},
\]
and a diffusion coefficient
\[
\kappa(t)=\left(\frac{\Gamma}{H}\right)^2\frac{D_z(t)^2}{t},
\]
which are then inserted into an advection-diffusion equation for tracer concentration [2604.09524]. In that formulation, gardening is not pure diffusion but the competition between downward burial and stochastic local mixing [2604.09524]. The model is validated against Apollo maturity profiles over \(1.4\times10^7\) to \(4.5\times10^8\) years and applied to \(^{60}\)Fe and \(^{244}\)Pu depth profiles [2604.09524]. Relative to the lateral overlap model of Hirabayashi et al., this is a vertical coarse-grained transport closure rather than a pointwise areal-thickness distribution [1802.03664; 2604.09524].

Analogous gardening concepts also appear outside the lunar case. On Vesta, photometric mapping with Dawn reveals a phase-curve slope parameter \(\nu\) that is interpreted as a proxy for regolith physical roughness [1702.00207]. High \(\nu\) occurs in ejecta of young craters such as Antonia and Cornelia, while low \(\nu\) correlates strongly with steep slopes and likely mass wasting [1702.00207]. Older crater ejecta are not photometrically extreme, leading to the inference that rough ejecta are smoothed over several tens of Myr, most likely by impact gardening [1702.00207]. This suggests that on Vesta gardening is used less as a thickness-production model than as a roughness-relaxation process.

For S-complex asteroids, a two-process weathering-plus-gardening model represents gardening as the conversion of weathered surface back into fresh surface on a timescale \(\tau_g\), with best-fit values \(\tau_w=960\pm160\ \mathrm{My}\) and \(\tau_g=2000\pm290\ \mathrm{My}\) [1004.3823]. The corresponding mean-field equation,
\[
\frac{dU}{dt}=\frac{1}{\tau_g}-\left(\frac{1}{\tau_g}+\frac{1}{\tau_w}\right)U,
\]
has a natural stochastic reading as a two-state patch process, although that interpretation is explicitly identified in the source synthesis as an inference rather than a statement of the paper itself [1004.3823]. The same work notes a discrepancy between the fitted gardening time and an impact-derived estimate of \(\tau_g\sim270\ \mathrm{My}\), proposing a “honeycomb” or armoring mechanism in which small impactors are inefficient resurfacers [1004.3823]. This controversy illustrates that effective gardening rates may depend strongly on unresolved assumptions about crater formation, ejecta efficiency, and surface structure.

A further limitation arises from constitutive state. The static stratification model of Schräpler et al. does not describe impact flux or overturn rates, but it provides a pressure-controlled compaction law,
\[
\Phi(\Sigma)=\Phi_2-\frac{\Phi_2-\Phi_1}{\exp\!\left(\frac{\log \Sigma-\log p_m}{\Delta}\right)+1},
\]
with \(p_m\propto r_p^{-4/3}\) for solid grains [1505.02923]. In a gardening context, this functions as a closure for post-disturbance porosity structure rather than a gardening law itself [1505.02923]. A plausible implication is that stochastic overturn and burial models require a separate constitutive description of how redeposited material compacts under self-weight.

Across these formulations, the phrase “stochastic model of regolith gardening” therefore denotes a family of related but nonidentical approaches. In the narrow sense established by Hirabayashi et al., it is an analytical description of the probability distribution of impact-generated regolith thickness produced by random crater emplacement, with breccia-lens infilling as the dominant source term [1802.03664]. In broader current usage, it can also refer to effective stochastic transport in depth [2604.09524], roughness evolution under impact reworking [1702.00207], or mean-field resurfacing of surface-state fractions [1004.3823]. The common element is that discrete impacts are treated as random events whose ensemble behavior governs the structure, thickness, or physical state of regolith on airless bodies.

Source: https://www.emergentmind.com/topics/stochastic-model-of-regolith-gardening