Distribution and Morphologies of Transverse Aeolian Ridges in ExoMars 2020 Rover Landing Site

Aeolian processes are believed to play a major role in the landscape evolution of Mars. Investigations on Martian aeolian landforms such as ripples, transverse aeolian ridges (TARs), and dunes, and aeolian sediment flux measurements are important to enhance our understanding of past and present wind regimes, the ongoing dust cycle, landscape evolution, and geochemistry. These aeolian bedforms are often comprised of loose sand and sharply undulating topography and thus pose a threat to mobility and maneuvers of Mars rovers. Here we present a first-hand account of the distribution, morphologies, and morphometrics of TARs in Oxia Planum, the recently selected ExoMars 2020 Rover landing site. The gridded mapping was performed for contiguous stretches of TARs within all the landing ellipses using 57 sub-meter high resolution imaging science experiment (HiRISE) scenes. We also provide the morphological descriptions for all types of TARs present within the landing ellipses. We use HiRISE digital terrain models (DTMs) along with the images to derive morphometric information for TARs in Oxia Planum. In general, the average areal TAR coverage was found to be 5.4% (±4.9% standard deviation), increasing from west to east within the landing ellipses. We report the average TAR morphometrics in the form of crest–ridge width (131.1 ± 106.2 m), down-wind TAR length (17.6 ± 10.1 m), wavelength (37.3 ± 11.6 m), plan view aspect ratio (7.1 ± 2.3), inter-bedform spacing (2.1 ± 1.1), slope (10.6◦ ± 6.1◦), predominant orientations (NE-SW and E-W), and height (1.2± 0.8 m). While simple TARs are predominant, we report other TAR morphologies such as forked TAR, wavy TAR with associated smaller secondary ripples, barchan-like TAR, networked TAR, and mini-TARs from the region. Our results can help in planning the rover traverses in terms of both safe passage and scientific returns favoring aeolian research, particularly improving our understanding of TARs.


