A New Spectral Transformation Approach and Quantitative Analysis for MarSCoDe Laser-Induced Breakdown Spectroscopy (LIBS) Data

: Zhurong rover successfully landed on the southern of Utopia Planet of Mars on 15 May 2021. One laser-induced breakdown spectroscopy (LIBS) system, the main payload of the Mars Surface Composition Detector (MarSCoDe), was installed on the Zhurong rover aimed to measure the elements and their abundance in Martian regolith. Now, there are three sets of LIBS system (ChemCam, SuperCam and MarSCoDe) working on Mars at difference landing sites with diverse geologic features. For Mars exploration, cross-validation is necessary to expand the model compatibility, test data validity, and get more available data of the same type payloads. Spectral transformation approach is the ﬁrst step and crucial for cross-validation of LIBS analysis model. Herein, a new 4-step spectral transformation approach was proposed to transform the LIBS spectra between three different LIBS systems (i.e., ChemCam, MarSCoDe, SDU-LIBS (recorded by self-built LIBS system)), whose data were partly different in spectral characteristics. Based on this approach, SDU-LIBS and MarSCoDe spectra data were transformed into ChemCam uniform and then the three kinds of LIBS data can have more similar spectral features and share one PLS (partial least squares) model for quantitative analysis. Our approach enables to make up the signal differences between different LIBS systems and gets acceptable quantitative analysis results of SDU-LIBS and MarSCoDe spectra using quantitative PLS model built by ChemCam calibration sample set. This work veriﬁed feasibility and availability of our approach for cross validation of different LIBS systems. Based on this method, MarSCoDe data were analyzed and got the preliminary satisfying results although no analysis model of laboratory replica payload was available under the existing conditions.


