2. Shallow Geological Structures of CE-4 Landing Site
Yutu-2 is the first rover to land on the far side of the Moon, and the LPR is capable of transmitting radar electromagnetic rader waves to the lunar surface [
8]. When electromagnetic waves propagate through the ground of the Moon, they are reflected at the interfaces of different materials because different materials have different dielectric constants [
13]. Therefore, the echo signals received by the receiving antenna of Yutu-2 are essentially the reflection signals generated at the interfaces when the electromagnetic waves pass through different media. After the raw data received by the LPR are processed using normalization and DC (direct component) removal and other processes [
28,
29]. By plotting these data on a 2-D plan view (B-scan), the location of the boundary of materials and the difference between the materials on two sides of the boundary can be observed in a more intuitive way, providing information support for the subsequent geological interpretation. The radar frequencies of the two channels are different, and so are the resolutions of them. Channel 1 is a low-frequency radar, and its effective penetration depth is roughly several hundred meters. Channel 2 is a high-frequency radar, and its effective penetration depth is roughly several tens of meters [
11,
13]. Therefore, when interpreting the shallow surface structure of the CE-4 landing zone, it is usually divided into high and low-frequency radar observations for respective interpretation.
2.1. High-Frequency LPR Radar Observation
According to the research of Xiao et al. [
30], the effective detection depth of the high-frequency radar carried by Yutu-2 is about 40 m. The signal beyond 40 m is close to the noise level and, thus, useless [
30]. Therefore, the observation data of the high-frequency radar can be used to discuss the stratification of geological structures from the surface to a 40 m depth. According to researchers, the geological subsurface observed by the radar of 40 m can be roughly divided into three geological units (
Figure 2), as in the following description:
The first layer is the shallow, fine-grained lunar regolith at a depth of 0–12 m below the surface [
8]. According to the radar data after tomographic reconstruction, most of the echo reflectivity in this layer is small, and the distribution of materials is relatively uniform. Therefore, it can be speculated that the material in this layer is mainly fine-grained lunar regolith particles less than 1 cm in diameter [
8]. According to the morphological method of the impact crater, the location of the low calcium pyroxene of the Von Kármán crater is about 8–13 m below the surface [
31], which is consistent with the study of Li et al. [
8]. Meanwhile, Li et al. [
8] believes that the materials in this layer are developed from the ejecta of craters like the Finsen, Von Kármán L, and L’ craters [
8]. Yuan et al. [
32] made a more detailed division in the division in this layer. Yuan et al. [
32] believes there is also a geological boundary 2.5 m below the lunar surface. Their research shows that the material between 0 and 2.5 m is from the von Kármán L and L ’crater ejecta (each crater contributes about 1 m in depth), while the materials between 2.5 m to about 12 m are from contemporary craters such as the Finsen [
32]. Furthermore, it is supported by the results obtained with optical and multispectral imagery [
33].
This layer was further subdivided into multiple geologic structures for a more refined study [
34]. Giannakis et al. [
34] suggests that ejecta from the Von Kármán L and L’ craters cover the entirety of the layer, while ejecta from the Finsen crater are further down. Based on data from Lu et al. [
35] and the chemical composition analysis from Huang et al. [
36], Giannakis et al. [
34] divided the subsurface 0–12 m in the detection area of the Yutu-2 rover into four more detailed sublayers (
Figure 2c). Among them, the material in the first two layers and the latter two layers originated from the ejecta of the Von Kármán L and the L’ craters, respectively.
Figure 2. (a) Radargram obtained within the first two lunar days. (b) Geological stratification diagram of the shallow surface layer in the detection area observed with the high-frequency radar. (c) Geological stratification from 0 to 12 m (HCP and LCP for high–calcium pyroxene and low–calcium pyroxene, respectively).
The second layer is a mixture layer ranging from 12 m to 24 m. The layer is characterized by a number of large rock fragments sandwiched between two layers of fine–grained lunar regolith. The Finsen crater is an Eratosthenian crater of about 3.2 Ga [
38,
39,
40], older than the von Kármán L and L’ impact crates. The Yutu-2 detection route is located within the secondary crater chain ejecta range of the Finsen crater, so the large numbers of rocks in this layer are interpreted as ejecta deposits from the Finsen crater. However, the formation of ejecting sediment beyond the ejecting blanket of the source impact crater is not only through carpeting but is also accompanied by a series of transformation processes such as shearing, mixing, and excavation [
41]. It makes the source of the material layer not necessarily homologous. The tiny–sized rock fragments inside the ejecta may have similar physical properties to the fine-grained lunar regolith. Therefore, it can be inferred that the material composition of the second layer consists of about three parts:
Using the model of the relation between the crater diameter and the distance of the crater to the CE-4 site, we can estimate an ejecta thickness of 7.2 to 15.5 m [
36]. Therefore, the layer can be further divided into two layers: 12–16 m and 16–24 m. The former came from the ejecta, the excavated material of the Finsen crater, and the paleo-regolith formed after the impact, while the latter was formed before the Finsen event.
The third layer is a 24–40 m thickness of accumulated ejecta. This layer is characterized by the alternating layers. These substances were made by ejecta from Alder [
36], Leibnitz [
38,
42] and other older craters [
43] before the Finsen event. Furthermore, it has a layer between 35 m and 40 m that is interpreted as the bedrock of basalt indicated by crater geomorphological evidence [
36,
44]. The signal noise and strong uncertainty are high below 40 m, so it is difficult to determine the accurate depth of Imbrium basalt and the more downward material composition [
8,
37].
In summary, although it can have a lot of nonuniqueness when interpreting the radar data [
45], the echo signals of the Yutu-2 high-frequency radar are relatively clear and clearly stratified, making the stratification of the detected area basically consistent with studies [
8,
30,
34,
36,
37,
46,
47]. However, due to different methods of processing the LPR data and different interpretations of the geological background of the landing zone, various researchers have different interpretations of the evolutionary process and the origin of different layers [
8,
30,
47]. Therefore, a more quantitative and accurate determination of the distribution of the dielectric constant of subsurface materials is urgently needed in the future.
2.2. Low-Frequency LPR Radar Observation
In the interpretation of low-frequency radar observations, different researchers have made various divisions from different aspects [
32,
37,
46,
48]. For the rough division explanation, we will mainly use the division pattern of Feng et al. [
49]. For a more visual and detailed look at the division, it is shown in the division graph (
Figure 3) of Feng et al. [
49]. Zhang et al. [
37] divided the 450 m below the lunar surface into three layers:
Figure 3. (a) Image result of low-frequency data. (b) Structure of the first 300 m of the Moon at the landing site.
-
The first layer (0~110 m) is the sedimentary layer ejected from other craters.
-
The second layer (110~240 m) is the Imbrium basalt layer.
-
The third layer (240~440 m) is ejecta from the Leibniz crater.
The first layer is the most obvious and uniform layer, but the source of the material is the most controversial. The reason for this may be that this layer contains high-frequency radar observations, but, on the one hand, different researchers have different geological backgrounds. On the another hand, interpreted high-frequency LPR observation obtained with different methods shows different clarity, leading to different opinions when inferring the origin of the material from different details. Therefore, although researchers agree with the division of the layer in general (i.e., the layer’s composition is the sedimentary ejecta from nearby craters), there are different opinions on the source of the material inferred from different details (see
Section 2.1). There are researchers who believe that part of the materials originates from the Schrödinger crater [
50]. So, there are different opinions on the interpretation of this layer in low-frequency radar data since the explanation of the formation of this layer by researchers is very similar to the explanation of the high-frequency LPR observation results.
The second layer is a low-reflectivity geological unit located from 110 m to 240 m. The reflectivity of this layer is significantly lower than the reflectivity of the first and third layers. The material in this layer is interpreted as the basalt from Imbrium. From the perspective of the volcanic event, according to the studies of Ivanov et al. [
51] and Huang et al. [
36], the Von Kármán crater was once or repeatedly covered by lava erupted from the Imbrium area, so it is reasonable to infer that the material in this layer is basalt. From the perspective of time, the Leibniz crater occurred about 3.3 Ga ago, and the basalt burst was about 3.1 Ga [
35], so it is reasonable that the molten basalt burst covered the ejecta of the Leibniz crater. In terms of thickness, according to the morphology of the crater, the coverage depth of the Imbrium-aged basalt in the CE-4 landing area is about 120 m, which is consistent with the thickness detected by the low-frequency LPR [
37,
52]. Zhang et al. [
46] divided this layer into more layers based on the data of the low-frequency LPR radar. They believe that the material from 110 to 130 m and 158 to 163 m below the Moon’s surface developed in the interval between the impact events of the 50~110 m, 130~158 m and 263~280 m layers, and the formation process is similar to the development process of the lunar regolith detected by the high-frequency LPR.
The third layer is a high-reflectivity geological unit located about 240 m to 440 m below the lunar surface. Through stratigraphic studies, the material in this layer may come from the ejecta of the Leibniz crater (∼245 km diameter) formed before about 3.3 Ga [
35,
53,
54]. Meanwhile, the results of numerical simulations infer that the thickness of the ejecta from the Leibniz crater is roughly 200 m thickness at the CE-4 landing site [
37]. It turns out to be a good agreement with the thickness of this layer detected by the low-frequency LPR radar. However, Zhang et al. [
46] considered that the material at 163~263 m below the lunar surface is Imbrium-aged basalt, which is an intermediate layer between lava flows. Meanwhile, there is a 17 m thick geological layer around 263~280 m below the lunar surface. It is assumed that this layer is ejecta material from the Alder crater and that the material in the 280~360 m below this is not from the Leibniz crater but is Imbrium-aged basalt. The possible reason for this variability in interpretation lies in the differences in the interpretation of the geological genesis by different researchers. Moreover, the constraints on the physical properties of the different layers below the lunar surface need to be further investigated.
2.3. Effectiveness of Low-Frequency Radar Data
Although high-frequency and low-frequency LPR data have helped us understand the Moon. However, there is a debate about the effectiveness of low-frequency LPR radar data. Pettinelli et al. [
55] considered that most of the low-frequency LPR radar data within 3000 ns are almost all invalid signals caused by instruments or other factors except those within 500–1800 ns and 2000–3000 ns. In their paper, Pettinelli et al. [
55] made average time amplitude(ATA) diagrams and spectral diagrams of the detection data obtained within 10 m of the low-frequency LPR carried by CE-3 and CE-4 and verified each other to point out problems. In terms of spectrograms, Pettinelli et al. [
55] found that the frequency of the signal within the time window of CE-3 and CE-4 within 4500 ns fluctuated from 0 to 100 MHz. However, almost all are 0–30 MHz low-frequency letter-based, with too much deviation from the designed frequency of the low-frequency radar (60 MHz). Moreover, the CE-3 signals within 3000 ns are composed of almost constant low-frequency signals, which means that the signal is noise. In terms of the ATA diagram, Pettinelli et al. [
55] highlighted two issues.
The first issue is that although the two different hemispheres of CE-3 and CE-4 landed on the Moon, the ATA diagram within 2000 ns is almost identical, with almost consistent peaks and troughs. Additionally, the data outside 2000 ns are very similar. The physical significance is that in the cross-section diagram made of two signals, both have almost equally smooth horizontal lines at similar depth positions at 2000 ns. The only difference between the two is that the cross-section of the CE-4 signal shows a little non-smooth area at the top surface.
The second issue is that, from the ATA diagram, the two signals have almost the same decay time. Both signals almost all reached the noise level at 3000 ns. In short, the ATA diagrams of the two sets of signals are very similar. In addition, Pettinelli et al. [
55] considered through the ATA diagram that the signal near 3700 ns was the only signal above the noise level, but it was also composed of signals at too-low frequencies.
In summary, Pettinelli et al. [
55] suggested that the signal in the first 2000 ns of the CE4 radar cross section is mainly composed of the time-invariant coherent noise of too-low frequencies. This noise appears as a horizontal band sequence in the cross-section. The signals after 3000 ns are either noise-level or invalid signals composed of too low-frequencies [
55].
Zhang et al. [
56] responded to Pettinelli et al. [
55]’s query. In response, Zhang et al. [
56] acknowledged that the ATA diagram of the signals of CE-3 and CE-4 were similar in the low-frequency signal region due to electromagnetic interference and noise, with the same overall decay time. However, they argue that differences in the ATA curves after 1100 ns are readily discernible. Meanwhile, they emphasize that the system disturbance and noise can be greatly reduced by subtracting the average trace from each curve because they are reproducible and recognizable. The signal after 1100 ns in the ATA plot that they made by subtracting the average trace (
Figure 4) also has a significant difference.
Figure 4. ATA diagram of the low-frequency LPR signal. (a) Comparison of the ATA plots between the low-frequency LPR data. The ATA curves were computed from the raw low-frequency LPR data acquired by CE-3 along its ∼114 m long track (blue line) and by CE-4 (red line). (b) The ATA plots after subtracting the background.
Furthermore, Zhang et al. [
56] responded to the issue raised by Pettinelli et al. [
55] about the same cross-section lines appearing at the same time and location in the cross-section plots made from the two signals. Zhang et al. [
56] agreed that this was an issue caused by the background of the system, but that this problem could be solved by choosing the right band-pass filtering. Zhang et al. [
56] made plots of the signal processed with either 10–80 MHz or 30–80 MHz band-pass filtering and concluded that the narrow and strong reflections within 1.4–1.8 µs and 2.0–2.2 µs have slopes between traces 1600 and 2600, indicative of subsurface boundaries rather than the invariant system noise. On top of that, there are more time periods where the signals can prove that it is a reflection rather than an over-interpretation.
Overall, Zhang et al. [
56] concluded that strong disturbances and noise can be largely reduced by subtracting the data averaged over several lunar days from each trace. The background can be significantly reduced using band-pass filtering, additionally enhancing the echoes from the lunar subsurface stratigraphy, especially for deep reflections after 1100 ns. Zhang et al. [
56] were able to respond well to some of the issues raised by Pettinelli et al. [
55]. However, other researchers have also given some warning of the pitfalls when interpreting radar data [
57]. In terms of importance, the effectiveness of the data from the low-frequency LPR radars is questionable, due to the fact that the noise problem is caused by the instrument itself. This would lead to difficult corroboration between the conclusions of structural layering within 110 m obtained by the high-frequency and low-frequency soundings respectively. On the other hand, the issues solved in the process of validating the radar data will also enable the academic community to understand the concerns of the radar.