Räumlich differenzierte Modellierung der Grundwasserströmung alluvialer Aquifere für mesoskalige Einzugsgebiete

Abstract

Physikalisch basierte Grundwassermodelle auf Einzugsgebietsebene sind als ein wesentlicher Bestandteil räumlich verteilter hydrologischer Modelle ein unverzichtbares Werkzeug der Global Change Forschung geworden. Dabei erfordert die Größe der Einzugsgebiete in der Regel eine sehr grobe Diskretisierung, um Gesamtmodelle mit akzeptablen Laufzeiten und Speicheranforderungen zu erhalten. Aber selbst bei fortschreitender Computerentwicklung verhindert die Datenlage eine genauere Auflösung, die der Heterogenität des Untergrundes gerecht wird. Als wesentliche Fragestellung folgt daraus, wie sich die Geometrie von Aquiferen auf einem groben Raster darstellen lässt und inwieweit diese Geometrie und die eingesetzten Parameter der letztendlich entwickelten Modelle mit den vorhandenen Messungen in Verbindung gebracht werden können. Die Beschreibung des Modellraumes (Aquifergeometrien, Flussläufe usw.) auf einem groben Raster wird besonders schwierig, wenn Strukturen abgebildet werden müssen, deren durchschnittliche Ausdehnung in der Größenordnung der Rasterweite liegt. Im Einzugsgebiet der Oberen Donau (78.000 km2 bis zum Pegel Passau-Achleiten), das ein Untersuchungsschwerpunkt des BMBF-Projektverbundes GLOWA ist, wurde ein Quadratkilometerraster als Grundlage der Modellierung gewählt. Dort existieren allerdings einige hydrologisch wichtige alluviale Aquifere, deren Ausdehnung im Bereich von einigen 100 Metern bis zu ein oder zwei Kilometern reicht. Diese Aquifere müssen abgebildet werden, da ohne sie der Wassertransfer aus den alpinen Bereichen zur Donau nicht gewährleistet ist. Für die Anpassung von hydrogeologisch wirksamen Grenzflächen, insbesondere der Aquiferbasis der alluvialen Grundwasserleiter, bietet sich eine Berechnung von klar definierten Gradienten auf diesen Grenzflächen an. Damit ist - analog zu den Fließrichtungen auf einem DGM - gewährleistet, dass das Grundwasser bis zum Auslass des Einzugsgebietes akkumuliert werden kann. Desweiteren können nicht plausible Modellergebnisse aufgrund von Interpolationsfehlern oder fehlerhafter Ausgangsdaten mit diesem Ansatz verhindert werden. Dieser neu entwickelte Ansatz wird in der vorliegenden Arbeit auf das Einzugsgebiet der Oberen Donau angewendet und getestet. Als numerischer Code wurde hierfür das Programm MODFLOW gewählt. Um den Einfluss der Diskretisierung auf die Modellergebnisse zu untersuchen, wurden für ein kleineres Einzugsgebiet (Einzugsgebiet der Ammer, 709 km2 bis zum Pegel Fischen am Ammersee) innerhalb des Gesamteinzugsgebiets zwei weitere Grundwassermodelle auf unterschiedlichen Skalen (100 bzw. 1000 m entwickelt. Das 100 m-Modell wurde als Referenzmodell kalibriert und validiert. Danach wurden diverse Upscaling-Techniken angewandt, um Durchlässigkeiten und Speicherkoeffizienten für das grobe Modell auf Basis des Referenzmodells zu optimieren. Letztendlich werden mit diesem Vergleich die Einsatzmöglichkeiten grob diskretisierter Modelle (z.B. in der Global Change Forschung) untersucht und diskutiert.

When assessing the effects of global change upon regional water resources, distributed hydrological models are valuable tools for understanding interactions between the hydrological cycle, the biosphere, land use or the change of land cover. Since these tools are supposed to be applied under changing conditions, they should be as physically-based as possible. One approach to achieve this is to separate the hydrological (and biological) cycle in several components and to represent these components by sets of partial differential equations interrelated by the concepts of continuity of mass and momentum. A spatial and temporal discretization must be chosen for solving these equations numerically. The major advantage of such models is that their parameters are comparable with measured data in the field or have at least a measurable meaning. Additionally, they provide a more detailed description by computing spatially distributed state variables, for example the groundwater level. These results can also be compared to field data. The greatest disadvantage of these models is the amount of parameters they generally need and which are often not known. To apply a physically-based approach on a regional (catchment) scale, the discretization has to be quite coarse in order to achieve reasonable calculation times of the models and manageable storage demand of the output data. At the same time, the coarser the discretization, the less physically-based are the determined parameters of the models. A current example of a distributed complex hydrological model built within the framework of global change research is the model DANUBIA of the program GLOWA-Danube, a research project funded by the German Federal Ministry of Education and Research (BMBF). The project aims at the development and utilization of a decision support system to investigate approaches of sustainable future water use. Therefore modelling the hydrological cycle under changing conditions is a decisive part of this project. The applied spatial discretisation of the overall model is a raster with a cellsize of 1000 m * 1000 m, a compromise between the participating scientific diciplines. For the groundwater component a finite difference scheme was chosen, since it is quite stable, fast and fits very well to a raster-based modelling approach. The research area of GLOWA-Danube, the Upper Danube Basin up to gauge Passau-Achleiten, covers 78.000 km2 in Germany, Austria and Switzerland. Following the identification of the important structures of the regional aquifer system and the creation of a hydrogeological conceptual model, it is important to implement this concept into a groundwater model such that a stable numerical solution of the model is attainable. The main problem in the Danube catchment is achieving a connected alluvial aquifer system which is able to receive the groundwater recharge in the mountain areas and which yields a reasonable base flow at existing gauging stations in the forelands. Due to the discrepancy between the finite difference cell size and the extent of the narrow, highly permeable aquifers, additional highly permeable cells have to be added in order to achieve a closed solution for groundwater flow using a finite difference scheme. In addition, it has to be ensured that each cell of the aquifer has at least one neighbouring cell (in the direction of groundwater flow) with a lower base to guarantee the connectivity of the aquifer. The connectivity can be ensured by upgrading the concept of the hydrological catchment drainage analysis. In this approach, the surface runoff is routed to the catchment gauge by changing the landsurface (the digital elevation model) in such a way that every point (raster cell) has a clearly defined flow direction. To reach this, all sinks in the digital elevation model have to be filled. Thereafter, the flow direction in flat areas must be determined. In this thesis, this concept is expanded to the groundwater runoff. The basis of the alluvial aquifers is adjusted in a similar way as the landsurface so that the groundwater flow can be accumulated to the catchment gauge and the connectivity of the groundwater flow is guaranteed. The landsurface is the upper boundary condition and a minimum aquifer thickness is required. Additionally, it is ensured that the river bottom lies between the aquifer top and bottom. With this approach, a consistent modelling of the geometry of the important components of runoff concentration is possible even with a coarse discretization. The results of this thesis demonstrate that modelling of important hydrogeological structures (here the alluvial aquifers) on a (necessarily) coarser scale than the extent of these structures is possible. The developed extended approach of the hydrological drainage analysis is easy to implement and the designed models are stable and fast and provide useful results for the basinwide analysis of groundwater flow.