Introduction
Aeolian processes have played a significant role in landscape evolution and development of surficial features on planets with considerable atmospheres such as Earth, Mars, and Venus [1]. For Mars, the substantial influx of orbiter, lander, and rover imagery and geophysical data with continuously improving resolutions in the past two decades have largely facilitated our understanding of the relict and active aeolian processes and landforms, proving aeolian processes to be one of the dominant surface-modifying agents on Mars [2]. Initially the coarser resolution images of Mariner 9 [3] and Viking Orbiters [4] during the 1970s revealed large, low albedo dune masses or large dark dunes (LDDs). The advent of higher-resolution remote sensors, such as the Mars Orbiter Camera (MOC) [5], High Resolution Stereo Camera [6], Context (CTX) Camera [7], and the high resolution imaging science experiment (HiRISE) [8], during the 1990s and 2000s have significantly aided Mars research [9] and made it possible to better study smaller aeolian landforms such as wind ripples, granule ripples, dust devils, yardangs, ventifacts, and transverse aeolian ridges (TARs). Even though the atmosphere of Mars is considerably thinner (i.e.,~1% of Earth's), the abundance of such landforms on Mars demonstrates strong atmosphere-lithosphere coupling over a large spatiotemporal domain [10]. Thus, these aeolian features and underlying interactions can tell us more about the climate and landscape evolutions on Mars and their presence highlights the regions with past or ongoing sediment transport [11]. In a majority of the cases, we are largely uncertain about the sediment volumes, their origin and traverse, extent, and their composition [11]. Due to the combination of these unknowns and the improvement in data resolutions and simulations, Mars aeolian research has been gaining momentum in the recent years. Several papers [12][13][14][15] over the past decades have reported on manual or semi-automated mapping of features like ripples and sand dunes on Mars at global and regional scales to help our understanding.
TARs are one of the predominant aeolian landforms found in all kinds of geologic settings in equatorial and mid-latitudes (mostly between +35 • and −55 • ) [11,16] on Mars. TARs have decameter-scale wavelengths, intermediate between the small centimeter-scale wind ripples and the large kilometer-scale dunes [17]. TARs were first named by Bourke et al. [18] and the precise information on their origin and sediment source, formation mechanisms, composition, distribution, activity, age, and superposition relationships are yet to be confirmed [19]. However, several interpretations [2,20] envisage them as granule ripples under a surface layer of mm-sized particles, too coarse to be moved by the present-day Martian winds [17]. There are several reported possible terrestrial analogs for Martian TARs, such as the megaripples in Argentina [21] and in the Lut Desert of Iran [19,22], megaripples at Wau-an-Namus, Libya [23], and reversing sand dunes in Idaho, USA [24]. TARs have generally brighter or similar albedo as the surrounding terrain and can occur singularly or as large contiguous stretches [10,22]. The albedo of bright TARs is as high as >0.21, nearly double the albedo of dark dunes [10]. They normally display simple, forked, sinuous, networked, and barchan-like patterns [10]. TARs are landforms with symmetric cross-sectional profiles normally displaying slopes of <28 • and maximum height-to-wavelength ratios of~0.13-0.17 [17,20]. However, in a recent global survey of their morphologies, TARs have been reported to display widely varying morphologies and uncorrelated distribution with latitudes and altitudes [17]. This survey [17] reports average TAR heights to be 2.6 m with a standard deviation of 1.4 m, while the tallest reported TAR reached a maximum height of 14 m. Another independent survey [25] reports typical dimensions of 88.5 m, 17.3 m, 1.3 m, and 25.8 m for length, width, height, and wavelength, respectively, of the mapped TARs. Geissler and Wilgus [17] further reported height-to-wavelength ratios to be weakly anti-correlated with TAR wavelengths, and no correlation between wavelengths and the local slope. Usually TARs in a region display similar orientations over large areas indicating their adaptation to the large-scale wind patterns irrespective of the local topography [26]. Granular megaripples are often considered as a terrestrial analog for TARs.
Oxia Planum, the selected ExoMars 2020 Rover landing site (Figures 1 and 2), was preferred over Mawrth Vallis [27] in the final ExoMars 2020 site selection meeting of European scientists and engineers held at the UK's National Space Centre in Leicester on 8 and 9 November 2018 due to the greater safety margins that the terrain in Oxia Planum offers for a successful landing and the rover's mobility. For ensuring the best scientific returns from the planned ExoMars mission, it is necessary to scrutinize the terrain and landforms in Oxia Planum at high resolutions to better guide future traverses and mobility of the rover from the point of landing. As explained above, aeolian landforms, particularly TARs, are not only interesting science targets [28], but they also pose a threat to landing and rover mobility [29] in terms of sand/dust depositions, direct physical damage to the lander or rover chassis and wheels, sinkage into the loose unconsolidated aeolian material [30,31], enhanced slippage [30,32], and permanent entrapment [33]. Sullivan et al. [28,29] provide detailed experimental background of Mars Exploration Rover (MER) Opportunity and MER Spirit where these rovers encountered TARs repeatedly in Meridiani Planum and Gusev Crater, respectively, because TARs were unavoidable during several traverse segments; but the scientific teams wisely used this as an opportunity to closely investigate TARs and other aeolian features. Similarly, the Mars Science Laboratory (MSL) Curiosity Rover during sols 526-538 of its operation at Dingo Gap Bedform in Gale Crater encountered a large 1 m-high megaripple and photographic, as well as spectroscopic cameras, onboard the rover were simultaneously used to study an intended wheel scuff [29]. The ExoMars 2020 Lander and Rover are also equipped with several instruments of this kind which can provide in-situ information to validate remote observations and help our understanding of TARs. For example, the lander has meteorological instruments such as HabitAbility: Brine Irradiation and Temperature (HABIT) and the meteorological package to provide meteorological context of the visual, spectroscopic, and radar observations made by the rover instruments, such as the Panoramic Camera (PanCam), Infrared Spectrometer for ExoMars (ISEM), Close-UP Imager (CLUPI), Water Ice and Subsurface Deposit Observation on Mars (WISDOM), Mars Multispectral Imager for Subsurface Studies (Ma_MISS), the visible and infrared imaging spectrometer MicrOmega, and the Raman spectrometer (RLS). In addition, the ExoMars 2020 is equipped with a drill that embeds Ma_MISS and can reach the subsurface to a maximum depth of 2 m. These instruments make the ExoMars 2020 rover an ideal platform to conduct aeolian research and promote a region-wide high-resolution mapping of TARs and other aeolian landforms. Thus, in the present study, we aim to map all the expected landing ellipses (3-σ) [34] covering a total area of 3211.07 km 2 in the Oxia Planum region for TARs ( Figure 2). We start with three main objectives to present an account of region-wide (1) distribution, (2) morphologies, and (3) morphometries of TARs in Oxia Planum. We first provide a brief introduction of Oxia Planum as the study site, followed by the methods used to reach the above-mentioned objectives. Our results can help plan the rover traverses in terms of both safe passage and scientific returns favoring aeolian research, particularly improving our understanding of TARs. Similar studies should be applied in the future to other potential regions of interest for landing missions on Mars, as TARs are not only relevant for mobility, but also for the safety assessment of the Entry, Descent, and Landing (EDL) phase.

Study Area
The Oxia Planum landing site (16 • to 19 • N, 337 to 332 • W) is located at the edge of an ancient basin called Chryse Planitia and is characterized by~4 billion-year-old sedimentary clay deposits. This region is marked by ancient highland cratered terrains highly eroded towards the highland-lowland boundary and falls in a wide basin at the outlet of the Cogoon Vallis System between Mawrth Vallis and Ares Vallis [35]. The geomorphological setting of the region and the composition of clay deposits indicate a possibility of them being lagoon or deltaic deposits near an ancient Martian ocean [35]. These layered exposures in the Oxia Planum region are rich in iron and magnesium phyllosilicates and possibly represent the south-western extension of the clay-rich Mawrth Vallis deposits [35].
For all the mapping and derived inferences within our study, the study area henceforth will refer to the region shown in Figure 1a [36] at 2 km/pixel resolution, the minimum and maximum TI values were 329.44 Jm −2 K −1 s −1/2 and 454.65 Jm −2 K −1 s −1/2 , respectively, with a mean TI of 395.52 Jm −2 K −1 s −1/2 and a standard deviation of 30.4 Jm −2 K −1 s −1/2 . The TI was found to be generally increasing from west-to-east ( Figure 1c) within the study area. We further analyzed TES-derived rock abundance [37] data of~7 km/pixel resolution for the study area and the % rock abundance was found to be varying between 10.8% and

Materials and Methods
The adopted methodology can be divided and described separately for the aforementioned three objectives as follows:

TAR Distribution
In addition to extensive stretches of TARs, Oxia Planum is also rich in so-called "mini-TARs", named so as these are about an order of magnitude smaller than most TARs; too small for their topography to be discerned using HiRISE digital elevation models (DEMs) [31]. However, these mini-TARs have length/height ratios similar to TARs suggesting them to be a possible variant of TARs [31]. In the present study, we performed the mapping at 1:2000 scale and did not consider these mini-TARs as they were not mappable at that scale. We have mapped TARs showing traditional decameter-scale wavelengths and all the results and discussion henceforth concern only them.
For the purpose of rover surface mobility, it is not the one-off TARs which are significantly dangerous as they can be avoided using even the rover mast cameras. The main problem arises when the rover encounters the contiguous stretches of TARs which can be actually hazardous for the rover mobility. Such closely located TARs are predominant in Oxia Planum. However, in our mapping, we identified and mapped both the singular TARs and the contiguous stretches of TARs in the selected grids. We used "Create Fishnet" tool of ArcGIS software and provided TES TI raster (Figure 1c) as the "template extent" dataset to derive corresponding grids of 2 × 2 km 2 dimensions ( Figure 2). As mentioned above, our study area is the region within the ellipse boundaries and therefore, to have a representative, balanced analysis of this wide region, we selected 66 grids uniformly distributed throughout the ellipses (Figure 2) and covering a total area of 264 km 2 for HiRISE-based TAR mapping. In total, 57 sub-meter resolution HiRISE scenes ( Figure 2) were

Materials and Methods
The adopted methodology can be divided and described separately for the aforementioned three objectives as follows:

TAR Distribution
In addition to extensive stretches of TARs, Oxia Planum is also rich in so-called "mini-TARs", named so as these are about an order of magnitude smaller than most TARs; too small for their topography to be discerned using HiRISE digital elevation models (DEMs) [31]. However, these mini-TARs have length/height ratios similar to TARs suggesting them to be a possible variant of TARs [31]. In the present study, we performed the mapping at 1:2000 scale and did not consider these mini-TARs as they were not mappable at that scale. We have mapped TARs showing traditional decameter-scale wavelengths and all the results and discussion henceforth concern only them.
For the purpose of rover surface mobility, it is not the one-off TARs which are significantly dangerous as they can be avoided using even the rover mast cameras. The main problem arises when the rover encounters the contiguous stretches of TARs which can be actually hazardous for the rover mobility. Such closely located TARs are predominant in Oxia Planum. However, in our mapping, we identified and mapped both the singular TARs and the contiguous stretches of TARs in the selected grids. We used "Create Fishnet" tool of ArcGIS software and provided TES TI raster (Figure 1c) as the "template extent" dataset to derive corresponding grids of 2 × 2 km 2 dimensions ( Figure 2). As mentioned above, our study area is the region within the ellipse boundaries and therefore, to have a representative, balanced analysis of this wide region, we selected 66 grids uniformly distributed throughout the ellipses ( Figure 2) and covering a total area of 264 km 2 for HiRISE-based TAR mapping.
In total, 57 sub-meter resolution HiRISE scenes ( Figure 2) were selected to form the composite that covered nearly all of the study area and assisted in TAR mapping. We used the physical account well-described in the literature [5,10,38] to identify TARs in the images. TARs characteristically form in groups of tens to thousands of individual units and can easily be characterized by their small, narrow, and transverse dimensions, and elongated form [10]. Slip faces are typically resolvable only in high-resolution HiRISE images, and the opposite slopes appear to be nearly symmetric [5,10,38].
Remote Sens. 2019, 11, x FOR PEER REVIEW 5 of 18 selected to form the composite that covered nearly all of the study area and assisted in TAR mapping. We used the physical account well-described in the literature [5,10,38] to identify TARs in the images. TARs characteristically form in groups of tens to thousands of individual units and can easily be characterized by their small, narrow, and transverse dimensions, and elongated form [10]. Slip faces are typically resolvable only in high-resolution HiRISE images, and the opposite slopes appear to be nearly symmetric [5,10,38]. Although our mapping grids were well-distributed across the study area and since TARs are known to usually form in groups [10], the expected results were sufficient to provide an acceptable first-hand assessment and information about TAR distribution across the region. However, the facts that TARs typically form in contiguous stretches and that such spatially correlated TARs were abundant in the study area also suggested that we further perform spatial interpolation using the information from the mapped grids to derive a spatially continuous TAR distribution map. For this, we selected 53 of the mapped grids (yellow grids in Figure 2) and used "Kriging" interpolator [39] in ArcGIS. Kriging is one of the best geostatistical interpolation techniques in soil science and geology [40] and is well-explored to interpolate Martian regolith parameters [27] as it considers both the degree of variation between known data points and the distance between them while estimating values in unknown areas. Kriging is based on the assumption that the distance and direction between sample points reflect a spatial correlation that can be further exploited to predict the variation in the surface [40]. We subsequently used the remaining uniformly distributed 13 grids Although our mapping grids were well-distributed across the study area and since TARs are known to usually form in groups [10], the expected results were sufficient to provide an acceptable first-hand assessment and information about TAR distribution across the region. However, the facts that TARs typically form in contiguous stretches and that such spatially correlated TARs were abundant in the study area also suggested that we further perform spatial interpolation using the information from the mapped grids to derive a spatially continuous TAR distribution map. For this, we selected 53 of the mapped grids (yellow grids in Figure 2) and used "Kriging" interpolator [39] in ArcGIS. Kriging is one of the best geostatistical interpolation techniques in soil science and geology [40] and is well-explored to interpolate Martian regolith parameters [27] as it considers both the degree of variation between known data points and the distance between them while estimating values in unknown areas. Kriging is based on the assumption that the distance and direction between sample points reflect a spatial correlation that can be further exploited to predict the variation in the surface [40]. We subsequently used the remaining uniformly distributed 13 grids (cyan grids in Figure 2) to evaluate the performance of the interpolation process. Upon quantifying the spatial uncertainty in the interpolation, all 66 grids were employed for deriving the final TAR distribution map for the study area.

Morphology
The morphological interpretations for the TARs in the study area were based on visually identifying the first-order crest-ridge plan view morphological forms (simple, forked, wavy, barchan-like, and networked) described in the literature [10]. We characterized the most common morphologies of TARs in the region.

Morphometry
We used all the three available HiRISE digital terrain models (DTMs) for the study area: (1) DTEEC_003195_1985_002694_1985_L01, (2) DTEEC_009880_1985_009735_1985_L01, (3) DTEEC_036925_1985_037558_1985_L01 and the corresponding images to perform high-resolution (~1 m/pixel) morphometry of several TARs. The HiRISE DTMs produced by the HiRISE team are reliable with the reported vertical precision of a few tenths of a meter as a result of a stereo matching precision of about 0.2 of a pixel [41].
We further searched for the inherent noise within the used HiRISE DTMs [42]. There are mainly four such artifacts as reported by the HiRISE team [43] and one needs to look into the terrain shaded relief to detect these: (1) boxes, (2) faceted areas, (3) charged coupled device (CCD) seams, and (4) manually interpolated areas. We scrutinized the three DTMs for such artifacts using terrain shaded relief views and were able to identify mainly two types of artifacts; faceted areas and CCD seams ( Figure 3). In fact, these artifacts cover large portions of the above-mentioned DEMs. We ensured that we sample and analyze only those TARs which are away from any of these artifacts. The "Hillshade" tool of ArcGIS was used to derive the shaded relief. Thus, we selected 150 TARs which were well-discernible in HiRISE DTMs, away from any artifacts, and displayed different sizes.
In the same way as Hugenholtz et al. [25], we also restricted our morphometric analysis to only 'simple' TARs as defined by Balme et al. [10]. These simple TARs have well-defined morphologies with symmetrical slopes, a single crestline, clearly marked edges, and low sinuosity [25], making the measurements easier and more accurate. We selected and digitized the crestlines of only those TARs which were clearly distinguishable in both the HiRISE image and its corresponding shaded relief (e.g., TARs indicated by the red arrow in Figure 3) to ensure that their topography was resolved in the DTM.
As defined by Balme et al. [10], TARs can be characterized by five morphometric parameters: (1) the crest-ridge width (W), (2) the down-wind TAR length (L), (3) the ridge-to-ridge mean spacing (λ), (4) plan view aspect ratio (a), where a = W/L, and (5) inter-bedform spacing (s), where s = λ/L (see Figure 4 for an illustration of these parameters in one specific example). L, as we take it, is the widest part of the TAR and should actually be recorded as the maximum "down-wind TAR length". Similarly, λ, as we measure it, is the average of distances between the consecutive TAR crestlines at approximately the middle of the TARs. In addition to calculating these parameters for the selected TARs, we also derived along-crestline and across-crestline elevation profiles for several of them. However, we do not limit ourselves only to these morphometric parameters, but we also derive slope, orientation, and heights (h) for the selected TAR samples. The Pythagoras' theorem was used for the h calculations and for the ease of understanding along with the transect analysis results, the description of height estimation is provided in Section 4.3. TAR orientations (N-S, NE-SW, NW-SE, or E-W) were analyzed based on the azimuthal position of TAR crest-ridge width.

DTEEC_003195_1985_002694_1985_L01,
(2) DTEEC_009880_1985_009735_1985_L01, (3) DTEEC_036925_1985_037558_1985_L01 and the corresponding images to perform high-resolution (~1 m/pixel) morphometry of several TARs. The HiRISE DTMs produced by the HiRISE team are reliable with the reported vertical precision of a few tenths of a meter as a result of a stereo matching precision of about 0.2 of a pixel [41].  Figure 2). The yellow ellipse highlights 'simple' TARs as defined by Balme et al. [10]. These simple TARs have well-defined morphologies with symmetrical slopes, a single crestline, clearly marked edges, and low sinuosity [25], making the measurements easier and more accurate. We selected and digitized the crestlines of only those TARs which were clearly distinguishable in both the HiRISE image and its corresponding shaded relief (e.g., TARs indicated by the red arrow in Figure 3) to ensure that their topography was resolved in the DTM.  Figure 4 for an illustration of these parameters in one specific example). L, as we take it, is the widest part of the TAR and should actually be recorded as the maximum "down-wind TAR length". Similarly, λ, as we measure it, is the average of distances between the consecutive TAR crestlines at

