Next Article in Journal
Impacts of the Microclimate of a Large Urban Park on Its Surrounding Built Environment in the Summertime
Next Article in Special Issue
Dependence of the Hydration of the Lunar Surface on the Concentrations of TiO2, Plagioclase, and Spinel
Previous Article in Journal
Phase Shift Migration with Modified Coherent Factor Algorithm for MIMO-SAR 3D Imaging in THz Band
Previous Article in Special Issue
High Resolution Apparent Thermal Inertia Mapping on Mars
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Uncertainty Introduced by Darkening Agents in the Lunar Regolith: An Unmixing Perspective

Image Analysis Group, TU Dortmund University, Otto-Hahn-Straße 4, 44227 Dortmund, Germany
*
Author to whom correspondence should be addressed.
Remote Sens. 2021, 13(22), 4702; https://0-doi-org.brum.beds.ac.uk/10.3390/rs13224702
Submission received: 24 September 2021 / Revised: 13 November 2021 / Accepted: 18 November 2021 / Published: 21 November 2021
(This article belongs to the Special Issue Cartography of the Solar System: Remote Sensing beyond Earth)

Abstract

:
On the Moon, in the near infrared wavelength range, spectral diagnostic features such as the 1- μ m and 2- μ m absorption bands can be used to estimate abundances of the constituent minerals. However, there are several factors that can darken the overall spectrum and dampen the absorption bands. Namely, (1) space weathering, (2) grain size, (3) porosity, and (4) mineral darkening agents such as ilmenite have similar effects on the measured spectrum. This makes spectral unmixing on the Moon a particularly challenging task. Here, we try to model the influence of space weathering and mineral darkening agents and infer the uncertainties introduced by these factors using a Markov Chain Monte Carlo method. Laboratory and synthetic mixtures can successfully be characterized by this approach. We find that the abundance of ilmenite, plagioclase, clino-pyroxenes and olivine cannot be inferred accurately without additional knowledge for very mature spectra. The Bayesian approach to spectral unmixing enables us to include prior knowledge in the problem without imposing hard constraints. Other data sources, such as gamma-ray spectroscopy, can contribute valuable information about the elemental abundances. We here find that setting a prior on TiO2 and Al2O3 can mitigate many of the uncertainties, but large uncertainties still remain for dark mature lunar spectra. This illustrates that spectral unmixing on the Moon is an ill posed problem and that probabilistic methods are important tools that provide information about the uncertainties, that, in turn, help to interpret the results and their reliability.

Graphical Abstract

1. Introduction

Reflectance spectroscopy can provide valuable insights into the composition of planetary surfaces. On the lunar surface, the mineral composition is dominated by plagioclase and pyroxenes [1], which show diagnostic absorption bands in the near-infrared wavelength range at 1- μ m and 2- μ m (e.g., [2]). These spectral features enable us to estimate the abundance of certain minerals or elements. Correlations between spectral parameters and spectra with known composition can be used to create large-scale maps of, for example, compounds of chemical elements, such as FeO or TiO2 (e.g., [3,4,5]). To determine the abundances of known materials, spectral unmixing can be performed (e.g., [6,7]). With higher spatial and spectral resolution data available in the last decade, such as data from the Moon Mineralogy Mapper (M3) instrument on-board the Chandrayaan-1 spacecraft [8], it is possible to investigate the composition of the Moon in more detail and to gain insight into the evolution and geology of the Moon.
Unmixing describes the process of determining the abundances of a set of constituent materials (endmembers) for a given spectrum employing radiative transfer modeling to describe the interaction of light with the different endmembers. Here, it must be differentiated between a surface where the endmembers build a macroscopic mixture, such that they are spatially separated as in a checkerboard pattern. This type of mixture can be described by a linear superposition of endmember reflectance spectra. The Moon, however, is covered by a layer of grains from crushed rocks, the so-called regolith. This constitutes a typical intimate mixture, where light interacts with the different particles by multiple scattering inside the medium. Intimate mixtures cannot be described by a linear mixture of reflectance spectra in the near infrared. It has, however, been shown (e.g., [6]) that the single scattering albedo (SSA) spectra can be unmixed linearly. Models such as the Hapke model [9] can be used to calculate the SSA spectra, thus transforming a non-linear problem into a linear unmixing problem.
The vast majority of unmixing approaches rely on classical optimization techniques (e.g., gradient descent or Levenberg Marquardt algorithm) to fit the parameters of the model, such as abundances or grain size, to match the measured spectrum. This has been done for several planetary bodies, such as on Mars (e.g., [10,11]) or on the Moon (e.g., [2,12,13]). The endmembers can for example be extracted from the image itself such that they span a simplex. Each spectrum in the image can then be calculated as a linear superposition of the edges of the simplex (endmembers) [6]. While this has the advantage of being applied easily, it has the disadvantage that the composition of these so called image endmembers might not be known. A different approach is to select endmember spectra from a library of laboratory samples of known composition [6]. This way the composition can be directly determined, but the endmembers must be selected carefully and some spectra might still lie outside of the resulting simplex.
Spectra obtained under lunar conditions, however, are systematically different compared to laboratory spectra of earth-analogs measured in the laboratory. Due to the lack of an atmosphere, the lunar regolith is subject to the continuous influence of the space environment. Micrometeorite bombardment and solar wind change the physical and chemical characteristics of the upper layers of the regolith [14,15]. These changes include (1) creation of submicroscopic iron particles (smFe0) [14,15], (2) selective comminution of larger grains into smaller grains, (3) melting of grains into agglutinates [15] and (4) an increase in porosity creating a ’fairy-castle’ structure [14]. The smFe0 particles can be grouped into two categories. Firstly, the smaller nanophase iron (npFe0) particles, which mostly redden the overall spectrum and dampen the absorption bands, but also darken. These particles are about 2–30 nm in diameter and accumulate in the rims of the mineral grains [15]. Secondly, the microphase iron (mpFe0) particles which darken, but do not lead to an increase in spectral slope [16]. All these changes significantly influence the optical properties of the lunar surface and make it substantially more difficult to distinguish between minerals. For example, Sunshine and Pieters [17] and Mustard and Pieters [2] found that, due to low spectral resolution and signal-to-noise ratio multiple solutions might result in similar least squares fits, but are physically implausible.
Therefore, unmixing lunar spectra either rely on image endmembers or mature samples returned by the Apollo missions. A catalog that covers a variety of compositions of the maria is the Lunar Soil Characterization Consortium (LSCC) catalog [18,19]. However, because the simplex is limited to only the main mare compositions and a few highland samples, this catalog cannot cover the entire mineralogy of the lunar surface. For example, the highlands are dominated by plagioclase (90%) [1] and the LSCC catalog cannot accurately describe the typical highland spectra [20]. Further, uncommon compositions such as the spinel-rich regions detected by Pieters et al. [21] cannot be explained by unmixing purely based on the LSCC catalog.
Usually, when unmixing a spectrum, only a single best-fit solution is obtained by using one of the optimization techniques mentioned above. Rommel et al. [13] used a similarity measure between the measured and the modeled spectrum to select the best combination between all possible combinations of endmembers. The differences between best fit solution for a particular combination of endmembers and the measured spectra for many combinations are, however, quite similar. Tiny changes in the spectrum could result in an entirely different combination to be selected as the best solution. Lapotre et al. [22] showed with a variety of experiments that the uncertainties of the endmember abundances and grain size are relatively large and many solutions are equally likely. On the Moon, these uncertainties are likely to be even higher, with the influence of space-weathering. In Rommel et al. [13] the best-fit solution often contained the ilmenite endmember, in contrast to the absence of ilmenite in the actual mixture. Mineral darkening agents such as ilmenite, which is abundant in some lunar maria, add another layer of complexity to the unmixing as they are almost featureless and, therefore, mostly darkened.
In this work, we develop a Bayesian approach to unmixing. In Bayesian inference (e.g., [23]), all parameters of the model are considered probabilistic variables. In this framework, in addition to the point estimates of the best fit, we also estimate the uncertainties of the parameters. Due to the darkening agents, which darken the spectra and dampen the absorption bands, the uncertainties of the mineral abundances are likely to be higher on the Moon than on other planetary bodies with an atmosphere (e.g., [22]). Additional prior information can be conveniently incorporated in the Bayesian framework. To reduce uncertainties of the mineral abundances, estimates of elemental abundances can be included such that the algorithm will prefer solutions that use endmember combinations with similar elemental compositions.

2. Methods

Generally, the methods employed in this work can be subdivided into two main frameworks. Firstly, the radiative transfer model describes the scattering and absorption behavior of the materials. This framework is built upon the Hapke model [9] and the space weathering model of Wohlfarth et al. [24]. Secondly, the Bayesian framework formalizes the probabilistic perspective of the problem. In the following section, we will describe the choices of the models and their parameters.

2.1. Radiative Transfer Models

2.1.1. Hapke Model

The semi-empirical Hapke model has become the quasi-standard for the analysis of particulate planetary surfaces. Over the years the original model [25] has been refined to also include opposition effects and better approximations for the phase function and Chandrasekhar’s H-function to model multiple scattering inside of the medium. The Modified Isotropic Multiple-Scattering Approximation (MIMSA) formulation of the model [26] is as follows:
r ( μ , μ 0 , g ) = ω 4 π μ 0 μ + μ 0 p ( g ) B S H ( g ) + M ( μ 0 , μ ) B C B ( g ) .
The single scattering albedo ω (SSA) can be calculated by inverting the equation for ω . The cosines of the incidence ( Θ i ) and emittance ( Θ e ) angles are denoted as μ 0 and μ , respectively. The phase angle g is the angle between the incident and emitted light rays. The phase function is denoted as p ( g ) , the opposition effect is subsumed in the terms for shadow hiding B S H and coherent backscattering B C B . We also include corrections for macroscopic surface roughness as described in Hapke [27], but for the laboratory mixtures used in this work, this effect is not relevant. The multiple scattering term M( μ 0 , μ ) models the scattering between the particles inside the medium and is described in detail in Hapke [26] and not repeated here.
In general, the phase function p(g) can be represented by a sum of Legendre polynomials P n ( cos ( g ) ) weighted with the material-specific Legendre coefficients b n :
p ( g ) = 1 + n = 1 b n P n ( cos ( g ) ) .
Specifically, the phase function p(g) of the soil can be approximated with the double Heney–Greenstein function [28]:
p ( g ) = 1 + c D H G 2 1 b D H G 2 1 2 b D H G cos ( g ) + b D H G 2 3 / 2 + 1 c D H G 2 1 b D H G 2 1 + 2 b D H G cos ( g ) + b D H G 2 3 / 2 ,
with the parameter c D H G being a measure for the amplitude of the backscattering lobe to the forward scattering lobe. c D H G > 0 denotes particles that are mainly backscattering. The parameter b D H G is a measure of the width of the lobe. For the lunar regolith, Warell [29] found the standard values for the parameter bDHG and cDHG to be 0.21 and 0.7, respectively. According to Yang et al. [30], ilmenite is more back-scattering while the other lunar minerals are mainly forward-scattering. We adopt the parameters for ilmenite from Yang et al. [30] and set bDHG = 0.19 and cDHG = 0.33 for the ilmenite endmembers. For the other endmembers we use the parameters from Warell [29] to be consistent with the previous work of, for example, Rommel et al. [13].
By expanding the Double Heney Greenstein function the Legendre coefficients can be derived from the parameters bDHG and cDHG as follows [26]:
b n = c D H G ( 2 n + 1 ) b D H G n .
The multiple scattering term M ( μ , μ 0 ) in Hapke [26] is derived from the Legendre coefficients of the single scattering phase function p(g). Using only the first 15 coefficients is sufficient, as the b n decreases with increasing n.
Under small phase angles, the surface appears brighter than would be expected. This is accounted for in the model by the correction terms BSH and BCB, denoting the shadow hiding and coherent backscatter effects, respectively [26]. The samples used in this work were acquired under a phase angle of 30°, therefore, we will assume that the opposition effect is negligible, if not mentioned otherwise.
This model is inverted to find the best fit albedo for the measured reflectance values. This procedure is applied to the endmember spectra that are then used according to the forward mixing model.

2.1.2. Forward Mixing Model

The single scattering albedos of the endmembers obtained by inverting the Hapke MIMSA model are subsequently used to calculate the SSA spectra of the mixture. According to Hapke [9] the single scattering albedo is defined as the scattering coefficients (S) divided by the extinction coefficient (E) of the mixture,
ω m i x = S E = i N i σ i Q S , i i N i σ i Q E , i ,
with N i as the number of particles of particle type i per unit volume and σ i = π a i 2 as the cross-sectional area, where a i denotes the radius of the corresponding particle type [9]. The scattering efficiency is Q S , i = Q E , i ω i , such that the mixing formula becomes:
ω m i x = i M i Q E , i ρ i D i ω i i M i Q E , i ρ i D i .
The bulk density of particle type i is denoted Mi and the average particle diameter is Di. Assuming that the endmembers have the same average grain size D i and the same extinction efficiencies Q E , i , the fraction can be reduced. The albedo of the mixture can then be expressed as a linear superposition of the endmember SSA spectra ω i weighted with the Hapke coefficients θ e m , i ,
ω m e d ( θ e m ) = i θ e m , i ω e m , i .
For a mixture of typical lunar minerals, the extinction efficiencies, grain size, and density are constant and approximately equal among the endmembers; therefore, these simplifications will be used to calculate the SSA spectrum of the mixed mineral endmembers. This spectrum is then subsequently used as the input to the space weathering framework to calculate the reflectance spectra of the mixture.
When the mixture directly results from the particular known endmembers used in the unmixing, the sum of the weights should be equal to one ( i θ e m , i = 1 ). In this work we are not enforcing this constraint, but prefer solutions where the sum is close to one. By not enforcing the sum-to-one constraint for the fractional abundances θ e m , i the sum can be seen as a “catch-all” parameter for all linear, wavelength independent, differences among the endmembers and compared to the examined sample. It is, however, physically implausible to have endmembers with a negative contribution. Therefore, all endmember weights must be larger than zero, that is, θ e m , i > 0 .
Apart from the mineral composition and maturity of the sample, other physical properties and the laboratory setups themselves also influence the measured spectra. In this work, we use samples and endmembers of the same grain size fraction to remove the influence of the grain size that has a strong effect on the overall spectra (e.g., [31]). The grain size also introduces large uncertainties for immature spectra [22] and would, therefore, influence the dependencies of the parameters we want to examine in this work.
The porosity or compactness [32] of a sample changes the brightness of a sample almost without introducing a change to the continuum slope. Therefore, this parameter introduces a scaling factor on the weights of the endmembers depending on the compactness of the sample.
Depending on the laboratory setup, other factors might influence the measured spectra. For example, the gain of the sensor might be different such that a linear offset is introduced, when comparing the measurements of the same mineral for different laboratories.

2.1.3. Space Weathering

Space weathering poses a major challenge for unmixing mineral abundances. In addition to the darkening and reddening, the diagnostic absorption features are obscured by the optical influence of the space weathering particles. This makes it harder to differentiate between the different minerals. While the band position remains constant the depth of the absorption that could normally be used to estimate the fraction of the constituent minerals is less clear.
Wohlfarth et al. [24] defined a framework based on ab-initio Mie modeling to characterize the effects of the space weathering particles. We include these characteristics in the unmixing framework. It is important to note that extinction efficiencies are not constant over the wavelength and are very different for the mineral endmembers and for each of the iron particles (see Figure 1a). As a result, the fraction in Equation (7) cannot be reduced. The SSA spectrum of the mixed mineral endmembers ( ω m e d ) according to Equation (7) is then used as the starting value for the mixture with the iron particles. It is here assumed that all endmembers have a similar grain size distribution, density and extinction efficiencies. The values are adapted from Wohlfarth et al. [24] and listed in Table 1. The mixture of the iron particles and the medium is then defined as:
ω t o t a l ( ω m e d ( θ e m ) , M n p F e , M m p F e ) = ( 1 M n p F e M m p F e ) Q E , m e d ρ m e d D m e d ω m e d + M n p F e Q E , n p F e ρ n p F e D n p F e ω n p F e + M m p F e Q E , m p F e ρ m p F e D m p F e ω m p F e ( 1 M n p F e M m p F e ) Q E , m e d ρ m e d D m e d + M n p F e Q E , n p F e ρ n p F e D n p F e + M m p F e Q E , m p F e ρ m p F e D m p F e .
The phase function of a mixture is also defined by the weighted sum of the phase functions of the components. The phase functions of the soil and the small npFe0 particles are independent of wavelength, but the larger mpFe0 particles have a wavelength dependent phase function (see also Figure 2), which is, therefore, given as a vector.
p t o t a l ( g ) = ( 1 M n p F e M m p F e ) Q E , m e d ω m e d ρ m e d D m e d p m e d + M n p F e Q E , n p F e ω n p F e ρ n p F e D n p F e p n p F e + M m p F e Q E , m p F e ω m p F e ρ m p F e D m p F e p m p F e ( 1 M n p F e M m p F e ) Q E , m e d ω m e d ρ m e d D m e d + M n p F e Q E , n p F e ω n p F e ρ n p F e D n p F e + M m p F e Q E , m p F e ω m p F e ρ m p F e D m p F e .
The phase functions are defined by the Legendre coefficients as in Equation (2). These b n are calculated based on the refractive indices, grain size, and the Legendre expansion of the Mie phase function from Fowler [33]. By equating the coefficients it can be shown that the b n of the mixture are also the weighted b n of the soil, npFe0, and mpFe0 particles. The Legendre coefficients of the mixture are then used to calculate the wavelength dependent phase function of the mixture and are additionally used in the multiple scattering term M( μ , μ 0 ). Finally, we obtain reflectance spectra of the mixture, according to Equation (1).

2.1.4. Endmember Catalogs

