Martian Ice Revealed by Modeling of Simple Terraced Crater Formation

Arcadia Planitia, a region in the northern midlatitudes of Mars, displays an uncommonly high abundance of simple craters with a concentric morphology, which is indicative of layering beneath the surface. Radar measurements suggest that the near surface layers could be made of excess water ice. In this study, we select two of these impact structures of similar size (D_c ~ 500 m), model their formation through iSALE shock physics code, and investigate the dependence of the final crater morphology on the material model parameters (cohesion and friction coefficient). Our parameter study shows that the intact and damaged cohesions of the nonporous ice play a fundamental role to obtain a good fit between our models and the topographic profiles taken from the digital terrain models in terms of crater diameter, crater wall inclination, and depth and size of the upper terrace. The central pit shape is instead controlled by the damaged friction coefficient of the basaltic crust, but it is mainly affected by projectile density and speed. Our results confirm that two layers of relatively pure water ice, each with different rheology and porosity, can explain the unique double-terraced morphology of impact craters in Arcadia Planitia. The low values of cohesion we find for the ice might point to snowfall as emplacement mechanism in the region. The different thicknesses of the ice layers in the two crater areas seem to suggest variations in ice deposition and/or evolution history across Arcadia Planitia.