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. One such approach, named HydroBlocks, accomplishes this task while maintaining computational efficiency via sub-grid tiles, or Hydrologic Response Units (HRUs), learned via a hierarchical clustering approach from available global high-resolution environmental data. However, until now, there has yet to be a macroscale river routing approach that is able to leverage HydroBlocks' approach to sub-grid heterogeneity, thus limiting the added value of field-scale land surface modeling in Earth System Models (e.g., riparian zone dynamics, irrigation from surface water, and interactive floodplains). This paper introduces a novel dynamic river routing scheme in HydroBlocks that is intertwined with the modeled field-scale land surface heterogeneity. The primary features of the routing scheme include:
To implement and test the novel approach, a 1.0-degree 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 existing tiling schemes within Earth system models.