Digital elevation models and derived flow directions are crucial to distributed hydrological modeling. As the spatial resolution of models is typically coarser than these data, we need methods to upscale flow direction data while preserving the river structure. We propose the Iterative Hydrography Upscaling (IHU) method and show it outperforms other often-applied methods. We publish the multi-resolution MERIT Hydro IHU hydrography dataset and the algorithm as part of the pyflwdir Python package.