The timeline of the early lunar bombardment remains unclear. The bombardment rate as a function of time is commonly modeled by three types of shapes: tail-end, sawtooth, and terminal cataclysm. Differently aged melt records the occurrence time of impact events and thus is crucial for constraining the timeline of the early lunar bombardment. Based on a spatially resolved numerical model, we simulate the evolving distribution of differently aged melt with a long-term impact mixing, where different shapes of impact rate function are considered. We compare the outcome of melt age distribution from different scenarios with the actual data from the lunar meteorites and the returned samples. The results suggest that, if the present data are representative of the melt age distribution on the Moon, the shape of the impact rate function is more likely comparable to the tail-end over the sawtooth and the terminal cataclysm, with the terminal cataclysm being least likely. In addition, using state-of-the-art U–Pb dating techniques, more abundant ancient basin melt is likely to be found in returned samples.
