Over the past decade, there has been appreciable progress towards modeling the water, energy, and carbon cycles at field scales (10–100 m) over continental to global extents in Earth system models (ESMs). One such approach, named HydroBlocks, accomplishes this task while maintaining computational efficiency via Hydrologic Response Units (HRUs), more commonly known as “tiles” in ESMs. In HydroBlocks, these HRUs are learned via a hierarchical clustering approach from available global high-resolution environmental data. However, until now there has yet to be a river routing approach that is able to leverage HydroBlocks' approach to modeling field-scale heterogeneity; bridging this gap will make it possible to more formally include riparian zone dynamics, irrigation from surface water, and interactive floodplains in the model. This paper introduces a novel dynamic river routing scheme in HydroBlocks that is intertwined with the modeled field-scale land surface heterogeneity. Each macroscale polygon (a generalization of the concept of macroscale grid cell) is assigned its own fine-scale river network that is derived from very high resolution (∼ 30 m) digital elevation models (DEMs); the inlet–outlet reaches of a domain’s macroscale polygons are then linked to assemble a full domain’s river network. The river dynamics are solved at the reach-level via the kinematic wave assumption of the Saint-Venant equations. Finally, a two-way coupling between each HRU and its corresponding fine-scale river reaches is established. To implement and test the novel approach, a 1.0° bounding box surrounding the Atmospheric Radiation and Measurement (ARM) Southern Great Plains (SGP) site in northern Oklahoma (United States) is used. The results show:
The resulting routing scheme provides an effective and efficient path forward to enable a two-way coupling between the high-resolution river networks and state-of-the-art tiling schemes in ESMs.