The impact melt that records the formation time of basins is essential for the understanding of the lunar bombardment history. To better understand melt distribution on the Moon, this study investigates mixing of melt by small impacts using a Monte Carlo numerical model. The obtained mixing behavior is then integrated into a larger scale model developed in previous work. While large impacts produce most of the melt volume in both the regolith and megaregolith, we find that the dominant source of melt near the surface is small impacts. Material in the top meter is affected mainly by impacts that form craters <5 km in diameter. In the uppermost 10 cm, melt with age <0.5 Ga is abundant; while as depth increases older melt is increasingly present. This may indicate that the excess of impact melt <0.5 Ga in lunar samples from the near surface is caused by the cumulative mixing of small impacts. A comparison of the age distribution of melt derived from craters of different sizes with that of impact glass constrains the size of spherule-forming impacts. Our model is consistent with observations if most impact glass spherules from the near surface are produced by <100 m craters and >100 m craters do not contribute abundant spherules. The distribution of the datable melt with depth is also analyzed, which is essential for future sampling missions. Excavated materials of young and large craters (>100 m on highlands; >10 km on maria) appear to be the most fruitful targets.
