Ecohydrological partitioning strongly influences the provision of water resources in terms of quantity, quality, location, and timing; thus, fundamentally affecting water-related ecosystem services (Brauman, 2015). In larger scale catchments (>100 km 2 ), ecohydrological functioning often exhibits strong spatio-temporal variability due to small-scale process heterogeneity and non-linearity, both of which result in often marked catchment-wide contrasts in landscape organization (Hunt et al., 2021;Skøien et al., 2003;Tetzlaff et al., 2011). Process-based models have been widely adopted for capturing such catchment functional heterogeneity, as well as its responses to environmental changes (Thirel et al., 2015). Efforts of the catchment modeling community have been progressively improving the physical realism of process conceptualization (Beven, 2002;Clark et al., 2017;Tetzlaff et al., 2008). As part of this trend, stable isotopes of water (i.e., 𝐴𝐴 𝐴𝐴 2 H and 𝐴𝐴 𝐴𝐴 18 O ), as ideal tracers, have received more attention and been increasingly integrated into process-based modeling . As complementary to the use of more traditional hydrometric data (Fenicia et al., 2008), isotopic tracer information has provided more direct evidence of the velocity of water particle transport and allowed inference of catchment-scale water age distributions (