Climate change on Mars and the formation of gullies, lobate debris aprons, and softened craters.

摘要

Recent data acquired from spacecraft missions has bettered our understanding of the nature and distribution of ice- and water-related features formed during recent periods of climate change on Mars. Numerical modeling of physical processes constrained by these newly acquired observations is an important tool with which hypotheses relating to the Martian climate can be tested. This work describes the development and implementation of a set of these models focused on the formation of a few young, ice- and water-related features. The subjects of this research are gullies, lobate debris aprons (LDAs), and craters with subdued topography known as "softened" craters. Flow of liquid water and ice over and/or within the Martian surface has been invoked in the formation of these features. Quantifying processes such as fluvial erosion and ice deformation using laboratory experiments is a Rosetta stone with which we can read into the climate history of Mars that is written on its surface.;We test the hypothesis that sediment transport on gully slopes occurs via fluvial transport processes by developing a numerical sediment transport model based on steep flume experiments performed by Smart [1984]. At 20° slopes, channels 1 m deep by 8m wide and 0.1 m deep by 3 m wide transport a sediment volume equal to the alcove volume of 6 x 105 m3 in 10 hours and 40 days, respectively, under constant flow conditions. Snowpack melting cannot produce the water discharge rates necessary for fluvial sediment transport, unless long-term (kyr) storage of the resulting meltwater occurs. If these volumes of water are discharged as groundwater, the required aquifer thicknesses and aquifer drawdown lengths would be unrealistically large for a single discharge event. More plausibly, the water volume required by the fluvial transport model could be discharged in ∼ 10 episodes for an aquifer 30 m thick, with a recurrence interval similar to that of Martian obliquity cycles (∼0.1 My).;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 the ice within lobate debris aprons based on their apparent age of 100 My. 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. 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 (≤ 4vol%), 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.;The presence of an extensive ice-rich layer in the near subsurface of the Martian regolith could result in viscous creep responsible for softening of craters at middle and high latitudes. The temperature of ground ice will vary both temporally and spatially due, respectively, to changes in Mars' obliquity and due to the slope effect on the effective angle of insolation. Numerical simulations of viscous creep indicate that these temperature variations cause the pole-facing slopes of craters to be systematically steeper than those of equator-facing slopes. Crater slopes should be most asymmetric between 25° and 400 latitude, depending on the thickness of the creeping layer. On the basis of the lack of any systematic slope asymmetry observed in the craters, we can place a conservative upper limit of ∼ 150 m on the thickness of the ice-rich creeping layer assuming a volumetric dust content of ≤ 70% and an exponentially decreasing regolith porosity with depth. If the creeping layer contains relatively clean ice, then the thickness of ice-rich material is limited to ∼ 100 m or less. The observations also suggest that the thickness of this creeping layer is reduced by ∼ 30% toward the equator. These results imply a global ice-rich regolith water volume of ∼ 10 7 km3, comparable to that proposed for a modest-sized northern plains ocean.