In this study we make use of 3-D hydraulic simulations to investigate the regional groundwater flow in the Upper Rhine Graben. The modeling is based on an existing detailed 3-D structural model covering the whole Upper Rhine Graben from the surface down to 14 km of depth. The overall goal of this study is to provide some quantitative analysis on the role of the hydraulic head topology in shaping the underground hydrodynamics by taking into account interactions with the heterogeneous subsurface sedimentary configuration of the basin system. Therefore, the main question addressed by this study can be summarized as follows: does the deep graben flow follow the topographic gradient and the flow direction of the river Rhine from the Alps northward to the northernmost area of the Upper Rhine Graben?
Our results demonstrate the presence of a regional subsurface flow in the
sedimentary rocks aligning from the graben flanks towards its center and in
the southern half of the graben from south to north. The graben-parallel
flow velocity is found to be about 1 order of magnitude lower than the velocity
predicted perpendicular to the main graben axis. Besides these general
trends, the modeling highlights local heterogeneities in the shallow 3-D
flow field. Those arise from the interaction between regional groundwater
flow and the heterogeneous sedimentary configuration. Within the Cenozoic
sediments forming the uppermost aquifer in the model, groundwater flows
are driven by imposed hydraulic gradients from recharge areas located at higher elevations in the Black Forest and Vosges Mountains to the discharge region at a lower elevation in the Rhine valley. The presence of a regional aquitard (Keuper) separating the shallow
and the deeper aquifer system (Muschelkalk, Buntsandstein, and Rotliegend) hinders hydraulic connection among the two aquifer systems. This
is exemplified by the development of a flow system in the deeper aquifers,
which shows a more continuous graben-parallel south–north direction.
Based on these results we can conclude that both the hydraulic head topology
and the level of structuration of the sedimentary sequence exert a 1st-order
role in shaping the regional flow system at depth. The regional model
predicts a heterogeneous flow system within the upper 4 km of the Upper
Rhine Valley, where flow velocities in the graben valley can reach up to 45 mm yr
The motivation behind this study is rooted in the search to quantify the intrinsic connection between the productivity of hydrothermal wells and the hydraulic conditions affecting the reservoir hosting the targeted resources. In order to achieve this goal, we demonstrate here that regional-scale (i.e., basin-wide) fluid flow models are a basic requirement prior to any exploitation strategy. Such models, integrating the regional geology and physics, are essential to understand the hydraulic conditions acting on the basin scale. The hydraulic conditions at that scale of the whole system provide the natural forcing exerted on the reservoirs (Frick et al., 2016; Noack et al., 2013).
The study area is the Upper Rhine Graben (URG hereafter), a still active rift developed during Eocene times by transtensional reactivation of the European Cenozoic Rift (Dèzes et al., 2004). With respect to the regional hydrodynamics, previous studies have identified the URG as a regional discharge valley. The main recharge is thought to occur along the graben flanks (Black Forest, Vosges Mountains and Odenwald) with groundwater discharging preferentially within the graben itself (Clauser and Villinger, 1990; Freymark et al., 2019; Guillou-Frottier et al., 2013; Lampe and Person, 2002; Pribnow and Clauser, 2000). Indeed, the deepest topographic point (also the lowest elevation of the hydraulic head) is found in the northern domain within the federal state of Hesse. Based on these observations, it is likely that the regional flow, as driven by hydraulic head gradients, is characterized by a dominant south-to-north component. This said, a physical-based quantification of this conceptual model has never been done. Therefore, in this study we aim to address the following basic questions: (a) what is the regional hydraulic configuration of the whole URG? (b) Does the deep fluid flow follow the shallow surface hydrodynamics as driven by the water table topology, or are those systems hydraulically disconnected? And to what degree?
To answer these questions we investigate the regional fluid flow within the URG, relying on a regional hydraulic model covering the entire URG. The structural model is based on an up-to-date detailed geological model (Freymark et al., 2017) limited here to 14 km of depth (Fig. 1), assuming this geological boundary is acting as a no-flow boundary. The finite-element model has been built within the commercial software FEFLOW© (Diersch, 2014); more information is given in Sect. 2. For more details about the structural modeling and the geological database used as input data, we refer to Freymark et al. (2017, 2015). It is worth mentioning that the input model integrates all available data as well as earlier modeling results and therefore provides a geologically consistent base for our numerical investigations on the shallow and deep fluid flow circulation dynamics with respect to their driving force. Given our scientific goal, and in order to minimize sources of uncertainties in the modeling studies, we neglect temperature, salinity, and mechanical effects, therefore solving only a hydraulic problem. The rationale behind our choice stems from the fact that it allows for the quantification of the pure pressure-driven regional fluid flow dynamics. These are in response to variations in hydraulic head gradients and variations in the subsurface hydraulic conductivity field.
As discussed in the Introduction, the input structural model integrates available seismic data, field and well observations, and published depth and thickness maps, and it assembles the information from previously published 3-D structural models (more detail in Freymark et al., 2017). The base of the model was limited to 14 km of depth assuming that the upper crust at that depth is hydraulically impermeable and that the fluid flow should be negligible at those depths (Fig. 1a). Data from groundwater measurements (BRGM, 2016; Herrmann, 2010; LUBW, 2016a, b) were integrated to determine an observation-based hydraulic head (imposed as a fixed upper boundary condition in the simulations), while the range of variations of hydraulic parameters was inferred from literature data (Table 1).
The modeling carried out is based on the assumption of fluid flow under steady-state conditions. The boundary conditions were set at all six borders: at the top and bottom and at the four vertical boundaries. The four vertical (south, west, north, and east) and the bottom model boundaries are closed to fluid flow (essential Neumann boundary condition). As an upper boundary condition a Dirichlet constant head is imposed based on hydraulic head measurements (black dots in Fig. 2). With this type of boundary condition the hydraulic head is fixed at each node at the top of the model, and inflow and outflow (recharge and discharge) are therefore allowed. In some areas of high topographic relief (e.g., Vosges Mountains, Jura Mountains, Alps, Black Forest, Hunsrück) the paucity of available data required us to rely on a certain degree of approximation by smoothing the topography at the highest elevations as a subdued replica of the groundwater surface (Fig. 2). In the Alps we fixed the groundwater surface to 100 m below the topography.
Applied upper Dirichlet boundary condition (BC) with input data for the interpolation in black (LUBW, 2016a, b; BRGM, 2016; Herrmann, 2010). Hydraulic head with respect to mean sea level.
Following our modeling assumption, the distribution of Darcy velocities
(
Physical input property hydraulic conductivity (
The numerical model comprises 11 geological units resolved by 68 numerical layers. The number of numerical layers has been imposed to ensure numerical stability (Fig. 3). The unstructured mesh consists of 6 736 896 tetrahedral elements. The model size is approximately 500 km in the south–north direction, 300 km in the west–east direction, and 14 km in vertical depth.
Numerical mesh illustrated as white lines:
We discuss the modeling outcomes based on specific profiles (their locations are shown by the green lines in Fig. 1): one cutting the graben parallel to its axis, while the other cross sections are perpendicular to the graben axis, with one located in the south, one in the center, and the last one in the north of the study domain.
A general feature common to all hydraulic head profiles is that the graben-perpendicular hydraulic head gradients (Fig. 4b–d) are stronger than those along the graben axis (Fig. 4a), mainly reflecting the topology of the groundwater head distribution used to impose the upper boundary condition (Fig. 2). Another observation that holds true for all investigated sections is the influence of the Keuper aquitard on the flow field; the latter is schematized in white in Fig. 4e. The flow in the deeper aquifers (pre-Keuper sequence) aligns according to the top of the crystalline basement. In contrast, the groundwater flow in the upper aquifers is decoupled from that observed at greater depth, being mainly driven by surface hydraulic head gradients. Hydraulic decoupling along the two regional aquifers is indicated by an inflection in the isolines of hydraulic heads in the places where the Keuper aquitard is present.
The hydraulic head distribution across the URG varies while moving from the southern to the northern domain. The Black Forest (groundwater head up to 1.2 km) to the east, the Jura Mountains (groundwater head up to 1.3 km) to the south, and the Vosges Mountains (groundwater head up to 1.25 km) to the west confine the groundwater flow around the southern URG. The hydraulic head gradient from the Vosges Mountains and the Black Forest to the center of the URG is approximately 2.5 m per 100 m. On the other hand, the gradient from the Jura Mountains to the lowest elevation point in the URG in the north amounts to only 0.3 m per 100 m (integrated over the whole length of the URG). Hydraulic head gradients perpendicular to the graben axis are about 7 times higher in magnitude in the southern part of the URG than the graben-parallel gradients. This is reflected by the modeled fluid flow, which shows a preferential component infiltrating from the graben shoulders into the graben itself. However, the magnitudes of hydraulic gradients vary along the graben. These gradients are highest in the southern domain (Fig. 4b), lowest in the central part (Fig. 4c), and moderate in the northern part of the URG (1.6 m per 100 m from the Odenwald recharge area to the graben center; Fig. 4d). The head gradient perpendicular to the graben axis in the central part of the URG is around 0.2–0.5 m per 100 m and 3 times higher than the parallel gradient, which is about 0.08 m per 100 km between Strasbourg and Heidelberg.
Simulated hydraulic head distribution
To illustrate the 3-D fluid flow across the study area we make use of
distinct maps on which we plot the
Going back to the main working hypothesis of this study, that is, whether or
not the regional fluid flow follows a preferential south-to-north axis over
the whole study area, we systematically discuss in the following the results
based on the
The northern part of the URG, between Mannheim and Frankfurt, is characterized by an even more complex flow regime lacking a dominant northward flow direction, though groundwater still ascends below the river Rhine. This requires a revision of the groundwater flow dynamics that, despite showing a south–north flow component, is locally overprinted by a small-scale circulation, leading to the occurrence of a heterogeneous flow field. In this region, the main flow component is provided by water recharging from the western and eastern graben shoulders (Figs. 4d and 5a). Moreover, the top of the basement in the northernmost part of the graben dips southward, thus imposing a north-to-south-oriented preferential direction to the flow.
Comparing the results obtained for the northern and the southern domains, we can identify hydrodynamics that share some characteristics, while being diverse overall. On a graben-wide scale, the recharge is from the graben shoulders to the graben center. This fluid flow component, perpendicular to the graben axis, is mainly present in the upper aquifer (post-Keuper sequence). It is driven by higher hydraulic gradients perpendicular to the URG than those parallel to the graben axis. However, the groundwater dynamics in the southern domain show a rather homogeneous character with consistent northward flow, which contrasts the one occurring further north. There, groundwater flow is controlled by the top of the basement deepening to the south, leading to a preferential north–south flow direction, which is therefore opposite to the dynamics imposed by topographic gradients.
Outside the URG area, predicted fluid flow patterns resemble the surface hydraulic head distribution, with infiltration in the Alpine region and pronounced fluid uprising in the bordering valleys (Fig. 5). Infiltration is also predicted to occur beneath the Black Forest and Vosges Mountains in response to elevated hydraulic gradients. For the latter, inflow occurs radially with a weak E–W component, while the inflow axis in the Black Forest is oriented SSW–NNE (Fig. 5).
Streamlines of simulated fluid flow with different colors for the
To better assess the reliability of the simulation results we cross-plot available data on mapped springs, here used as indicators for discharge areas. Therefore, we compare the location of mineral and thermal springs (green in Fig. 5c) and artesian conditions (yellow in Fig. 5c) with the model prediction (upflows and downflows). Based on the fit between model and hydraulic data we can distinguish two main areas: (1) domains in which the model predicts uprising water where springs or artesian conditions are observed (areas 1 to 4 in Fig. 5c, f) and (2) domains in which the model only partly predicts the observed springs or artesian conditions (areas 5 to 7 in Fig. 5c, f).
Along the Western Border Fault of the URG the model can reproduce the
artesian conditions (springs in area number 1 in Fig. 5c). Recharge
occurs, though with moderate rates of approximately 1
The misfit found for areas 5 to 7 also likely relates to uprising fluid
flow, mainly fault and fracture controlled. Here we propose that a better
fit could be achieved by resolving these features either as part of the mesh
(discrete) or as an equivalent continuum. Likewise, the springs in area number
6 of the Black Forest are not reproduced. This is because the Black Forest
acts as a recharge area in this regional-scale model, for which local-scale
features are not resolved. In area 7 groundwater infiltrates from the
surface following a northwest direction, parallel to the hydraulic
gradients, from the northernmost area of the Black Forest into the area
between Heidelberg and Karlsruhe. The observation that we cannot reproduce
the artesian conditions in this area can be taken as proof
that the model could not resolve some important local structural features in
these domains. However, previous modeling studies (Cherubini et
al., 2014; Freymark et al., 2019) have demonstrated that even
large faults implemented as features that are more permeable than the
geological units have only limited lateral effects on the flow field and do
not change the regional, pressure-driven flow pattern. In particular, the
recent modeling work of Freymark et al. (2019) deals with the influence of
the border faults of the URG on the general fluid flow and heat transport.
They could show that the effects of the main border faults are less
important than hydraulic head and permeability contrasts between the sediments
and basement. The aperture and permeability of the faults have no
significant influence on the overall simulated flow and the thermal field.
In order to check the plausibility of the model results concerning recharge and
discharge rates, we calculated the latter for the upper model boundary. The
results show that the majority of values lies in a reasonable range of
Based on the results discussed above, we end up with a revised conceptual model of the regional flow field within the URG as exemplified in Fig. 6. Higher hydraulic head gradients perpendicular to the graben than parallel to the graben axis result in a predominant flow from the graben shoulders to the center of the graben. This characteristic is applicable to the total extent of the URG and is consistent with previous studies (e.g., Freymark et al., 2019) already demonstrating that a flow regime from a recharging eastern boundary fault to a discharging western boundary fault (Clauser and Villinger, 1990) is not supported by 3-D models considering 1st-order structural characteristics. In addition, the larger spatial extent of our model compared to the model of Freymark et al. (2019) allows for the detection of additional influencing factors. Accordingly, the flow from south to north is dominant in the southern area of the graben because of the high head elevation of the Jura Mountains south of the URG. In the central part of the URG, the upflow axis shifts to the eastern area of the graben. There, the influence of the top of the Variscan basement (deepest area of the basement) plays an important role in the deep and shallow aquifer systems. In the north, the deep fluid flow is not following the surface topography northward (as the river Rhine) and the question posed at the beginning, if the deep fluid flow in the URG valley follows the topography from south to north, has to be revised. In the north of the graben, the aquitard is not fully separating the deep aquifers from the shallow aquifer, though the fluid flow is a complex result of deep and shallow hydrodynamics.
Schematic drawing of regional fluid flow in the URG based on the hydraulic model results.
Based on the results and discussion the following major conclusions can be
drawn (summarized in Fig. 6).
The general fluid flow direction is dependent on the hydraulic head
gradients and the geometric configuration of the permeable units. The high
hydraulic gradients at the graben shoulders and low gradients in the Rhine
Valley in concert with the northward-inclined base of the Cenozoic sediment
fill cause a flow direction perpendicular to the flow direction of
the river Rhine and perpendicular to the graben axis, from the graben
shoulders towards its center. Local flow heterogeneities in the graben can be structurally linked to the
topology of the top basement configuration of the Variscan basement. Two separated flow systems evolve in the Upper Rhine Graben that show
different degrees of decoupling through the local presence or absence of the
Keuper aquitard.
Ongoing activities focus on potential modifications of the flow field if full thermo-hydraulic coupling is considered and on the role of the tectonic structures on the regional flow component, which might be overprinted locally as identified by the clustering of localized hydrothermal surface manifestations (springs as per Fig. 5c).
The data used to support the findings of this study are available from the corresponding author upon request.
NK and JB prepared the concept for this study. NK and MF built the model. NK prepared the paper with contributions from all coauthors. MS-W, JB, and MC supervised this work.
The authors declare that they have no conflict of interest.
This article is part of the special issue “European Geosciences Union General Assembly 2019, EGU Division Energy, Resources & Environment (ERE)”. It is a result of the EGU General Assembly 2019, Vienna, Austria, 7–12 April 2019.
We are grateful to the Federal Ministry for Economic Affairs and Industry for funding this research (grant number: 0325944B). Moreover, we are thankful for the cooperation with the Technical University of Darmstadt within the project Hessen 3D 2.0 and to HLNUG for providing hydrological data. For the sponsored MIKE Powered by DHI license files, we kindly acknowledge DHI WASY. We want to thank the anonymous reviewers for their contributions, which helped improve the quality of this paper.
The article processing charges for this open-access publication were covered by a Research Centre of the Helmholtz Association.
This paper was edited by Michael Kühn and reviewed by two anonymous referees.