Simulation and optimization of Ma_MISS surveys
Authors: Rossi Lorenzo , Ferrari Marco , De Sanctis Maria Cristina , Frigeri Alessandro , De Angelis Simone , Costa Nicole , Altieri Francesca , Formisano Michelangelo , Ammannito Eleonora 
Published: September 2022

&lt;p&gt;&lt;strong&gt;Introduction&amp;#160;&lt;/strong&gt;&lt;br /&gt;We developed a simple model to simulate Ma_MISS data acquisition and visualization and to support the development of acquisition strategies.&lt;/p&gt;
&lt;p&gt;&lt;strong&gt;The Ma_MISS instrument&amp;#160;&lt;br /&gt;&lt;/strong&gt;Ma_MISS (Mars Multispectral Imager for Subsurface Studies) is the miniaturized VNIR spectrometer embedded in the drill system of the ExoMars rover&lt;sup&gt;[1]&lt;/sup&gt;. Ma_MISS will perform spectral reflectance measurements inside holes drilled into the surface of Mars up to a depth of 2 m. It will characterize the mineralogy and stratigraphy of the borehole walls and will allow the in-situ study of the subsurface environment, before any samples are collected and extracted. Taking advantage of the finely controllable drill tool motion, Ma_MISS can map the borehole walls. In addition to single point acquisitions, Ma_MISS can also be programmed to acquire &amp;#8220;Ring&amp;#8221; and &amp;#8220;Column&amp;#8221; scans, where the spectrometric acquisitions are interleaved with step rotations or translations of the drill tool.&lt;/p&gt;
&lt;p&gt;&lt;strong&gt;Borehole stratigraphy model&amp;#160;&lt;/strong&gt;&lt;br /&gt;To easily simulate Ma_MISS surveys, we first developed a simple model of a borehole. As further described in [2], six different simulant samples were prepared. Reflectance spectra were acquired on each of the samples with DAVIS-MOT, the new Ma_MISS laboratory model&lt;sup&gt;[3]&lt;/sup&gt;. The reflectance spectra of these samples are shown in Figure 1. Six different regions were drawn as polylines on a 2D projection of the cylindrical borehole walls. Out of the six simulants, a different one was assigned to each region. As shown in Figure 2, this model was made to represents a 250 mm deep section of a borehole, exhibiting a varied stratigraphy with 5 different layers and a vertical vein.&lt;/p&gt;
&lt;p&gt;&lt;img src=&quot;&quot; alt=&quot;&quot; width=&quot;500&quot; height=&quot;393&quot; /&gt; &lt;br /&gt;Figure 1: Reflectance spectra of samples.&lt;/p&gt;
&lt;p&gt;&lt;img src=&quot;&quot; alt=&quot;&quot; width=&quot;480&quot; height=&quot;380&quot; /&gt;&lt;/p&gt;
&lt;p&gt;Figure 2: Borehole stratigraphy model&lt;/p&gt;
&lt;p&gt;&lt;strong&gt;Survey simulation&amp;#160;&lt;br /&gt;&lt;/strong&gt;The borehole model can be used to quickly simulate a Ma_MISS survey. Once the acquisition strategy of the survey is defined, the coordinates of the points where spectrometric acquisitions are to be simulated are computed. Small random errors can be added to the coordinates to simulate the accuracy of the drill tool actuators and sensors. For each point of this set, we generate a spectrum that simulates the real acquisition. Each simulated spectrum is generated from the measured spectra of one of the six simulants, based on its region in the stratigraphy model, with the addition of random shifts, slopes, and noise to simulate out-of-focus acquisitions and instrumental noise.&lt;/p&gt;
&lt;p&gt;&lt;strong&gt;Development of acquisition strategies&amp;#160;&lt;br /&gt;&lt;/strong&gt;The simulated data can be used to develop acquisition strategies, aimed at maximizing the information gained while limiting the time required to complete a survey of a drilled borehole. The time required to complete the survey depends on both the spectrometer acquisition settings and the drill tool motion. For example, the duration of a drill movement depends on angular and vertical step sizes (hereafter referred to as &amp;#916;&amp;#952; and &amp;#916;z respectively), with vertical translations generally taking longer than rotations. Multiple combination of these parameters can be tried out on the same stratigraphy model to evaluate which strategy can provide the best compromise between spatial resolution and time required.&lt;/p&gt;
&lt;p&gt;&lt;strong&gt;Simulated data examples&lt;br /&gt;&lt;/strong&gt;For the 4&lt;sup&gt;th&lt;/sup&gt; RSOWG (Rover Science Operations Working Group) Simulation, held in November 2021, we provided data from a simulated Ma_MISS acquisition survey. The simulated Ma_MISS survey (&amp;#8220;survey A&amp;#8221;) consisted of:&lt;br /&gt;&amp;#160;-6 rings of 360 points each, with &amp;#916;&amp;#952;=1&amp;#176; and displaced vertically by &amp;#160;&amp;#916;z=50mm;&lt;br /&gt;&amp;#160;-8 columns made up of 126 points each, taken with &amp;#916;&amp;#952;=45&amp;#176; and a vertical step of &amp;#916;z=2mm.&lt;br /&gt;This amounts to a total of 3168 individual spectra. A view of the simulated data is shown in Figure 3 and Figure 5A. This acquisition strategy provides a high angular resolution in the rings with &amp;#916;&amp;#952;=1&amp;#176; (which can be useful to identify small vertical features or individual grains) and a good vertical resolution in the column with &amp;#916;z=2mm (useful to reconstruct the stratigraphy of the subsurface). &lt;br /&gt;Another survey (&amp;#8220;survey B&amp;#8221;) that was simulated with the same model is shown in Figure 4 and Figure 5B-C. In this case the acquisition was made up of 35 rings of 72 points each, with &amp;#916;&amp;#952;=5&amp;#176; and &amp;#916;z=6.25 mm, for a total of 2952 spectra. This strategy offers worse spatial resolution in both the vertical direction and along the circumference but can nonetheless provide a more regular overall view of the borehole stratigraphy.&lt;/p&gt;
&lt;p&gt;&lt;img src=&quot;&quot; alt=&quot;&quot; width=&quot;550&quot; height=&quot;390&quot; /&gt;&lt;br /&gt;Figure 3: Simulated survey A&lt;/p&gt;
&lt;p&gt;&lt;img src=&quot;&quot; alt=&quot;&quot; width=&quot;550&quot; height=&quot;403&quot; /&gt;&lt;/p&gt;
&lt;p&gt;Figure 4: Simulated survey B&lt;/p&gt;
&lt;p&gt;&lt;strong&gt;Data visualization&amp;#160;&lt;br /&gt;&lt;/strong&gt;In Figure 3 and 4, each dot represents a simulated spectrum. In this case the colour of each dot was computed from its spectrum such that it looks somewhat like it would appear to the human eye. The same colours were interpolated across the surface (with a scattered data bicubic interpolation method) to obtain views like those shown in Figure 5A-B. Other useful visualizations of these spectra can be obtained by highlighting other spectral features, e.g., band depths, band ratios or specific IR bands. An example is shown in Figure 5C, where the 1.95 &amp;#181;m band depth is computed on the simulated spectra of survey B and its value is interpolated across the cylindrical surface.&lt;br /&gt;&lt;img src=&quot;&quot; alt=&quot;&quot; width=&quot;450&quot; height=&quot;500&quot; /&gt;&lt;/p&gt;
&lt;p&gt;Figure 5: Visualization of simulated data on the borehole surface&lt;/p&gt;
&lt;p&gt;&lt;strong&gt;Future work&amp;#160;&lt;br /&gt;&lt;/strong&gt;The developed simulation model can be easily extended to other borehole stratigraphy models, with different geometries and spectra. It will also be complemented with tools to estimate the time needed to perform each simulated survey, which will make it a useful tool for the optimization of acquisition strategies. The developed data visualization tools can also be used with real, measured data. We are planning to verify the simulated data on real surveys. In fact, we can perform a complete survey with the Ma_MISS laboratory breadboard on an actual drilled sample. The acquired data will be visualized and used to reconstruct the sample stratigraphy.&lt;/p&gt;
&lt;p&gt;&lt;strong&gt;Acknowledgments&amp;#160;&lt;br /&gt;&lt;/strong&gt;This work is supported by the Italian Space Agency (ASI) grant ASI-INAF n. 2017-412-H.0.&lt;/p&gt;
&lt;p&gt;&lt;strong&gt;References&lt;br /&gt;&lt;/strong&gt;[1] De Sanctis M.C. et al. (2017): Astrobiology, 17, 6-7&lt;br /&gt;[2] Ferrari M. et al. (2022): LPSC53, Abstract#1339&lt;br /&gt;[3] De Angelis S. et al. (2022): LPSC53, Abstract#1796&lt;/p&gt;