In order to generate spectra similar to lunar conditions, we need to select endmembers representative of the general mineralogy of the Moon. We are using only mineral endmembers and are not trying to model the influence of agglutinate glasses. This glassy fraction usually has a similar composition to the mineral fraction [34], as it is created endogenously. The absorption band depth is also not affected by glasses, but by a large amount of smFe0 in the agglutinates [35]. Agglutinates here are only relevant for the unmixing of the LSCC catalog in Section 3.3, but this should be kept in mind when evaluating the results.
The highlands are dominated by plagioclase (around 90 wt.%) with the remaining fraction mainly consisting of pyroxenes and minor phases of olivine and spinels [1]. The composition of the maria is more diverse. About 30–40 wt.% are pyroxenes and plagioclase still forms a major part of the composition with up to 60 wt.% of the mineral composition. Olivine contributes up to 10 wt.% and the opaque mineral ilmenite is abundant in some mare areas with up to 20 wt.% at maximum [1].
Among these groups of minerals are also variations depending on the elemental composition. For example, Sun and Lucey [36] showed that by generating albedo spectra based on the magnesium number and grain size the abundance estimation of olivine and pyroxenes can be improved. We will use fixed albedo spectra as endmembers that are chosen to be a good representation of a wide variety of possible true endmembers while retaining the most distinct endmembers.
In general, we want to have one representative of each of the major minerals on the Moon. Pyroxenes can be subdivided into low calcium pyroxenes (LCP), mainly orthopyroxenes (OPX), and high calcium pyroxenes (HCP), usually clinopyroxenes (CPX). The LCP pyroxenes exhibit an absorption 0.9- μ m and 1.9- μ m [37,38]. The HCP pyroxenes can be further subdivided into two spectrally distinct groups. Type A, with an absorption at 0.9- μ m and at 1.15- μ m, but no absorption around 2- μ m [37] similar to olivine. Type B, exhibits absorption bands at 1.05- μ m and 2.35- μ m [37]. The absorption band positions, however, are not fixed, but vary with Mg2+, Fe2+, and Ca2+ content [39,40]. For the unmixing, the distinction between these three types of pyroxenes (HCP, LCP-A, and LCP-B) is, however, beneficial, because we want to minimize the linear dependencies between the endmember spectra.

2.2. Bayesian Inference

From experience of using unmixing algorithms based on an exhaustive search of all possible endmember combinations (e.g., [13]), we know that several possible combinations can produce similar results. The solution is, therefore, not unique. There might be multiple local minima and some endmembers that are similar, when they become mature, are even gradually interchangeable, without significantly reducing the quality of the fit. Classical unmixing techniques, that are based on finding the best fit in the least squares sense, ignore these uncertainties and only take a single solution, which by a somewhat arbitrary similarity measure best describes the measured spectrum. By exploring the parameter space more thoroughly we can obtain uncertainties and can also learn more about the dependencies among the model parameters.
The process of estimating the probability distributions of parameters of a model based on a set of data and additional assumptions or knowledge about the data is called Bayesian inference (see e.g., [23]). All parameters of the model are described as probability distributions. The goal, therefore, is to estimate the probability densities of the parameters Θ given some data D. This is called the posterior distribution p ( Θ D ) . According to Bayes’ rule the posterior density can be described as:
p ( Θ D ) = p ( Θ ) p ( D Θ ) p ( D ) .
The data are not changing for one observation, therefore, we can assume that p ( D ) is constant. Consequently, we can define the unnormalized posterior density [23] as:
p ( Θ D ) p ( Θ ) p ( D Θ ) .
The probability of the parameters p ( Θ ) is the prior distribution and can introduce already known information about the distribution of the parameters into the procedure or it can be uninformative. The probability of the data given a set of parameters Θ is denoted as the likelihood ( p ( D Θ ) ). The likelihood is evaluated as a function of Θ and gives us an estimate on how likely the data have originated from the current parameters.
This unnormalized posterior distribution can be calculated analytically if the probability distributions are chosen in a specific way, which might be limiting in practice. With the advancement of computer technology and the introduction of Markov Chain Monte Carlo (MCMC) methods [41], such as the Metropolis-algorithm [42], it is possible to estimate complex probability distributions numerically. The algorithm used in this work to sample from the posterior density distribution is described in the following section.

2.2.1. Sampling Algorithm

The goal of the Monte Carlo simulation is to approximate the posterior distribution. Most of the sampling algorithms as well as our approach are based on the Metropolis algorithm [42], which was subsequently generalized to non-symmetrical proposal distributions and is then called Metropolis–Hastings algorithm [43]. This algorithm describes a random walk through the posterior space such that for a large number of samples, the samples converge to the posterior distribution. Generally, at each step of the Markov Chain a proposal ( Θ t π ( Θ t 1 ) ) is randomly generated based on the current location in parameter space. This proposal is then either accepted or rejected and the current parameters are consequently used for the next proposal. If the proposal is more probable than the current sample, then the proposal is always accepted; if the proposal is less likely than the current sample, it is accepted with a certain probability. Therefore, less probable solutions are also sometimes accepted.
The random walk behavior of the Metropolis algorithm can, however, lead to slow convergence because many proposals are rejected (e.g., [23]) and correlated parameters can make the algorithm much less efficient [44]. Modified versions of the Metropolis algorithm improve the proposals by, for example, estimating the covariance matrix on the previous examples and scaling the proposal distribution accordingly (e.g., [45]) or by additionally proposing an alternative step that is chosen based on the rejection of the first proposal [46]. However, it still remains that these methods tend to be less efficient for non-linearly correlated parameters, due to their inherent random walk behavior [44].
Hamiltonian Monte Carlo (HMC) borrows ideas from Hamiltonian dynamics [47] and introduces a ’momentum’ variable. As in Hamiltonian dynamics, the momentum variable and the position are building a dynamic system. The entire sampling process is subdivided into simulations of L leapfrog steps each of length ϵ . To calculate the momentum at a position in posterior space we need to calculate the derivative of the posterior density for the local environment. With the advances in automatic differentiation in tools like, for example, Theano [48], which is used in the pymc3 framework [49], it is possible to quickly calculate the derivatives of more complex functions.
The two parameters of HMC (L and ϵ ), however, must be tuned by hand and strongly influence the performance of the algorithm. The parameter L determines the length of one simulation of Hamiltonian dynamics and, if chosen too large, the leapfrog steps tend to do a U-turn and return to the starting location [23]. The No U-Turn Sampler (NUTS) [44] is an extension to HMC that automatically stops one simulation of Hamiltonian dynamics when the leapfrog step would not increase the distance to the starting location and, therefore, removes the need to tune the L parameter by hand [44]. This way, the NUTS sampler can effectively explore the posterior space.
This algorithm is likely the best choice for our problem, because we can effectively sample from a posterior distribution where we expect some parameters to be correlated and the target distribution might be multimodal.

2.2.2. Likelihood

The unmixing problem according to Section 2.1 is that the abundances ( θ e m ) of the endmembers and the space weathering components npFe0 and mpFe0 ( θ s m F e ) with known albedo spectra, extinction efficiencies and phase functions must be determined. They are the parameters we estimate with the help of MCMC. The known data that is used to calculate the reflectance spectra are denoted as X . The likelihood/goodness of fit of the current parameter estimates is evaluated with both the reflectance spectrum ( r m i x ) and the continuum removed spectrum ( r c r ) to emphasize the absorption band depths over the absolute fit. The continuum removed spectrum is calculated in a straightforward way such that certain fixed wavelengths are chosen as the points of the continuum (see, Figure 3). Usually, the continuum is defined by the convex hull of the spectrum (e.g., [50]). Calculating the continuum in this way is very time consuming considering that the function and its derivative must be evaluated thousands of times. The convex hull is not a linear function such that it is inefficient to implement it in Theano because the derivatives cannot be determined easily. The indices of the local maxima would be dependent on the current spectrum. Because for each of the current parameters of the Markov Chain the continuum is different and depends on the slope, it cannot be computed prior to the sampling. Both the method using the convex hull and the method using fixed wavelengths used in this work are able to characterize the absorption bands usually found on the Moon.
Therefore, the concatenation of the spectrum and the continuum removed spectrum represents the data S. The output of the model is, therefore, a function of the parameters θ e m and θ s m F e such that:
S m i x ( θ e m , θ s m F e ) = [ r m i x , r c r ] .
The likelihood is then evaluated by building the sum over all log likelihoods or the product of all probability density function values at the measured data ( S ), given a normal distribution around the modeled data and some variance σ 2 , which is also considered a parameter of the model and, therefore, is also sampled.
p ( S θ e m , θ s m F e , X ) j = 1 2 K N ( S j μ = S m i x , j , σ 2 ζ j ) .
K denotes the number of channels. If we know that the model’s accuracy is wavelength dependent, it is reasonable to include this information in the likelihood function. Therefore, we introduce a wavelength dependent scaling factor ( ζ ) that can be used to emphasize, or de-emphasize certain wavelength channels.

2.2.3. Prior

Bayesian inference offers a convenient way to include a-priori knowledge in the sampling process. As described in Equation (11), the posterior is proportional to the likelihood and the prior p ( Θ ) . For multiple priors, the joint probability is given by the product of the individual probability. Because our parameters are usually limited to the interval of [0, 1], a good choice for the priors on our parameters is the Beta distribution, which is defined by two shape parameters α and β as follows [51]:
B e t a ( x α , β ) = Γ ( α + β ) Γ ( α ) Γ ( β ) x α 1 ( 1 x ) β 1 ,
with Γ denoting the gamma-function. Choosing α and β controls the shape of the Beta distribution. Some examples are shown in Figure 4. For example, for a distribution with α = 2 and β = 5 , the mode of the distribution is 0.2, but for a distribution with α = 1 and β = 3 , the mode is zero and the probability decreases for increasing values. For values outside the interval [0, 1], the probability is zero.
The prior distribution on the endmember abundances can be set as follows, with the shape parameters being chosen depending on the application:
p ( θ e m ) i = 1 N B e t a ( θ e m , i α = α e m , i , β = β e m , i ) .
The space weathering components, namely npFe0 and mpFe0, form only a small fraction of the composition. Therefore, we can choose the prior similar to α = 1 and β = 3 in Figure 4 for each abundance of npFe0 and mpFe0 ( θ s m F e )
p ( θ s m F e ) i = 1 2 B e t a ( θ s m F e , i α = α s m F e , i , β = β s m F e , i ) .
The shape parameters will be chosen such that small values are preferred by the sampler.
We do not enforce the sum to one constraint to account for differences in the laboratory setups of the endmembers as well as porosity and other influences. Nonetheless, it is sensible to set a prior on the sum of the endmember weights θ e m :
p ( θ e m ) N ( θ e m μ = 1.0 , σ s u m 2 ) .
This prior is centered around one such that solutions where the sum of the endmembers is close to one are favored. The variance σ s u m 2 is chosen depending on the differences among the endmembers and the investigated samples.
The elemental abundances can be relatively well characterized remotely with, for example, gamma ray spectroscopy (e.g., [52]). Spectral parameters in the UV (e.g., [3,4]) or in the NIR (e.g., [5]) can also be used to estimate the elemental abundances, but should be used with care. For example, according to the TiO2 maps of Sato et al. [4] the swirl Reiner Gamma is deprived of TiO2 compared to its surroundings, which is physically implausible, as the optical differences are likely due to a difference in maturity and/or compaction [53]. Setting a prior, however, does not rule out solutions that do not confirm with prior knowledge, as long as the probability of the prior does not become zero. If, for example, the spectrum fits better for a high ilmenite concentration, this solution will be found to be more likely.
The RELAB database provides elemental abundances for many of the samples (e.g., [37]) and the elemental composition of the endmembers from Rommel et al. [13] were determined by electron microprobe analysis. We here use endmembers with well-defined elemental abundances ( E e m , M × N matrix, for M elemental abundances and N endmembers) and use the estimated elemental abundances of the samples ( e s a m p l e , M × 1 ) as prior information. If the samples are mixed directly from the endmembers with weights c ( N × 1 ) the actual elemental abundances can be calculated as
e s a m p l e = E e m c .
For realistic mixtures, this relationship produces highly non-unique solutions for the mineral weights because, for example, the MgO abundance could be attributed to a pyroxene or an olivine. For each sample the elemental abundances are calculated based on the proposed parameters of the mineral abundances θ e m according to:
θ e l e m = E e m θ e m i = 1 n θ e m , i .
As the beta distribution should be centered around the estimated elemental abundance we want to enforce the mode of that distribution to be equal to the estimated elemental abundance. The mode of the Beta distribution for α > 1 and β > 1 is given by:
mode ( x ) = α 1 α + β 2 = ^ e s a m p l e .
By increasing α ( α 1 ), the variance of the distribution decreases. Therefore, we choose a α e l e m approximately according to the desired variance of that particular element and calculate β e l e m according to Equation (20) such that:
β e l e m = α e l e m 1 e s a m p l e + 2 α e l e m .
Consequently, we define the prior for the elemental abundances to be:
p ( θ e m ) i = 1 M B e t a ( θ e l e m , i α = α e l e m , i , β = β e l e m , i ) .
This distribution then effectively acts as a prior on the weights of the mineral endmembers, but there are multiple combinations that fit equally well.

3. Results

The measured spectrum of a mixture is influenced by several factors, namely (1) composition (2), space weathering, (3) grain size, and (4) porosity. To disentangle the effects, we designed several experiments. The grain size is constant for all experiments and is expected to be known. Firstly, we calculate mature synthetic spectra to remove effects not described by the mixing model itself, except for noise. This can also be interpreted as a general validation of the sampling approach. Secondly, to remove the influence of space weathering and to test the general viability of the mixing model we used laboratory mixtures of known endmember spectra. Thirdly, we apply the full framework to unmix the LSCC spectra with endmember spectra from the RELAB catalog and include space weathering particles. These LSCC spectra can be seen as a good representation of different maturity levels and the general mineralogy of the lunar maria.

3.1. Synthetic Mixtures

