On the non-uniqueness of the hydro-geomorphic responses in a zero-order catchment with respect to soil moisture

Abstract

This study advances mechanistic interpretation of predictability challenges in hydro-geomorphology related to the role of soil moisture spatial variability. Using model formulations describing the physics of overland flow, variably saturated subsurface flow, and erosion and sediment transport, this study explores (1) why a basin with the same mean soil moisture can exhibit distinctly different spatial moisture distributions, (2) whether these varying distributions lead to non-unique hydro-geomorphic responses, and (3) what controls non-uniqueness in relation to the response type. Two sets of numerical experiments are carried out with two physically-based models, HYDRUS and tRIBS+VEGGIE+FEaST, and their outputs are analyzed with respect to pre-storm moisture state. The results demonstrate that distinct spatial moisture distributions for the same mean wetness arise because near-surface soil moisture dynamics exhibit different degrees of coupling with deeper-soil moisture and the process of subsurface drainage. The consequences of such variations are different depending on the type of hydrological response. Specifically, if the predominant runoff response is of infiltration excess type, the degree of non-uniqueness is related to the spatial distribution of near-surface moisture. If runoff is governed by subsurface stormflow, the extent of deep moisture contributing area and its “readiness to drain” determine the response characteristics. Because the processes of erosion and sediment transport superimpose additional controls over factors governing runoff generation and overland flow, non-uniqueness of the geomorphic response can be highly dampened or enhanced. The explanation is sediment composed by multi-size particles can alternate states of mobilization or surface shielding and the transient behavior is inherently intertwined with the availability of mobile particles. We conclude that complex nonlinear dynamics of hydro-geomorphic processes are inherent expressions of physical interactions. As complete knowledge of watershed properties, states, or forcings will always present the ultimate, if ever resolvable, challenge, deterministic predictability will remain handicapped. Coupling of uncertainty quantification methods and space-time physics-based approaches will need to evolve to facilitate mechanistic interpretations and informed practical applications.