Abstract. This study explores the suitability of a single hillslope as a parsimonious representation of a catchment in a physically based model. We test this hypothesis by picturing two distinctly different catchments in perceptual models and translating these pictures into parametric setups of 2-D physically based hillslope models. The model parametrizations are based on a comprehensive field data set, expert knowledge and process-based reasoning. Evaluation against streamflow data highlights that both models predicted the annual pattern of streamflow generation as well as the hydrographs acceptably. However, a look beyond performance measures revealed deficiencies in streamflow simulations during the summer season and during individual rainfall-runoff events as well as a mismatch between observed and simulated soil water dynamics. Some of these shortcomings can be related to our perception of the systems and to the chosen hydrological model, while others point to limitations of the representative hillslope concept itself. Nevertheless, our results confirm that representative hillslope models are a suitable tool to assess the importance of different data sources as well as to challenge our perception of the dominant hydrological processes we want to represent therein. Consequently, these models are a promising step forward in the search for the optimal representation of catchments in physically based models.

Abstract. According to Dooge (1986) intermediate-scale catchments are systems of organized complexity, being too organized and yet too small to be characterized on a statistical/conceptual basis, but too large and too heterogeneous to be characterized in a deterministic manner. A key requirement for building structurally adequate models precisely for this intermediate scale is a better understanding of how different forms of spatial organization affect storage and release of water and energy. Here, we propose that a combination of the concept of hydrological response units (HRUs) and thermodynamics offers several helpful and partly novel perspectives for gaining this improved understanding. Our key idea is to define functional similarity based on similarity of the terrestrial controls of gradients and resistance terms controlling the land surface energy balance, rainfall runoff transformation, and groundwater storage and release. This might imply that functional similarity with respect to these specific forms of water release emerges at different scales, namely the small field scale, the hillslope, and the catchment scale. We thus propose three different types of "functional units" -specialized HRUs, so to speak -which behave similarly with respect to one specific form of water release and with a characteristic extent equal to one of those three scale levels. We furthermore discuss an experimental strategy based on exemplary learning and replicate experiments to identify and delineate these functional units, and as a promising strategy for characterizing the interplay and organization of water and energy fluxes across scales. We believe the thermodynamic perspective to be well suited to unmask equifinality as inherent in the equations governing water, momentum, and energy fluxes: this is because several combinations of gradients and resistance terms yield the same mass or energy flux and the terrestrial controls of gradients and resistance terms are largely independent. We propose that structurally adequate models at this scale should consequently disentangle driving gradients and resistance terms, because this optionally allows Published by Copernicus Publications on behalf of the European Geosciences Union. E. Zehe et al.: HESS Opinions: Thermodynamic reinterpretation of the HRU conceptequifinality to be partly reduced by including available observations, e.g., on driving gradients. Most importantly, the thermodynamic perspective yields an energy-centered perspective on rainfall-runoff transformation and evapotranspiration, including fundamental limits for energy fluxes associated with these processes. This might additionally reduce equifinality and opens up opportunities for testing thermodynamic optimality principles within independent predictions of rainfall-runoff or land surface energy exchange. This is pivotal to finding out whether or not spatial organization in catchments is in accordance with a fundamental organizing principle.

Abstract. The study deals with the identification and characterization of rapid subsurface flow structures through pedo- and geo-physical measurements and irrigation experiments at the point, plot and hillslope scale. Our investigation of flow-relevant structures and hydrological responses refers to the general interplay of form and function, respectively. To obtain a holistic picture of the subsurface, a large set of different laboratory, exploratory and experimental methods was used at the different scales. For exploration these methods included drilled soil core profiles, in situ measurements of infiltration capacity and saturated hydraulic conductivity, and laboratory analyses of soil water retention and saturated hydraulic conductivity. The irrigation experiments at the plot scale were monitored through a combination of dye tracer, salt tracer, soil moisture dynamics, and 3-D time-lapse ground penetrating radar (GPR) methods. At the hillslope scale the subsurface was explored by a 3-D GPR survey. A natural storm event and an irrigation experiment were monitored by a dense network of soil moisture observations and a cascade of 2-D time-lapse GPR trenches. We show that the shift between activated and non-activated state of the flow paths is needed to distinguish structures from overall heterogeneity. Pedo-physical analyses of point-scale samples are the basis for sub-scale structure inference. At the plot and hillslope scale 3-D and 2-D time-lapse GPR applications are successfully employed as non-invasive means to image subsurface response patterns and to identify flow-relevant paths. Tracer recovery and soil water responses from irrigation experiments deliver a consistent estimate of response velocities. The combined observation of form and function under active conditions provides the means to localize and characterize the structures (this study) and the hydrological processes (companion study Angermann et al., 2017, this issue).

