A Global High-Resolution Hydrological Model to Simulate the Dynamics of Surface Liquid Reservoirs: Application on Mars

Abstract. Surface runoff shapes planetary landscapes, but global hydrological models often lack the resolution and flexibility to simulate dynamic surface water bodies beyond Earth. Recent studies of Mars have revealed abundant geological and mineralogical evidence for past surface water, including valley networks, crater lakes, deltas and possible ocean margins dating from late Noachian to early Hesperian times. These features suggest that early Mars experienced periods allowing liquid water stability, runoff and sediment transport. To investigate where surface water could accumulate and how it may have been redistributed, we developed a global high-resolution (km-scale) surface hydrological model. The model uses a pre-computed hydrological database that maps topographic depressions, their spillover points, hierarchical connections between basins, and lake volume-area-elevation relationships. This database approach greatly accelerates simulations by avoiding repeated geomorphic processing. The model dynamically forms, grows, merges and dries lakes and putative seas without prescribing fixed coastlines, by transferring water volumes between depressions according to their storage capacities and overflow rules. We explore model behavior over the present-day Mars' topography measured by MOLA (Mars Orbiter Laser Altimeter) topography for a range of evaporation rates and total water inventories expressed as Global Equivalent Layer (GEL). Simulations are iterated to steady state under the assumption that precipitation balances evaporation plus overflow. The model outputs the extent and depth of surface water bodies and identifies main drainage pathways using overflow fluxes as runoff indicators. Results show a transition toward a contiguous northern ocean between low (1–10 m) GEL values and increasing concentration of water in northern lowlands and major impact basins at higher GEL. We discuss the model's limitations, including its dependence on topography and the absence of subsurface flows, and propose future improvements. This framework provides a quantitative tool to link preserved geomorphology with plausible past hydrological states. Future work will couple the model with a 3D global climate model into a Planetary Evolution Model (PEM) to study transient water redistribution and climate-hydrology feedbacks.