A global map of gullied hillslopes on Mars

The distribution of gullies on Mars has been used to support different modes of formation involving CO2, H2O or entirely dry processes. We mapped the extent of gullied hillslope using a global mosaic of ConTeXt camera (CTX) at 6 m/pixel and High Resolution Imaging Science Experiment (HiRISE) browse products at ~3 m/pixel. Our new globalcatalog of gullied hillslopes allows for quantitative, area-based analysis of gullies both at global and regional scale, which was not possible with point-based or more localized previous datasets. We analyzed aspect and gradient information derived from MOLA topography for every mapped gullied hillslope. We have confirmed the trends identified in previous works, which are that the mean gradient of gullied hillslopes is lower than 30°, the angle of repose of dry material on Mars, and that the preferred aspect of gullied hillslopes is latitude-dependent, with a shift between pole-facing to equator-facing happening around 40° in both hemispheres. We established a hierarchy of factors that can explain the global distribution of gullied hillslopes: Latitude is the strictest control because gullies are constrained between 26°-83°S and 28°-76°N. The availability of steep slopes is a strong control over gully distribution, and we observed that gullied hillslopes aspect reflects the orientation of large regional reliefs, which implies that hillslope aspect is less of a factor than the hillslope gradient in explaining gully distribution. We analyzed the correlation of gully distribution with a subsurface water ice consistency map, and we found that gullies are preferentially distributed on areas with some ground ice rather than no ice at all. We identified for the first time a bias toward an east-facing trend for southern hemisphere gullies poleward of ~40° of latitude. We also linked our gullied hillslopes dataset with a global crater morphometry dataset: 63% of the gullied area on Mars is found in craters, and we find that the area covered by gullies in a given crater is positively correlated with the depth and diameter ratio of this crater. Our observations point toward a complex interplay between local slope, insolation, and thermal properties of the substrate to explain the spatial distribution of the current population of gullies. Finally, we have assessed the very high resolution (sub-meter) imagery coverage of gullies to be 36% and that the repeated coverage reaches 21% of gullies, hence this would need expanding to better understand the spatial trends in active gullies.