“…Some works noted the existence of weak heterogeneity of the lunar regolith [11,[33][34][35]. However, to date, no numerical simulation has been conducted specifically on the weak heterogeneity of the lunar regolith, probably due to three aspects: 1. lack of methods describing heterogeneous media; 2. technical difficulty in constructing random models; 3. excessive computational cost of numerical simulations [16] compared with traditional methods. In this study, we assumed that the weak heterogeneity of the lunar regolith obeys the self-organization model [36], which is basically a stochastic approach for dealing with scattering waves through randomly inhomogeneous structures and is well known in the seismology field; then, we used the finite-difference time-domain method [16] to carry out a series of numerical simulations using various models; next, we analyzed the influence of different parameters on the weak reflection characteristics; finally, we constrained the possible range of each controlling parameter through grid searching techniques by matching the numerical simulation results with the Yutu-2 LPR observations, in terms of waveforms, envelopes, and reflection patterns.…”