ALBERT

Description: The objective of this work is the simulation of saturated and unsaturated flow in a karstified aquifer using a double continuum approach. The HydroGeoSphere code (Therrien et al., 2006) is employed to simulate spring discharge with the Richards equations and van Genuchten parameters to represent flow in the (1) fractured matrix and (2) conduit continuum coupled by a linear exchange term. Rapid vertical small-scale flow processes in the unsaturated conduit continuum are accounted for by applying recharge boundary conditions at the bottom of the saturated model domain. An extensive sensitivity analysis is performed on single parameters as well as parameter combinations. The transient hydraulic response of the karst spring is strongly controlled by the matrix porosity as well as the van Genuchten parameters of the unsaturated matrix, which determine the head dependent inter-continuum water transfer when the conduits are draining the matrix. Sensitivities of parameter combinations partially reveal a non-linear dependence over the parameter space. This can be observed for parameters not belonging to the same continuum as well as combinations, which involve the exchange parameter, showing that results of the double continuum model may depict a certain degree of ambiguity. The application of van Genuchten parameters for simulation of unsaturated flow in karst systems is critically discussed.

Description: The deduction by conventional means of qualitative and quantitative information about groundwater discharge into lakes is complicated. Nevertheless, at least for semi-arid regions with limited surface water availability, this information is crucial to ensure future water availability for drinking and irrigation purposes. Overcoming this lack of discharge information, we present a satellite-based multi-temporal sea-surface-temperature (SST) approach. It exploits the occurrence of thermal anomalies to outline groundwater discharge locations using the example of the Dead Sea. Based on a set of 19 Landsat Enhanced Thematic Mapper (ETM+) images 6.2 (high gain), recorded between 2000 and 2002, we developed a novel approach which includes (i) an objective exclusion of surface-runoff-influenced data which would otherwise lead to erroneous results and (ii) a temporal SST variability analysis based on six statistical measures amplifying thermal anomalies caused by groundwater. After excluding data influenced by surface runoff, we concluded that spatial anomaly patterns of the standard deviation and range of the SST data series spatially fit best to in situ observed discharge locations and, hence, are most suitable for detecting groundwater discharge sites.

Description: Karst aquifers are characterized by highly conductive conduit flow paths embedded in a less conductive fissured and fractured matrix, resulting in strong permeability contrasts with structured heterogeneity and anisotropy. Groundwater storage occurs predominantly in the fissured matrix. Hence, most mathematical karst models assume quasi-steady-state flow in conduits neglecting conduit-associated drainable storage (CADS). The concept of CADS considers storage volumes, where karst water is not part of the active flow system but hydraulically connected to conduits (for example karstic voids and large fractures). The disregard of conduit storage can be inappropriate when direct water abstraction from karst conduits occurs, e.g., large-scale pumping. In such cases, CADS may be relevant. Furthermore, the typical fixed-head boundary condition at the karst outlet can be inadequate for water abstraction scenarios because unhampered water inflow is possible. The objective of this work is to analyze the significance of CADS and flow-limited boundary conditions on the hydraulic behavior of karst aquifers in water abstraction scenarios. To this end, the numerical discrete-continuum model MODFLOW-2005 Conduit Flow Process Mode 1 (CFPM1) is enhanced to account for CADS. Additionally, a fixed-head limited-flow (FHLQ) boundary condition is added that limits inflow from constant head boundaries to a user-defined threshold. The effects and the proper functioning of these modifications are demonstrated by simplified model studies. Both enhancements, CADS and FHLQ boundary, are shown to be useful for water abstraction scenarios within karst aquifers. An idealized representation of a large-scale pumping test in a karst conduit is used to demonstrate that the enhanced CFPM1 is able to adequately represent water abstraction processes in both the conduits and the matrix of real karst systems, as illustrated by its application to the Cent Fonts karst system.