Introduction
Mars is by far the most understood terrestrial planet other than Earth. NASA and other space agencies have conducted multiple missions to Mars over the past 58 years, including orbital and in situ probes. As China's first independent Mars mission, Tianwen-1 implemented orbital and in situ detection at one time [1][2][3][4]. The Zhurong rover successfully landed on the south of Utopia Planitia of Mars on 15 May 2021 where it is believed to have experienced mud volcano, water/ice, and other sulfate mineral formed by geologic activities associated with water [5][6][7][8]. The understanding of Martian geology will deepen by Zhurong rover exploration.
As one of the most important scientific instruments on Zhurong Rover, Mars Surface Composition Detector (MarSCoDe) adopted LIBS to obtain the elemental compositions and their abundance for surface materials on Mars [1,9]. LIBS can quickly and conveniently detect the elemental compositions of the target at sub-millimeter scale in a real time and remote way without sample preparation. Luminous plasma will be produced when LIBS focuses the laser on the target. The elements in the sample can be excited and their emission lines can be acquired and analyzed. LIBS is a powerful tool in planetary exploration [10][11][12][13][14][15][16] which has been applied for Mars explorations such as ChemCam aboard Curiosity rover at Gale Crate [17][18][19] and SuperCam on Perseverance rover at Jezero Crater (19 February 2021) [14,15]. ChemCam revealed the possible biohabitability in the past and has a systematic understanding of the geological history and geological features of the landing site of Curiosity in Gale crater by LIBS detection [20][21][22][23][24][25][26][27][28].
LIBS spectra can be interpreted by using a variety of multiple linear analysis (MLA) methods, for example, PLS is a popular method for modeling relations between sets of intuitive variables by means of latent variables and be used to construct a linear multivariate regression model [29,30]. It is an effective way which has been successfully applied to interpret in situ ChemCam LIBS spectra and obtained derived values similar as the actual values of elemental abundance [17][18][19].
Recently, Zhurong rover has traveled some distances from landing site to south and obtained 32 in situ LIBS spectra from dunes, soils, and rocks. It will continue to move southward to a regarded ancient geologically active area [31]. The LIBS data from Zhurong will continue to deepen our understanding of Mars geological evolution. However, the harsh environment, especially the various temperatures, will cause a certain degree of spectral variation. A PSO (particle swarm optimization) approach was used to alleviate the situation [32]. Temperature influences are still present in the LIBS data in a certain degree, and affect the further analysis. At present, MarSCoDe lacks an available prediction model because testing conditions of MarSCoDe replica is still lacking for model building.
However, for planetary exploration, it is important for data cross-validation, which is necessary to expand the model compatibility, test data validity, and get more available data for conjoint analysis. Therefore, spectral transformation approach as the first stage for cross-validation of LIBS analysis model is crucial for unifying the data format and getting ready for quantitative analysis.
Spectral analysis technology is a powerful tool for complex target analysis. Due to different spectral response and working parameters, the spectra data of different instruments are little difference in spectral response characteristics (such as spectral intensity, spectral resolution) which restrict the further analysis and need to compensate by data transformation approach. Similarly, "calibration transfer" [33][34][35] has been reported to treat near infrared (NIR) spectroscopy data aiming to build a quantitative analysis model to analyze the data of other NIR system. Large amounts of NIR data from two instruments, much training time, and several spectral processing methods were needed to train their models. However, a similar approach has not yet been reported on LIBS. Moreover, those approaches are hard to implement in the case of a small amount of data from MarSCoDe.
The LIBS system (SDU-LIBS) [11], similar to MarSCoDe with different paraments, and a Mars chamber [36] have been built at Shandong University. A Martian Analogues Library (MAL) [37] and the corresponding LIBS database recorded by SDU-LIBS have also been established aiming to interpret the Martian in situ spectra [11]. In this study, spectral data acquired by SDU-LIBS were used for analysis.
In this paper, a 4-step spectral transformation approach was proposed to transform MarSCoDe and SDU-LIBS [11,36]) spectra to ChemCam spectral format, regarded as standard data format, to make those three sources of data have the same spectral data format (step size (the distance between two adjacent pixels), spectral band range, spectral intensity unit, et al.). At same time, necessary spectral normalization was also performed to partly remove the spectral variations effect due to temperature variations of Martian external environment and laser power variation during target excitation. The 4-step approach was done to unify the spectral data formats of ChemCam, MarSCoDe, and SDU-LIBS without training. The small amount of MarSCoDe spectral data does not prevent the preformation of our 4-step approach. After data were transformed by our proposed 4-step approach, SDU-LIBS and MarSCoDe LIBS spectra can be quantitatively analyzed using the PLS model built by ChemCam calibration database. The elemental abundant analytical ability of our approach was verified by quantitative analysis of SDU-LIBS samples which have definite element abundance. Then the MarSCoDe spectra were quantitatively analyzed by the validated PLS model to interpret the regolith components detected by Zhurong rover.

Samples and Data
At the beginning of this work, only 32 in situ spectra and 44 on-board calibration targets spectra from MarSCoDe are available. Every LIBS spectrum is the average of 60 shoots. A total of 45 high SNR (signal-to-noise ratio) spectra (26 in situ spectra and 19 on-board calibration target spectra) were analyzed in this work.
ChemCam calibration samples were excited with 250 consecutive laser pulse on the surface of each sample, and a total of five average spectra form five different locations were recorded on one target.
A total of 11 identical ChemCam calibration samples (listed in Table 1) were specially selected and recorded by SDU-LIBS system aiming to verify the effectiveness of our proposed approach in this work. Three spectra (the average of 10 shots, integration time of 10 ms per shot) were obtained from three different locations using SDU-LIBS system [11] under simulated Mars conditions (CO 2 , 700 Pa). Wavelength calibration and spectral radiance calibration of SDU-LIBS have been done before SDU-LIBS data recordation. FeO T refers to total iron, including both ferric and ferrous state [19].
Routine data pre-processing was also performed in this work before the application of our approach. Wavelet transform was used to remove the white noise caused by dark currents [18,19]. The continuum background, was identified by the local minimum values and spline interpolation [18], and then removed from original spectra.
Because the LIBS spectral bands of MarSCoDe and ChemCam do not match well, only the same spectral bands of those two systems with a higher SNR (246.8 nm < λ < 388.5 nm in the UV, 387.9 nm < λ < 469.1 in the VIO and 515.6nm < λ < 849.1 nm in the VNIR) [18] were selected for further spectral procession. Because the instruments have different step size, the points recorded by the three systems are not aligned with each other. For example, the small difference in one selected Na(I) peak position of three systems also was demonstrated in Figure 1. The emission line of Na(I) in MarSCoDe spectra was located at 588.6 nm to 588.8 nm, while at 589.1 nm in ChemCam spectra and 589.7 nm in SDU-LIBS. The differences were eliminated by our approach and are described in the next section.

