SHARAD Data Analysis with High Resolution Digital Terrain Models
Authors: Desage Léopold , Herique Alain , Kofman Wlodek , Zine Sonia 
Published: September 2022

&lt;p&gt;&lt;strong&gt;Abstract&lt;/strong&gt;&lt;/p&gt;
&lt;p&gt;The SHAllow RADar (SHARAD) is a synthetic aperture radar (SAR) onboard Mars Reconnaissance Orbiter,&amp;#160; in Martian orbit since 2006 [1]. With its 20MHz frequency, it can fathom the first hundreds of meters of the subsurface with a range resolution of about 10m in typical Martian materials. In order to detect subsurface interfaces with radars, we need to eliminate the echoes coming from the surface : as the design of most radars in orbit gives them a large antenna lobe, off-nadir echoes could arrive at the same delay as a subsurface nadir reflection. The classical method to remove this so called &amp;#8220;clutter&amp;#8221; consists in comparing the radar signal to simulations of the surface echoes using Digital Terrain Models (DTMs). Our aim is to study the first tens of meters of the Martian subsurface with SHARAD. To do so, we need high resolution DTMs because in theory, the higher we get in resolution, the more detailed the simulation will be and the clearer the discrimination between surface and subsurface features will be. We will present results on high resolution simulations performed with SPRATS, our coherent simulator. We will show that while being high resolution, models obtained by photogrammetry sometimes contain artifacts that can be misleading for radar data interpretation.&amp;#160;&amp;#160;&lt;/p&gt;
&lt;p&gt;&amp;#160;&lt;/p&gt;
&lt;p&gt;&lt;strong&gt;Simulations with SPRATS : first study with MOLA DTMs&lt;/strong&gt;&lt;/p&gt;
&lt;p&gt;SPRATS is a toolset developed at IPAG that allows to perform both coherent radar simulations of surfaces and 3D SAR processing of them [2], [3]. Those capabilities enable the simulation of the actual signal sensed by SHARAD, with the same processing applied to it, in order to get as close as possible from the instrument results. It &amp;#160;allows for a direct power comparison and thus finer analysis. To begin, we performed simulations with DTMs generated with MOLA, a laser altimeter [4]. While having a relatively low resolution, the nature of the acquisition method give those DTMs a high accuracy and precision, resulting in very low artifacts. To study deep reflectors in areas of relatively low rugosity (i.e. the northern plains [5]), those models are sufficient. But the low resolution is a limiting factor when studying close subsurface, or simply to reproduce surface roughness effects on the radar signal. To improve the simulations, higher resolution&amp;#160; models are necessary.&lt;/p&gt;
&lt;p&gt;&amp;#160;&lt;/p&gt;
&lt;p&gt;&lt;strong&gt;HRSC models, higher resolution but sometime containing artifacts&lt;/strong&gt;&lt;/p&gt;
&lt;p&gt;With a resolution of 50 to 100m, HRSC DTMs [6] yield better results in simulating smaller details. It allows to confirm or discard reflectors identified with MOLA [7]. However, these models are acquired by photogrammetry, a technique that is an estimation of the surface topography, compared to MOLA which is a direct measurement. Photogrammetry introduces artifacts that are not easy to estimate, because they depend on the actual topography. We will show a comparative study of simulations with HRSC and MOLA models on a region of interest located in Terra Cimmeria, following a study made by [8]. The amplitude of the artifacts on the HRSC models is too high to study the first tens of meters of the subsurface with SHARAD. Following the idea of getting as close as possible to the actual SHARAD data, we need models that describe the surface at a resolution better that the radar&amp;#8217;s wavelength.&lt;/p&gt;
&lt;p&gt;&amp;#160;&lt;/p&gt;
&lt;p&gt;&lt;strong&gt;Simulations with models at wavelength-scale resolution (CTX)&lt;/strong&gt;&lt;/p&gt;
&lt;p&gt;We will present a comparative study of simulations using CTX models &amp;#8212; with a resolution of 12m &amp;#8212; and HRSC models. We will also show that the scale of the artifacts on these DTMs being below the SHARAD&amp;#8217;s wavelength, CTX DTMs yield near perfect surface echoes simulation,&amp;#160; allowing for a fine detail comparative analysis of the SHARAD data. However, given their acquisition method [9], CTX DTMs have a relatively poor surface coverage compared to HRSC, so we used photoclinometry with CTX images on lower resolution models [10] to keep the high resolution information. Comparing high resolution simulations using these models to SHARAD data allowed to highlight small-scale artifacts on the CTX DTMs, as they introduce noise in the radargram.&lt;/p&gt;
&lt;p&gt;&amp;#160;&lt;/p&gt;
&lt;p&gt;&lt;strong&gt;Conclusions&lt;/strong&gt;&lt;/p&gt;
&lt;p&gt;This study shows that wavelength-scale or smaller artifacts on DTMs are needed to perform shallow subsurface analysis of SHARAD data.&amp;#160; It also showed that high resolution models acquired by photogrammetry are prone to artifacts, which can perturb the simulated signal. This artifacts issue can prove to be helpful for DTM quality estimation, especially for missions where no laser altimeter is present to validate the altimetry measurements.&lt;/p&gt;
&lt;p&gt;&amp;#160;&lt;/p&gt;
&lt;p&gt;&lt;strong&gt;References&lt;/strong&gt;&lt;/p&gt;
&lt;p&gt;[1] R. Seu &lt;em&gt;et al.&lt;/em&gt;, &amp;#171;&amp;#160;SHARAD sounding radar on the Mars Reconnaissance Orbiter&amp;#160;&amp;#187;, 2007, doi: 10.1029/2006JE002745.&lt;/p&gt;
&lt;p&gt;[2] Y. Berquin, A. Herique, W. Kofman, et E. Heggy, &amp;#171;&amp;#160;Computing low-frequency radar surface echoes for planetary radar using Huygens-Fresnel&amp;#8217;s principle: COMPUTING RADAR SURFACE ECHOES&amp;#160;&amp;#187;, oct. 2015, doi: 10.1002/2015RS005714.&lt;/p&gt;
&lt;p&gt;[3] J.-F. Nouvel, A. Herique, W. Kofman, et A. Safaeinili, &amp;#171;&amp;#160;Radar signal simulation: Surface modeling with the Facet Method: RADAR SIGNAL SIMULATION&amp;#160;&amp;#187; , febr. 2004, doi: 10.1029/2003RS002903.&lt;/p&gt;
&lt;p&gt;[4] M. T. Zuber &lt;em&gt;et al.&lt;/em&gt;, &amp;#171;&amp;#160;The Mars Observer laser altimeter investigation&amp;#160;&amp;#187;, 1992, doi: 10.1029/92JE00341.&lt;/p&gt;
&lt;p&gt;[5] C. M. Stuurman &lt;em&gt;et al.&lt;/em&gt;, &amp;#171;&amp;#160;SHARAD detection and characterization of subsurface water ice deposits in Utopia Planitia, Mars: SHARAD DETECTION OF ICE UTOPIA PLANITIA&amp;#160;&amp;#187;, sept. 2016, doi: 10.1002/2016GL070138.&lt;/p&gt;
&lt;p&gt;[6] G. Neukum et R. Jaumann, &amp;#171;&amp;#160;HRSC: the High Resolution Stereo Camera of Mars Express&amp;#160;&amp;#187;, 2004.&lt;/p&gt;
&lt;p&gt;[7] C. W. Cook &lt;em&gt;et al.&lt;/em&gt;, &amp;#171;&amp;#160;Sparse subsurface radar reflectors in Hellas Planitia, Mars&amp;#160;&amp;#187;, sept. 2020, doi: 10.1016/j.icarus.2020.113847.&lt;/p&gt;
&lt;p&gt;[8] S. Adeli, E. Hauber, G. G. Michael, P. Fawdon, I. B. Smith, et R. Jaumann, &amp;#171;&amp;#160;Geomorphological Evidence of Localized Stagnant Ice Deposits in Terra Cimmeria, Mars&amp;#160;&amp;#187;, juin 2019, doi: 10.1029/2018JE005772.&lt;/p&gt;
&lt;p&gt;[9] M. C. Malin &lt;em&gt;et al.&lt;/em&gt;, &amp;#171;&amp;#160;Context Camera Investigation on board the Mars Reconnaissance Orbiter&amp;#160;&amp;#187;, 2007, doi: 10.1029/2006JE002808.&lt;/p&gt;
&lt;p&gt;[10] S. Doute et C. Jiang, &amp;#171;&amp;#160;Small-Scale Topographical Characterization of the Martian Surface With In-Orbit Imagery&amp;#160;&amp;#187;, janv. 2020, doi: 10.1109/TGRS.2019.2937172.&lt;/p&gt;