Description: We describe a semi-automatic method to objectively and reproducibly extract lineaments based on the global one arc-second ASTER GDEM. The combined method of linear filtering and object-based classification ensures a high degree of accuracy resulting in a lineament map. Subsequently lineaments are differentiated into geological and morphological lineaments to assign a probable origin and hence a hydro-geological significance. In the western catchment area of the Dead Sea (Israel) the orientation and location of the differentiated lineaments are compared to characteristics of known structural features. The authors demonstrate that a strong correlation between lineaments and structural features exist, being either influenced by the Syrian Arc paleostress field or the Dead Sea stress field or by both. Subsequently, we analyse the distances between lineaments and wells thereby creating an assessment criterion concerning the hydraulic significance of detected lineaments. Derived from this analysis the authors suggest that the statistic analysis of lineaments allows a delineation of flow-paths and thus significant information for groundwater analysis. We validate the flow-path delineation by comparison with existing groundwater model results based on well data.

Description: Due to their heterogeneous nature, karst aquifers pose a major challenge for hydrogeological investigations. Important procedures like the delineation of catchment areas for springs are hindered by the unknown locations and hydraulic properties of highly conductive karstic zones. In this work numerical modeling was employed as a tool in delineating catchment areas of several springs within a karst area in southwestern Germany. For this purpose, different distributive modeling approaches were implemented in the Finite Element simulation software Comsol Multiphysics®. The investigation focuses on the question to which degree the effect of karstification has to be taken into account for accurately simulating the hydraulic head distribution and the observed spring discharges. The results reveal that the representation of heterogeneities has a large influence on the delineation of the catchment areas. Not only the location of highly conductive elements but also their geometries play a major role for the resulting hydraulic head distribution and thus for catchment area delineation. The size distribution of the karst conduits derived from the numerical models agrees with knowledge from karst genesis. It was thus shown that numerical modeling is a useful tool for catchment delineation in karst aquifers based on results from different field observations.

Description: Information on groundwater discharge over large spatial scales are essential for groundwater management particularly in (semi-) arid regions. If discharge areas are known, direct measurements over larger spatial scales are complicated to obtain by conventional means, why thermal remote sensing is increasingly applied to localize and quantify groundwater discharge. In this context, mostly unconsidered is (i) the influence of surface-runoff that can negatively affect groundwater focused studies and (ii) the representativeness of remotely sensed groundwater discharge based on single thermal images, against the background of discharge intermittency. Addressing these issues we apply a multi-temporal SST data approach based on 19 Landsat ETM+ band 6.2 (high gain) data from the year 2000 until 2002 at the example of the (semi-)arid vicinity of Dead Sea. To be independent of auxiliary rain data we develop a novel approach to identify surface-runoff influenced images solely using image statistics. Compared to foregoing rain events the result reveals a general influence-time of at least two days, but also that a simple time-difference criterion to exclude possible surface-runoff is not advisable. In the second part of the study we evaluate the significance of six statistical measures calculated on a per-pixel basis on the remaining 12 surface-runoff uninfluenced sea-surface-temperature (SST) data using in situ discharge measurements of the Israel Hydrological Service (IHS). We found that the spatial patterns of the standard deviation and range on the SST data series best fit to the IHS observed discharge locations and hence are suitable for detecting groundwater discharge areas.