Abstract. The phrase form and function was established in architecture and biology and refers to the idea that form and functionality are closely correlated, influence each other, and co-evolve. We suggest transferring this idea to hydrological systems to separate and analyze their two main characteristics: their form, which is equivalent to the spatial structure and static properties, and their function, equivalent to internal responses and hydrological behavior. While this approach is not particularly new to hydrological field research, we want to employ this concept to explicitly pursue the question of what information is most advantageous to understand a hydrological system. We applied this concept to subsurface flow within a hillslope, with a methodological focus on function: we conducted observations during a natural storm event and followed this with a hillslope-scale irrigation experiment. The results are used to infer hydrological processes of the monitored system. Based on these findings, the explanatory power and conclusiveness of the data are discussed. The measurements included basic hydrological monitoring methods, like piezometers, soil moisture, and discharge measurements. These were accompanied by isotope sampling and a novel application of 2-D time-lapse GPR (ground-penetrating radar). The main finding regarding the processes in the hillslope was that preferential flow paths were established quickly, despite unsaturated conditions. These flow paths also caused a detectable signal in the catchment response following a natural rainfall event, showing that these processes are relevant also at the catchment scale. Thus, we conclude that response observations (dynamics and patterns, i.e., indicators of function) were well suited to describing processes at the observational scale. Especially the use of 2-D time-lapse GPR measurements, providing detailed subsurface response patterns, as well as the combination of stream-centered and hillslope-centered approaches, allowed us to link processes and put them in a larger context. Transfer to other scales beyond observational scale and generalizations, however, rely on the knowledge of structures (form) and remain speculative. The complementary approach with a methodological focus on form (i.e., structure exploration) is presented and discussed in the companion paper by Jackisch et al. (2017).

No abstract

Analysis of time-lapse ground-penetrating radar (GPR) data can provide information regarding subsurface hydrological processes, such as preferential flow. However, the analysis of time-lapse data is often limited by data quality; for example, for noisy input data, the interpretation of difference images is often difficult. Motivated by modern image-processing tools, we have developed two robust GPR attributes, which allow us to distinguish amplitude (contrast similarity) and time-shift (structural similarity) variations related to differences between individual time-lapse GPR data sets. We tested and evaluated our attributes using synthetic data of different complexity. Afterward, we applied them to a field data example, in which subsurface flow was induced by an artificial rainfall event. For all examples, we identified our structural similarity attribute to be a robust measure for highlighting time-lapse changes also in data with low signal-to-noise ratios. We determined that our new attribute-based workflow is a promising tool to analyze time-lapse GPR data, especially for imaging subsurface hydrological processes.

We present a 2D/3D topographic migration scheme for ground‐penetrating radar (GPR) data which is able to account for variable velocities by using the root mean square (rms) velocity approximation. We test our migration scheme using a synthetic 2D example and compare our migrated image to the results obtained using common GPR migration approaches. Furthermore, we apply it to 2D and 3D field data. These examples are recorded across common subsurface settings including surface topography and variations in the GPR subsurface velocity field caused by a shallow ground water table. In such field settings, our migration strategy provides well focused images of common‐offset GPR data without the need for a detailed interval velocity model. The synthetic and field examples demonstrate that our topographic migration scheme allows for accurate GPR imaging in the presence of variations in surface topography and subsurface velocity.

Ice complex deposits are characteristic, ice-rich formations in northern East Siberia and represent an important part in the arctic carbon pool. Recently, these late Quaternary deposits are the objective of numerous investigations typically relying on outcrop and borehole data. Many of these studies can benefit from a 3D structural model of the subsurface for upscaling their observations or for constraining estimations of inventories, such as the local carbon stock. We have addressed this problem of structural imaging by 3D groundpenetrating radar (GPR), which, in permafrost studies, has been primarily used for 2D profiling. We have used a 3D kinematic GPR surveying strategy at a field site located in the New Siberian Archipelago on top of an ice complex. After applying a 3D GPR processing sequence, we were able to trace two horizons at depths below 20 m. Taking available borehole and outcrop data into account, we have interpreted these two features as interfaces of major lithologic units and derived a 3D cryostratigraphic model of the subsurface. Our data example demonstrated that a 3D surveying and processing strategy was crucial at our field site and showed the potential of 3D GPR to image geologic structures in complex ice-rich permafrost landscapes.