Spectral Transformation Approach
MarSCoDe shots the laser on the Martian regolith/rock in two directions and collected one set of LIBS spectra at every direction of every sampling site during Zhurong rover exploration respectively. At the same time, some important information such as the rover location, detector temperature were also recorded. All the wavelengths of LIBS data were calibrated by PSO method [32] during every LIBS data collection. However, wavelength drift [1] was still observed which should be corrected to get more stable LIBS data.

Spectral Transformation Approach
MarSCoDe shots the laser on the Martian regolith/rock in two directions and collected one set of LIBS spectra at every direction of every sampling site during Zhurong rover exploration respectively. At the same time, some important information such as the rover location, detector temperature were also recorded. All the wavelengths of LIBS data were calibrated by PSO method [32] during every LIBS data collection. However, wavelength drift [1] was still observed which should be corrected to get more stable LIBS data.
At present, an accessible way is worth trying to build a quantitative analysis model by using ChemCam data as calibration dataset instead of MarSCoDe calibration dataset which has not been built until now. All spectra must be consistent in one format which is the highest priority to perform before next analysis. To this end, our proposed spectral transformation approach, made up of 4 step (Figure 2), was carried out to accomplish this transformation as described briefly below. At present, an accessible way is worth trying to build a quantitative analysis model by using ChemCam data as calibration dataset instead of MarSCoDe calibration dataset which has not been built until now. All spectra must be consistent in one format which is the highest priority to perform before next analysis. To this end, our proposed spectral transformation approach, made up of 4 step (Figure 2), was carried out to accomplish this transformation as described briefly below. (1) Step 1 Intensity Unit Transformation The physical unit of ChemCam LIBS spectra (photons per second [10]) and MarSCoDe LIBS spectra (Watt [1]) original intensity were different which must be unified first. Based on Equation (1), their intensity unit can be consistent with each other.
RIAfterTrans i and RIMarSCoDe i are the intensity before and after transformation of the pixel i, h is the Planck constant, c is the lightspeed, and λi is the wavelength of the RIMarSCoDe i. The physical unit of ChemCam LIBS spectra (photons per second [10]) and MarSCoDe LIBS spectra (Watt [1]) original intensity were different which must be unified first. Based on Equation (1), their intensity unit can be consistent with each other.
RI AfterTrans i and RI MarSCoDe i are the intensity before and after transformation of the pixel i, h is the Planck constant, c is the lightspeed, and λ i is the wavelength of the RI MarSCoDe i .
(2) Step 2 Wavelength Recalibration The wavelength recalibration aimed to calibrate/unify the peak position of MarSCoDe to ChemCam spectra. Some obvious peaks of major elements without the self-absorption (listed in Table 2) were selected to perform the wavelength recalibration. The wavelength recalibration functions of every spectrometer were fitted by using the wavelength position of ChemCam spectra and the pixel positions of the selected peaks in each MarSCoDe spectra, and then applied to each MarSCoDe spectra to accomplish the wavelength recalibration. The spectral assignment based on NIST Atomic Spectra Database [38]. (3) Step 3 Interpolation After the intensity and wavelength of MarSCoDe were unified, the difference in spectral step size remained unsolved. As shown in Figure 3a, the wavelength of ChemCam spectra (the red verticals along the X axis) were added to the MarSCoDe spectra (the black verticals along the X axis), and then the peak intensity of ChemCam spectral wavelength (the red points) were linearly interpolated based on MarSCoDe peak profile. Then the wavelength and associated intensity which were not present in ChemCam spectra (the black verticals and the black points) were removed from the MarSCoDe spectra after interpolation. Then a new dataset retained most information of raw MarSCoDe spectra in the same data format as ChemCam spectra was obtained. The peak shapes have also been retained. However, as shown in Figure 3b, the peak top might be lost in some situation, named "peak top loss", because the wavelength at peak top of two kinds of LIBS data occasionally mismatch, which was only observed in the narrower peaks after interpolation. black verticals and the black points) were removed from the MarSCoDe spectra after interpolation. Then a new dataset retained most information of raw MarSCoDe spectra in the same data format as ChemCam spectra was obtained. The peak shapes have also been retained. However, as shown in Figure 3b, the peak top might be lost in some situation, named "peak top loss", because the wavelength at peak top of two kinds of LIBS data occasionally mismatch, which was only observed in the narrower peaks after interpolation.