Description: Karst aquifers are characterized by highly conductive conduit flow paths embedded in a less conductive fissured and fractured matrix resulting in strong permeability contrasts with structured heterogeneity and anisotropy. Groundwater storage occurs predominantly in the fissured matrix. Hence, most karst models assume quasi steady-state flow in conduits neglecting conduit associated drainable storage (CADS). The concept of CADS considers storage volumes, where karst water is not part of the active flow system but rather hydraulically connected to conduits (for example karstic voids and large fractures). The disregard of conduit storage can be inappropriate when direct water abstraction from karst conduits occurs, e.g. large scale pumping. In such cases, CADS may be relevant. Furthermore, the typical fixed head boundary condition at the karst outlet can be inadequate for water abstraction scenarios because unhampered water inflow is possible. The objective of this paper is to analyze the significance of CADS and flow-limited boundary conditions on the hydraulic behavior of karst aquifers in water abstraction scenarios. To this end, the numerical hybrid model MODFLOW-2005 Conduit Flow Process Mode 1 (CFPM1) is enhanced to account for CADS. Additionally, a fixed-head limited-flow (FHLQ) boundary condition is added that limits inflow from constant head boundaries to a user-defined threshold. The affect and proper functioning of these modifications is demonstrated by simplified model studies. Both enhancements, CAD storage and the FHLQ boundary, are shown to be useful for water abstraction scenarios within karst aquifers. An idealized representation of a large-scale pumping test in a karst conduit is used to demonstrate that the enhanced CFPM1 is potentially able to adequately represent water abstraction processes in both the conduits and the matrix of real karst systems.

Description: Assessing the hydraulic parameters of karst aquifers is a challenge due to their high degree of heterogeneity. The unknown parameter field generally leads to a high ambiguity for flow and transport calibration in numerical models of karst aquifers. In this study, a distributive numerical model was built for the simulation of groundwater flow and solute transport in a highly heterogeneous karst aquifer in south western Germany. Therefore, an interface for the simulation of solute transport in one-dimensional pipes was implemented into the software Comsol Multiphysics® and coupled to the three-dimensional solute transport interface for continuum domains. For reducing model ambiguity, the simulation was matched for steady-state conditions to the hydraulic head distribution in the model area, the spring discharge of several springs and the transport velocities of two tracer tests. Furthermore, other measured parameters such as the hydraulic conductivity of the fissured matrix and the maximal karst conduit volume were available for model calibration. Parameter studies were performed for several karst conduit geometries to analyse the influence of the respective geometric and hydraulic parameters and develop a calibration approach in a large-scale heterogeneous karst system. Results show that it is not only possible to derive a consistent flow and transport model for a 150 km2 karst area, but that the combined use of groundwater flow and transport parameters greatly reduces model ambiguity. The approach provides basic information about the conduit network not accessible for direct geometric measurements. The conduit network volume for the main karst spring in the study area could be narrowed down to approximately 100 000 m3.

Description: The coastDat data sets were produced to give a consistent and homogeneous database mainly for assessing weather statistics and long-term changes for Europe, especially in data sparse regions. A sequence of numerical models was employed to reconstruct all aspects of marine climate (such as storms, waves, surges etc.) over many decades. Here, we describe the atmospheric part of coastDat2 (Geyer and Rockel, 2013, doi:10.1594/WDCC/coastDat-2_COSMO-CLM). It consists of a regional climate reconstruction for entire Europe, including Baltic and North Sea and parts of the Atlantic. The simulation was done for 1948 to 2012 with a regional climate model and a horizontal grid size of 0.22° in rotated coordinates. Global reanalysis data were used as forcing and spectral nudging was applied. To meet the demands on the coastDat data set about 70 variables are stored hourly.

Description: The coastDat data sets were produced to give a consistent and homogeneous database mainly for assessing weather statistics and climate changes since 1948, e.g., in frequencies of extremes for Europe, especially in data sparse regions. A sequence of numerical models was employed to reconstruct all aspects of marine climate (such as storms, waves, surges, etc.) over many decades. The acronym coastDat stands for the set of consistent ocean and atmospheric data, where the atmospheric data where used as forcing for the reconstruction of the sea state. Here, we describe the atmospheric part of coastDat2 (Geyer and Rockel, 2013; doi:10.1594/WDCC/coastDat-2_COSMO-CLM). It consists of a regional climate reconstruction for the entire European continent, including the Baltic Sea and North Sea and parts of the Atlantic. The simulation was done for 1948 to 2012 with the regional climate model COSMO-CLM (CCLM) and a horizontal grid size of 0.22 degree in rotated coordinates. Global reanalysis data of NCEP1 were used as forcing and spectral nudging was applied. To meet the demands on the coastDat data set about 70 variables are stored hourly.