Constraints on martian lobate debris apron evolution and rheology from numerical modeling of ice flow

Radar observations in the Deuteronilus Mensae region by Mars Reconnaissance Orbiter have constrained the thickness and dust concentration found within mid-latitude ice deposits, providing an opportunity to more accurately estimate the rheology of ice responsible for the formation of lobate debris aprons based on their apparent age of similar to 100 Myr. We developed a numerical model simulating ice flow under martian conditions using results from ice deformation experiments, theory of ice grain growth based on terrestrial ice cores, and observational constraints from radar profiles and laser altimetry. By varying the ice grain size, the ice temperature, the subsurface slope, and the initial ice volume we determine the combination of parameters that best reproduce the observed LDA lengths and thicknesses over a period of time comparable to the apparent ages of LDA surfaces (90-300 Myr). We find that an ice temperature of 205 K, an ice grain size of 5 mm, and a flat subsurface slope give reasonable ages for many LDAs in the northern mid-latitudes of Mars. Assuming that the ice grain size is limited by the grain boundary pinning effect of incorporated dust, these results limit the dust volume concentration to less than 4%. However, assuming all LDAs were emplaced by a single event, we find that there is no single combination of grain size, temperature, and subsurface slope which can give realistic ages for all LDAs, suggesting that some or all of these variables are spatially heterogeneous. Based on our model we conclude that the majority of northern mid-latitude LDAs are composed of clean (<= 4 vol%), coarse (>= 1 mm) grained ice, but regional differences in either the amount of dust mixed in with the ice, or in the presence of a basal slope below the LDA ice must be invoked. Alternatively, the ice temperature and/or timing of ice deposition may vary significantly between different mid-latitude regions. Either eventuality can be tested with future observations. (C) 2011 Elsevier Inc. All rights reserved.