(4) Step 4 Energy Density Correction
Laser focusing energy density on the surface of target is one of the main influencing factors for the variation of intensity and SNR of LIBS emission lines [39]. The energy density will change with various distances from target during MarSCoDe shooting. The laser energy of MarSCoDe and ChemCam is also different. Therefore, the energy density correction is necessary. Norite, the same LIBS flight calibration samples of MarSCoDe and ChemCam, is the best candidate for energy density correction between two instruments. To achieve the goal, the ratios of MarSCoDe (3 spectra, after above 3 steps processing) and ChemCam (41 spectra) in situ spectra of norite are calculated. Then, the average of these spectral ratios (a total of 123), a single vector containing pixel-by-pixel correction factors, was applied to the MarSCoDe data for energy correction [19].
After all the above four steps, the MarSCoDe spectra were transformed into Chem-Cam data format. Now, the most important question is how to evaluate the effectiveness of our method. MarSCoDe LIBS data recorded on Mars do not have the element abundance values to assess our approach which also need to be quantitatively analyzed. Under this situation, the SDU-LIBS spectra that have the definite main element abundance values were used to test the effectiveness of our proposed approach which were transformed into ChemCam data format using the same 4-step approach as the MarSCoDe.

(4) Step 4 Energy Density Correction
Laser focusing energy density on the surface of target is one of the main influencing factors for the variation of intensity and SNR of LIBS emission lines [39]. The energy density will change with various distances from target during MarSCoDe shooting. The laser energy of MarSCoDe and ChemCam is also different. Therefore, the energy density correction is necessary. Norite, the same LIBS flight calibration samples of MarSCoDe and ChemCam, is the best candidate for energy density correction between two instruments. To achieve the goal, the ratios of MarSCoDe (3 spectra, after above 3 steps processing) and ChemCam (41 spectra) in situ spectra of norite are calculated. Then, the average of these spectral ratios (a total of 123), a single vector containing pixel-by-pixel correction factors, was applied to the MarSCoDe data for energy correction [19].
After all the above four steps, the MarSCoDe spectra were transformed into ChemCam data format. Now, the most important question is how to evaluate the effectiveness of our method. MarSCoDe LIBS data recorded on Mars do not have the element abundance values to assess our approach which also need to be quantitatively analyzed. Under this situation, the SDU-LIBS spectra that have the definite main element abundance values were used to test the effectiveness of our proposed approach which were transformed into ChemCam data format using the same 4-step approach as the MarSCoDe.

