AbstractThe continually evolving large ice sheets present in the Northern Hemisphere during the last glacial cycle caused significant changes to river pathways both through directly blocking rivers and through glacial isostatic adjustment. Studies have shown these river pathway changes had a significant impact on the ocean circulation through changing the pattern of freshwater discharge into the oceans. A coupled Earth system model (ESM) simulation of the last glacial cycle thus requires a hydrological discharge model that uses a set of river pathways that evolve with Earth's changing orography while being able to reproduce the known present-day river network given the present-day orography. Here, we present a method for dynamically modelling river pathways that meets such requirements by applying predefined corrections to an evolving fine-scale orography (accounting for the changing ice sheets and isostatic rebound) each time the river directions are recalculated. The corrected orography thus produced is then used to create a set of fine-scale river pathways and these are then upscaled to a coarser scale on which an existing present-day hydrological discharge model within the JSBACH land surface model simulates the river flow. Tests show that this procedure reproduces the known present-day river network to a sufficient degree of accuracy and is able to simulate plausible paleo-river networks. It has also been shown this procedure can be run successfully multiple times as part of a transient coupled climate model simulation.