This experiment is designed to estimate the uncertainties of the mixtures when all influences not covered by the model are removed. Spectra are generated with the mixing model including additional Gaussian noise with a standard deviation of σ n = 0.0008 . This corresponds to a signal-to-noise ratio (SNR) between approximately 170 and 400 (with SNR = μ r , c h σ n ) depending on the channel.
We selected the endmembers according to the categorizations above (see Section 2.1.4). A representative of each of the classes LCP, HCP Type A, HCP Type B, Olivine, Plagioclase, and Ilmenite is selected from the RELAB catalog (https://pds-geosciences.wustl.edu/spectrallibrary/default.htm, accessed on 1 November 2020) and the SSA values are obtained by inverting the Hapke model and displayed in Figure 5. All endmembers are samples of the 0–45 μ m size fraction. The elemental abundances of the endmembers are listed in Table A1.
The generated mixtures are based on a mare composition of the endmembers of 10 wt.% hypersthene, 40 wt.% augite, 5 wt.% olivine, 30 wt.% plagioclase, and 15 wt.% diopside. The ilmenite abundance is varied between 0 wt.% and 16 wt.% and the remaining fraction is normalized, such that the endmember weights again sum up to unity. The abundances of npFe0 and mpFe0 is varied between 0 wt.% and 2.2 wt.%, and 0 wt.% and 1.2 wt.%, respectively. Ultimately, 48 synthetic mixtures (mix0–mix47) are created. The detailed results for all mixtures are listed in Table A3 and Table A4.
Generally, all endmembers are measured under similar conditions. According to these conditions the priors have to be chosen. In our case, we do not want to directly employ knowledge about mineral abundances ( θ e m ). Therefore, we choose an uninformative prior, that is:
p ( θ e m ) i = 1 N B e t a ( θ e m , i α = 1 , β = 1 ) .
Both shape parameters of the Beta distribution are set to one for the prior distribution, such that the prior is a uniform distribution between zero and one. Because negative weights are physically implausible, these samples are always rejected.
The mixtures are synthetically generated, such that we set σ s u m = 0.02 to prefer solutions where the sum of the fractional abundances is close to one. The experiment is carried out once with no prior on the elemental abundances and once with an elemental prior on the abundances of TiO2 and Al2O3 (see also Table 2). For theoretical abundances below 1 wt.% the variance becomes very low when using the method mentioned above. Therefore, we set α e l e m = 1 and β e l e m = 10 for the shape of the prior, making the mode zero and the PDF is given in Figure 6.
In general, we set the shape parameters for npFe0 and mpFe0 to α s m F e = 1 and β s m F e = 30 . These parameters were tuned by hand such that 80% of the probability density is distributed in the interval [0 wt.%, 5 wt.%]. Higher values than 5% are physically implausible. Trang and Lucey [54] found about 2 wt.% of smFe at the maximum for the Moon, excluding the polar regions. The model of Trang and Lucey [54], first versions of which were introduced in Trang et al. [55] and Lucey and Riner [16], differs from the model of Wohlfarth et al. [24], for example, Trang et al. [55] do not consider the influence of the phase function. Therefore, the predicted concentration of npFe0 and mpFe0 is not necessarily the same, leading us to set a relatively uninformative prior on each of the smFe abundances. The prior on the abundance of smFe particles is consistent throughout this work and is not changed between experiments.
Generally, the mean predicted solution does not necessarily correspond to the actual mixture coefficients used, but the theoretical values are always part of the solution. Classical unmixing approaches are limited to just one solution and small changes may lead to different optimal solutions, which are also not necessarily equal to the ground truth values.
Ilmenite shows, compared to the mafic minerals, an almost featureless spectrum in the NIR wavelength range. Because of its low albedo, the abundance of ilmenite also dampens the absorption bands and darkens the overall spectrum. This is, to some extent, similar to the spectral influence of mpFe0 particles. If the absolute reflectance is poorly defined, so is the ilmenite abundance. Consequently, depending on the prior setup and the sum constraint the uncertainties of ilmenite are different. If no prior on the elemental or mineral abundances is set, the abundance of ilmenite is relatively poorly defined. For the mixtures that, in theory, do not contain ilmenite, a majority of the accepted mixtures contain significant amounts of ilmenite (see Figure 7). A combination of ilmenite and other endmembers can, therefore, equally explain the measured spectrum. This effect can be mitigated when a prior is used (see Figure 8). For the fresh mixture (mix0), without any smFe particles, the uncertainty of the spectrum is small (see Figure 9). The mineral uncertainties are small, but for the other mixtures, the uncertainties significantly increase. The majority of posterior density of plagioclase (PLG) is generally lower than the actual abundance, and for olivine (OLV) and/or diopside (CPXA) the majority of posterior density is below the actual value. This suggests that a mixture of olivine/diopside and ilmenite produces very similar results to plagioclase. Olivine and diopside have similar absorptions at about 1- μ m, making them interchangeable to some extent when the maturity increases. For example, for mix15 there are two distinct modes with either high concentration of diopside or olivine. The mineral uncertainties of mix5 are large, but the spectrum is clearly defined. This highlights that there are several solutions that produce nearly identical spectra. Sometimes, however, the modal abundance values do not sum up to one, because, for example, just one of the distributions is skewed. This can lead to mode predicted spectra outside of the confidence band. The modal abundance values were likely never sampled in this particular combination.
Ilmenite is the only endmember that contains significant amounts of TiO2. Including prior knowledge about the abundance of TiO2 can, therefore, favor solutions that are close to the true ilmenite value. In Figure 8 the violin plots of the acceptable solutions considering the TiO2 and Al2O3 prior are displayed. The abundance of ilmenite is then much more clearly defined. The uncertainties of the other minerals also mostly decrease, while the uncertainties of the spectra remain similar.
The results for the most mature spectrum with the highest abundance of ilmenite are shown in Figure 10. With the prior, the uncertainties of all endmembers can be reduced. The confidence plots remain very similar. Figure 11 shows the distribution of the differences between theoretical values and mean predicted abundance. By including the prior the differences can be reduced. The mean of olivine is usually above the true value, while the mean of plagioclase and clinopyroxene is usually below the theoretical value. The theoretical abundance of olivine is relatively small compared to the other minerals, consequently, the mean is usually higher. The uncertainties of olivine and clinopyroxenes are high for both with and without the prior. This can be seen in Figure 12. Including prior knowledge about the elements Al2O3 and TiO2 the uncertainties of all mineral abundances are mitigated. This is despite the fact that these elements are almost depleted in the pyroxenes and olivine.
The abundance of npFe0 and mpFe0 is usually well defined, as changes to the abundance of npFe0 and mpFe0 can have a strong influence on the modeled spectrum. In Figure 7, Figure 8 and Figure 10 it can be seen that the predicted npFe0 and mpFe0 do not change significantly with the introduction of a prior. In Table 3 the correlation coefficients between true smFe abundance and the standard deviation of the minerals is listed. For all minerals it can be seen that the correlation coefficients are positive. Therefore, the higher the abundance of smFe is, the larger the uncertainties become. The uncertainty of ilmenite is, however, uncorrelated to the smFe abundance, even if no prior knowledge is used. Because ilmenite does not show a prominent absorption band such a feature cannot be obscured by the smFe particles. This is visible in Figure 13. While the overall uncertainties decrease with a prior on the elemental abundances, the trend that the uncertainties increase with increasing maturity is visible for both with and without a prior.
Generally, it can be seen that the ilmenite abundance is difficult to determine reliably without prior information, for both immature and mature spectra. For mature spectra, the uncertainties of all mineral abundances become very large, even for an almost perfect case, where influences which are not modeled, are removed. This illustrates that unmixing on the Moon is an ill-posed problem, and additional knowledge about the composition is necessary. The addition of a prior to indirectly include information about the plagioclase and ilmenite abundance is effective in reducing some of the uncertainties.

3.2. Fresh Laboratory Mixtures

This experiment is set up to unmix real laboratory mixtures (see also Rommel et al. [13]) from laboratory endmembers. Thus, we select the same endmembers as in Rommel et al. [13] but omit one of the two labradorites (plagioclases) as the two are very similar, and we want to select endmembers that are spectrally distinct. We also omit the diallagite endmember, to reduce the catalog to the same size as for the other experiments and this endmember is only present in one mixture, which is, consequently, also omitted here. The amphibole endmember pargasite (PRG) was selected in Rommel et al. [13] to test whether uncommon minerals can also be detected by the algorithm. We keep this endmember for our analysis even though it is not present on the Moon.
For the inversion of the Hapke model, the parameters are chosen in the same way as in Rommel et al. [13] for comparability. Therefore, the term for the shadow hiding opposition effect is not set to zero, but B S H 0 = 3.1 and h S H = 0.11 are adopted from Rommel et al. [13]. The resulting SSA spectra are shown in Figure 14.
The elemental abundances of the endmembers were determined in Rommel et al. [13] based on electron microprobe analysis in the laboratory (see Table A5) and the modes and shape parameters are determined according to Section 2.2.3. The results were obtained by using no elemental priors on the one hand and using a prior for the TiO2 abundance on the other hand (see first two rows of Table 2).
The true values of the abundance of submicroscopic iron particles in this experiment is zero. In order to test whether fresh and mature spectra can be reliably distinguished from each other, we still employ the same prior for smFe as for the other experiments.
All endmembers and mixtures were measured in the same laboratory and the mixtures are directly mixed from the endmembers; therefore, the weights of the endmembers should sum up to a value close to one. We set a prior on the sum according to Equation (17) with a standard deviation of σ s u m = 0.02 .
The detailed results of the unmixing without a prior on the elemental abundances and with a prior on the TiO2 abundance are listed in Table A7 and Table A8, respectively. Generally, the mean predicted results are very similar to those of Rommel et al. [13] considering that the endmembers labradorite A and diallagite were omitted. Even though we did not use a strong prior on the abundance of iron particles, the predicted abundances are almost exclusively below 0.02 wt.% or a fraction smaller than 2 × 10 4 .
For the version without a prior on the TiO2 abundance, ilmenite is erroneously predicted in the mixtures of pargasite and plagioclase. For example, Figure 15 shows the violin plots of all similarly well-fitting solutions. It can be seen that ilmenite is predicted to be present, and a solution without ilmenite is not part of the uncertainty. Therefore, a solution with ilmenite fits better than one without ilmenite. If we, however, prefer solutions where the contribution of ilmenite is close to zero, we obtain a better result and the errors in pargasite and plagioclase abundance are reduced. In Figure 16 the confidence plots are plotted. Here it can be seen that the mean and mode predicted spectra fit quite well for both prior setups. By introducing the prior for the TiO2 abundance, the uncertainties of the fit increase ( σ = 0.006 compared to σ = 0.003 for no prior). Therefore, slightly less well fitting solutions are accepted if the ilmenite abundance confirms better with the prior knowledge.
In the mixtures containing plagioclase and ferrosilite, plagioclase is confused with pargasite if no prior knowledge is included. The two endmembers are spectrally very similar and are both bright (see Figure 14). In Figure 17, an example violin plot is shown, where the confusion of plagioclase and pargasite can be seen in Figure 17a. Small contributions of ilmenite can already dampen the absorption band of the pargasite endmember, such that the similarity to plagioclase increases. The mean predicted abundance of ilmenite is quite low with 0.48 wt.% for the version without a prior, but still leads to a significantly different solution. The confidence plots in Figure 18 look nearly identical, except that for the TiO2 prior to the uncertainty of the 0.9 μ m absorption band slightly increases. This suggests that both solutions, with and without ilmenite, represent regions of high posterior density. It is, however, very difficult for the sampler to traverse from one region to the other, as the two are far apart and due to the prior on the sum it is unlikely for a limited number of samples that both regions are sampled. However, with a relatively uninformative prior as in Figure 6 this ambiguity can be removed, such that the solutions with no ilmenite are preferred. This trend illustrated by the two examples above, can also be seen in Figure 19. The differences between mean predicted solution and true abundance for the endmembers augite, olivine, and ferrosilite are not improved by including prior knowledge about the TiO2 abundance. This is likely because of inaccuracies of the data and the simple mixing model and also a problem in Rommel et al. [13]. The inherent simplifications of the Hapke model makes it such that the fractional abundances are not necessarily equal to their actual weight fraction (e.g., [56]). Similar to a regularization term in classical optimizations, we can use the prior to favor certain solutions, without imposing hard constraints. Generally, without prior knowledge the abundance of plagioclase is usually underestimated, while the abundance of ilmenite and pargasite is overestimated. The uncertainties are mostly small, both with and without TiO2 prior (see Figure 20).

3.3. LSCC Samples

In the final experiment, we unmix spectra from the LSCC catalog with laboratory mineral endmembers and the npFe0 and mpFe0 particles. The LSCC catalog contains mature samples that were returned from the lunar surface by the Apollo missions and is, therefore, a good representation of lunar mineralogy. We use the <45 μ m size fraction spectra for the unmixing, and because there is no data on the mineral abundances of this size fraction we use the mineral abundances from the 20–45 μ m size fraction, which correspond to a large portion of the weight, as a reference for evaluating the results. The endmembers are the same samples from the RELAB library listed in Section 3.1. These endmembers do not correspond to the actual constituent minerals and there are additional minerals that are not covered by the endmember catalog, for example, spinel. Additionally, the endmember catalog does not contain agglutinates in comparison to the LSCC samples with large fractions of agglutinates. The ground truth abundances are obtained from the tables of Taylor et al. [18,19]. The abundances of the minerals were grouped into clinopyroxene, orthopyroxene, olivine, plagioclase, and ilmenite. The abundances were then normalized to the mineral fraction because we are not considering glasses. We have seen that, without a prior on the elemental abundances, the uncertainties become very large for mature spectra. Therefore, we will use the same prior once on TiO2 and Al2O3 abundances, which in the previous experiments reduced the uncertainties effectively for plagioclase and ilmenite, but also for many of the other minerals. Additionally, we will employ a prior on all elements with a relevant abundance on the Moon, which might contribute to distinguishing between minerals (see also Table 2). The elemental abundances of the LSCC samples are listed in Table A9. The spectra were measured under laboratory conditions, but the mixtures are much more complex compared to the laboratory mixtures from Section 3.2. We are, therefore, not enforcing the sum-to-one constraint as strongly, as we set σ s u m = 0.1 . Because plagioclase and ilmenite have both almost featureless flat spectra, the uncertainties between these two increase significantly the more relaxed the prior on the sum is. Therefore, the prior knowledge about TiO2 and Al2O3 abundance becomes even more essential.
Detailed results for the unmixing of all LSCC spectra are listed in Table A10 and Table A11 for the two experiments with different priors. Generally, the uncertainties are similar between both versions, but for some samples the priors on the other elements help in reducing the uncertainties. This can be seen in Figure 21. Therefore, we will focus on the Ti/Al prior version.
In some cases, the plagioclase abundance is underestimated; such an example can be seen in Figure 22. For the violin plots, the two clinopyroxenes are grouped together. The spectrum fits quite well and the true plagioclase abundance is not part of the uncertainty. Therefore, a prior cannot mitigate the differences between mean predicted and theoretical abundance. Instead, olivine and orthopyroxene abundances are overestimated. A similar problem can be seen for 12001 and 15041. These samples are the ones with the highest fraction of agglutinic glasses (12001: 56.2 wt%, and 15041: 51.3 wt%, and 10084: 53.9 wt%) among the mare samples. For all other samples, the predicted abundances are close to the true values.
Figure 23 shows the typical results for a mature mare spectrum. Most of the abundances are relatively well defined with the prior on the TiO2 and Al2O3 abundance. Compared to sample 10084 in Figure 22, there are significantly fewer npFe0 and mpFe0 particles predicted to be present. As sample 14260 contains less FeO (9.65 wt.% compared to 14.8 wt.% in 10084) it is likely that the smFe abundance saturates more quickly, such that additional smFe cannot easily be created as quickly anymore. The ilmenite abundance is comparatively small, which might also be a contributing factor to the uncertainties being relatively small. Similarly, for the mature highland spectrum in Figure 24, the iron abundance in the highlands is generally much smaller so that less smFe can be created. According to the Is/FeO parameter, this is also a very mature spectrum and the absorption bands are weakly pronounced. However, the smFe abundance is small compared to the mare samples. In Figure 25, the results for sample 79221, with both the Ti/Al prior and also with the prior on all major elemental abundances, is shown. This sample is rich in ilmenite and also relatively mature, therefore the uncertainties are quite high. In such cases, the results can be improved when additional elemental abundance information is available. For most other LSCC samples the Ti/Al prior is, however, sufficient to constrain the procedure. The npFe0 and mpFe0 abundance is not much effected by the change in prior knowledge.
In Figure 26 the correlations between Is/FeO or Is and mean predicted npFe0 or mpFe0 abundance are displayed. The abundances are only weakly correlated with the Is/FeO maturity index. In contrast the correlation between the ferromagnetic resonance Is and the abundance is very strong. Even for the small dataset this correlation is significant. Because the Is is a measure on the amount of iron particles in the size fraction between 4 nm and 33 nm [57], this suggests that the abundance of npFe0 can be obtained relatively reliably even with a limited endmember catalog. This difference between Is/FeO and Is also supports the findings of Trang and Lucey [54] that the saturation limit of npFe0 particles is dependent on the FeO abundance. The saturation limit increases with increasing FeO abundance. Figure 27 shows the relationship between npFe0 and mpFe0 abundance in the LSCC samples. Usually, the factor between the abundance of npFe0 and mpFe0 abundance is approximately constant, where the amount of npFe0 is typically twice as high as the amount of mpFe.

4. Discussion

Generally, the three experiments show that a probabilistic perspective provides important insights into the unmixing process overall. When trying to unmix lunar spectra many factors have to be considered to obtain reliable results. A non-probabilistic approach needs to make many assumptions about the composition and the maturity of the surface to obtain plausible results. Even if these assumptions are accurate only a single best fit solution is obtained. All experiments shown in this work illustrate that there is no unique solution and that a variety of combinations can produce similar results and, thus, a number of different interpretations are possible.
The synthetic mixtures in Section 3.1 were used to estimate the influence of the darkening agents on the overall uncertainties of the mineral abundances. The results showed that with increasing smFe abundance the uncertainties increased. Ilmenite or TiO2 also contributed to uncertainties in the data, but there appears to be no direct relationship between true ilmenite abundance and uncertainties. However, introducing prior knowledge about the TiO2 and Al2O3 abundance can significantly reduce the uncertainties of all mineral abundances as well as the errors between mean predicted and true mineral abundances and is, therefore, suitable for the application of the approach to real lunar spectra. For very mature spectra the uncertainties become very large. This illustrates the need for a probabilistic approach to the problem of unmixing spectra for mineral and smFe particle abundances. The abundance of the smFe particles could, however, still be reliably detected for any maturity level. The uncertainties of the minerals all increased with increasing smFe concentration, except for ilmenite. This might be explained by the fact that ilmenite with its featureless spectrum is already difficult to determine accurately if the absolute brightness of the sample is not well defined. Ilmenite, therefore, is not easily distinguished by any absorption bands, which become more difficult to quantify with increasing maturity.
Ambiguities between spectrally similar endmembers such as plagioclase, pargasite, and ilmenite in Section 3.2 can lead to larger errors (up to 57 wt.%) in the best fit solution of the Hapke model and the measured spectra. When these ambiguities are mitigated with using prior information, the errors are usually well below 15 wt.%. All models have to include some simplifications such that they cannot perfectly represent reality. These problems inherent to the mixing model itself usually cannot be predicted by the MCMC sampler. Ambiguities and uncertainties due to the choice of model parameters on the other hand can be quantified and may indicate problems with the choice of endmembers and/or assumptions made about the data. In some cases where plagioclase and pargasite were confused, two local minima far apart from each other in parameter space were not found by the sampler, even though the mean σ of the confidence plots is very similar. Generally, the fresh laboratory spectra could easily be distinguished as such. The abundance of smFe is predicted to be effectively absent.
The final experiment of unmixing the LSCC spectra with RELAB mineral endmembers demonstrated that this approach can be applied to the lunar surface. When the abundance of agglutinitic glasses is high, olivine and plagioclase abundances tend to be poorly defined. The reconstruction is very accurate, but the uncertainties of the mineral abundances are high for mature spectra. This shows that classical optimization techniques are prone to fail on the task of unmixing the spectra with accurate abundances. The abundances of the smFe particles are, however, well defined. The results align with the findings of Trang and Lucey [54] that the npFe and mpFe abundance for mature spectra depends on the FeO abundance. Furthermore, the predicted abundance of npFe and mpFe particles correlates well with the Is parameter of the LSCC samples. For npFe, this can be interpreted to be a sign that the algorithm predicts plausible results. The npFe and mpFe abundances are also strongly correlated, such that the correlation of mpFe might be attributed to this nearly constant factor. According to Lucey and Riner [16] mpFe sized particles are not present in the samples investigated, but the larger particle size is necessary to model the effects of space weathering. Wohlfarth et al. [24] further suggested that clusters or layers of npFe particles that effectively act as a larger mpFe particle are responsible for the darkening. This effect can also be simulated with light scattering theory as in Arnaut et al. [58]. This would suggest that npFe and mpFe (or clusters of npFe) are both created by the same processes. The cluster/layer theory is consistent with the results of the LSCC unmixing.
All endmembers and investigated mixtures have the same average grain size, which usually is not the case for remotely measured datasets. For the Martian surface, Lapotre et al. [22] and Lapotre et al. [59] investigated the influence of the grain size on the uncertainties and showed that the solution to this inverse problem is highly non-unique. Integrating the grain size into the existing framework would add another layer of complexity to the model, making it even more challenging to obtain meaningful results. Therefore, prior knowledge about the grain size would likely be necessary to constrain the procedure. If the grain size varied significantly between endmembers, systematical errors with respect to their relative fraction would be introduced. Some simplifications to the model, however, have to be made for this highly complex inverse problem to be solvable. The effects of the grain size in combination with space weathering components and the other agents that darken the overall spectrum, discussed in this work, can be investigated in possible future research.
The laboratory samples examined in this work all have a relatively high SNR in comparison to, for example, the Moon Mineralogy Mapper data. When the SNR decreases, the standard deviation of the likelihood function is likely to increase, such that more solutions might be considered acceptable. This, in turn, most likely also increases the uncertainties of the fractional abundances. Nonetheless, we can indirectly measure the influence of the SNR based on the average differences between model prediction and observed spectra. In our work, which follows a Bayesian approach, the variance of the differences is a part of the model and is thus also estimated (cf. Equation (13)). However, averaging over neighboring pixels, where the composition is similar, improves the SNR. Due to the high computational effort to sample from these complex distributions, the number of pixels that can be investigated is limited. Therefore, clustering the data and applying the unmixing to the centroids also mitigates the effects of a low SNR and makes computing the abundances over large areas feasible (e.g., [60,61]).

5. Conclusions

Unmixing can be used to obtain quantitative information about the composition of remote surfaces. Building on the Hapke model (e.g., [9]) and the space weathering model of Wohlfarth et al. [24], we have developed a combined unmixing model. The parameters of the model and their corresponding uncertainties were obtained by Bayesian inference with an MCMC sampler. The results showed that, especially on the Moon, there is no unique solution that would perfectly describe the measured spectra. The absorption bands at around 1- μ m and 2- μ m are the main features for estimating the abundances. These features are obscured by submicroscopic iron particles, making it harder to quantify the corresponding mineral abundances. Some mineral combinations even become gradually interchangeable, only slightly changing the modeled spectrum. Therefore, additional knowledge is necessary to constrain the procedure. We find that information about the elemental abundances, which can be obtained from different independent data sources, is well suited to helping to distinguish between different minerals and we prefer solutions that correspond better with this prior knowledge. Inherent simplifications of the mixing model can lead to systematic errors. In such cases, the uncertainties might still be small, therefore, it is difficult to predict such problems with the model itself if no ground truth is available. Overall, the uncertainties can be accurately estimated from the MCMC sampler. For very mature spectra, these uncertainties become large, suggesting that a classical optimization approach is destined to fail at consistently predicting accurate results. The abundance of submicroscopic iron particles, responsible for the optical effects of maturity, could reliably be determined with the proposed approach. The results of the unmixing on LSCC spectra showed that the probabilistic unmixing approach is applicable to realistic lunar spectra. The knowledge about the uncertainties inherent to the problem are, however, essential for the interpretation of the results obtained for lunar spectra.

Author Contributions

Conceptualization, M.H., T.W. and C.W.; methodology, M.H., T.W. and K.W.; software, M.H.; validation, M.H.; formal analysis, M.H.; investigation, M.H.; data curation, M.H.; writing—original draft preparation, M.H.; writing—review and editing, T.W., C.W. and K.W.; visualization, M.H.; supervision, C.W. All authors have read and agreed to the published version of the manuscript.

Funding

Thorsten Wilhelm and Christian Wöhler acknowledge funding from the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) project number 269661170.

