** The goal of the Subsurface Water Ice Mapping (SWIM) project is to provide a set of mapping products using existing spacecraft data that delineate subsurface ice in the mid-latitudes of Mars**. We aim to identify and map indicators of possible subsurface ice in each data set and use a combination of all data sets to assess the likelihood of ice being present in shallow (< 5 m depth) and deep (> 5 m depth) zones.

Not all existing ice will be identified, since the current data sets impose limits on lateral and vertical resolution as well as sensing depth (Figure 1). In broad terms, the available data sets examine four zones: 1) the surface itself (i.e., to a few microns depth) for imagery and elevation data, 2) within the upper ~1 m for thermal and neutron spectrometers, 3) within the upper ~5 m for radar reflections from the surface, and 4) between ~15 m and several 100 m depths for radar reflections from subsurface interfaces. Geomorphologic indicators of ice span all four zones. Apart from transient ice after a recent impact, there is no ice exposed at the surface in the areas of interest (0 to 60° latitude).

The interpretation of ice presence from these data sets is to a degree subjective and does not lend itself to precise calculation of probabilities. However, some quantification of our confidence in our identification and mapping of ice is warranted and highly desirable for planning future landing sites that rely on the presence of ice for resources or scientific studies. To that end, we have developed the ** SWIM Equation**, produced in the spirit of the famous

We begin by defining * ice consistency* values C

Term | Data Sets | Ice Depth (m) |
---|---|---|

C_{I} |
All | < 5 and > 5 |

C_{N} |
Consistency of neutron-detected hydrogen with shallow ice | < 1 |

C_{T} |
Consistency of thermal behavior with shallow ice | < 1 |

C_{G} |
Consistency of geomorphology with shallow and deep ice | All |

C_{RS} |
Consistency of radar surface echoes with shallow ice | < 5 |

C_{RD} |
Consistency of radar dielectric properties with deep ice | > 5 |

In assessing the composite ice consistency C_{I} of all data sets, we started with a “hands-off” approach wherein we chose not to apply any weightings or separation by sensing depths during the __first season of SWIM (2018-2019):__

C_{I} = ( C_{N} + C_{T} + C_{G} + C_{RS} + C_{RD} ) / 5

Our northern hemisphere products available since 2019 use this introductory formulation of the SWIM Equation. In our __second season of SWIM (2019-2020)__, we have introduced a new approach wherein we provide separate equations for three depth zones: 0-1 m, which is dominated by neutron and thermal spectrometer data; 1-5 m, which is dominated by shallow geomorphic and radar surface-return data; and >5 m, which is dominated by deeper geomorphic and radar subsurface dielectric data. To enable this approach, we divided the mapped geomorphic landforms into two groups corresponding to ice presence above and below 5 m depths, producing terms C_{GS} and C_{GD}. In addition, we introduce weighting factors s_{Y}, the shallowness of method Y, which is determined by dividing the depth of interest for each equation by the sensing depth of each method. Thus, for the first depth zone, we have:

C_{I}[0-1m] = ( s_{N}C_{N} + s_{T}C_{T} + s_{GS}C_{GS} + s_{RS}C_{RS} ) / (s_{N} + s_{T} + s_{GS} + s_{RS} )

= ( C_{N} + C_{T} + 0.2*C_{GS} + 0.2*C_{RS} ) / 2.4

where:

s_{N} = 1 m / 1 m = 1

s_{T} = 1 m / 1 m = 1

s_{GS} = 1 m / 5 m = 0.2

s_{RS} = 1 m / 5 m = 0.2

For the second depth zone, we have:

C_{I}[1-5m] = ( s_{GS}C_{GS} + s_{RS}C_{RS} + s_{RD}C_{RD} ) / ( s_{GS} + s_{RS} + s_{RD})

= ( C_{GS} + C_{RS} + 0.3*C_{RD} ) / 2.3

where:

s_{GS} = 5 m / 5 m = 1

s_{RS} = 5 m / 5 m = 1

s_{RD} = 5 m / 15 m = 0.3

For the third depth zpne, we have:

C_{I}[>5m] = ( s_{GD}C_{GD} + s_{RD}C_{RD} ) / ( s_{GD} + s_{RD} )

= ( C_{GD} + C_{RD} ) / 2.0

where:

s_{GS} = 5 m / 5 m = 1

s_{RS} = 5 m / 5 m = 1

s_{RD} = 5 m / 15 m = 0.3

This formulation of the SWIM equation is by no means the only approach possible, and other methods may be brought to bear, depending on one’s goals.

* Example 1:* From a thermal perspective (represented by the term C

*Measurement Caveat:* A buried high thermal inertia detection is also consistent with the presence of rock. Hence the need for multiple data sets.

* Example 2:* In the case of radar surface returns (C

*Measurement Caveat:* A low power return is also consistent with the presence of non-ice low density material such as thick (> 1 m) accumulations of dust.

Taking the two examples above, we see that a positive value for either C_{T} or C_{RS} in isolation is not a unique detection of ice by itself. The power behind the SWIM Equation is that it integrates multiple measurements to build confidence in the detection of ice. Within this in mind, we see that for a given locale, a buried high thermal inertia measurement __in addition__ to a low radar power return are inconsistent with the occurrence of either buried rock or thick dust layers, and instead they argue strongly for the presence of ice.