Partial Least Squares (PLS) Regression Model
After the above transformation, the LIBS spectra of MarSCoDe can quantitatively analyze the major elements abundance (SiO 2 , TiO 2 , Al 2 O 3 , FeO T , MgO, CaO, Na 2 O, K 2 O) by ChemCam analysis model. All spectra of ChemCam calibration set were transformed by z-score normalization [40] before building the model. As shown in Table 3, among those data, 80 spectra were randomly selected from the dataset as test set and the rest were used to train this model. Earth-to-Mars correction has also been applied to ChemCam spectra before PLS modeling and testing because of the slight differences between the LIBS data which were recorded on Mars and in the laboratory (under simulated Mars conditions) due the unexpected impacts (such as variations of environmental or the working status) on Mars during LIBS shooting [18,19]. The RMSE (root mean squared error) was used as the criterion for PCs (PLS latent variable) selection to ensure the accuracy of the model [17,40]. In Section 3, the reliability was tested by comparing the derived values of the SDU-LIBS spectra from the PLS model with actual chemical composition value.

Reliability Validation of Spectral Transformation Approach
Before quantitively analyzing MarSCoDe data, the effectiveness of 4 step transformation approach (Figure 2) was tested by SDU-LIBS spectra.
The working spectral bands and their step size of ChemCam, MarSCoDe, and SDU-LIBS system are listed in Table 4. In UV and VIS bands, the step size of SDU-LIBS were similar as MarSCoDe and bigger than ChemCam. All of three instruments have similar step size in VNIR band. Therefore, the 4-step transformation approach can also be applied to the SDU-LIBS spectra. The SDU-LIBS spectra were treated as test set to verify the effectiveness of the transformation approach. As described above, more peak top loss will be observed in some narrow peaks after interpolation performance. Moreover, the degree of information loss before and after interpolation performance of MarSCoDe will be tested and discussed. The selected peaks listed in Table 2 rather than full spectrum were used to comparatively analyze the degree of information loss. Then, the PLS model described in Section 2.3 was used to derive the abundance of major elements of transformed SDU-LIBS data and compared with their actual values. Table 4. The step size of 3 LIBS systems in various working spectral bands.

Instrument
Step Size (nm)

Effect of Spectral Transformation Approach
As an example, the original spectrum of sample NO.1 in Table 1 before transformation and the spectrum after transformation is shown in Figure 4 to demonstrate the effect of our proposed spectral transformation approach.

Effect of Spectral Transformation Approach
As an example, the original spectrum of sample NO.1 in Table 1 before transformation and the spectrum after transformation is shown in Figure 4 to demonstrate the effect of our proposed spectral transformation approach.  Table 2 were used to recalibrate the wavelength of SDU-LIBS spectra. The result is shown in Table 5. As shown in Figure 4 and Table 5, most of the peaks have calibrated to the desired position.  Table 2 were used to recalibrate the wavelength of SDU-LIBS spectra. The result is shown in Table 5. As shown in Figure 4 and Table 5, most of the peaks have calibrated to the desired position. The differences in intensity of LIBS data before and after the 4-step transformation were compared in Figure 4, which are derived from the step of interpolation because the peak shape and the spectral information will be lost (as mentioned in step 3) by interpolation. In this case, the above differences were assessed by the RMSE of original intensity and the relative intensity of O at 777 nm [41]. The RMSE results are shown in Table 6, the smaller RMSE values indicated the smaller information loss. The RMSE values are about 0.01 for original peak intensity and from 0.04 to 0.08 for relative intensity of O. The latter are a little bigger because of the greater peak top loss. Both RMSE of SDU-LIBS and ChemCam in Table 6 show that original SDU-LIBS spectra are little different in response to ChemCam calibration spectra as mentioned above. Table 6. The RMSE comparison of before and after interpolation (column A is intensity RMSE between after and before interpolation; column B is relative intensity RMSE between after and before interpolation; column C is intensity RMSE between SDU-LIBS before interpolation and ChemCam; column D is relative intensity RMSE between SDU-LIBS before interpolation and ChemCam).

No.
A Similar to MarSCoDe, SDU-LIBS data were corrected by using the energy correction parameter matrix which was calculated by randomly choosing one sample from 11 samples (listed in Table 1). As analyzed above, the energy density correction parameter matrix of SDU-LIBS has lower precision than MarSCoDe and ChemCam because of different response of various LIBS system. It will influence the quantitative analysis accuracy.
After Earth-to-Mars correction, SDU-LIBS data were used to PLS model for quantitative analysis. Similarly, the actual values of FeO T in ChemCam standard database are also mostly distributed in the smaller and narrower range of 0-10 wt.%, which causes an overfitted model and a lower R 2 of FeO T . In addition, the RMSEC and RMSEP of these models in Figure 5 are similar, and mostly points are tightly distributed around the 1:1 line. Therefore, the PLS model built in this study is available.