Results and Discussion
The following sections summarize the key results of our study within the predefined objectives and we discuss the results in view of the previous studies.

TAR Distribution
We present the gridded TAR mapping results in Figure 5. The % TAR area within each grid varies between 0.3% and 46.3% with a mean value of 5.78% and a standard deviation of 9.13% (Figure 5a). The result of Kriging interpolation using the first 53 mapped grids ( Figure 2) is shown in Figure 5b. Two distinct nodes of high TAR density in the eastern margin of the study area are clearly observed (Figure 5b). The % TAR area varies between 0% and 44.354% with a mean value of 5.66% and a standard deviation of 6.97% (Figure 5b). The remaining 13 mapped grids ( Figure 2) were selected as the validation dataset to evaluate the performance of the interpolator. Owing to the uniform distribution of the used grids for interpolation, the observed uncertainty (i.e., actual % area, interpolated % area) in interpolation was significantly less with a mean value of 3.27% and a standard deviation of 12.47% of the grid area. We further used these evaluation grids to interpolate the uncertainty across the study area and observed that for~98.4% of the terrain the uncertainty was within ±5% of the grid area (Figure 5c). We observed the highest uncertainty of~44% for one of the eastern-most grids and uncertainty values in the range of 5%-20% for the surrounding grids (Figure 5c). With such good interpolation results, we repeated the spatial interpolation using all the 66 mapped grids to derive the final TAR distribution map for the study area (Figure 5d). The % TAR area in the final TAR distribution map varied between 0.3% and 46.3% and the mean % TAR area for the region was 5.39% with a standard deviation of 4.87%. Remote Sens. 2019, 11, x FOR PEER REVIEW 9 of 18 As visible in Figure 5d, the TAR density and areal coverage is significantly higher towards the eastern margin of the study area as compared to the rest of the region. In fact, 63.5% of the region shows % TAR coverage of <5% and 80.32% of the study area displays <10% TAR coverage. There are very few conference abstracts [34,44] and one MPhil thesis [45], which have presented some results on partial gridded mapping of TARs in the region and our results are in good agreement with those studies. For example, Bridges et al. [44] presented TAR mapping results for 180 km 2 (~5.6% of the study area) at 1:2000 scale (the same scale as selected by us) and reported an average TAR coverage of 4.1% with a very high (≤43%) TAR coverage towards the eastern margin. Our results for the entire study area (3211.07 km 2 ) also show a mean TAR coverage of 5.39% with highest mean areal coverage of ~43% in the eastern margin. Bridges et al. [34] presented the updated results of TAR mapping for the study area and reported a mean TAR coverage of 4.4% for the sampled and mapped grids. Henson [45], for several more of the grids in the study area, reported an average TAR coverage of 4.9% ± 5.9%. Henson [45] also reported the highest TAR coverage of 42.7% (close to our estimate of 46.3%) in the east and the lowest overall coverage of 0.4% (close to 0.3% as we report). Our TAR estimates are in good agreement with the ones provided by other researchers which further suggest the effectiveness of our mapping and spatial interpolation and provide a holistic spatially continuous distribution of TAR across the region. In summary, Figure 5d is a reliable estimate of the distribution of TARs within the ExoMars landing ellipses in Oxia Planum, which can be used for surface operation planning and EDL assessment. It can also be used to compare with global and meso-scalar atmospheric circulation models to understand the role of present-day aeolian processes on the distribution and orientation of these landforms. We further tried to compare the TI values with TAR areal coverage for the respective grids, and as reported by previous studies for other regions on Mars [17,26], we too did not observe any direct correlation between the two. On an average, for <20% of the grid coverage by TARs, the observed TI was 393.94 Jm -2 K -1 s -1/2 while for As visible in Figure 5d, the TAR density and areal coverage is significantly higher towards the eastern margin of the study area as compared to the rest of the region. In fact, 63.5% of the region shows % TAR coverage of <5% and 80.32% of the study area displays <10% TAR coverage. There are very few conference abstracts [34,44] and one MPhil thesis [45], which have presented some results on partial gridded mapping of TARs in the region and our results are in good agreement with those studies. For example, Bridges et al. [44] presented TAR mapping results for 180 km 2 (~5.6% of the study area) at 1:2000 scale (the same scale as selected by us) and reported an average TAR coverage of 4.1% with a very high (≤43%) TAR coverage towards the eastern margin. Our results for the entire study area (3211.07 km 2 ) also show a mean TAR coverage of 5.39% with highest mean areal coverage of~43% in the eastern margin. Bridges et al. [34] presented the updated results of TAR mapping for the study area and reported a mean TAR coverage of 4.4% for the sampled and mapped grids. Henson [45], for several more of the grids in the study area, reported an average TAR coverage of 4.9% ± 5.9%. Henson [45] also reported the highest TAR coverage of 42.7% (close to our estimate of 46.3%) in the east and the lowest overall coverage of 0.4% (close to 0.3% as we report). Our TAR estimates are in good agreement with the ones provided by other researchers which further suggest the effectiveness of our mapping and spatial interpolation and provide a holistic spatially continuous distribution of TAR across the region. In summary, Figure 5d is a reliable estimate of the distribution of TARs within the ExoMars landing ellipses in Oxia Planum, which can be used for surface operation planning and EDL assessment. It can also be used to compare with global and meso-scalar atmospheric circulation models to understand the role of present-day aeolian processes on the distribution and orientation of these landforms. We further tried to compare the TI values with TAR areal coverage for the respective grids, and as reported by previous studies for other regions on Mars [17,26], we too did not observe any direct correlation between the two. On an average, for <20% of the grid coverage by TARs, the observed TI was 393.94 Jm −2 K −1 s −1/2 while for >20% TAR coverage, the observed TI was 424.35 Jm −2 K −1 s −1/2 . One interesting observation that we make is regarding the comparison between the TAR coverage ( Figure 5d) and the rock abundance maps (Figure 1d). While, the eastern margin of the study area is abundant in TARs with lesser rock abundance, it is the middle-to-western part of the study area, which while having lesser TAR coverage, displays higher than average rock abundance with maximum values reaching up to 35.4% (Figure 1d). This means that the entire Oxia Planum region within the landing ellipses can pose serious challenges during the ExoMars 2020 landing and maneuvers as the presence of dense TARs and high rock abundance have nearly mutually exclusive spatial distribution and dangerous levels throughout the study area.

Morphology
While performing the morphological observations, we were able to identify most of the first-order crest-ridge plan view morphologies of TARs as explained by Balme et al. [10], including the mini-TARs [31] in the study area ( Figure 6). We have not focused on several other reported morphologies such as "raked pattern TARs" [46] and "dust drifts" [16]. The most prevalent morphology is simple TAR, abundantly present throughout the study area. This is followed by the forked TARs, which although significantly lesser in number than the simple TARs, are as abundant on the crater floors as outside the craters. The wavy, barchan-like, and networked morphologies are rare in the study area and are mainly present on the crater floors. Sometimes all or a combination of these five morphologies can co-exist within the same TAR deposit [10]. For example, in the study area, the forked and wavy morphologies with associated smaller secondary ripples are found to be co-existing mostly with simple morphologies. The barchan-like morphologies are found encompassing other morphologies within several of the crater floors. The mini-TARs are abundant in two regions, the eastern margin as well as the south-west region, with usually NE-SW orientations. Figure 6f can be used for comparison purposes between TARs and mini-TARs, which are about an order of magnitude smaller [31]. Balme et al. [31] have explained these mini-TARs to be potentially more hazardous for a rover as they are located in dense clusters.
We also observed several examples of TAR transitions as reported by Balme et al. [10] for the Meridiani region. Figure 7a highlights one case where barchan-like TARs coalesce to form joint and simple TARs through merging of their crest ridges. These barchan-like TARs usually display crest-ridge widths (W) of~20-50 m and down-wind TAR lengths (L) of~10-15 m. The simple TARs which they form can display average W of~500 m. However, this transition is rare in the study area. Spatial transitions from simple, forked, and wavy TARs into networked TARs are relatively common in the study area (Figure 7b-d). In Figure 7b, we observe gradual transitions from simple to forked to wavy, and ultimately to networked morphologies. Figure 7c presents a case where secondary ridge crests forming approximately perpendicular to the largest simple TAR crest ridges display a transition to networked TARs. These TARs with secondary ridge crests are also called "feathered" TARs by Bourke et al. [18]. The W of these feathered TARs in the study area usually vary between 200-400 m. Figure 7d presents an interesting case of transition from simple to networked morphology as a possible combination of downslope winds and steep topography. Figure 7d also shows an elongated depression as a result of two merged craters. The shadows on the left suggest the topographic gradient with higher elevations in the smaller crater and gradual deepening towards the larger crater. Moving downslope from the border of the merged craters to the center of the larger crater, the simple TARs start showing transition and change in orientation along the marked red arrow (Figure 7d). We also observed several examples of TAR transitions as reported by Balme et al. [10] for the Meridiani region. Figure 7a highlights one case where barchan-like TARs coalesce to form joint and simple TARs through merging of their crest ridges. These barchan-like TARs usually display crest-ridge widths (W) of ~20-50 m and down-wind TAR lengths (L) of ~10-15 m. The simple TARs which they form can display average W of ~500 m. However, this transition is rare in the study area. Spatial transitions from simple, forked, and wavy TARs into networked TARs are relatively common in the study area (Figure 7b-d). In Figure 7b, we observe gradual transitions from simple to forked to wavy, and ultimately to networked morphologies. Figure 7c presents a case where secondary ridge crests forming approximately perpendicular to the largest simple TAR crest ridges display a transition to networked TARs. These TARs with secondary ridge crests are also called shows an elongated depression as a result of two merged craters. The shadows on the left suggest the topographic gradient with higher elevations in the smaller crater and gradual deepening towards the larger crater. Moving downslope from the border of the merged craters to the center of the larger crater, the simple TARs start showing transition and change in orientation along the marked red arrow (Figure 7d).

Morphometry
To our knowledge, we did not find any research paper which exclusively provides published accounts of TAR morphometrics in Oxia Planum. In addition to the conventional TAR morphometric parameters (W, L, λ, a, and s) as defined by Balme et al. [10], we also derived slope, orientation, and heights for the selected TAR samples. However, before we discuss the results, we would reiterate that the selected TAR samples for this analysis are representative of only simple TAR morphologies in the study area. Figure 8 presents two of the transects which were used to derive height (h) of 50 TARs with clearly distinct crests and troughs. The Pythagoras' theorem was used for the h calculations (Figure 8a). BF (i.e., L), BD, and EC in Figure 8a can directly be derived using the planimetric measurements on the HiRISE images. Similarly, elevation differencing

Morphometry
To our knowledge, we did not find any research paper which exclusively provides published accounts of TAR morphometrics in Oxia Planum. In addition to the conventional TAR morphometric parameters (W, L, λ, a, and s) as defined by Balme et al. [10], we also derived slope, orientation, and heights for the selected TAR samples. However, before we discuss the results, we would reiterate that the selected TAR samples for this analysis are representative of only simple TAR morphologies in the study area. Figure 8 presents two of the transects which were used to derive height (h) of 50 TARs with clearly distinct crests and troughs. The Pythagoras' theorem was used for the h calculations ( Figure 8a). BF (i.e., L), BD, and EC in Figure 8a can directly be derived using the planimetric measurements on the HiRISE images. Similarly, elevation differencing between the two points on the DTMs can directly provide AE, CF and AD. Afterwards, using the Equations (1) to (6), we can derive h. We have five right-angled triangles in Figure 8a to derive the required height (h); ∆CFB, ∆ADB, ∆AEC, ∆AGB, and ∆AGC. From these right-angled triangles, we can derive the following equations: By equating Equations (4) and (5) we get; Now, we can calculate the equivalent of BG from Equation (6) using values of AB, AC, and BC derived from Equations (1)-(3). When we put this value of BG in Equation (5), we can derive the value of h. AB = (BD 2 + AD 2 ) 1/2 , h = (AC 2 -GC 2 ) 1/2 = (AC 2 -(BC -BG) 2 ) 1/2 , h = (AB 2 -BG 2 ) 1/2 , By equating Equations (4) and (5) we get; (AB 2 -BG 2 ) = (AC 2 -(BC -BG) 2 ) => AB 2 -AC 2 + BC 2 = 2BC.BG => BG = (AB 2 -AC 2 + BC 2 ) / 2BC, Now, we can calculate the equivalent of BG from Equation (6) using values of AB, AC, and BC derived from Equations (1), (2), and (3). When we put this value of BG in Equation (5), we can derive the value of h.  The study area is full of small TARs with λ = 2-3 m or L = 2-3 m [45]. Even in our samples which were chosen based on their clear visibility in HiRISE DTM shaded reliefs, we were able to characterize a TAR with L = 3.7 m (Table 1). Similarly, our other estimations regarding various TAR morphometrics are in good agreement with the ones reported for other regions by previous researchers.  [11] reported a range of~3-10 • , while Shockey and Zimbelman [47] and Zimbelman [20] reported the average slopes to be 14.5 • and~15 • , respectively. All these dimensions and morphometrics are of a similar order of magnitude and nature as observed by us in Oxia Planum (Table 1) for a localized and smaller TAR population. We found a wide variability of TAR dimensions across the three DTMs within the study area with up to a factor of~7-10 difference in TAR morphometrics as also reported by Geissler and Wilgus [17]. Overall, we report the average TAR morphometrics (±standard deviation) of W = 131.  Table 1). Usually the longer TARs are wider too in most of the instances, but we did not observe any prevalent or statistically significant relationships of the W and L with the h dimensions of TARs in the study area as also evident from the transect profiles in Figure 8.

Conclusions
The present research provides first-hand detailed information about the distribution, morphologies, and morphometrics of TARs in the ExoMars 2020 landing site within Oxia Planum. TARs are widespread aeolian landforms on contemporary Mars, but we are yet to be certain of their formation mechanism, age and temporal evolution, composition and geophysics, and their relationship with global sedimentation and climatic evolution. Thus, although the presence of TARs in the study area with a particularly high density in the entire eastern margin of the landing ellipses poses a significant potential threat to the ExoMars 2020 mission, it also equips us with a unique opportunity to investigate these enigmatic landforms at unprecedented spectral resolutions and corroborate the results with precise meteorological parameters. The rover platform is equipped with infrared spectrometers and multispectral imagers, and the provision of a drill with an embedded multispectral camera provides us the opportunity to perform profiling starting from the tougher top layer of TARs and reaching to the supposedly loosely-packed subsurface. Considering that the drill can reach up to a maximum depth of 2 m and the average TAR height reported by us is 1.25 m, we can actually aim towards the complete profiling of a given TAR in terms of its composition and sedimentation. Although the TARs in Oxia Planum are similar in their dimensions to the TARs reported from other regions, we identify a strong NE-SW and E-W orientation for a majority of them. The presence of dense populations of mini-TARs in the region will be an additional challenge for the rover maneuvers as the maximum height of these mini-TARs can reach up to 25 cm [31] and they cannot be characterized using HiRISE DTMs. One of the most useful products of future studies could be a map of all TAR crestlines in the study area with their orientation represented as a rose plot. In the present work, the crestline analysis was part of our morphometric analyses which we performed only for selected 150 TARs, which were clearly represented in the DTMs. Performing such a directional analysis for all the TAR crestlines in the study area will take a considerable effort. Additionally, given the highly variable morphologies in this vast area, we will have to first fix some set of rules to perform such a directional analysis and it can be considered as one of the future scopes of the present work. In summary, our results can prove to be useful for ExoMars surface operation planning and EDL assessment. Furthermore, such studies are of interest to compare with global and meso-scalar atmospheric circulation models to understand the role of present-day aeolian processes on the distribution and orientation of TARs. Future observations from the ExoMars Rover and surface platform will provide ground truth to these remote sensing and wind modeling analyses.