Data Availability Statement

RELAB database (https://pds-geosciences.wustl.edu/spectrallibrary/default.htm, accessed on 1 November 2021). LSCC catalog of lunar returned samples (https://pgi.utk.edu/lunar-soil-characterization-consortium-lscc-data/, accessed on 12 October 2021). Laboratory mineral and mixture spectra from Rommel et al. [13].

Acknowledgments

We want to thank the RELAB and LSCC teams for making the data to the catalog of spectra available publicly. We also want to thank Daniela Rommel for all the hard work, measuring the spectra and preparing the mixtures in Section 3.2. And finally the authors wish to thank the Department of Geochemistry of the University of Göttingen for providing the mineral samples examined in Section 3.2.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Appendix A.1. Synthetic Mixtures

Table A1. Elemental abundances of endmember catalog for experiments in Section 3.1 and Section 3.3. The sample spectra and elemental abundances are taken from the RELAB library (https://pds-geosciences.wustl.edu/spectrallibrary/default.htm, accessed on 1 November 2020).
Table A1. Elemental abundances of endmember catalog for experiments in Section 3.1 and Section 3.3. The sample spectra and elemental abundances are taken from the RELAB library (https://pds-geosciences.wustl.edu/spectrallibrary/default.htm, accessed on 1 November 2020).
SiO2TiO2Al2O3Cr2O3MgOCaOMnOFeONa2OK2OP2O5Fe2O3
AG-TJM-00954.090.161.230.7526.791.520.4915.220.050.050.00.00
AG-TJM-01050.730.748.730.0016.6515.820.135.371.270.000.01.08
PO-EAC-05640.420.000.030.1348.250.190.1511.110.000.000.00.00
PL-EAC-02954.850.0627.710.020.0010.970.010.005.150.400.00.46
PP-ALS-10548.411.055.290.0312.3722.150.266.150.360.000.03.79
SC-EAC-0340.2047.610.010.020.010.000.0245.430.000.000.00.00
Table A2. Origin of endmember samples in Section 3.1 and Section 3.3. All samples except for the AG-TJM-009 Hypersthene endmember, which is of meteoritic origin, originate from Earth.
Table A2. Origin of endmember samples in Section 3.1 and Section 3.3. All samples except for the AG-TJM-009 Hypersthene endmember, which is of meteoritic origin, originate from Earth.
Sample IDMineralOriginAbbreviation
AG-TJM-009HyperstheneJohnstown MeteoriteOPX
AG-TJM-010AugiteKakanui, New ZealandCPXB
PO-EAC-056OlivineGreen sand beach near S Point, HIOLV
PL-EAC-029PlagioclaseYlijarvi, Ylamaaa, Kimi, FinlandPLG
PP-ALS-105DiopsideTrail Creek, Grand Co., COCPXA
SC-EAC-034IlmeniteTelemark, NorwayILM
Table A3. Detailed results for the unmixing of the synthetically generated by using the forward mixing model. Only a uniform prior on the elemental and mineral abundances was set.
Table A3. Detailed results for the unmixing of the synthetically generated by using the forward mixing model. Only a uniform prior on the elemental and mineral abundances was set.
mix0mix1mix2
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX9.979.960.0410.009.809.820.1410.0010.7610.660.2710.00
CPXB39.9540.000.1740.0040.1940.330.8940.0037.3137.751.5840.00
OLV4.584.440.365.008.298.351.845.0018.2819.203.195.00
PLG29.8529.930.4130.0029.1128.451.7130.0013.8916.753.0530.00
CPXA15.0215.380.5015.0012.6012.432.1315.0018.9613.993.5615.00
ILM0.170.830.630.000.521.501.210.002.033.602.210.00
npFe0.000.000.000.000.000.000.000.000.000.010.010.00
mpFe0.000.000.000.000.400.400.000.400.780.790.010.80
SUM99.55100.54 100.00100.51100.87 100.00101.23101.96 100.00
mix3mix4mix5
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX10.079.910.3410.0010.3410.460.1910.009.879.870.2410.00
CPXB35.5039.321.6940.0038.5639.691.3040.0038.7238.861.6040.00
OLV40.6432.955.255.006.372.482.055.007.988.404.565.00
PLG1.6012.864.5330.0025.9628.151.9730.0027.4927.862.8330.00
CPXA7.462.732.1715.0018.3618.762.9215.0013.5613.215.0315.00
ILM9.804.192.490.000.241.621.400.004.304.612.420.00
npFe0.070.050.010.000.570.570.010.600.600.610.020.60
mpFe1.181.160.021.200.010.010.000.000.400.400.010.40
SUM105.08101.95 100.0099.82101.15 100.00101.93102.81 100.00
mix6mix7mix8
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX8.698.630.3310.0010.5210.680.5810.0011.4311.640.3910.00
CPXB46.1545.511.5940.0033.2139.952.9140.0033.8131.752.4440.00
OLV18.7414.985.365.0010.3711.967.175.0011.082.412.635.00
PLG26.5626.744.5130.0021.1029.305.3130.0019.8518.352.7630.00
CPXA0.512.262.0715.0023.737.815.9715.0022.6335.255.7815.00
ILM1.764.322.500.001.631.901.550.000.691.641.380.00
npFe0.680.690.020.600.550.570.030.601.341.320.021.40
mpFe0.750.760.020.801.281.290.031.200.010.010.000.00
SUM102.41102.45 100.00100.56101.60 100.0099.48101.04 100.00
mix9mix10mix11
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX10.6810.460.4910.009.719.440.3710.009.778.860.5410.00
CPXB39.8036.922.9340.0043.5944.301.6040.0039.6446.412.7640.00
OLV7.629.756.685.005.395.743.675.008.749.255.805.00
PLG22.4622.453.8130.0034.5936.192.8930.0020.7827.354.1130.00
CPXA4.6219.448.0215.007.013.592.9815.0020.167.795.4715.00
ILM0.491.741.460.000.132.071.590.001.781.541.310.00
npFe1.401.400.041.401.441.460.021.401.381.440.031.40
mpFe0.370.370.020.400.810.790.020.801.181.220.031.20
SUM85.67100.77 100.00100.41101.33 100.00100.88101.19 100.00
mix12mix13mix14
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX10.3910.090.5010.0011.5411.440.4810.006.486.710.6010.00
CPXB38.6238.962.8240.0037.7336.572.2440.0037.7239.381.6440.00
OLV11.319.795.695.0033.0835.296.595.001.922.241.725.00
PLG19.8819.743.1730.006.577.233.5430.0045.9447.812.0430.00
CPXA23.4420.657.9515.0012.517.605.6715.006.132.571.8915.00
ILM1.691.811.490.001.561.931.560.002.101.771.430.00
npFe2.152.150.042.202.222.230.032.202.252.290.042.20
mpFe0.000.010.010.000.320.320.020.400.730.740.020.80
SUM105.32101.04 100.00102.98100.06 100.00100.29100.48 100.00
mix15mix16mix17
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX10.2911.090.7410.009.179.160.079.209.079.070.189.20
CPXB38.1635.753.9940.0036.8936.790.3636.8038.8738.891.1336.80
OLV24.3811.318.605.004.925.390.994.600.715.902.804.60
PLG18.2318.964.3630.0026.7326.321.1827.6027.0828.032.4127.60
CPXA0.5921.9910.5615.0013.7113.661.5113.8017.809.523.2413.80
ILM1.961.531.310.009.989.871.988.0015.519.712.568.00
npFe2.172.170.062.200.000.000.000.000.010.010.010.00
mpFe1.181.180.031.200.000.000.000.000.410.410.010.40
SUM93.61100.63 100.00101.40101.20 100.00109.03101.12 100.00
mix18mix19mix20
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX9.049.100.309.2010.3710.390.339.208.698.680.209.20
CPXB37.7640.481.7536.8036.8937.341.6436.8035.9535.971.2436.80
OLV7.868.053.454.6012.2214.883.524.602.262.771.964.60
PLG21.7326.483.4927.6021.0525.063.4027.6030.0529.811.8527.60
CPXA16.189.283.7713.809.784.852.9713.8015.5914.843.5113.80
ILM5.824.082.388.007.956.382.678.008.067.692.658.00
npFe0.000.010.010.000.000.010.010.000.600.600.010.60
mpFe0.780.790.010.801.191.210.021.200.000.000.000.00
SUM98.3997.46 100.0098.2598.89 100.00100.6099.77 100.00
mix21mix22mix23
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX9.759.520.309.2010.5910.380.419.209.868.920.519.20
CPXB33.8736.211.7536.8032.0133.402.6936.8034.5940.232.3336.80
OLV7.8415.356.164.6012.453.793.244.6013.349.925.964.60
PLG24.2723.854.1027.6018.8823.853.9927.6019.3328.955.1727.60
CPXA13.797.675.1113.8018.5121.586.1113.8013.034.223.5913.80
ILM7.955.712.728.005.915.002.678.0012.187.873.068.00
npFe0.580.610.020.600.580.570.020.600.640.650.020.60
mpFe0.400.390.010.400.770.790.020.801.171.220.031.20
SUM97.4798.31 100.0098.3597.99 100.00102.33100.11 100.00
mix24mix25mix26
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX8.768.840.349.207.977.870.459.208.518.310.579.20
CPXB32.6336.471.9136.8028.4333.902.6536.8036.4538.132.5236.80
OLV6.8710.074.204.601.577.635.544.6015.526.575.014.60
PLG23.4828.492.9527.6022.7927.964.1727.6033.0334.104.2027.60
CPXA19.238.475.4313.8031.6315.137.0013.800.476.235.0113.80
ILM9.696.862.818.002.914.792.658.003.555.382.758.00
npFe1.411.410.021.401.421.410.031.401.421.430.031.40
mpFe0.000.010.010.000.330.330.020.400.830.850.020.80
SUM100.6799.20 100.0095.3097.28 100.0097.5498.72 100.00
mix27mix28mix29
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX10.5810.400.739.209.859.220.519.208.967.920.649.20
CPXB34.7435.523.7736.8033.6336.952.2836.8031.5837.272.9336.80
OLV17.5418.779.944.600.001.881.804.608.914.143.354.60
PLG21.8718.426.8827.6033.2438.103.1627.6029.9835.604.1727.60
CPXA9.7810.647.8713.806.557.374.3513.8013.888.375.7313.80
ILM3.614.752.768.006.403.352.268.004.025.152.808.00
npFe1.351.340.041.402.172.190.032.202.122.210.042.20
mpFe1.321.320.051.200.000.010.000.000.460.440.020.40
SUM98.1198.50 100.0089.6896.88 100.0097.3398.43 100.00
mix30mix31mix32
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX7.817.530.629.209.098.670.779.208.398.420.118.40
CPXB36.2139.002.6536.8035.8937.862.6636.8034.7434.780.5333.60
OLV2.693.312.814.6011.934.453.484.603.723.531.274.20
PLG31.2335.794.0327.6027.8435.143.8727.6027.1227.421.4725.20
CPXA15.277.145.0813.808.125.794.2413.8010.1310.021.7712.60
ILM5.714.992.628.002.775.022.728.0015.8615.642.2816.00
npFe2.232.250.042.202.282.270.052.200.000.000.000.00
mpFe0.720.760.030.800.991.030.041.200.000.000.000.00
SUM98.9297.76 100.0095.6396.94 100.0099.9699.81 100.00
mix33mix34mix35
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX8.888.810.208.408.658.310.418.409.128.350.548.40
CPXB31.1330.991.1433.6028.7333.863.0533.6029.9734.013.1133.60
OLV7.957.582.044.2010.495.003.824.201.732.411.874.20
PLG21.1020.712.6825.2013.9522.965.3425.2028.4831.943.1325.20
CPXA15.7315.242.6312.6023.2015.106.0912.6016.198.715.2212.60
ILM17.8017.592.7216.0011.5312.022.8416.0012.2711.802.4716.00
npFe0.000.000.000.000.010.020.010.000.010.020.010.00
mpFe0.390.390.010.400.710.750.020.801.171.180.011.20
SUM102.59100.92 100.0096.5597.25 100.0097.7697.23 100.00
mix36mix37mix38
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX8.928.830.238.408.298.180.338.408.878.700.358.40
CPXB34.3832.561.2633.6033.2532.422.0633.6029.2632.761.8333.60
OLV0.071.050.954.205.096.784.664.2015.6517.316.274.20
PLG27.4828.192.0725.2024.4224.113.8325.2014.2621.095.2625.20
CPXA10.3814.402.7812.600.9112.185.6712.6013.643.743.2912.60
ILM13.9713.482.7416.0016.3316.272.9516.0021.7416.693.1016.00
npFe0.590.590.010.600.630.620.020.600.600.610.020.60
mpFe0.000.000.000.000.370.370.010.400.800.820.020.80
SUM95.2198.51 100.0088.2999.94 100.00103.42100.30 100.00
mix39mix40mix41
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX8.018.030.498.409.699.680.458.4012.2811.470.608.40
CPXB31.4831.332.6233.6028.6129.512.6633.6021.4627.393.0733.60
OLV6.5915.457.724.208.433.162.584.2021.6525.189.204.20
PLG21.6321.886.0325.2020.9022.753.7625.203.738.455.3825.20
CPXA7.936.935.5912.6017.6620.386.3112.6025.1512.358.6212.60
ILM17.6917.023.2216.0012.5412.153.0316.0014.2912.783.0816.00
npFe0.630.630.030.601.351.350.031.401.321.350.041.40
mpFe1.221.220.031.200.000.010.010.000.400.400.030.40
SUM93.34100.64 100.0097.8397.63 100.0098.5797.63 100.00
mix42mix43mix44
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX7.787.450.508.407.166.170.698.409.689.220.608.40
CPXB37.1738.831.6133.6027.7332.312.3533.6028.1231.304.2533.60
OLV1.452.772.214.201.734.393.394.200.771.321.274.20
PLG29.7532.752.8625.2035.0636.523.6025.2019.3222.585.3125.20
CPXA0.542.842.5012.6011.963.983.2512.6027.0920.938.7112.60
ILM14.5714.323.0116.0013.8914.963.0016.0011.8111.682.4016.00
npFe1.471.460.031.401.451.520.041.402.062.100.062.20
mpFe0.770.780.020.801.051.040.031.200.000.010.010.00
SUM91.2798.95 100.0097.5398.33 100.0096.7997.04 100.00
mix45mix46mix47
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX8.278.680.568.406.657.770.688.408.118.370.838.40
CPXB29.7330.842.7433.6034.5534.152.3933.6033.1433.913.5333.60
OLV10.2616.037.454.2012.514.384.144.2017.419.246.994.20
PLG21.5221.325.2625.2024.3434.423.7525.2022.0826.875.7125.20
CPXA2.558.106.4412.605.014.513.9012.600.136.115.1912.60
ILM14.0213.423.2416.0016.3912.243.0516.0015.1912.323.2816.00
npFe2.212.190.042.202.182.230.052.202.192.220.062.20
mpFe0.440.430.030.400.800.790.030.801.081.110.051.20
SUM86.3598.39 100.0099.4497.47 100.0096.0796.83 100.00
Table A4. Detailed results for the unmixing of the synthetically generated by using the forward mixing model. A prior on the TiO2 and the Al2O3 was set with α e l e m = 40 . All results are displayed in wt.%.
Table A4. Detailed results for the unmixing of the synthetically generated by using the forward mixing model. A prior on the TiO2 and the Al2O3 was set with α e l e m = 40 . All results are displayed in wt.%.
mix0mix1mix2
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX9.989.960.0410.009.829.820.1410.0010.4610.590.2610.00
CPXB39.9940.040.1640.0040.7240.510.8840.0037.6539.101.5040.00
OLV4.624.650.285.008.808.811.785.0018.9318.692.895.00
PLG30.2930.310.2130.0029.4129.281.5030.0019.1020.832.4830.00
CPXA14.9714.980.3215.0012.5411.431.9915.0016.0810.463.1915.00
ILM0.010.170.130.000.050.190.140.000.150.230.160.00
npFe0.000.000.000.000.000.000.000.000.000.010.010.00
mpFe0.000.000.000.000.400.400.000.400.800.800.010.80
SUM99.86100.12 100.00101.33100.04 100.00102.3699.90 100.00
mix3mix4mix5
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX9.709.810.3410.0010.2210.450.2010.009.879.870.2310.00
CPXB40.0840.801.5340.0038.5540.261.3140.0040.1539.751.4640.00
OLV19.9024.584.535.004.492.241.605.009.8310.864.625.00
PLG21.4021.223.3830.0030.6629.871.6530.0030.7730.552.4330.00
CPXA8.553.122.3415.0020.6517.283.1215.008.748.754.4315.00
ILM0.270.310.180.000.020.150.130.000.010.250.170.00
npFe0.020.030.010.000.570.580.010.600.620.620.020.60
mpFe1.191.190.021.200.020.020.000.000.400.400.010.40
SUM99.9099.84 100.00104.59100.24 100.0099.36100.03 100.00
mix6mix7mix8
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX8.748.710.3210.0010.5710.600.5210.0011.2611.440.3910.00
CPXB46.5345.721.5140.0039.6840.782.3940.0033.6433.522.2940.00
OLV12.4412.664.605.008.1611.205.915.006.232.281.855.00
PLG30.3030.603.5130.0031.5531.714.1430.0020.4121.462.5230.00
CPXA1.322.031.8015.0015.526.074.4115.0028.4231.324.8615.00
ILM0.270.280.180.000.010.260.180.000.000.110.100.00
npFe0.680.680.010.600.550.570.020.601.341.330.021.40
mpFe0.770.770.020.801.311.300.031.200.000.010.010.00
SUM99.6099.99 100.00105.49100.63 100.0099.96100.14 100.00
mix9mix10mix11
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX9.9210.300.4510.009.959.510.3610.008.778.770.4710.00
CPXB39.8938.582.5440.0042.8643.971.5640.0047.1847.212.1940.00
OLV19.9910.136.045.000.476.673.905.0018.309.475.325.00
PLG23.5225.363.1630.0034.6636.592.8130.0025.2029.053.2430.00
CPXA5.8615.296.9915.001.893.142.5615.000.655.724.5815.00
ILM0.050.200.170.000.390.270.170.000.040.190.150.00
npFe1.451.410.031.401.451.450.021.401.441.440.031.40
mpFe0.350.370.020.400.790.790.020.801.231.220.021.20
SUM99.2399.85 100.0090.23100.15 100.00100.14100.40 100.00
mix12mix13mix14
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX9.899.860.4710.0011.5911.360.4810.008.006.910.6710.00
CPXB41.2240.942.4940.0036.2738.912.1940.0035.5238.951.8040.00
OLV16.0111.355.425.0021.6124.655.875.001.892.732.155.00
PLG22.9222.602.7530.0012.4514.572.8130.0043.6747.632.2730.00
CPXA14.7715.136.9615.0017.259.785.8015.009.643.002.4515.00
ILM0.040.170.140.000.230.240.160.000.320.310.180.00
npFe2.152.170.032.202.212.230.032.202.242.280.042.20
mpFe0.000.010.010.000.340.340.020.400.710.730.030.80
SUM104.85100.04 100.0099.4099.50 100.0099.0499.54 100.00
mix15mix16mix17
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX10.7410.710.6810.009.169.160.079.209.109.070.189.20
CPXB38.6438.873.3540.0036.7636.820.3436.8038.8038.931.1436.80
OLV30.8411.697.305.005.545.550.894.605.976.022.664.60
PLG18.0623.213.6230.0026.5626.630.8427.6027.8928.342.0827.60
CPXA1.2415.308.4015.0013.2213.381.2513.809.589.243.0113.80
ILM0.040.180.140.009.009.321.378.009.839.181.518.00
npFe2.182.190.052.200.000.000.000.000.000.010.010.00
mpFe1.201.190.031.200.000.000.000.000.410.410.010.40
SUM99.5699.96 100.00100.25100.85 100.00101.17100.79 100.00
mix18mix19mix20
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX8.979.070.309.2010.3410.310.349.208.588.680.209.20
CPXB39.5440.161.7136.8037.2937.601.7136.8037.7235.961.2236.80
OLV7.047.043.254.6013.7813.893.354.606.292.641.864.60
PLG24.8024.462.9727.6020.7824.973.0127.6030.2729.361.2927.60
CPXA3.2211.423.6713.809.015.162.9113.809.2815.163.3413.80
ILM8.057.441.338.007.488.091.418.007.318.471.458.00
npFe0.000.010.010.000.000.010.010.000.620.600.010.60
mpFe0.790.790.010.801.211.220.021.200.000.000.000.00
SUM91.6399.59 100.0098.69100.02 100.0099.45100.27 100.00
mix21mix22mix23
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX9.459.460.319.2010.3410.320.409.209.048.910.499.20
CPXB36.6436.451.8236.8031.9633.402.4436.8036.5340.282.1636.80
OLV13.4013.255.714.607.172.622.024.605.849.985.344.60
PLG24.6023.823.1827.6022.5222.892.9927.6027.5528.403.9927.60
CPXA4.868.895.4613.8020.3022.875.1613.8012.384.343.5813.80
ILM9.017.871.398.007.277.481.378.0010.208.751.528.00
npFe0.590.610.020.600.550.560.020.600.640.650.020.60
mpFe0.400.400.010.400.790.790.010.801.201.220.031.20
SUM97.9699.74 100.0099.5699.58 100.00101.53100.66 100.00
mix24mix25mix26
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX8.928.780.339.207.787.850.449.208.308.310.579.20
CPXB35.2436.751.7736.8036.2833.652.5736.8035.1437.552.7536.80
OLV10.0310.363.894.6017.547.065.104.6019.537.524.994.60
PLG26.9927.792.3727.6020.2326.113.2227.6027.9031.063.8827.60
CPXA14.378.104.9113.8014.9516.776.8913.800.887.745.6013.80
ILM8.448.381.478.006.997.631.408.007.188.241.488.00
npFe1.421.410.021.401.431.410.031.401.431.420.031.40
mpFe0.000.010.010.000.330.330.010.400.820.840.020.80
SUM103.98100.15 100.00103.7799.07 100.0098.92100.42 100.00
mix27mix28mix29
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX10.1710.160.739.209.929.360.559.207.627.940.659.20
CPXB36.8037.103.4636.8033.2735.562.7036.8034.5836.303.0736.80
OLV7.1812.136.964.600.031.761.664.600.855.143.964.60
PLG22.4621.374.1327.6030.2933.733.3027.6027.3032.253.6827.60
CPXA26.8512.288.0013.8018.3711.525.4913.8012.1710.436.2413.80
ILM7.667.471.388.006.797.431.398.008.748.131.478.00
npFe1.351.330.041.402.122.170.042.202.192.210.042.20
mpFe1.351.340.041.200.000.000.000.000.440.430.020.40
SUM111.12100.50 100.0098.6699.37 100.0091.25100.19 100.00
mix30mix31mix32
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX8.227.500.649.208.788.730.819.208.438.420.118.40
CPXB34.7238.722.6336.8031.5536.893.1536.8034.8834.750.5333.60
OLV2.553.933.034.607.315.654.134.603.153.341.194.20
PLG29.2433.343.3227.6023.9431.733.8327.6027.1627.071.3125.20
CPXA17.008.034.9613.8017.807.615.3113.8011.2310.361.6312.60
ILM10.288.181.478.009.508.081.468.0016.2916.291.8716.00
npFe2.212.250.042.202.182.260.052.200.000.000.000.00
mpFe0.740.750.030.801.011.020.041.200.000.000.000.00
SUM102.0299.70 100.0098.8898.68 100.00101.14100.23 100.00
mix33mix34mix35
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX8.808.780.208.408.588.240.398.408.208.130.398.40
CPXB31.3231.251.0833.6031.6834.372.4033.6031.7934.732.0933.60
OLV7.487.331.934.203.793.652.874.207.683.332.644.20
PLG20.9021.422.3325.2018.8922.733.3525.2025.2229.443.0725.20
CPXA15.0214.692.3012.6020.5115.345.5012.6012.488.714.4012.60
ILM17.2617.272.1016.0016.8414.592.0116.0010.6514.552.2316.00
npFe0.000.000.000.000.010.020.010.000.030.030.010.00
mpFe0.390.390.010.400.750.760.010.801.131.170.021.20
SUM100.78100.73 100.00100.2998.93 100.0096.0298.89 100.00
mix36mix37mix38
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX8.748.820.228.407.988.170.328.408.648.640.358.40
CPXB32.5832.551.2133.6032.0832.492.0233.6033.8933.401.6033.60
OLV0.160.980.874.2014.845.944.034.2012.9314.694.964.20
PLG29.4027.141.7025.2024.1424.353.0425.2023.2023.263.6325.20
CPXA8.8514.892.5912.6017.1412.605.4812.600.883.773.1912.60
ILM14.5415.232.0216.0017.7116.682.1116.0016.6116.572.1616.00
npFe0.590.590.010.600.620.620.020.600.610.600.020.60
mpFe0.000.000.000.000.380.370.010.400.830.830.020.80
SUM94.2899.62 100.00113.88100.23 100.0096.15100.34 100.00
mix39mix40mix41
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX8.948.130.598.409.439.610.448.4011.0111.400.608.40
CPXB28.4731.692.4533.6029.5029.792.5033.6029.7929.473.1433.60
OLV19.0112.955.564.201.853.032.554.206.8511.396.904.20
PLG22.7924.613.5725.2022.3021.883.0625.2015.6216.223.6625.20
CPXA5.216.644.3312.6021.8820.385.8912.6014.3016.618.8212.60
ILM13.9916.282.2716.0013.9814.332.0716.0014.1113.702.0316.00
npFe0.610.620.020.601.341.360.031.401.341.340.041.40
mpFe1.211.230.031.200.000.010.010.000.440.450.020.40
SUM98.42100.30 100.0098.9599.02 100.0091.6998.78 100.00
mix42mix43mix44
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX7.947.420.528.407.276.400.778.408.818.720.508.40
CPXB37.4238.801.7033.6027.3531.093.0033.6035.1434.992.7133.60
OLV4.273.422.584.202.925.293.464.200.481.851.714.20
PLG27.4431.162.6125.2030.1133.353.4825.2018.7524.893.4225.20
CPXA2.323.052.4012.6014.106.284.8712.6016.4314.235.7412.60
ILM16.4616.302.1916.0018.4617.292.0316.0013.8113.842.0316.00
npFe1.451.470.031.401.481.510.041.402.162.160.042.20
mpFe0.770.770.020.801.041.040.031.200.000.010.010.00
SUM95.86100.15 100.00100.2199.71 100.0093.4198.53 100.00
mix45mix46mix47
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX8.838.590.548.407.407.720.698.408.438.260.818.40
CPXB30.1631.742.3633.6030.1233.812.6033.6030.2833.863.1233.60
OLV7.1013.916.054.200.755.233.934.204.599.465.674.20
PLG21.8622.243.5925.2028.0931.503.5725.2021.3825.014.0825.20
CPXA16.908.085.9212.608.945.684.3012.6019.246.885.6212.60
ILM14.0914.872.0816.0015.9015.362.1316.0012.7214.822.1116.00
npFe2.132.200.042.202.272.230.052.202.212.210.062.20
mpFe0.460.440.020.400.780.780.030.801.051.100.041.20
SUM98.9599.42 100.0091.1999.31 100.0096.6698.28 100.00

Appendix A.2. Fresh Laboratory Samples

Table A5. Elemental abundances of endmember catalog for Section 3.2 taken from Rommel et al. [13].
Table A5. Elemental abundances of endmember catalog for Section 3.2 taken from Rommel et al. [13].
SiO2TiO2Al2O3Cr2O3MgOCaOMnOFeONa2OK2OP2O5Fe2O3
Pargasite43.801.0615.630.0318.8512.390.050.252.841.720.000.28
Augite50.101.264.810.6015.2222.700.082.410.490.020.002.35
Olivine37.700.070.040.0335.400.160.3926.490.030.020.000.23
Plagioclase54.700.0028.800.000.0011.220.000.284.780.600.000.00
Ferrosilite47.900.050.320.035.180.881.9244.700.030.030.000.20
Ilmenite0.0250.400.020.126.150.000.2634.100.000.000.007.90
Table A6. Origins of endmember catalog for Section 3.2 taken from Rommel et al. [13]. All samples were provided by the University of Göttingen and originate from Earth.
Table A6. Origins of endmember catalog for Section 3.2 taken from Rommel et al. [13]. All samples were provided by the University of Göttingen and originate from Earth.
MineralOriginAbbreviation
PargasiteMerelani Hill, TanzaniaPRG
AugitePaškopole, Czech RepublicAUG
OlivineChinaOLV
PlagioclaseNordingrå, SwedenPLG
FerrosiliteMansjöberg, SwedenFS
IlmeniteRogaland, NorwayILM
Table A7. Detailed results for the unmixing of fresh laboratory spectra from Rommel et al. [13]. The results were obtained without the use of an elemental prior.
Table A7. Detailed results for the unmixing of fresh laboratory spectra from Rommel et al. [13]. The results were obtained without the use of an elemental prior.
mix0mix1mix2
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
PRG25.3725.580.609.2036.7436.950.5918.5046.1946.240.4728.10
AUG0.020.140.110.000.030.110.090.000.000.100.080.00
OLV0.000.050.050.000.000.050.040.000.000.040.040.00
PLG74.5974.130.7090.8062.9462.830.6781.5053.7353.880.5771.90
FS0.280.280.090.000.150.140.080.000.080.100.060.00
ILM14.4913.911.010.0013.6214.040.980.0011.6011.470.840.00
npFe0.000.000.00 0.000.000.00 0.000.000.00
mpFe0.000.000.00 0.000.000.00 0.000.000.00
SUM114.75114.09 100.00113.48114.12 100.00111.61111.83 100.00
mix3mix4mix5
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
PRG78.2878.190.3268.0083.1182.980.4278.5090.1090.040.3689.10
AUG0.080.140.080.000.120.130.080.000.010.040.040.00
OLV0.000.060.050.000.000.050.050.000.000.030.030.00
PLG22.5222.360.3932.0017.5517.550.4821.5010.6710.670.3810.90
FS0.020.060.040.000.030.070.050.000.030.070.050.00
ILM6.766.740.600.009.289.320.740.007.727.800.640.00
npFe0.000.000.00 0.000.000.00 0.000.000.00
mpFe0.000.000.00 0.000.000.00 0.000.000.00
SUM107.67107.55 100.00110.08110.10 100.00108.53108.65 100.00
mix6mix7mix8
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
PRG18.7618.851.710.0030.0230.113.690.0044.4146.994.150.00
AUG0.000.080.080.000.000.190.190.000.120.300.290.00
OLV0.530.620.340.003.593.800.890.006.106.821.100.00
PLG70.9871.211.8992.8043.1543.114.1985.1018.9614.114.7976.90
FS11.9311.880.247.2025.7425.640.5614.9033.4932.960.6823.10
ILM0.060.610.570.000.040.680.630.000.150.820.790.00
npFe0.000.000.00 0.000.000.00 0.000.000.00
mpFe0.000.000.00 0.000.000.00 0.000.000.00
SUM102.27103.25 100.00102.54103.52 100.00103.24101.99 100.00
mix9mix10mix11
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
PRG33.2738.701.840.0016.8921.221.760.0011.8811.520.820.00
AUG0.220.720.720.004.181.271.100.000.090.570.500.00
OLV5.664.771.110.000.732.211.130.000.040.500.400.00
PLG7.452.091.8058.804.321.581.5938.000.090.720.6526.30
FS57.0557.110.7641.2080.5580.170.8762.0093.3493.390.5273.70
ILM0.570.640.610.000.290.480.460.000.000.430.410.00
npFe0.000.000.00 0.000.010.01 0.000.010.01
mpFe0.000.000.00 0.000.000.00 0.000.000.00
SUM104.23104.03 100.00106.97106.93 100.00105.44107.12 100.00
mix12mix13mix14
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
PRG3.853.731.150.003.897.303.790.0053.5255.491.0956.50
AUG3.602.811.460.0018.793.513.2254.100.440.560.540.00
OLV1.021.050.750.006.446.813.890.000.190.340.320.00
PLG0.070.610.5813.700.693.102.720.000.230.800.740.00
FS99.9799.790.2086.3073.1683.113.1945.9050.3149.680.8643.50
ILM0.010.360.360.000.010.850.800.000.080.510.490.00
npFe0.000.000.00 0.000.030.03 0.000.010.01
mpFe0.000.000.00 0.000.010.01 0.000.000.00
SUM108.52108.35 100.00102.97104.68 100.00104.77107.38 100.00
mix15mix16
ModeMeanStdTheoModeMeanStdTheo
PRG0.131.401.250.0031.8832.0416.290.00
AUG0.321.161.140.002.187.306.500.00
OLV20.7828.542.7553.607.477.106.050.00
PLG0.031.131.050.0012.3215.6713.110.00
FS76.6276.172.6646.402.695.924.5854.60
ILM0.050.460.450.0024.8824.8511.5245.40
npFe0.000.010.01 4.876.821.85
mpFe0.000.010.01 0.000.060.06
SUM97.93108.85 100.0081.4292.88 100.00
Table A8. Detailed results for the unmixing of fresh laboratory spectra from Rommel et al. [13].
Table A8. Detailed results for the unmixing of fresh laboratory spectra from Rommel et al. [13].
mix0mix1mix2
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
PRG21.7521.670.909.2033.1533.630.9518.5044.3044.070.7328.10
AUG0.000.040.040.000.000.040.040.000.000.030.030.00
OLV0.000.020.020.000.000.020.020.000.000.020.020.00
PLG81.8481.820.9190.8069.8369.840.9781.5059.0559.060.7571.90
FS0.000.040.040.000.000.030.030.000.000.020.020.00
ILM0.000.010.010.000.000.030.030.000.020.070.060.00
npFe0.010.010.00 0.010.010.00 0.010.010.00
mpFe0.000.000.00 0.000.000.00 0.000.000.00
SUM103.59103.60 100.00102.98103.60 100.00103.37103.26 100.00
mix3mix4mix5
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
PRG77.4377.620.3668.0082.2982.360.5678.5090.0589.820.4689.10
AUG0.000.020.020.000.000.030.030.000.000.020.020.00
OLV0.000.010.010.000.000.020.020.000.000.010.010.00
PLG24.5924.500.4032.0020.2020.370.6021.5012.2612.400.5210.90
FS0.000.010.010.000.000.020.010.000.000.010.010.00
ILM1.181.130.420.001.161.070.440.001.871.880.590.00
npFe0.000.000.00 0.000.000.00 0.000.000.00
mpFe0.000.000.00 0.000.000.00 0.000.000.00
SUM103.20103.31 100.00103.66103.86 100.00104.19104.15 100.00
mix6mix7mix8
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
PRG0.000.040.040.000.000.070.070.000.030.130.120.00
AUG0.000.020.020.000.000.050.050.000.000.080.080.00
OLV0.000.070.070.000.010.260.230.000.090.420.370.00
PLG90.5590.520.1892.8075.7875.950.3685.1065.6165.870.5376.90
FS13.0012.970.217.2027.5328.270.3914.9037.4137.460.5523.10
ILM0.000.000.000.000.000.000.000.000.000.000.000.00
npFe0.000.000.00 0.000.000.00 0.000.000.00
mpFe0.000.000.00 0.000.000.00 0.000.000.00
SUM103.55103.63 100.00103.32104.60 100.00103.15103.96 100.00
mix9mix10mix11
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
PRG0.000.240.220.000.040.410.360.000.600.650.510.00
AUG0.010.150.140.000.030.240.220.000.070.220.210.00
OLV0.210.570.500.000.170.810.680.000.040.330.300.00
PLG44.4343.970.7458.8020.8624.010.8638.0010.2011.870.6826.30
FS60.4860.540.7741.2079.3681.980.8262.0094.1394.140.5473.70
ILM0.000.000.000.000.000.010.010.000.000.010.010.00
npFe0.000.010.01 0.000.010.01 0.000.010.01
mpFe0.000.000.00 0.000.000.00 0.000.000.00
SUM105.12105.47 100.00100.45107.45 100.00105.04107.22 100.00
mix12mix13mix14
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
PRG0.740.790.620.008.077.413.770.0053.3455.501.1156.50
AUG0.160.480.410.000.683.633.1754.100.430.550.540.00
OLV2.352.480.870.003.136.653.830.000.050.340.330.00
PLG2.674.261.0413.700.672.992.610.002.590.860.800.00
FS99.9999.820.1886.3084.3983.103.0745.9049.7049.720.8543.50
ILM0.000.010.010.000.891.000.230.000.100.130.110.00
npFe0.000.010.00 0.010.030.03 0.000.010.01
mpFe0.000.000.00 0.000.010.01 0.000.000.00
SUM105.91107.82 100.0097.84104.78 100.00106.22107.10 100.00
mix15mix16
ModeMeanStdTheoModeMeanStdTheo
PRG0.190.570.510.0015.1120.7811.020.00
AUG0.090.450.410.000.866.575.860.00
OLV25.4129.632.5653.605.088.536.560.00
PLG0.191.341.200.005.3910.368.500.00
FS76.0776.512.6646.406.436.344.5054.60
ILM0.000.020.010.0042.0641.835.4945.40
npFe0.000.010.01 4.095.651.63
mpFe0.000.010.01 0.000.040.05
SUM101.94108.51 100.0074.9394.42 100.00

Appendix A.3. LSCC Samples

Table A9. Elemental abundances used for the LSCC samples taken from Taylor et al. [18] and Taylor et al. [19]. The elemental and mineral abundances of the 20–45 μ m size fraction were used for the priors, but the 0–45 μ m size fraction spectra were used for the unmixing. The 0–45 μ m size fraction mineral abundances were not available.
Table A9. Elemental abundances used for the LSCC samples taken from Taylor et al. [18] and Taylor et al. [19]. The elemental and mineral abundances of the 20–45 μ m size fraction were used for the priors, but the 0–45 μ m size fraction spectra were used for the unmixing. The 0–45 μ m size fraction mineral abundances were not available.
SiO2TiO2Al2O3Cr2O3MgOCaOMnOFeONa2OK2OP2O5Fe2O3
1008441.38.3012.000.308.4611.600.2115.500.390.120.120.17
1200145.33.2011.000.4110.609.830.2116.900.390.210.230.11
1203046.13.7410.500.409.949.090.2317.600.410.260.190.12
1414147.21.9615.000.2611.0010.100.1511.600.590.470.260.07
1416347.12.0015.400.2311.0010.200.1511.500.570.410.210.08
1425947.11.9915.800.2410.7010.500.1511.000.600.430.260.09
1426047.41.8616.300.2210.4010.700.1410.700.600.440.220.10
1504146.12.0312.500.3911.209.910.2015.200.360.160.190.11
1507145.82.3312.400.4311.409.810.2115.600.360.140.150.10
6114144.50.5826.100.116.5615.200.085.150.460.100.060.05
6122144.50.5627.200.095.4515.900.064.620.460.070.050.04
6223144.50.5825.700.116.5915.300.095.310.420.090.070.08
6480144.60.6326.500.106.0915.600.084.820.440.120.060.10
6746144.40.4427.300.095.1116.100.074.930.410.050.030.07
6748144.70.4926.700.095.9815.600.085.190.450.060.050.04
7018140.78.1111.500.4310.1010.300.2216.000.350.080.060.16
7106139.29.489.330.4810.809.580.2318.500.340.070.040.17
7150138.410.709.940.469.979.940.2417.800.350.070.070.17
7922140.57.3811.600.4010.9010.300.2215.800.380.090.060.17
Table A10. Detailed results for the unmixing of LSCC spectra with the catalog of RELAB endmembers. A prior on the TiO2 and the Al2O3 was set with α e l e m = 40 .
Table A10. Detailed results for the unmixing of LSCC spectra with the catalog of RELAB endmembers. A prior on the TiO2 and the Al2O3 was set with α e l e m = 40 .
100841200112030
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX11.5010.730.470.8023.1422.920.611.9922.8322.880.826.48
CPX40.0935.271.8536.6946.2045.631.7043.7542.3337.702.2350.40
OLV7.1111.704.193.287.022.311.908.249.024.593.557.24
PLG12.9411.942.5039.347.1511.682.0932.4612.3821.343.9925.77
OTH 4.92 7.27 5.73
ILM6.5212.012.1514.986.015.501.036.308.047.881.284.38
npFe1.861.940.04 1.561.510.03 0.350.350.02
mpFe0.690.680.03 0.790.840.02 0.330.340.02
SUM78.1681.66 100.0089.5288.03 100.0094.6094.39 100.00
141411416314259
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX17.0316.990.5113.9311.8711.860.4816.3412.5512.670.4920.03
CPX6.757.711.6322.534.776.361.6824.335.529.771.6229.09
OLV14.1714.455.817.3622.7721.996.146.0314.465.123.796.22
PLG54.0954.115.0848.9647.6048.625.1147.5155.1060.573.2238.16
OTH 3.72 3.77 2.98
ILM4.324.290.733.503.804.190.712.014.084.090.693.52
npFe0.160.160.01 0.580.590.02 1.071.070.02
mpFe0.110.110.01 0.200.200.02 0.350.360.02
SUM96.3697.56 100.0090.8293.03 100.0091.7192.21 100.00
142601504115071
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX14.8315.010.5113.9823.3123.260.668.1323.3423.120.656.66
CPX9.2712.131.8226.8834.1640.742.6140.3536.2639.392.1238.97
OLV1.674.163.156.2713.244.733.347.113.414.423.268.06
PLG56.0957.353.0446.5913.3816.892.8433.4116.2820.432.7037.42
OTH 3.58 8.41 4.96
ILM4.403.650.652.693.183.190.642.593.434.050.783.93
npFe1.111.130.02 1.601.590.04 1.331.340.03
mpFe0.400.400.02 0.800.830.03 0.710.720.02
SUM86.2692.30 100.0087.2788.81 100.0082.7391.42 100.00
611416122162231
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX4.634.540.283.394.464.460.174.184.254.300.284.23
CPX0.541.810.545.390.381.180.336.220.791.590.466.19
OLV6.464.212.374.2418.1117.152.005.5010.416.202.835.90
PLG80.9884.452.3085.7775.0974.801.8282.8378.5581.982.6182.45
OTH 0.61 0.42 0.61
ILM1.171.110.200.611.081.140.200.851.181.100.180.61
npFe0.430.420.01 0.060.060.00 0.490.470.01
mpFe0.150.150.01 0.030.030.00 0.140.150.01
SUM93.7796.11 100.0099.1298.73 100.0095.1895.17 100.00
648016746167481
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX4.264.200.264.405.865.900.183.975.035.060.144.09
CPX0.101.830.525.403.594.410.525.800.460.970.295.06
OLV7.213.142.043.905.817.892.233.3611.7012.801.805.55
PLG80.7885.122.0985.2181.3481.081.9786.3379.3779.241.6384.88
OTH 0.43 0.13 0.28
ILM1.151.200.210.650.790.810.160.400.970.990.180.14
npFe0.440.430.01 0.120.110.00 0.150.150.00
mpFe0.120.120.01 0.030.030.00 0.040.040.00
SUM93.5095.49 100.0097.40100.10 100.0097.5299.05 100.00
701817106171501
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX13.6613.640.373.250.010.380.342.3312.4012.800.482.62
CPX27.3025.511.6130.630.012.160.4939.1631.6231.591.2835.99
OLV10.697.643.897.7517.9711.646.677.900.711.871.546.56
PLG19.9324.472.9436.4054.8961.006.3528.1418.3318.952.2230.05
OTH 2.80 1.42 2.37
ILM18.7019.813.0419.1736.2637.295.8421.0630.7927.044.0422.40
npFe1.771.770.03 1.031.030.04 1.781.770.03
mpFe0.650.660.02 0.120.140.03 0.640.640.02
SUM90.2891.08 100.00109.14112.47 100.0093.8592.24 100.00
79221
ModeMeanStdTheo
OPX10.7610.760.443.48
CPX11.6212.801.8227.70
OLV24.7123.596.3711.35
PLG24.6427.313.0939.97
OTH 0.24
ILM16.4916.342.5917.27
npFe1.821.810.04
mpFe0.600.600.02
SUM88.2290.80 100.00
Table A11. Detailed results for the unmixing of LSCC spectra with the catalog of RELAB endmembers. For all relevant elemental abundances a prior was introduced.
Table A11. Detailed results for the unmixing of LSCC spectra with the catalog of RELAB endmembers. For all relevant elemental abundances a prior was introduced.
100841200112030
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX11.1110.900.470.8023.0322.880.611.9923.4722.880.876.48
CPX44.6142.222.4436.6947.3445.811.5843.7542.6737.882.1450.40
OLV1.322.551.933.281.741.331.138.240.851.801.467.24
PLG12.0513.442.7439.3410.8612.031.6132.4619.7023.192.2725.77
OTH 4.92 7.27 5.73
ILM15.0314.582.2214.985.516.961.376.307.4410.181.844.38
npFe1.931.930.05 1.511.510.03 0.330.350.02
mpFe0.680.710.02 0.860.840.02 0.350.350.02
SUM84.1383.69 100.0088.4889.02 100.0094.1395.95 100.00
141411416314259
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX17.0016.940.4813.9311.4411.760.4616.3413.0112.750.5020.03
CPX14.268.501.4322.5311.138.421.5824.3318.669.381.6929.09
OLV7.1512.813.417.3616.9216.653.466.034.449.803.106.22
PLG55.3154.753.1548.9649.3451.883.2347.5151.9555.642.8738.16
OTH 3.72 3.77 2.98
ILM4.685.140.953.504.605.020.932.015.104.990.903.52
npFe0.160.160.01 0.580.590.02 1.061.070.02
mpFe0.120.110.01 0.210.210.01 0.330.340.02
SUM98.4098.13 100.0093.4393.72 100.0093.1692.56 100.00
142601504115071
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX15.8815.140.5513.9823.3623.360.678.1323.2823.190.646.66
CPX19.9412.822.2726.8848.1742.502.9440.3546.7140.362.2838.97
OLV6.427.512.986.275.712.852.247.113.922.952.188.06
PLG45.6852.693.6046.598.0716.753.0233.4117.0320.632.3337.42
OTH 3.58 8.41 4.96
ILM4.294.380.852.693.753.760.822.594.754.920.973.93
npFe1.031.120.04 1.551.580.04 1.341.340.03
mpFe0.390.390.02 0.830.830.03 0.710.720.02
SUM92.2192.54 100.0089.0789.21 100.0095.6992.05 100.00
611416122162231
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX4.804.740.263.394.474.330.174.184.524.500.304.23
CPX0.701.180.355.393.482.110.426.220.011.350.356.19
OLV10.8110.481.764.2413.9914.021.755.5012.5311.421.905.90
PLG77.3978.331.7785.7775.0577.151.5182.8375.8876.571.9382.45
OTH 0.61 0.42 0.61
ILM1.301.320.260.611.131.260.260.851.171.280.250.61
npFe0.440.440.01 0.060.060.00 0.480.480.01
mpFe0.130.130.01 0.030.030.00 0.140.140.01
SUM95.0096.05 100.0098.1198.88 100.0094.1195.12 100.00
648016746167481
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX4.754.420.264.405.955.930.163.975.065.050.144.09
CPX4.891.420.365.404.264.390.475.800.391.090.295.06
OLV8.189.091.633.909.778.291.543.3612.1912.591.465.55
PLG75.7779.081.7285.2181.2980.611.3786.3379.2579.281.3284.88
OTH 0.43 0.13 0.28
ILM1.091.420.280.650.840.950.210.401.091.130.230.14
npFe0.440.440.01 0.120.120.00 0.150.150.00
mpFe0.110.110.01 0.030.030.00 0.040.040.00
SUM94.6895.45 100.00102.11100.17 100.0097.9899.15 100.00
701817106171501
ModeMeanStdTheoModeMeanStdTheoModeMeanStdTheo
OPX13.6213.710.383.250.130.400.372.3312.4312.770.472.62
CPX29.0728.561.8230.637.992.630.5039.1634.3632.031.2535.99
OLV1.714.282.567.7512.9418.693.067.900.131.461.226.56
PLG19.8324.092.7836.4053.9954.663.5328.1417.3118.642.0630.05
OTH 2.80 1.42 2.37
ILM19.0821.833.1119.1735.5232.994.3921.0628.9427.943.6122.40
npFe1.771.760.03 1.001.030.04 1.741.770.03
mpFe0.670.660.02 0.130.110.02 0.660.640.02
SUM83.3192.47 100.00110.56109.37 100.0093.1792.84 100.00
79221
ModeMeanStdTheo
OPX11.2611.100.403.48
CPX23.1123.641.5627.70
OLV10.7611.373.7511.35
PLG27.5727.742.8839.97
OTH 0.24
ILM19.8719.002.8417.27
npFe1.761.760.03
mpFe0.630.620.02
SUM92.5692.85 100.00

References

  1. Papike, J.; Taylor, L.; Simon, S. Lunar minerals. In Lunar Sourcebook; a User’s Guide to the Moon; Heiken, G., Vaniman, D., French, B., Eds.; Cambridge University Press: Cambridge, UK, 1991; pp. 121–182. [Google Scholar]
  2. Mustard, J.F.; Pieters, C.M. Quantitative abundance estimates from bidirectional reflectance measurements. J. Geophys. Res.: Solid Earth 1987, 92, E617–E626. [Google Scholar] [CrossRef] [Green Version]
  3. Lucey, P.G.; Blewett, D.T.; Jolliff, B.L. Lunar iron and titanium abundance algorithms based on final processing of Clementine ultraviolet-visible images. J. Geophys. Res. Planets 2000, 105, 20297–20305. [Google Scholar] [CrossRef]
  4. Sato, H.; Robinson, M.S.; Lawrence, S.J.; Denevi, B.W.; Hapke, B.; Jolliff, B.L.; Hiesinger, H. Lunar mare TiO2 abundances estimated from UV/Vis reflectance. Icarus 2017, 296, 216–238. [Google Scholar] [CrossRef]
  5. Bhatt, M.; Wöhler, C.; Grumpe, A.; Hasebe, N.; Naito, M. Global mapping of lunar refractory elements: Multivariate regression vs. machine learning. Astron. Astrophys. 2019, 627, A155. [Google Scholar] [CrossRef]
  6. Keshava, N.; Mustard, J.F. Spectral unmixing. IEEE Signal Process. Mag. 2002, 19, 44–57. [Google Scholar] [CrossRef]
  7. Bishop, J.L.; Bell, J.F., III; Moersch, J.E. Remote Compositional Analysis: Techniques for Understanding Spectroscopy, Mineralogy, and Geochemistry of Planetary Surfaces; Cambridge University Press: Cambridge, UK, 2019; Volume 24. [Google Scholar]
  8. Pieters, C.M.; Boardman, J.; Buratti, B.; Chatterjee, A.; Clark, R.; Glavich, T.; Green, R.; Head, J.; Isaacson, P.; Malaret, E.; et al. The Moon Mineralogy Mapper (M3) on Chandrayaan-1. Curr. Sci. 2009, 96, 500–505. [Google Scholar]
  9. Hapke, B. Theory of Reflectance and Emittance Spectroscopy; Cambridge University Press: Cambridge, UK, 2012. [Google Scholar]
  10. Poulet, F.; Carter, J.; Bishop, J.; Loizeau, D.; Murchie, S. Mineral abundances at the final four curiosity study sites and implications for their formation. Icarus 2014, 231, 65–76. [Google Scholar] [CrossRef]
  11. Liu, Y.; Glotch, T.D.; Scudder, N.A.; Kraner, M.L.; Condus, T.; Arvidson, R.E.; Guinness, E.A.; Wolff, M.J.; Smith, M.D. End-member identification and spectral mixture analysis of CRISM hyperspectral data: A case study on southwest Melas Chasma, Mars. J. Geophys. Res.: Planets 2016, 121, 2004–2036. [Google Scholar] [CrossRef]
  12. Li, L. Quantifying lunar soil composition with partial least squares modeling of reflectance. Adv. Space Res. 2008, 42, 267–274. [Google Scholar] [CrossRef]
  13. Rommel, D.; Grumpe, A.; Felder, M.P.; Wöhler, C.; Mall, U.; Kronz, A. Automatic endmember selection and nonlinear spectral unmixing of Lunar analog minerals. Icarus 2017, 284, 126–149. [Google Scholar] [CrossRef]
  14. Hapke, B. Space weathering from Mercury to the asteroid belt. J. Geophys. Res.: Planets 2001, 106, 10039–10073. [Google Scholar] [CrossRef]
  15. Pieters, C.M.; Noble, S.K. Space weathering on airless bodies. J. Geophys. Res.: Planets 2016, 121, 1865–1884. [Google Scholar] [CrossRef] [PubMed]
  16. Lucey, P.G.; Riner, M.A. The optical effects of small iron particles that darken but do not redden: Evidence of intense space weathering on Mercury. Icarus 2011, 212, 451–462. [Google Scholar] [CrossRef]
  17. Sunshine, J.M.; Pieters, C.M. Determining the composition of olivine from reflectance spectroscopy. J. Geophys. Res. Planets 1998, 103, 13675–13688. [Google Scholar] [CrossRef]
  18. Taylor, L.A.; Pieters, C.; Keller, L.P.; Morris, R.V.; McKAY, D.S.; Patchen, A.; Wentworth, S. The effects of space weathering on Apollo 17 mare soils: Petrographie and chemical characterization. Meteorit. Planet. Sci. 2001, 36, 285–299. [Google Scholar] [CrossRef]
  19. Taylor, L.A.; Pieters, C.; Patchen, A.; Taylor, D.H.S.; Morris, R.V.; Keller, L.P.; McKay, D.S. Mineralogical and chemical characterization of lunar highland soils: Insights into the space weathering of soils on airless bodies. J. Geophys. Res. Planets 2010, 115. [Google Scholar] [CrossRef]
  20. Hess, M.; Wilhelm, T.; Arnaut, M.; Wöhler, C. Uncertainty-Based Unmixing of Space Weathered Lunar Spectra. In Proceedings of the Lunar and Planetary Science Conference, Virtual Meeting. 15–19 March 2021; Number 2548, p. 1048. [Google Scholar]
  21. Pieters, C.M.; Besse, S.; Boardman, J.; Buratti, B.; Cheek, L.; Clark, R.N.; Combe, J.P.; Dhingra, D.; Goswami, J.N.; Green, R.O.; et al. Mg-spinel lithology: A new rock type on the lunar farside. J. Geophys. Res. Planets 2011, 116, E6. [Google Scholar] [CrossRef] [Green Version]
  22. Lapotre, M.G.; Ehlmann, B.L.; Minson, S.E. A probabilistic approach to remote compositional analysis of planetary surfaces. J. Geophys. Res. Planets 2017, 122, 983–1009. [Google Scholar] [CrossRef] [Green Version]
  23. Gelman, A.; Carlin, J.B.; Stern, H.S.; Dunson, D.B.; Vehtari, A.; Rubin, D.B. Bayesian Data Analysis; CRC Press: Boca Raton, FL, USA, 2013. [Google Scholar]
  24. Wohlfarth, K.S.; Wöhler, C.; Grumpe, A. Space Weathering and Lunar OH/H2O Insights from Ab Initio Mie Modeling of Submicroscopic Iron. Astron. J. 2019, 158, 80. [Google Scholar] [CrossRef]
  25. Hapke, B. Bidirectional reflectance spectroscopy: 1. Theory. J. Geophys. Res. Solid Earth 1981, 86, 3039–3054. [Google Scholar] [CrossRef]
  26. Hapke, B. Bidirectional Reflectance Spectroscopy: 5. The Coherent Backscatter Opposition Effect and Anisotropic Scattering. Icarus 2002, 157, 523–534. [Google Scholar] [CrossRef] [Green Version]
  27. Hapke, B. Bidirectional reflectance spectroscopy: 3. Correction for macroscopic roughness. Icarus 1984, 59, 41–59. [Google Scholar] [CrossRef] [Green Version]
  28. Kattawar, G.W. A three-parameter analytic phase function for multiple scattering calculations. J. Quant. Spectrosc. Radiat. Transf. 1975, 15, 839–849. [Google Scholar] [CrossRef] [Green Version]
  29. Warell, J. Properties of the Hermean regolith: IV. Photometric parameters of Mercury and the Moon contrasted with Hapke modelling. Icarus 2004, 167, 271–286. [Google Scholar] [CrossRef]
  30. Yang, Y.; Li, S.; Milliken, R.E.; Zhang, H.; Robertson, K.; Hiroi, T. Phase Functions of Typical Lunar Surface Minerals Derived for the Hapke Model and Implications for Visible to Near-Infrared Spectral Unmixing. J. Geophys. Res. Planets 2019, 124, 31–60. [Google Scholar] [CrossRef]
  31. Howari, F.M.; Acbas, G.; Nazzal, Y.; AlAydaroos, F. Hapke-based computational method to enable unmixing of hyperspectral data of common salts. Chem. Cent. J. 2018, 12, 1–15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Hapke, B. Bidirectional reflectance spectroscopy: 6. Effects of porosity. Icarus 2008, 195, 918–926. [Google Scholar] [CrossRef] [Green Version]
  33. Fowler, B.W. Expansion of Mie-theory phase functions in series of Legendre polynomials. JOSA 1983, 73, 19–22. [Google Scholar] [CrossRef]
  34. Baker, A.; Jolliff, B.; Yasanayake, C.; Denevi, B.; Lawrence, S. Lunar Agglutinate Glass Compositions and Comparison to Soil Grain-Size Fractions. In Proceedings of the Lunar and Planetary Science Conference, Virtual Meeting. 16–20 March 2020; Number 2326, p. 1321. [Google Scholar]
  35. Denevi, B.; Yasanayake, C.; Jolliff, B.; Lawrence, S.; Hiroi, T.; Martin, A. The Spectral Properties of Lunar Agglutinates. In Proceedings of the Lunar and Planetary Science Conference, Virtual Meeting. 15–19 March 2021; Number 2548, p. 2368. [Google Scholar]
  36. Sun, L.; Lucey, P.G. Unmixing Mineral Abundance and Mg# With Radiative Transfer Theory: Modeling and Applications. J. Geophys. Res. Planets 2021, 126, e2020JE006691. [Google Scholar]
  37. Cloutis, E.A. Pyroxene reflectance spectra: Minor absorption bands and effects of elemental substitutions. J. Geophys. Res. Planets 2002, 107, 6:1–6:12. [Google Scholar] [CrossRef] [Green Version]
  38. Burns, R.G.; Burns, R.G. Mineralogical Applications of Crystal Field Theory; Number 5; Cambridge University Press: Cambridge, UK, 1993. [Google Scholar]
  39. Klima, R.L.; Pieters, C.M.; Dyar, M.D. Spectroscopy of synthetic Mg-Fe pyroxenes I: Spin-allowed and spin-forbidden crystal field bands in the visible and near-infrared. Meteorit. Planet. Sci. 2007, 42, 235–253. [Google Scholar] [CrossRef]
  40. Klima, R.L.; Dyar, M.D.; Pieters, C.M. Near-infrared spectra of clinopyroxenes: Effects of calcium content and crystal structure. Meteorit. Planet. Sci. 2011, 46, 379–395. [Google Scholar] [CrossRef]
  41. Neal, R.M. Probabilistic Inference Using Markov Chain Monte Carlo Methods; Department of Computer Science, University of Toronto Toronto: Ontario, ON, Canada, 1993. [Google Scholar]
  42. Metropolis, N.; Rosenbluth, A.W.; Rosenbluth, M.N.; Teller, A.H.; Teller, E. Equation of state calculations by fast computing machines. J. Chem. Phys. 1953, 21, 1087–1092. [Google Scholar] [CrossRef] [Green Version]
  43. Hastings, W.K. Monte Carlo sampling methods using Markov chains and their applications. Biometrika 1970, 57, 97–109. [Google Scholar] [CrossRef]
  44. Hoffman, M.D.; Gelman, A. The No-U-Turn sampler: Adaptively setting path lengths in Hamiltonian Monte Carlo. J. Mach. Learn. Res. 2014, 15, 1593–1623. [Google Scholar]
  45. Haario, H.; Saksman, E.; Tamminen, J. An adaptive Metropolis algorithm. Bernoulli 2001, 7, 223–242. [Google Scholar] [CrossRef] [Green Version]
  46. Haario, H.; Laine, M.; Mira, A.; Saksman, E. DRAM: Efficient adaptive MCMC. Stat. Comput. 2006, 16, 339–354. [Google Scholar] [CrossRef]
  47. Neal, R.M. MCMC using Hamiltonian dynamics. In Handbook of Markov Chain Monte Carlo; Chapman and Hall/CRC: London, UK, 2011; Volume 2, Ch. 5. [Google Scholar]
  48. Bastien, F.; Lamblin, P.; Pascanu, R.; Bergstra, J.; Goodfellow, I.; Bergeron, A.; Bouchard, N.; Warde-Farley, D.; Bengio, Y. Theano: New features and speed improvements. arXiv 2012, arXiv:1211.5590. [Google Scholar]
  49. Salvatier, J.; Wiecki, T.V.; Fonnesbeck, C. Probabilistic programming in Python using PyMC3. Peerj Comput. Sci. 2016, 2, e55. [Google Scholar] [CrossRef] [Green Version]
  50. Fu, Z.; Robles-Kelly, A.; Caelli, T.; Tan, R.T. On automatic absorption detection for imaging spectroscopy: A comparative study. IEEE Trans. Geosci. Remote Sens. 2007, 45, 3827–3844. [Google Scholar] [CrossRef]
  51. Bishop, C.M. Pattern Recognition and Machine Learning; Springer: New York, NY, USA, 2006. [Google Scholar]
  52. Lawrence, D.; Feldman, W.; Barraclough, B.; Binder, A.; Elphic, R.; Maurice, S.; Thomsen, D. Global elemental maps of the Moon: The Lunar Prospector gamma-ray spectrometer. Science 1998, 281, 1484–1489. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Hess, M.; Wöhler, C.; Bhatt, M.; Berezhnoy, A.; Grumpe, A.; Wohlfarth, K.; Bhardwaj, A.; Shevchenko, V. Processes governing the VIS/NIR spectral reflectance behavior of lunar swirls. Astron. Astrophys. 2020, 639, A12. [Google Scholar] [CrossRef]
  54. Trang, D.; Lucey, P.G. Improved space weathering maps of the lunar surface through radiative transfer modeling of Kaguya multiband imager data. Icarus 2019, 321, 307–323. [Google Scholar] [CrossRef]
  55. Trang, D.; Lucey, P.G.; Izenberg, N.R. Radiative transfer modeling of MESSENGER VIRS spectra: Detection and mapping of submicroscopic iron and carbon. Icarus 2017, 293, 206–217. [Google Scholar] [CrossRef]
  56. Grumpe, A.; Mengewein, N.; Rommel, D.; Mall, U.; Wöhler, C. Interpreting spectral unmixing coefficients: From spectral weights to mass fractions. Icarus 2018, 299, 1–14. [Google Scholar] [CrossRef]
  57. Housley, R.; Cirlin, E.; Goldberg, I.; Crowe, H. Ferromagnetic resonance studies of lunar core stratigraphy. In Proceedings of the Lunar and Planetary Science Conference Proceedings, Houston, TX, USA, 15–19 March 1976; Volume 7, pp. 13–26. [Google Scholar]
  58. Arnaut, M.; Wohlfarth, K.; Wöhler, C. The interaction between multiple nanophase iron particles changes the slope of lunar reflectance spectra. In Proceedings of the European Planetary Science Congress, Virtual Meeting. 13–24 September 2021. [Google Scholar] [CrossRef]
  59. Lapotre, M.G.; Ehlmann, B.; Minson, S.E.; Arvidson, R.; Ayoub, F.; Fraeman, A.; Ewing, R.; Bridges, N. Compositional variations in sands of the Bagnold Dunes, Gale Crater, Mars, from visible-shortwave infrared spectroscopy and comparison with ground truth from the Curiosity rover. J. Geophys. Res. Planets 2017, 122, 2489–2509. [Google Scholar] [CrossRef]
  60. Ren, Z.; Sun, L.; Zhai, Q. Improved k-means and spectral matching for hyperspectral mineral mapping. Int. J. Appl. Earth Obs. Geoinf. 2020, 91, 102154. [Google Scholar] [CrossRef]
  61. Arnaut, M.; Wohlfarth, K.; Wöhler, C. A Study of Hyperspectral Clustering—How To Speed Up Spectral Unmixing Of Moon Mineralogy Mapper Data. In Proceedings of the Lunar and Planetary Science Conference, Virtual Meeting. 16–20 March 2020; p. 3008. [Google Scholar]
Figure 1. Extinction efficiencies and SSA spectra of the space weathering components npFe0 and mpFe0 calculated according to Wohlfarth et al. [24]. (a) Extinction efficiencies Q E ; (b) SSA.
Figure 1. Extinction efficiencies and SSA spectra of the space weathering components npFe0 and mpFe0 calculated according to Wohlfarth et al. [24]. (a) Extinction efficiencies Q E ; (b) SSA.
Remotesensing 13 04702 g001
Figure 2. Phase function values for 30° phase angle. The Soil and the npFe0 particles phase functions are independent of the wavelength in contrast to the mpFe0 particles.
Figure 2. Phase function values for 30° phase angle. The Soil and the npFe0 particles phase functions are independent of the wavelength in contrast to the mpFe0 particles.
Remotesensing 13 04702 g002
Figure 3. Example calculation of the continuum removed spectrum for endmember AG-TJM-009 from the RELAB database space weathered with 1.5 wt.% npFe0 and 0.8 wt.% mpFe. The continuum is calculated as the linear interpolation between the reflectance factor at 730.48 nm, 1578.86 nm and 2497.11 nm. These channels usually lie outside the absorption bands. (a) Reflectance spectrum; (b) continuum removed spectrum.
Figure 3. Example calculation of the continuum removed spectrum for endmember AG-TJM-009 from the RELAB database space weathered with 1.5 wt.% npFe0 and 0.8 wt.% mpFe. The continuum is calculated as the linear interpolation between the reflectance factor at 730.48 nm, 1578.86 nm and 2497.11 nm. These channels usually lie outside the absorption bands. (a) Reflectance spectrum; (b) continuum removed spectrum.
Remotesensing 13 04702 g003
Figure 4. Probability density function (PDF) of the Beta distribution for some examples of shape parameters α and β .
Figure 4. Probability density function (PDF) of the Beta distribution for some examples of shape parameters α and β .
Remotesensing 13 04702 g004
Figure 5. The endmembers for the experiments in Section 3.1 and Section 3.3. All the endmembers have a grain size of 0–45 μ m. The RELAB IDs are: OPX is the hypersthene AG-TJM-009; CPX Type B is the augite AG-TJM-010; OLV is the olivine PO-EAC-056; PLG is the plagioclase PL-EAC-029; CPX Type A is the pyroxene diopside PP-ALS-105; ILM is the ilmenite SC-EAC-034.
Figure 5. The endmembers for the experiments in Section 3.1 and Section 3.3. All the endmembers have a grain size of 0–45 μ m. The RELAB IDs are: OPX is the hypersthene AG-TJM-009; CPX Type B is the augite AG-TJM-010; OLV is the olivine PO-EAC-056; PLG is the plagioclase PL-EAC-029; CPX Type A is the pyroxene diopside PP-ALS-105; ILM is the ilmenite SC-EAC-034.
Remotesensing 13 04702 g005
Figure 6. The probability distribution function of the prior on the elemental abundances for a theoretical abundance below 1 wt.%.
Figure 6. The probability distribution function of the prior on the elemental abundances for a theoretical abundance below 1 wt.%.
Remotesensing 13 04702 g006
Figure 7. Violin plots for synthetic mixtures 0, 5, and 15 for the MCMC sampling without a prior on the elemental abundances. These mixtures do not contain ilmenite, but the posterior distribution contains many solutions with significant amounts of ilmenite. Yellow diamonds mark the true abundances that were used for mixing the endmembers. (a) Without elemental prior (mix0); (b) Without elemental prior (mix0); (c) Without elemental prior (mix5); (d) Without elemental prior (mix5); (e) Without elemental prior (mix15); (f) Without elemental prior (mix15).
Figure 7. Violin plots for synthetic mixtures 0, 5, and 15 for the MCMC sampling without a prior on the elemental abundances. These mixtures do not contain ilmenite, but the posterior distribution contains many solutions with significant amounts of ilmenite. Yellow diamonds mark the true abundances that were used for mixing the endmembers. (a) Without elemental prior (mix0); (b) Without elemental prior (mix0); (c) Without elemental prior (mix5); (d) Without elemental prior (mix5); (e) Without elemental prior (mix15); (f) Without elemental prior (mix15).
Remotesensing 13 04702 g007
Figure 8. Violin plots for synthetic mixtures 0, 5, and 15 for the MCMC sampling including a prior on the elemental abundances. Compared to the version without a prior the ilmenite abundance could effectively be constrained. This also slightly reduced the uncertainties associated with the other minerals. Yellow diamonds mark the true abundances that were used for mixing the endmembers. (a) With TiO2 and Al2O3 prior (mix0); (b) With TiO2 and Al2O3 prior (mix0); (c) With TiO2 and Al2O3 prior (mix5); (d) With TiO2 and Al2O3 prior (mix5); (e) With TiO2 and Al2O3 prior (mix15); (f) With TiO2 and Al2O3 prior (mix15).
Figure 8. Violin plots for synthetic mixtures 0, 5, and 15 for the MCMC sampling including a prior on the elemental abundances. Compared to the version without a prior the ilmenite abundance could effectively be constrained. This also slightly reduced the uncertainties associated with the other minerals. Yellow diamonds mark the true abundances that were used for mixing the endmembers. (a) With TiO2 and Al2O3 prior (mix0); (b) With TiO2 and Al2O3 prior (mix0); (c) With TiO2 and Al2O3 prior (mix5); (d) With TiO2 and Al2O3 prior (mix5); (e) With TiO2 and Al2O3 prior (mix15); (f) With TiO2 and Al2O3 prior (mix15).
Remotesensing 13 04702 g008
Figure 9. Confidence plots for synthetic mixtures 0, 5, and 15. In the left column, the plots of the sampling without an elemental prior are displayed. For the confidence plots in the right column, prior information about the elements TiO2 and Al2O3 was included. (a) Without elemental prior (mix0); (b) With TiO2 and Al2O3 prior (mix0); (c) Without elemental prior (mix5); (d) With TiO2 and Al2O3 prior (mix5); (e) Without elemental prior (mix15); (f) With TiO2 and Al2O3 prior (mix15).
Figure 9. Confidence plots for synthetic mixtures 0, 5, and 15. In the left column, the plots of the sampling without an elemental prior are displayed. For the confidence plots in the right column, prior information about the elements TiO2 and Al2O3 was included. (a) Without elemental prior (mix0); (b) With TiO2 and Al2O3 prior (mix0); (c) Without elemental prior (mix5); (d) With TiO2 and Al2O3 prior (mix5); (e) Without elemental prior (mix15); (f) With TiO2 and Al2O3 prior (mix15).
Remotesensing 13 04702 g009
Figure 10. Results for the most mature spectrum, which also has the highest ilmenite concentration. The spectrum becomes very dark. The abundance of npFe0 and mpFe0 is well defined either with or without the prior information. The uncertainties of the mineral abundances decrease when including the prior. The yellow diamonds in the violin plots mark the true abundances that were used for mixing the endmembers. (a) Without elemental prior (mix47); (b) Without elemental prior (mix47); (c) With TiO2 and Al2O3 prior (mix47); (d) With TiO2 and Al2O3 prior (mix47); (e) Without elemental prior (mix47); (f) With TiO2 and Al2O3 prior (mix47).
Figure 10. Results for the most mature spectrum, which also has the highest ilmenite concentration. The spectrum becomes very dark. The abundance of npFe0 and mpFe0 is well defined either with or without the prior information. The uncertainties of the mineral abundances decrease when including the prior. The yellow diamonds in the violin plots mark the true abundances that were used for mixing the endmembers. (a) Without elemental prior (mix47); (b) Without elemental prior (mix47); (c) With TiO2 and Al2O3 prior (mix47); (d) With TiO2 and Al2O3 prior (mix47); (e) Without elemental prior (mix47); (f) With TiO2 and Al2O3 prior (mix47).
Remotesensing 13 04702 g010
Figure 11. Differences between mean predicted solution and true abundance. The differences in ilmenite and plagioclase abundance can be reduced with a prior on the TiO2 and Al2O3 abundance. Each circle represents one mixture. (a) Without prior; (b) TiO2 and Al2O3 prior.
Figure 11. Differences between mean predicted solution and true abundance. The differences in ilmenite and plagioclase abundance can be reduced with a prior on the TiO2 and Al2O3 abundance. Each circle represents one mixture. (a) Without prior; (b) TiO2 and Al2O3 prior.
Remotesensing 13 04702 g011
Figure 12. Uncertainties associated with each of the mineral abundances. With a prior all uncertainties, except for the orthopyroxene, which is already quite low, can be reduced. Each circle represents one mixture. (a) Without prior; (b) TiO2 and Al2O3 prior.
Figure 12. Uncertainties associated with each of the mineral abundances. With a prior all uncertainties, except for the orthopyroxene, which is already quite low, can be reduced. Each circle represents one mixture. (a) Without prior; (b) TiO2 and Al2O3 prior.
Remotesensing 13 04702 g012
Figure 13. Correlation between true smFe abundance and the uncertainties (standard deviations) of the plagioclase abundance. It can be seen that the uncertainties generally decrease with the elemental prior. (a) Without prior; (b) TiO2 and Al2O3 prior.
Figure 13. Correlation between true smFe abundance and the uncertainties (standard deviations) of the plagioclase abundance. It can be seen that the uncertainties generally decrease with the elemental prior. (a) Without prior; (b) TiO2 and Al2O3 prior.
Remotesensing 13 04702 g013
Figure 14. The endmembers for the experiment in Section 3.2 taken from Rommel et al. [13]. PRG: pargasite, AUG: augite, OLV: olivine, PLG: plagioclase (labradorite B), FS: ferrosilite, ILM: hemo-ilmenite.
Figure 14. The endmembers for the experiment in Section 3.2 taken from Rommel et al. [13]. PRG: pargasite, AUG: augite, OLV: olivine, PLG: plagioclase (labradorite B), FS: ferrosilite, ILM: hemo-ilmenite.
Remotesensing 13 04702 g014
Figure 15. Violin plots for mix0 of the fresh laboratory mixtures. Yellow diamonds mark the theoretical abundances that were used for mixing the endmembers. (a) Without prior; (b) TiO2 prior.
Figure 15. Violin plots for mix0 of the fresh laboratory mixtures. Yellow diamonds mark the theoretical abundances that were used for mixing the endmembers. (a) Without prior; (b) TiO2 prior.
Remotesensing 13 04702 g015
Figure 16. Confidence for mix0 of the fresh laboratory mixtures. Uncertainty of the spectrum slightly increases with the introduction of a prior on the TiO2 abundance. (a) Without prior; (b) TiO2 prior.
Figure 16. Confidence for mix0 of the fresh laboratory mixtures. Uncertainty of the spectrum slightly increases with the introduction of a prior on the TiO2 abundance. (a) Without prior; (b) TiO2 prior.
Remotesensing 13 04702 g016
Figure 17. Violin plots for mix10 of the fresh laboratory mixtures. Yellow diamonds mark the theoretical abundances that were used for mixing the endmembers. (a) Without prior; (b) TiO2 prior.
Figure 17. Violin plots for mix10 of the fresh laboratory mixtures. Yellow diamonds mark the theoretical abundances that were used for mixing the endmembers. (a) Without prior; (b) TiO2 prior.
Remotesensing 13 04702 g017
Figure 18. Confidence for mix10 of the fresh laboratory mixtures. The prediction of the 0.9 μ m absorption is slightly less certain for the version with a prior on TiO2 abundance. Mean predicted spectra are almost identical, even though the weights are quite different. For the version without an elemental prior the sampled σ is 0.011 and 0.012 for the version with a TiO2 prior. (a) Without prior; (b) TiO2 prior.
Figure 18. Confidence for mix10 of the fresh laboratory mixtures. The prediction of the 0.9 μ m absorption is slightly less certain for the version with a prior on TiO2 abundance. Mean predicted spectra are almost identical, even though the weights are quite different. For the version without an elemental prior the sampled σ is 0.011 and 0.012 for the version with a TiO2 prior. (a) Without prior; (b) TiO2 prior.
Remotesensing 13 04702 g018
Figure 19. Prior on the element TiO2 mainly present in ilmenite, significantly reduces the differences in predicted plagioclase and also pargasite abundance and theoretical abundances. These spectra of the three endmembers are relatively flat and do not show significant absorption bands. (a) Without prior; (b) TiO2 prior.
Figure 19. Prior on the element TiO2 mainly present in ilmenite, significantly reduces the differences in predicted plagioclase and also pargasite abundance and theoretical abundances. These spectra of the three endmembers are relatively flat and do not show significant absorption bands. (a) Without prior; (b) TiO2 prior.
Remotesensing 13 04702 g019
Figure 20. The uncertainties of the respective endmembers are in general relatively low. Only the mixture containing ilmenite and ferrosilite is highly uncertain. For that particular sample, no good solution has been found considering the prior on the sum of the endmember weights. (a) Without prior; (b) TiO2 prior.
Figure 20. The uncertainties of the respective endmembers are in general relatively low. Only the mixture containing ilmenite and ferrosilite is highly uncertain. For that particular sample, no good solution has been found considering the prior on the sum of the endmember weights. (a) Without prior; (b) TiO2 prior.
Remotesensing 13 04702 g020
Figure 21. Uncertainties associated with each of the mineral abundances for the unmixing of the LSCC spectra. (a) TiO2 and Al2O3 prior; (b) All prior.
Figure 21. Uncertainties associated with each of the mineral abundances for the unmixing of the LSCC spectra. (a) TiO2 and Al2O3 prior; (b) All prior.
Remotesensing 13 04702 g021
Figure 22. Results for LSCC sample 10084 with Ti/Al prior. Yellow diamonds mark the estimated theoretical mineral fractional abundances. Plagioclase is underestimated, while olivine is generally overestimated. This is a mature spectrum with an Is/FeO value of 88, therefore, also the estimated npFe0 and mpFe0 abundance is high. The spectrum is generally very dark. (a) Endmember abundances (10084); (b) smFe abundances (10084); (c) Confidence plot (10084).
Figure 22. Results for LSCC sample 10084 with Ti/Al prior. Yellow diamonds mark the estimated theoretical mineral fractional abundances. Plagioclase is underestimated, while olivine is generally overestimated. This is a mature spectrum with an Is/FeO value of 88, therefore, also the estimated npFe0 and mpFe0 abundance is high. The spectrum is generally very dark. (a) Endmember abundances (10084); (b) smFe abundances (10084); (c) Confidence plot (10084).
Remotesensing 13 04702 g022
Figure 23. Results for LSCC sample 14260 with Ti/Al prior. This mature mare spectrum (Is/FeO = 93.3) can be well characterized. Yellow diamonds mark the estimated theoretical mineral fractional abundances. Except for a small difference in theoretical and predicted clinopyroxene abundance, the fit is good. (a) Endmember abundances (14260); (b) smFe abundances (14260); (c) Confidence plot (14260).
Figure 23. Results for LSCC sample 14260 with Ti/Al prior. This mature mare spectrum (Is/FeO = 93.3) can be well characterized. Yellow diamonds mark the estimated theoretical mineral fractional abundances. Except for a small difference in theoretical and predicted clinopyroxene abundance, the fit is good. (a) Endmember abundances (14260); (b) smFe abundances (14260); (c) Confidence plot (14260).
Remotesensing 13 04702 g023
Figure 24. Results for LSCC sample 62231 with Ti/Al prior. This LSCC sample is a mature (Is/FeO = 116.7) highland spectrum. There are almost no absorption bands visible. Yellow diamonds mark the true mineral fractional abundances. The mean predicted values are close to the theoretical mineral abundances. (a) Endmember abundances (62231); (b) smFe abundances (62231); (c) smFe abundances (62231).
Figure 24. Results for LSCC sample 62231 with Ti/Al prior. This LSCC sample is a mature (Is/FeO = 116.7) highland spectrum. There are almost no absorption bands visible. Yellow diamonds mark the true mineral fractional abundances. The mean predicted values are close to the theoretical mineral abundances. (a) Endmember abundances (62231); (b) smFe abundances (62231); (c) smFe abundances (62231).
Remotesensing 13 04702 g024
Figure 25. Results for LSCC sample 79221 both with the prior on TiO2 and Al2O3 abundance and a prior on all relevant elements. Including prior information about the other elements as well can help in mitigating uncertainties between clinopyroxene and olivine. This is also a mature mare spectrum with an Is/FeO value of 91.0. Yellow diamonds mark the theoretical mineral fractional abundances. (a) Endmember abundances with Ti/Al prior (79221); (b) smFe abundances with Ti/Al prior (79221); (c) Endmember abundances all elements prior (79221); (d) smFe abundances all elements prior (79221); (e) Confidence Ti/Al prior (79221); (f) Confidence all elements prior (79221).
Figure 25. Results for LSCC sample 79221 both with the prior on TiO2 and Al2O3 abundance and a prior on all relevant elements. Including prior information about the other elements as well can help in mitigating uncertainties between clinopyroxene and olivine. This is also a mature mare spectrum with an Is/FeO value of 91.0. Yellow diamonds mark the theoretical mineral fractional abundances. (a) Endmember abundances with Ti/Al prior (79221); (b) smFe abundances with Ti/Al prior (79221); (c) Endmember abundances all elements prior (79221); (d) smFe abundances all elements prior (79221); (e) Confidence Ti/Al prior (79221); (f) Confidence all elements prior (79221).
Remotesensing 13 04702 g025aRemotesensing 13 04702 g025b
Figure 26. Relationship of npFe0 and mpFe0 abundance to the Is/FeO maturity parameter and the ferromagnetic resonance (FMR) Is. The FMR parameter is sensitive to small iron particles in the size range between about 4 nm and 33 nm [57]. The correlation between npFe0 and Is is clearly visible as well as for the larger mpFe0 particles. (a) npFe0, r = 0.39; (b) mpFe0, r = 0.34; (c) npFe0, r = 0.88; (d) mpFe0, r = 0.87.
Figure 26. Relationship of npFe0 and mpFe0 abundance to the Is/FeO maturity parameter and the ferromagnetic resonance (FMR) Is. The FMR parameter is sensitive to small iron particles in the size range between about 4 nm and 33 nm [57]. The correlation between npFe0 and Is is clearly visible as well as for the larger mpFe0 particles. (a) npFe0, r = 0.39; (b) mpFe0, r = 0.34; (c) npFe0, r = 0.88; (d) mpFe0, r = 0.87.
Remotesensing 13 04702 g026
Figure 27. The correlation between the mean predicted npFe0 and mpFe0 abundance is r = 0.90. Therefore, with increasing npFe0 abundance typically also the mpFe0 abundance increases. The factor npFe mpFe approximately amounts to 2.
Figure 27. The correlation between the mean predicted npFe0 and mpFe0 abundance is r = 0.90. Therefore, with increasing npFe0 abundance typically also the mpFe0 abundance increases. The factor npFe mpFe approximately amounts to 2.
Remotesensing 13 04702 g027
Table 1. Material parameters.
Table 1. Material parameters.
ParameternpFe0mpFe0Medium
ρ [g/cm3]7.5687.5683.5
D [nm]1015030  μ m
Q E ( λ ) Figure 1aFigure 1a2.1
Table 2. Prior setups for the elemental abundances used during this work. The shape parameter β is then calculated according to Equation (21). If α e l e m > 1 and the theoretical elemental abundance is below 1 wt.%, the shape parameters will be set to α e l e m = 1 and β e l e m = 10 .
Table 2. Prior setups for the elemental abundances used during this work. The shape parameter β is then calculated according to Equation (21). If α e l e m > 1 and the theoretical elemental abundance is below 1 wt.%, the shape parameters will be set to α e l e m = 1 and β e l e m = 10 .
SiO2TiO2Al2O3Cr2O3MgOCaOMnOFeONa2OK2OP2O5SO3
None α e l e m 111111111111
TiO2  α e l e m 1401111111111
Ti/Al α e l e m 14040111111111
All α e l e m 303030130301301111
Table 3. Pearson r coefficients for the relationship between theoretical smFe abundance and standard deviation of the respective endmembers. Except for ilmenite with an elemental prior all correlation coefficients are positive. Therefore, with increasing smFe abundance, the uncertainties of the mineral abundances increase.
Table 3. Pearson r coefficients for the relationship between theoretical smFe abundance and standard deviation of the respective endmembers. Except for ilmenite with an elemental prior all correlation coefficients are positive. Therefore, with increasing smFe abundance, the uncertainties of the mineral abundances increase.
OPXCPXBOLVPLGCPXAILM
rNone0.930.710.470.540.480.13
rTi/Al0.940.810.560.700.580.02
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hess, M.; Wilhelm, T.; Wöhler, C.; Wohlfarth, K. Uncertainty Introduced by Darkening Agents in the Lunar Regolith: An Unmixing Perspective. Remote Sens. 2021, 13, 4702. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13224702

AMA Style

Hess M, Wilhelm T, Wöhler C, Wohlfarth K. Uncertainty Introduced by Darkening Agents in the Lunar Regolith: An Unmixing Perspective. Remote Sensing. 2021; 13(22):4702. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13224702

Chicago/Turabian Style

Hess, Marcel, Thorsten Wilhelm, Christian Wöhler, and Kay Wohlfarth. 2021. "Uncertainty Introduced by Darkening Agents in the Lunar Regolith: An Unmixing Perspective" Remote Sensing 13, no. 22: 4702. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13224702

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop