Knowledge of groundwater flow is of high relevance for groundwater management and the planning of different subsurface utilizations such as deep geothermal facilities. While numerical models can help to understand the hydrodynamics of the targeted reservoir, their predictive capabilities are limited by the assumptions made in their set up. Among others, the choice of appropriate hydraulic boundary conditions, adopted to represent the regional to local flow dynamics in the simulation run, is of crucial importance for the final modelling result. In this publication we present the hydrogeological models to obtain results to quantify how and to which degree different upper hydraulic boundary conditions and vertical cross boundary fluid movement influence the calculated deep fluid conditions Therefore, we take the central Upper Rhine Graben area as a natural laboratory. The presented three models are set up with different sets of boundary conditions. The Reference Model uses the topography as upper hydraulic pressure surface of 0 kPa. In model S1, a subdued replica of the topography, which was built on the base of hydraulic head measurements is applied as the upper hydraulic boundary condition and in model S2 vertical cross boundary flow is implemented.

Based on our results, we illustrate in the landing paper that for the Upper Rhine Graben specific characteristics of the flow field are robust and insensitive to the choice of imposed hydraulic boundary conditions, while specific local characteristics are more sensitive. Accordingly, these robust features characterizing the first order groundwater flow dynamics in the Upper Rhine Graben include: (i) a regional groundwater flow component descending from the graben shoulders to rise at its centre; (ii) infiltration of fluids across the graben shoulders, which locally rise along the main border faults; (iii) the presence of heterogeneous hydraulic potentials at the rift shoulders. The configuration of the adopted boundary conditions influence primarily calculated flow velocities and the absolute position of the upflow axis within the graben sediments. In addition, the choice of upper hydraulic boundary conditions exerts a direct control on the evolving local flow dynamics, with the degree of influence gradually decreasing with increasing depth.

With respect to regional flow modelling of basin hosted, deep water resources, the main conclusions derived from this study are: (i) the often considered water table as an exact replica of the model topography (Reference Model) likely introduces a source of error in the simulations in regional hydraulic modelling approaches. Here, we show that these errors can be minimized by making use of a water table as upper boundary condition derived from available hydraulic head measurements (model S1). If the study area is part of a supra-regional flow system - like the central Upper Rhine Graben is part of the whole Upper Rhine Graben - the in- and outflow across vertical boundaries need to be considered (model S2).