Small craters of the lunar maria are observed to be in a state of equilibrium, in which the rate of production of new craters is, on average, equal to the rate of destruction of old craters. Crater counts of multiple lunar terrains over decades consistently show that the equilibrium cumulative size-frequency distribution (SFD) per unit area of small craters of radius >r is proportional r^(-2), and that the total crater density is a few percent of so-called geometric saturation, which is the maximum theoretical packing density of circular features. While it has long been known that the primary crater destruction mechanism for these small craters is steady diffusive degradation, there are few quantitative constraints on the processes that determine the degradation rate of meter to kilometer scale lunar surface features. Here we combine analytical modeling with a Monte Carlo landscape evolution code known as the Cratered Terrain Evolution Model to place constraints on which processes control the observed equilibrium size-frequency distribution for small craters. We find that the impacts by small distal ejecta fragments, distributed in spatially heterogeneous rays, is the largest contributor to the diffusive degradation which controls the equilibrium SFD of small craters. Other degradation or crater removal mechanisms, such cookie cutting, ejecta burial, seismic shaking, and micrometeoroid bombardment, likely contribute very little to the diffusive topographic degradation of the lunar maria at the meter scale and larger.

Impact cratering is likely a primary agent of regolith generation on airless bodies. Regolith production via impact cratering has long been a key topic of study since the Apollo era. The evolution of regolith due to impact cratering, however, is not well understood. A better formulation is needed to help quantify the formation mechanism and timescale of regolith evolution. Here, we propose an analytically derived stochastic model that describes the evolution of regolith generated by small, simple craters. We account for ejecta blanketing as well as regolith infilling of the transient crater cavity. Our results show that the regolith infilling plays a key role in producing regolith. Our model demonstrates that, because of the stochastic nature of impact cratering, the regolith thickness varies laterally, which is consistent with earlier work. We apply this analytical model to the regolith evolution at the Apollo 15 site. The regolith thickness is computed considering the observed crater size-frequency distribution of small, simple lunar craters (< 381 m in radius for ejecta blanketing and < 100 m in radius for the regolith infilling). Allowing for some amount of regolith coming from the outside of the area, our result is consistent with an empirical result from the Apollo 15 seismic experiment. Finally, we find that the timescale of regolith growth is longer than that of crater equilibrium, implying that even if crater equilibrium is observed on a cratered surface, it is likely the regolith thickness is still evolving due to additional impact craters.