Effect of PLS Model Prediction
For SDU-LIBS spectra after 4-step transformation, their derived values of PLS model are shown in Figure 6. As analyzed above, the elements abundance derived values are influenced by the low precision energy density correction parameter matrix (detailed in Section 3.1), especially the elements whose emission lines local at narrow range such as SiO 2 (RMSEP 12.33), MgO (RMSEP 10.44), and Al 2 O 3 (RMSEP 6.22), the R 2 of the latter two elements is even negative. The reason is that the less and narrow located emission lines reduced the energy density correction effect. As shown in Equation (1), the intensity of emission line in near infrared band will be enlarged after the intensity unit transformation. The emission lines of K in spectra of SDU-LIBS are mainly located in VNIR band, respectively, whose peak intensity had inevitably been influenced by intensity unit transformation. Therefore, the derived values of K 2 O are higher than actual values. The emission lines of Ti, Fe, Ca, Na are distributed in a wider range which has little influence on the information loss due to interpolation and the energy density correction. The RMSEP of SiO 2 , TiO 2 , CaO, Na 2    The derived values from PLS model for SDU-LIBS spectra after 4-step spectral transformation approach is acceptable, which implied that most information is retained after the transformation. The quantitative analysis results of MarSCoDe data should be much better than SDU-LIBS model because of the better energy density parameter and more information retention.

Application of Our Approach to MarSCoDe LIBS Data
After the effect of this 4-step approach was tested by SDU-LIBS data, this approach was applied to the MarSCoDe spectra. Similarly, RMSE was used to contrastively analyze the MarSCoDe spectra data after and before interpolation (shown in Table 7). Expect for a few outliers, the RMSE of intensity and relative intensity is stable in the range of 0.01 to 1.5 which is reasonable for MarSCoDe data. Table 7. The RMSE of MarSCoDe before and after interpolation. (a) The RMSE of in situ Martian spectra. (b) The RMSE of on-board calibration target spectra (column A is intensity RMSE between the MarSCoDe spectra of before and after interpolation; column B is relative intensity RMSE between the MarSCoDe spectra of before and after interpolation; column C is intensity RMSE between the MarSCoDe spectra of before and after interpolation; column D is relative intensity RMSE between the MarSCoDe spectra of before and after interpolation).  For energy density correction of MarSCoDe data, one on-board calibration target of Zhurong rover, norite, was used to get correction parameter. The norite spectra recorded by MarSCoDe and ChemCam are much similar, so the energy density correction parameter is more effective. The spectral transformation effect of MarSCoDe spectra is better than that of SDU-LIBS spectra. Table 8 shows the derived value average of norite against its actual value, which is reasonable within the error except a lower FeO T and a higher K 2 O abundant. This can be explained that the emission lines of Fe and K in the in situ spectra of MarSCoDe are mainly located in UV band and VNIR band, respectively, and their peak intensity had been more influenced by intensity unit transformation. The CaO-derived value average of the gypsum and olivine on-board calibration target are higher than other on-board calibration targets. The situation corresponds to the actual situation. Table 9 is the elemental abundancederived values (average derived values of 26 in situ spectra) at Tianwen-1 landing site acquired by MarSCoDe and GRS (Gamma Ray Spectrometer) on board the 2001 Mars Odyssey Mission [42]. The SiO 2 -derived values of MarSCoDe are consistent with GRS and the difference of derived values FeO T and K 2 O due to the change of influence of intensity unit transformation. Another explanation is that GRS probes the average of a wide range of regions and MarSCoDe detects the value of a point in situ, thus the existence of difference is inevitable. The elemental abundance derived values of MarSCoDe spectra using our transformation approach is relatively accurate and reasonable except FeO T .   Figure 7 is the total alkali-silica plot of the MarSCoDe in situ spectra derived values. Because of the higher derived value of K 2 O, the truly total alkali axis coordinates might be lower~1 unit than that shown in Figure 7. Silica derived value lied between Martian soil average and Martian upper crust, and tends to the SiO 2 abundance of Martian upper crust. The silica-derived value average is 49.09 which indicated that Tianwen-1 landing is a kind of basic rock which fill in the in situ detection gap from basic rock to intermediate rock, and may account for the transition component from the Mars upper crust to finegrained components. As mentioned above, the range of alkali derived values from 1 to 4 and in fact the average is about from 2 to 3. It suggests that alkali of Tianwen-1 landing site is similar as other landing sites in Figure 7 and indicates the near alkali level of the northern lowlands. The dataset of MarSCoDe distribute under the grey curve (the alkalinesubalkaline boundary) which implied MarSCoDe targets may be subalkaline composition different from results of Opportunity, Spirit, and Curiosity exploration.

