Large-domain, spatially contiguous hydrological and land surface models are important tools for managing our water supplies. Hydrological information on the continental or global scale is needed to handle new emerging international and global water management challenges, which include topics like water allocation in international, national, and large river basins, operational flood forecasting services, global water security or the influence of climate extremes on water resources (Archfield et al., 2015). These applications are particularly challenging in areas without hydrologic measurements, which includes a majority of basins worldwide that are effectively ungauged (Hrachowitz et al., 2013). This results in the need for further development in large-domain hydrological modeling to simulate water fluxes and states in both gauged and ungauged basins in different climates in a spatially consistent manner (Rakovec et al., 2019).In 1982, Jim Dooge stated that "the parameterization of hydrologic processes to the grid-scale of general circulation models is a problem that has not been tackled, let alone solved" (Dooge, 1982) and shortly after that Leavesley et al. (1983) concluded that optimization of distributed parameters of hydrological models is an "ill-posed" problem due to the large number of degrees of freedom. Since then, model parameterization is still one of the major unsolved problems in hydrology (Blöschl et al., 2019). One way to potentially solve this problem is to relate hydrological model parameters/structures to landscape properties (e.g., K.