Eric Zechner

Variable-density flow in heterogeneous porous media - laboratory experiments and numerical simulations.

Konz, M., Ackerer, P., Younes, A., Huggenberger, P., Zechner, E., 2009a. 2D Stable Layered Laboratory-scale Experiments for Testing Density-coupled Flow Models. Water Resources Research, 45. doi:10.1029/2008WR007118., a series of laboratory-scale 2D tank experiments were conducted and accurately simulated for density driven flow problems on homogeneous porous media. In the present work, we extended the numerical and experimental studies to heterogeneous problems. The heterogeneous porous medium was constructed with a low permeability zone in the centre of the tank and had well-defined parameters and boundary conditions. Concentration distributions were measured in high resolution using a photometric method and an image analysis technique. The numerical model used for the simulations was based on efficient advanced approximations for both spatial and temporal discretizations. The Method Of Lines (MOL) was used to allow higher-order temporal discretization. Three different boundary conditions, corresponding to different localizations of the inflow and the outflow openings at the opposite edges of the tank, were applied to investigate different flow scenarios in the heterogeneous porous medium flow tank. Simulation results of all three density coupled experiments revealed a density-dependent behavior of dispersion. Thus, a reduction of dispersivites was required to obtain a good matching of the experimental data. The high quality of the experiments enabled a detailed testing of numerical variable-density flow codes under heterogeneous conditions. Therefore, the experiments were considered to be reliable benchmark tests.

On the effects of subsurface parameters on evaporite dissolution (Switzerland)

Uncontrolled subsurface evaporite dissolution could lead to hazards such as land subsidence. Observed subsidences in a study area of Northwestern Switzerland were mainly due to subsurface dissolution (subrosion) of evaporites such as halite and gypsum. A set of 2D density driven flow simulations were evaluated along 1000 m long and 150 m deep 2D cross sections within the study area that is characterized by tectonic horst and graben structures. The simulations were conducted to study the effect of the different subsurface parameters that could affect the dissolution process. The heterogeneity of normal faults and its impact on the dissolution of evaporites is studied by considering several permeable faults that include non-permeable areas. The mixed finite element method (MFE) is used to solve the flow equation, coupled with the multipoint flux approximation (MPFA) and the discontinuous Galerkin method (DG) to solve the diffusion and the advection parts of the transport equation. Results show that the number of faults above the lower aquifer that contains the salt layer is considered as the most important factor that affects the dissolution compared to the other investigated parameters of thickness of the zone above the halite formation, a dynamic conductivity of the lower aquifer, and varying boundary conditions in the upper aquifer.

The role of tectonic structures and density-driven groundwater flow for salt karst formation

Groundwater circulation in evaporite bearing horizons and the resulting evolution of karst frequently causes geotechnical problems such as land subsidence or collapses. Even comparably small subsidence rates can significantly affect sensitive urban infrastructure. Dissolution of deep-seated evaporites occurs when groundwater, undersaturated with respect to highly soluble minerals, enters into contact with an evaporite body. Groundwater can interact from all directions in space. Laboratory dissolution experiments were carried out to study the effect of freshwater in contact with rock salt, accounting for the effect of dipping interfaces as an abstraction of inclined geological formations. The dissolution experiments were carried out both in 2D and 3D. The approach allows studying dissolution in a normal fault juxtaposing a developing intrastratal salt karst against a salt layer or against a poorly soluble layer. Modeling experiments were conducted to study the effect of variable density flow in inclined aquifers bordering halite formations and the effect of different fault geometries. The numerical simulations of density driven flow were evaluated along 1,000 m long and 150 m deep 2D cross sections to study the effect of tectonic horst and graben structures. Results show that a large fault zone consisting of a set of smaller hydraulically higher conductive faults is a major factor for dissolution.

Examples and Case Studies

The presented examples and case studies illustrate specific applications of adaptive management of water resources in the region of Basel, Northwestern Switzerland. Such concepts together with the setup of tools and process-oriented experiments allow testing hypotheses. The applied methods facilitated us to fill several gaps of knowledge of subsurface processes. The examples focus on questions with practical as well as research. Most topics are relevant for urban areas and the sustainable use of subsurface resources in general.