Effects of soil heterogeneity on martian ground-ice stability and orbital estimates of ice table depth

Data from the Mars Odyssey Gamma-Ray Spectrometer (GRS) instrument suite and results from numerical simulations of subsurface ground-ice stability have been used to estimate the depth of martian ground-ice. Geographic correlation between these estimates is remarkable; the relative ice table depth distributions also agree well. However, GRS-based estimates of ice table depth are generally deeper than predictions based on ground-ice stability simulations. This discrepancy may be related to heterogeneities in the martian surface such as rocks, dust, and albedo variations. We develop a multi-dimensional numerical model of ground-ice stability in a heterogeneous martian subsurface and use it to place the first quantitative constraints on the response of the ice table to meter-scale heterogeneities. We find that heterogeneities produce significant undulations/topography in the ice table at horizontal length scales of a few meters. Decimeter scale rocks create localized areas of deep ice, producing a vertical depression of 10-60 cm in the ice table over a horizontal range of 1-2 rock radii. Decimeter scale dust lenses produce locally shallow ice; however the magnitude of the vertical deflection of the ice table is small (1-4 cm). The effects of decimeter scale albedo variations on the ice table are nearly negligible, although albedo very weakly enhances the effects of dark rocks and bright dust on the ice table. Additionally, we investigate the role played by rocks in estimates of ice table depth based on orbital data. Surface rocks can account for more than half of the discrepancy between ice table depths inferred from GRS data and those predicted by theoretical ice-stability simulations that utilize thermophysical observations. Our results have considerable relevance to the up-coming Mars Scout Mission, Phoenix, because they indicate that the uncertainty in the ice table depth of a given region is greater than differences between current depth estimates. Likewise, small-scale depth variability due to heterogeneities at the eventual landing site is potentially greater than differences between current depth estimates. (c) 2006 Elsevier Inc. All rights reserved.