Discussion
Because of the higher derived value of K2O, the truly total alkali axis coordinates might be lower ~1 unit than that shown in Figure 7. Silica derived value lied between Martian soil average and Martian upper crust, and tends to the SiO2 abundance of Martian upper crust. The silica-derived value average is 49.09 which indicated that Tianwen-1 landing is a kind of basic rock which fill in the in situ detection gap from basic rock to intermediate rock, and may account for the transition component from the Mars upper crust to finegrained components. As mentioned above, the range of alkali derived values from 1 to 4 and in fact the average is about from 2 to 3. It suggests that alkali of Tianwen-1 landing site is similar as other landing sites in Figure 7 and indicates the near alkali level of the northern lowlands. The dataset of MarSCoDe distribute under the grey curve (the alkaline-subalkaline boundary) which implied MarSCoDe targets may be subalkaline composition different from results of Opportunity, Spirit, and Curiosity exploration. Figure 7. Total alkali-silica plot of the results of the quantitative analysis from MarSCoDe in situ spectra (modified from [43]). The grey curve indicates the alkaline-subalkaline boundary [44]. The Martian soil average from [45]. The hexagon represents the Mars upper crust composition and composition averages of landing sites. Data from landers/rovers are used only as a reference to the composition of Martian regolith/rock from small-scale detection of the landing site. Figure 7. Total alkali-silica plot of the results of the quantitative analysis from MarSCoDe in situ spectra (modified from [43]). The grey curve indicates the alkaline-subalkaline boundary [44]. The Martian soil average from [45]. The hexagon represents the Mars upper crust composition and composition averages of landing sites. Data from landers/rovers are used only as a reference to the composition of Martian regolith/rock from small-scale detection of the landing site.

Conclusions
A spectral transformation approach was proposed for data cross-calibration between different LIBS systems such as, ChemCam, MarSCoDe, and SDU-LIBS data which include 4-step transformation. Our method enables the spectra of MarSCoDe to be analyzed using the ChemCam PLS model. SDU-LIBS spectra were used to validate the reliability of this approach and training a PLS model based on ChemCam database. The final results of validation show our 4-step approach can retain the original spectral information to a large level and obtain reliable PLS quantitative analysis results.
Based on the results of the quantitative analysis from MarSCoDe data, our result shows that Tianwen-1 landing site was located in a region with basic rock. The alkali of Tianwen-1 landing site is similar to other landing sites. Analysis of Tianwen-1 landing site from MarSCoDe LIBS spectra will deepen our understanding of Martian geology. Further analysis in the future is required for obtaining more detailed information.
The effect of this 4-step approach can be optimized by compensating the emission lines of specific elements and adjusting the correction parameter matrix in the future. The derived values will be more robust with more accurate PLS model. The approach can transform LIBS spectra of different formats to similar ones. This approach can combine specified format LIBS spectra such as MarSCoDe spectra with other LIBS database to analyze and find out more information in the future.