Although these changes will profoundly influence groundwater recharge and discharge in mountainous environments (Hayashi, 2019), they have largely been ignored so far (Somers et al., 2019). Since surface water and groundwater resources are closely coupled, an improved understanding of surface water-groundwater interactions is highly relevant for a sustainable water governance as well as for water-dependent ecosystems in mountainous regions (e.g., Holman, 2006;Krause et al., 2014;Schilling et al., 2020).Within the last two decades, studies on river-aquifer exchange dynamics have substantially improved the understanding of the drivers (e.g., river discharge) and controls (e.g., riverbed hydraulic conductivity) of water exchange patterns and their impact on biogeochemical cycling of solutes (e.g., reviews by Boano et al., 2014;Brunner et al., 2017; Lewandowski et al., 2019 and references therein). Particularly, the continued recognition and investigation of riverbed dynamics as key controls on river-aquifer exchange have brought substantial scientific progress in the field of surface water-groundwater interactions (e.g., Mutiti & Levy, 2010;Tang et al., 2018). However, the spatiotemporal dynamics of surface water-groundwater interactions still remain elusive, mainly due to a lack of high-resolution field data (Barthel & Banzhaf, 2016;

Integrated modeling of coupled surface-subsurface flow and ensuing role in diverse Earth system processes is of current research interest to characterize nonlinear rainfall-runoff response and also to understand land surface energy balances, biogeochemical processes, geomorphological dynamics, etc. A growing number of complex models have been developed for water-related research, and many of these are made available to the Earth science community. However, relatively few resources have been made accessible to the potentially large group of Earth science and engineering users. New users have to invest an extraordinary effort to study the models. To provide a stimulating experience focusing on the learning curve of integrated modeling of coupled surface-subsurface flow, we describe use cases of an open source model, the Penn State Integrated Hydrologic Model, PIHM. New users were guided through data processing and model application by reproducing a numerical benchmark problem and a real-world watershed simulation. Specifically, we document the PIHM application and its computational workflow to enable intuitive understanding of coupled surface-subsurface flow processes. In addition, we describe the user experience as important evidence of the significance of reusability. The interaction shows that documentation of data, software, and computational workflow in research papers is a promising method to foster open scientific collaboration and reuse. This study demonstrates how open science practice in research papers would promote the utility of open source software. Addressing such open science practice in publications would promote the utility of journal papers. Further, popularization of such practice will require coordination among research communities, funding agencies, and journals.

Abstract. Essentially all hydrogeological processes are strongly influenced by the subsurface spatial heterogeneity and the temporal variation of environmental conditions, hydraulic properties, and solute concentrations. This spatial and temporal variability needs to be considered when studying hydrogeological processes in order to employ adequate mechanistic models or perform upscaling. The scale at which a hydrogeological system should be characterized in terms of its spatial heterogeneity and temporal dynamics depends on the studied process and it is not always necessary to consider the full complexity. In this paper, we identify a series of hydrogeological processes for which an approach coupling the monitoring of spatial and temporal variability, including 4D imaging, is often necessary: (1) groundwater fluxes that control (2) solute transport, mixing and reaction processes, (3) vadose zone dynamics, and (4) surface-subsurface water interaction occurring at the interface between different subsurface compartments. We first identify the main challenges related to the coupling of spatial and temporal fluctuations for these processes. Then, we highlight some recent innovations that have led to significant breakthroughs in this domain. We finally discuss how spatial and temporal fluctuations affect our ability to accurately model them and predict their behavior. We thus advocate a more systematic characterization of the dynamic nature of subsurface processes, and the harmonization of open databases to store hydrogeological data sets in their four-dimensional components, for answering emerging scientific question and addressing key societal issues.

Abstract. Essentially all hydrogeological processes are strongly influenced by the subsurface spatial heterogeneity and the temporal variation of environmental conditions, hydraulic properties, and solute concentrations. This spatial and temporal variability generally leads to effective behaviors and emerging phenomena that cannot be predicted from conventional approaches based on homogeneous assumptions and models. However, it is not always clear when, why, how, and at what scale the 4D (3D + time) nature of the subsurface needs to be considered in hydrogeological monitoring, modeling, and applications. In this paper, we discuss the interest and potential for the monitoring and characterization of spatial and temporal variability, including 4D imaging, in a series of hydrogeological processes: (1) groundwater fluxes, (2) solute transport and reaction, (3) vadose zone dynamics, and (4) surface–subsurface water interactions. We first identify the main challenges related to the coupling of spatial and temporal fluctuations for these processes. We then highlight recent innovations that have led to significant breakthroughs in high-resolution space–time imaging and modeling the characterization, monitoring, and modeling of these spatial and temporal fluctuations. We finally propose a classification of processes and applications at different scales according to their need and potential for high-resolution space–time imaging. We thus advocate a more systematic characterization of the dynamic and 3D nature of the subsurface for a series of critical processes and emerging applications. This calls for the validation of 4D imaging techniques at highly instrumented observatories and the harmonization of open databases to share hydrogeological data sets in their 4D components.