1. Introduction
Our current knowledge of the lunar crust suggests that lateral variations in composition are crucial for identifying the structure and distribution of materials on the Moon [
1,
2]. Based on the composition of surface materials, the lunar crust can be categorized into distinct geological regions [
3]. However, these materials also possess unique physical properties that have a notable impact on the spatial differences in near-surface composition. One such property is the dielectric constant, which measures how easily a material can be polarized in response to an electric field. The dielectric constant is a complex parameter: (a) the real part quantifies the electric polarizability of a dielectric medium, whereas (b) the imaginary part quantifies a dielectric material’s inherent dissipation of electromagnetic (EM) energy. A high dielectric constant indicates that the material can store a significant amount of electrical energy, while a low dielectric constant implies that the material is less able to store electrical energy.
Studying the dielectric properties of the lunar regolith is important for several reasons: (a) Material distribution: The dielectric constant reveals important clues about the distribution of materials within the regolith. Different materials have distinct dielectric properties, allowing us to identify and map their presence across the lunar surface. This knowledge is useful for gaining a better understanding of the Moon’s composition and its geological history. (b) Unraveling surface features: Analyzing the dielectric constant helps us understand the origin and evolution of lunar surface features. Variations in dielectric properties indicate changes in material composition, porosity, or moisture content, providing insights into the processes that have shaped the Moon’s surface over time. (c) Volatile detection: As the dielectric constant, to a first order, is used as a density proxy, it can also provide important information about the volatiles trapped in the lunar regolith. Volatiles are essential resources for future lunar missions, offering potential sources of water for life support and fuel production. By quantifying dielectric properties, we can improve our ability to locate and precisely characterize these valuable resources.
In the 1970s, the Apollo and Luna missions collectively retrieved 387.1 kg and 0.3211 kg of lunar samples, respectively. Subsequent laboratory analyses, conducted in studies such as [
4,
5,
6,
7], encompassed measurements of density, complex permittivity, and chemical composition. The outcomes of these investigations revealed significant correlations. Specifically, the real component of the dielectric constant in lunar samples demonstrated dependence on bulk density, while the imaginary component was found to be associated with both bulk density and the abundance of ilmenite (FeO + TiO
2). Moreover, China’s Chang’E-5 mission achieved a successful retrieval of 1.731 kg of lunar samples from the Oceanus Procellarum, as documented in studies by Qian et al. [
8] and Li et al. [
9]. While in situ sampling from the lunar surface offers the most direct insights into the regolith, it is important to note that these samples are presently gathered exclusively at 10 landing sites located in the central part of the Moon’s nearside. Consequently, the information derived from these samples only provides a comprehensive understanding of the lunar regolith within the specific regions where sampling has occurred, as discussed in [
10].
Remotely sensed radar data provides the potential to characterize the dielectric properties of the near surface and subsurface of the Moon’s regolith on a global scale. This information is not only valuable for scientific purposes but also essential for supporting both human and robotic exploration endeavors. Earth-based radar studies have been used to investigate the physical nature of the lunar near surface for resource characterization, landing site hazards, lunar geologic history, and evaluation of physical properties (e.g., [
11,
12]). However, these data have limitations in accessing the lunar farside and in the geometries available for observing the lunar nearside. In this context, the Miniature Radio Frequency (Mini-RF) instrument on the NASA Lunar Reconnaissance Orbiter (LRO) offers global coverage of the lunar surface [
13,
14].
Measuring the dielectric constant of the lunar surface from remotely sensed radar data is not a straightforward process. Radar backscatter is sensitive to physical properties such as roughness, dielectric constant, regolith thickness, FeO + TiO
2 content, and buried rock abundance [
15]. Moreover, radar frequency and incidence angle additionally impact the backscatter. Radiative transfer (RT) equations can be used to express the magnitude of radar backscatter as a function of physical properties and radar parameters, which can form the basis of a forward model [
15]. While a forward model can establish a highly complex and nonlinear relationship between radar signal and surface properties, we additionally require an inversion approach along with forward modeling to quantify the dielectric constant from radar observations.
Currently, only a few studies have attempted to estimate the dielectric constant of lunar regolith through the inversion of radar data. Campbell et al. [
16,
17] conducted Earth-based radar observations to investigate the scattering properties of planetary surfaces. They developed an inversion model for the dielectric constant based on the co-polarization ratio of radar echoes. Their model assumes negligible cross-polarization components and that the radar backscatter from the subsurface is significantly stronger than that from the surface. As a result, this model is only applicable to flat lunar surface regions and often overestimates the dielectric constant [
18]. On the other hand, Kumar et al. [
19] developed a dielectric constant inversion technique based on measured LH and LV returns of Mini-RF. Here, a ratio of Fresnel reflection coefficient is derived using a three-component decomposition method, which is then numerically solved on a pixel-per-pixel basis to estimate the real part of the dielectric constant. A limitation of this approach is the numerical solver itself, which does not always converge to a numeric value. Specifically, the decomposition method used to separate different scattering processes from the radar data can be problematic. The separation is not always reliable, as it depends heavily on the accuracy of the model assumptions and the numerical solver. This can lead to results that are not directly correlated with the physical properties of interest, such as roughness or the dielectric constant, due to the inherent difficulties in quantifying these processes accurately [
20].
Given the complex forward relationship of radar signals with physical properties, machine learning, especially artificial neural network (ANN), has demonstrated its utility as an emerging inversion approach in Earth-based retrieval studies [
21]. Shukla [
22] and Shukla et al. [
23] employed the Integral Equation Model (IEM), which only includes surface scattering, as a forward model and multi-layer perceptron (MLP) ANN for inversion and applied it to lunar regolith. Another study utilized the IEM modeling technique in two-dimensional and three-dimensional simulations using the Chandrayaan-2 L-band DFSAR data for permanently shadowed regions (PSRs) with ANN as the inversion [
24]. The penetration depth of L- and S-band suggests that the scattering contribution may not be likely from the top surface alone, thereby making the IEM less usable for the inversion approach. Furthermore, Gao et al. [
25] applied the extended Bragg (X-Bragg) model to Mini-RF data and inverted the dielectric constant using a minimization approach. The proposed model, however, is limited to regolith fines in lunar PSRs and mid-latitude regions. None of the studies have fully integrated the Mini-RF data for global and polar perspectives of the dielectric constant.
We present a new inversion model for dielectric constant that uses spatially semi-controlled global and controlled polar Mini-RF S-band mosaics [
2]. Here, we develop a two-layer lunar regolith model with five basic processes through direct scattering from the top surface and bedrock, diffuse scattering from buried inclusions, and scattering from the interaction between bedrock and buried inclusions. We parameterize the IEM for rough surface and derive the scattering contribution from a low dielectric layer of Rayleigh spheres with irregular regolith-bedrock boundaries. As part of the inversion, we develop a novel deep learning-based model to retrieve dielectric constant values. The robustness and reliability of the inversion are evaluated by comparing the retrieved dielectric constant with lunar sample drive core measurements [
10]. Our approach addresses the limitations of previous models, which either involved surface scattering in the forward model or could only be used to invert the local dielectric constant for regolith fines and smooth mid-latitude regions.
2. Data and Methods
Mini-RF is a side-looking synthetic aperture radar (SAR) instrument that transmits circular polarization and receives coherently on orthogonal linear polarizations (H and V) at 12.6 cm (S-band) or 4.2 cm (X/C-band) wavelengths [
26]. The data is acquired in either ‘baseline’ mode with 150 m spatial resolution or in ‘zoom’ mode with 30 m spatial resolution. In our work, we focus on S-band zoom data products, which are sensitive to scatterers in the size range of >~0.1–1.26 m and penetrate to similar vertical depths [
26]. Mini-RF data can be represented by the Stokes parameters
[
27]:
Here,
is the complex voltage in the subscripted polarization. Model inputs for the horizontal (
) and vertical (
) components are derived from Mini-RF Stokes parameter-controlled polar [
28] and semi-controlled global [
2] products. For this work, the global and polar products were down-sampled to spatial resolutions of 64 pixels/degree and 512 pixels/degree respectively. A lower resolution was used for the global product because it was semi-controlled, and a higher resolution for the polar products because they were fully controlled. The polar products included separate east-looking and west-looking data. These data were defined as consisting of individual Mini-RF data strips that had azimuthal look directions between 0
–180
and 180
–360
, respectively. Both datasets were included in the analysis to reduce potential look direction bias in the inverted results for the polar products. The proposed methodology is illustrated in
Figure 1.
2.1. Impact of Physical Properties on Radar Backscatter: A Radiative Transfer Modeling Approach
Radiative transfer modeling is an effective approach to understanding the interaction between electromagnetic waves and the surface in the microwave region, thereby gaining insights into surface characteristics. Backscattering models provide a foundation for evaluating the electromagnetic (EM) response and the factors that influence it. Simulating the EM response of the Moon is a complex task that requires a comprehensive understanding of the EM properties of materials at a desired frequency range. While lunar samples have been extensively tested at frequencies below 1 GHz [
10], challenges arise when dealing with frequencies > 1 GHz. This poses a time-consuming endeavor that demands significant effort. To address this, we attempted to replicate the EM response using a widely-used model commonly employed in soil studies [
21].
For our analysis, we employ a two-layer lunar regolith model that considers the regolith to be a homogeneous fine-grained layer with buried inclusions between the top surface and bottom subsurface of well-defined roughness and dielectric constant, as shown in
Figure 2. We represent five basic processes: (a) rough surface scattering from the top regolith, (b) volume scattering from buried inclusions, (c) subsurface scattering from bedrock, (d) subsurface-volume scattering from the interaction of the radar signal with bedrock first followed by buried inclusions, and (e) volume-subsurface scattering from the radar interaction with buried inclusions first followed by the bedrock.
To model these five scattering processes, we parameterize the Integral Equation Model (IEM) for the rough surface and derive the scattering contribution from a low dielectric layer of Rayleigh spheres with irregular regolith-bedrock boundaries [
29,
30]. The Rayleigh parameters are extracted for estimating the absorption, scattering and transmission losses of the radar signal when interacting with the regolith at different depths. In this regard, the total radar backscatter from our two-layer lunar regolith model can be expressed as
The rough surface scattering component from the top regolith is given by
where
Here,
is a wavenumber and hence, dependent on the incident radar wavelength,
is the vertical component of roughness, represented by root mean squared (RMS) height,
is the incidence angle of the sensor,
is the Kirchhoff field coefficient, and
represents the complementary field coefficient for either HH or VV polarization channels [
30]. The horizontal component of roughness or autocorrelation length is expressed in the form of an exponential surface correlation function, i.e., by
. The exponential function is found to be well-suited for representing the natural surface conditions on both the Earth’s surface and other planetary settings [
15].
The second term is associated with the subsurface bedrock, where the radar signal is scattered from its direct incidence on the bedrock [
15]. The radar signal, in such a scenario, is attenuated twice while passing from the top surface layer to the bottom and vice versa.
Here, the depth of the regolith layer is denoted by
,
,
, and
, where
is the scattering angle,
is the transmission angle, and
is the incidence angle, respectively (as shown in
Figure 2). The extinction coefficient is represented by
and, hence, the optical depth is
. The power Fresnel transmission coefficient is denoted by
, where
and
decide the respective angle of the radar signal as it travels in the regolith medium [
29].
The third term is the volume scattering component from the buried inclusions, which may be either rocks or water ice.
Here, the term is the Rayleigh phase function for horizontal and vertical polarizations, and is the radar albedo, represented by the ratio of volume scattering coefficient and extinction coefficient .
The subsurface-volume and volume-subsurface terms are derived by considering the reciprocity condition, which necessarily makes both terms the same. In this scattering regime, the bedrock exhibits a significant coherent reflectivity, given by
[
29]. Here,
is the radar wavenumber in the layer and
is the RMS height of bedrock.
The Rayleigh phase function,
is equal to
. Since the above term is equal to
, we sum this term two times with the other scattering components in the total backscatter equation. For a detailed derivation of individual scattering contributions, we refer to Fung and Chen [
29]. Note that our model does not consider coherent backscatter opposition effects and multiple scattering events among the rocks.
2.1.1. Parameterization of Radiative Transfer Model
In our simulation, the autocorrelation length is assumed to be the same as the wavelength of incident radar waves (i.e., 12.6 cm) and is used with an exponential surface correlation function [
15]. This assumption is due to the lack of surface roughness information comparable to the Mini-RF wavelength [
15]. The variation in roughness is defined by the range of RMS height between 0 cm (very smooth) and 5 cm (very rough). Based on the Apollo and Luna samples, we have taken the FeO+TiO
2 wt% from 0 to 30 with a step of 0.5 [
15]. Bulk density is assumed to be between 0.75 and 3.32 g/cm
3 [
10]. One important parameter representing the quantitative presence of buried inclusions in the regolith is the volume fraction. We describe the limits of volume fraction between 0 and 0.1 with a step of 0.001. This parameter is then used in Lichtenecker’s mixing rule to compute the effective dielectric constant of regolith with buried inclusions, such as silicate rocks [
15]. We model the buried inclusion as a Rayleigh spherical scatterer with a radius varying from 0.5 cm to 5 cm [
15]. Most of the scatterers are non-spherical. However, when they are randomly oriented and distributed within a layer, they can be approximated as a collection of spherical scatterers [
29]. Regolith thickness varies from 4 m to 15 m in the model. The incidence angle is further ranged between 0
and 80
, thereby capturing all the incidence angles of current state-of-the-art orbital S-band radar systems in lunar orbit [
15]. In
Table 1, the range of all the above parameters is considered either from the Apollo science data or from literature-based surveys [
10,
15,
22]. Using these parameters, we parameterize our radiative transfer (RT) model to simulate the total backscatter for horizontal and vertical polarizations. We then prepare a comprehensive lookup table based on the RT simulations, which is further used for training our novel deep learning-based inversion model.
2.1.2. Sensitivity Analysis of Radar Backscatter
Different lunar materials and surfaces have different dielectric constants, which affect how radar waves interact with them. Moreover, the incidence angle at which radar signals strike the surface affects the scattering properties. For example, understanding how radar backscatter changes with incidence angle can help differentiate between scattering from a rough surface and volume scattering due to buried inclusions in the regolith. To investigate the impact of incidence angle and dielectric constant on radar backscatter, we examine two scenarios with our two-layer model, as shown in
Table 2. The first scenario is for sensitivity to incidence angle, while the second corresponds to the dielectric constant. For both scenarios, we consider the regolith thickness to be 5 m, the RMS height as 1 cm, and the radius of the buried rock inclusions as 1 cm with 25% volume in the regolith layer [
15]. For the first scenario, we use a dielectric constant of
for the regolith layer (corresponding to a bulk density of 1.525 g/cm
3 and 7 wt% ilmenite content) and
for buried rock (corresponding to solid rock with a density of 3.2 g/cm
3 and 5 wt% FeO + TiO
2). The incidence angle ranges between 0
and 80
. For the second scenario, a 49
incidence angle of Mini-RF is used, while the real part of the dielectric constant is varied from 2 to 10 with a fixed imaginary part of 0.003. Since we are using both east- and west-looking polar mosaics, it is necessary to examine how the incidence angle influences the radar backscatter and thus induces a bias in the prediction. Such an analysis could be important to understand the backscattering behavior of the lunar regolith before inverting the dielectric constant from Mini-RF (see
Section 3.1).
2.2. Retrieval of Dielectric Constant from Mini-RF SAR—A Deep Learning-Based Inversion Model
Deep learning algorithms are powerful tools for modeling complex systems that involve multiple variables and exhibit nonlinear behaviors [
31]. When dealing with inversion problems, which possess both challenging characteristics, leveraging deep learning-based models becomes a practical choice to achieve highly accurate representations. These models are inspired by the structure and functioning of the human brain [
31]. They consist of interconnected computational units called neurons. These neurons are organized into layers, with each layer responsible for specific computations. In essence, neural networks are mathematical models that can capture intricate patterns and relationships in data. The fundamental building block of a neural network is the perceptron, a simplified model of a biological neuron [
31]. Perceptrons take input data, apply weights to these inputs, and then sum them up. The result is passed through an activation function to produce an output. This output is used as input for subsequent layers of the network, forming a multi-layer perceptron neural network. The architecture is inherently robust to noise [
31]. The presence of multiple hidden layers in deep neural networks enhances the complex feature learning process of the network even in the presence of noise, reducing the need for rigorous data preprocessing.
Inversion problems often involve deducing the properties of a system from observed data. These problems are known for their complexity and nonlinearity, making them challenging to solve using traditional approaches. Deep learning-based models excel in handling such complexities. They can automatically discover relationships between input and output data, capturing nonlinearities and multivariate interactions that may be challenging to specify using traditional mathematical models. However, there is a caveat with deep learning: it typically requires a substantial amount of training data to generate accurate representations. This is because neural networks need numerous examples to learn and generalize from. We use the lookup table from our RT model simulations to train the deep learning model; therefore, the size of the training data is not a concern in our case. By utilizing pre-existing simulations as training data, we can leverage the advantages of deep learning without the need for massive in situ datasets.
For our setup, we fine-tuned a multi-layer perceptron neural network with five hidden layers containing 128, 64, 32, 8, and 4 nodes, respectively, to achieve a reliable and highly accurate model of inversion. The number of nodes was chosen based on a grid search. The activation function for all the layers is “ReLU”, as it is known for robust performance compared to other choices in Earth-based retrieval studies [
21,
22] and for filtering out irrelevant variations introduced by noise [
21]. As part of the input-output configurations, we have four inputs, namely incident wavelength, total radar backscatter at horizontal and vertical polarization (
,
), and incidence angle (
), and two outputs, namely real part of dielectric constant (
) and imaginary part of dielectric constant (
). Before training the model, the dataset was shuffled and then split into two parts: 20% was kept as validation data, and the remaining 80% was used for training. Moreover, during the training process, K-Fold shuffling with five folds was applied. Cross-validation with five folds aims to estimate how a model will perform and generalize on unseen, noisy data, striking a balance between two aspects, i.e., reducing both variance and bias. If we use too few folds (e.g., 2 or 3), the model’s performance estimates can be highly variable, leading to unreliable assessments of its generalization capability (high variance). On the other hand, if we use too many folds (e.g., 10 or more), each fold’s test set becomes very small, which can introduce bias into the performance estimates because the test data might not be representative of the overall distribution. The inherent noise in lunar radar data is expected to be within the noise-handling capacity of multi-layer perceptrons [
21,
22]. Ultimately, the evaluation of the inversion model is performed by computing the statistical variables of mean squared error (
) and the coefficient of determination (
).
The testing of our model is performed on the Mini-RF global and polar mosaic (east and west) of horizontal and vertical components, with an incident wavelength of 12.6 cm and an incidence angle of 49
. The model output consists of dielectric constant maps, with scales set at 512 per degree for polar regions and 64 per degree for global regions (see
Section 3.3). We use horizontal and vertical components because the transformation of radar polarization information into hybrid-pol Stokes parameters follows the same mathematical principles regardless of the initial scattering matrix representation [
32].
2.3. Validation with Apollo Samples
To validate our inversion model, we conducted tests using 7.5 m
15 m spatial scale Mini-RF data of the Apollo 11, 14, and 16 landing sites. Although the Mini-RF data strips cover the surroundings of Apollo 12, 15, and 17 landing sites, the exact locations of these landing sites are not covered, limiting our capability to perform validation for those sites. Since the Mini-RF resolution is coarser compared to the point-scale sample acquisition from Apollo landing sites, we considered an average of the magnitude of dielectric constant of 65 samples within the Apollo landing sites from Carrier III et al. [
10]. There are a total of 12, 17, and 36 measurements from the Apollo 11, 14 and 16 landing sites, respectively. We then generated dielectric constant images for the respective sites from Mini-RF data using our inversion model. To facilitate comparison with dielectric constant measurements at the sample scale, we employed a radial averaging approach. The circle’s radius started at 7.5 m, corresponding to the minimum Mini-RF spatial scale, and extended from 10 m to 40 m (encompassing an area of 5 km
2) with increments of 2 m. We ultimately visualize the variations in dielectric constant within each site, due to multiple increments, in the form of violin plots (see
Section 3.2). This approach enhances confidence in comparing the inverted dielectric constant obtained from Mini-RF with laboratory-measured dielectric constants.
4. Discussion and Implications
Mini-RF S-band observations have been used to derive the dielectric properties of the Moon at global and polar scales. The predictions from our inversion model are in line with laboratory-measured dielectric constants of Apollo samples, proving its robustness and capability in discerning lunar materials based on dielectric grounds. The use of the east- and west-looking Mini-RF mosaics for polar regions improves dielectric prediction compared to using only one look direction. In this way, we reduce the uncertainty in our maps caused by Mini-RF’s look direction. For the first time, we now understand the spatial variability of dielectric constant across the lunar surface, which could provide new insights into the vertical physical characteristics of the regolith. The observed variability in dielectric constant can mainly be attributed to differences in mineralogical content, bulk density, and the population of excavated regolith materials.
S-band radar exhibits sensitivity to scatterers at centimeter-to-meter scales. Notably distinct from typical composition of lunar regolith, certain near-surface terrains showcase unique scatterer size distributions. Among these terrains are craters and pyroclastic deposits. When focusing on young, fresh craters, radar data of the lunar surface prominently highlight their walls, floors, and ejecta blankets [
39]. This is due to a substantial proportion of scatterers in the centimeter-to-meter range, which are not only deposited on the surface but also mixed within the lunar regolith (a contribution from impact melt is also possible). As a result, the dielectric constant for these craters surpasses that of the adjacent terrain due to increased effective bulk density. Moreover, the notable variations in dielectric constant within the craters could likely be caused by progressive weathering of regolith material, thereby shifting the particle size distribution towards a range that falls beyond the sensitivity of S-band radar. We thus foresee in the future to explore the relationship between the dielectric constant and crater age, similar to Bell et al. [
40], Fassett et al. [
41], and Nypaver et al. [
42].
One potential use of our dielectric maps is in understanding the nature, abundance, and distribution of surficial water ice in the lunar PSRs. According to Heggy et al. [
43], regolith fines on the floors of permanently shadowed craters smaller than 5 km in diameter are optimal targets for the unambiguous detection of water-ice enrichment using S-band radar observations. Moreover, regions with high loss tangents could indicate significant EM wave absorption at S-band frequency by the regolith, which may help determine the proportion of water-ice present in the regolith. By understanding the dielectric properties of these regions from our maps at the polar scale and their correlation with crater characteristics, we can enhance our ability to identify and study water ice in these intriguing lunar environments. The remarkably low dielectric constant and loss tangents of lunar materials, for example, in the crater floor of Rozhdestvenskiy and Cabeus or in the PKT regions, could indicate their high transparency to electromagnetic energy at S-band frequency. This means that radar waves could easily penetrate lunar soils to greater depths.
Our dielectric maps, when merged with scattering properties and other new LRO observations, could prove important in actively connecting our knowledge about the physical properties and composition of lunar surface materials. In addition, it is worth comparing the dielectric constant with the maps of H from Diviner, as both are notionally density proxies. This comparison could be made by functionally relating them and identifying whether mineralogy is the main residual. Such a combined use of measurements could refine previous interpretations of lunar terrane divisions (such as PKT, FHT, and SPA). This could further our knowledge of the evolution of physical properties and their connection with space weathering and small- or large-scale impact modification processes. At the same time, there are several telemetry methods for inverting dielectric constant from SELENE Lunar Radar Sounder (LRS) and microwave brightness temperature that could be used as an additional cross-validation source for our dielectric maps. According to Hongo et al. [
44], LRS provides detailed subsurface information but requires careful interpretation of radar wave interactions. Microwave brightness temperature measurements offer broad spatial coverage but may be influenced by surface temperature and other factors. We thus emphasize the need to integrate multiple radar data sources to achieve a comprehensive understanding of the dielectric properties of the lunar surface.
5. Conclusions
Mini-RF S-band observations provide a unique opportunity to map the physical properties of lunar regolith. In this study, we present semi-controlled global and fully controlled polar S-band maps of lunar dielectric constant at 64 and 512 pixels per degree spatial resolution, respectively. These maps could be used to understand the distribution of materials within different stratigraphic layers at multiple depths. Moreover, the variability in dielectric constant can help evaluate the radar response, thereby inferring the spatial distribution of lunar features. For this, we developed a novel deep learning inversion architecture for dielectric retrievals from Mini-RF S-band data. The model is trained using data simulated from a physics-based radiative transfer model, parameterized with realistic values of physical properties from Apollo science data under Mini-RF radar configuration. We tested the trained model on S-band Mini-RF global and polar mosaics. The predictions from the model are in strong agreement with field measurements at Apollo sites. Mare basalts in the PKT exhibited low value of dielectric constant, mainly due to high ilmenite content and relatively low bulk density. We notice some localized hotspots within PKT, i.e., in the crater ejecta and pyroclastic deposits. On the contrary, there is a clear distinction in dielectric constant between PKT and highland materials of FHT. One of the reasons is the composition by itself: highland materials have low ilmenite content. It is observed that the variation in dielectric constant is influenced by changes in ilmenite content, which could easily mask the effect of density on loss tangent. Moreover, it is difficult to differentiate the SPA and FHT using dielectric contrast. We further notice intriguing and complex dielectric properties in the PSRs of the lunar south and north poles, providing an opportunity to quantify the fractional volume of water ice content in the polar regolith. In this regard, our findings highlight the capability of the inversion model to capture subtle differences in dielectric constant due to surface composition and buried volumes discernible at Mini-RF radar wavelengths. From the sensitivity study, we found that for S-band, radar backscatter decreases with incidence angle but increases with dielectric constant at Mini-RF incidence angles. Our dielectric constant maps can thus be utilized to understand the lunar surface processes at global and polar scales, which could play an important role in ARTEMIS mission planning.
Currently, the inversion model is also capable of producing dielectric maps of polar regions at the highest available spatial resolution of 2048 pixels per degree. The future scope of our study is to calibrate the model for local incidence angles. We expect the model to perform better in capturing the finest details of the regolith, when provided with local incidence angle information. We also aim to use advanced radiative transfer theories for better representation of lunar regolith. For instance, the use of the dense radiative transfer model (DMRT) [
45] could be an optimal option for quantifying radar backscatter as a function of physical properties. Moreover, the advanced integral equation model (AIEM) [
46] is valid over large scales of roughness and performs better in simulating rough surface scattering processes than traditional physics-based models, such as IEM. This could provide an opportunity to include a wide range of realistic roughness values while training the model. We also intend to improve our deep learning model in terms of input-output configuration, i.e., to invert both surface roughness and dielectric constant. Such an output could help us separate the contribution of roughness from volume scattering processes, thereby quantifying the water ice content from SAR observations.
We note that there might be moisture contamination in Apollo samples [
4,
5,
6,
7], which could lead to unreliable dielectric constant measurements. However, in Strangway and Olhoeft [
7], the effect of moisture levels on the dielectric constant is found to be most significant at lower frequencies, i.e., <10 kHz, wherein the loss tangent increases drastically by almost a factor of 20. We thus recommend analyzing more returned lunar samples from future lunar missions with great care, ensuring negligible effect of moisture on the results. This will also greatly help the planetary remote sensing community validate inversion frameworks with great confidence. As an alternative, lunar meteorites and sample-returned rocks could also be used as proxies for lunar regolith (particularly regolith breccias or Apollo/Luna regolith powders). In fact, new techniques could be useful for studying the porosity and mechanical properties of meteorites, which might be used as simulants for future laboratory studies (e.g., [
47,
48]). In addition, using a tomographic technique, we could identify the main minerals, density, and moisture levels in each sample (even on the subsurface) that impact the dielectric constant and thus could be useful in testing our model. At the same time, conducting experiments under varying conditions and with different techniques can present a challenge for direct comparison. Given that we know the relationship between dielectric constant, density, moisture levels, and temperature, there exist normalization methods that could provide a unified measurement value [
4,
5,
6,
7]. Once the data is normalized, statistical methods such as regression analysis or multivariate analysis can be employed to compare the results across different conditions. These techniques are standard in the field and provide a robust framework for analyzing and comparing the dielectric properties of lunar samples. Furthermore, to mitigate this, we foresee the combined use of dielectric constant maps augmented with other LRO, Chandrayaan-2 products, and compositional information to infer the evolution of lunar features due to weathering and impact cratering processes.