In der vorliegenden Dissertation sollten die hydrogeochemischen Prozesse herausgearbeitet werden, die für die Wasserqualitätsänderung während eines ASR Experiments in Bolivar, Südaustralien und während der Versickerung in einem künstlichen Grundwasseranreicherungsbecken in Berlin von Bedeutung waren. Reaktive Stofftransportmodellierung des ASR Experiments in Bolivar, Südaustralien zeigte, dass die hydrochemischen Veränderungen in der direkten Umgebung des Injektionsbrunnens während der Speicherphase nur durch rapide Änderungen der Sauerstoff- und Nitrat reduzierenden Bakterienmasse erklärt werden können. Die hydrochemischen Veränderungen in größerer Distanz zum Injektionsbrunnen wurden überwiegend durch Ionenaustauschprozesse und Kalzitlösung verursacht. Geochemische und hydraulische Messungen unter einem Sickerbecken in Berlin zeigten, dass die beobachteten geochemischen Änderungen im Sickerwasser mit den periodisch auftretenden wassergesättigten/wasserungesättigen Bedingungen unter dem Becken einhergehen. Während der ungesättigten Periode wird Luft unter das Becken gezogen und führt zur plötzlichen Reoxidierung von bereits in der gesättigten Periode gebildeten Eisensulfiden und zur beschleunigten Mineralisation von sedimentärem organischem Kohlenstoff. Reaktive Stofftransportmodellierung auf größerer Skale zeigte, dass allein die saisonalen Temperaturunterschiede im Infiltrationswasser für die beobachtete zeitliche und räumliche Dynamik der Redoxzonen im weiteren Abstrom des Sickerbeckens verantwortlich sind. Das Abbauverhalten der Arzneimittelsubstanz Phenazon hängt ausschließlich von der Verfügbarkeit von gelöstem Sauerstoff und damit indirekt von der Wassertemperatur im Aquifer ab. In der vorliegenden Arbeit wird deutlich, dass ein adäquates Verständnis der wasserqualitätsändernden Prozesse in künstlichen Anreicherungsystemen nur dann erreicht werden kann wenn Strömung, Transport und reaktive Prozesse, im Feld als auch in der Modellierung, simultan betrachtet werden. / In this thesis, three major studies were carried out in order to understand the key factors controlling the water quality changes that occurred during a reclaimed water Aquifer Storage and Recovery (ASR) experiment at Bolivar, South Australia and during ponded infiltration in Berlin, Germany. Multi-component reactive transport modelling of the ASR experiment suggested that during the storage phase, dynamic changes in bacterial mass have a significant influence on the local geochemistry in the vicinity of the injection well. Water quality changes further away from the injection well were mainly driven by ion exchange and calcite dissolution. Geochemical and hydraulic measurements below an artificial recharge pond in Berlin, Germany, showed that the observed dynamic changes of the hydrochemistry within the seepage water are strongly linked to the periodic saturated/unsaturated hydraulic conditions below the pond. During unsaturated conditions, atmospheric oxygen penetrates from the pond margins to the centre below the pond, leading to (i) a sudden re-oxidation of sulphide minerals that have formed previously during saturated conditions and (ii) an enhanced mineralisation of sedimentary particulate organic carbon. Reactive transport modelling showed that at larger scale, seasonal temperature changes of the infiltration water are the key control for the observed temporal and spatial redox dynamics further downstream the recharge pond. Moreover, the degradation behaviour of the pharmaceutically residue phenazone solely depends on the availability of dissolved oxygen, and thus indirectly on the water temperature within the aquifer. Overall this thesis shows that a sound understanding and analysis of the key processes affecting the water quality changes during artificial recharge of groundwater could only be achieved when flow, transport and reactive processes are considered simultaneously, both in the field and during modelling.

Die Bodenzone, in der Literatur vielfach auch Wurzelzone, Aerationszone oder ungesättigte Zone genannt, ist geprägt durch variabel-wassergesättigte Verhältnisse und nimmt in vielen Disziplinen eine wichtige Rolle ein. Aus Sicht des Schutzguts Grundwasser stellt sie eine Schutz- und Pufferzone vor oberirdischen Umwelteinflüssen dar, in der eindringende oder eingebrachte Schadstoffe durch die dort ablaufenden Transport-, Abbau- und Sorptionsprozesse retardiert, teilweise bis vollständig abgebaut oder in andere Stoffe umgesetzt werden können, und somit eine Verunreinigung des Grundwassers verhindern kann. Um potenzielle Gefährdungen des Grundwassers anhand einer Altlast oder eines Schadensfalls abschätzen zu können, ist in Deutschland eine Sickerwasserprognose nach dem Bundesbodenschutzgesetz und der Bundesbodenschutzverordnung vorgeschrieben. Hierbei übernimmt die ungesättigte Zone die Funktion des Quell- und Transportterms für den Schadstoff. Der Quellterm dient der Beschreibung des zeitlichen Austragsverhaltens von Schadstoffen aus der Schadstoffquelle mit dem Sickerwasser, der Transportterm beschreibt den Wirkungspfad im Boden von der Geländeoberkante bis zur Grundwasseroberfläche. Die Anforderungen und Aufgaben des vom BMBF geförderten Forschungsvorhabens „Prognose des Schadstoffeintrags in das Grundwasser mit dem Sickerwasser“ (SiWaP) motivierten die Entwicklung des Programms PCSiWaPro. Innerhalb des Vorhabens sollte die Möglichkeit geschaffen werden, mit geringem Aufwand eine modellgestützte Sickerwasserprognose unter Berücksichtigung der Forschungsergebnisse aus SiWaP durchführen zu können. Kommerziell verfügbare Software blieb dabei außen vor, da die Implementierung eigener Prozesse, Datenbanken und Parameter damit nicht möglich ist. Gleichzeitig war eine komplexe Betrachtung der ablaufenden Prozesse erforderlich sowie die Dokumentation der Ein- und Ausgabedaten für eine entsprechende Nachweispflicht. Dies führte zur Entwicklung einer grafischen Benutzeroberfläche (GUI) mit einem Assistenten, der den Anwender in 5 sequenziell ablaufenden Schritten zu einem physikalisch begründeten Ergebnis führt (Protokoll). Alle notwendigen Eingaben werden dazu mit sinnvollen Werten vorbelegt und bei Änderung durch den Nutzer auf Plausibilität geprüft. Gleichzeitig sollte die Funktionalität nicht auf die Möglichkeiten des Assistenten beschränkt bleiben und dem erfahrenen Modellierer alle Optionen der numerischen Simulation bereitstellen. Die Dokumentation der Ein- und Ausgabedaten wird dabei durch die Verwendung von Datenbanken sichergestellt. Für den Einsatz in Ingenieurbüros, Behörden oder auch international war die GUI mehrsprachig zu implementieren. Diese Anforderungen begründeten die Entwicklung eines Simulationssystems, um den Wasserhaushalt und Stofftransport in ungesättigten Böden auch unter komplexen Bedingungen berechnen zu können. Das aus dem zuvor genannten BMBF-Verbundvorhaben SiWaP entstandene Programm PCSiWaPro war wesentlicher Bestandteil nachfolgender Forschungsvorhaben, deren Ergebnisse in die weitere Entwicklung des Programms einflossen und dessen Anwendungsgebiete außerhalb der Sickerwasserprognose erweiterten. So sind erforderliche Eingangsdaten wie bodenhydraulische und Stofftransportparameter oft mit Unsicherheiten behaftet oder können nur in Wertebereichen gefasst werden. Um derartige Unschärfen auch in den Berechnungsergebnissen von numerischen Simulationen ausweisen zu können, wurde die Fuzzy-Set-Theorie verwendet, die eine Zuordnung der Unsicherheiten über sogenannte α-Schnitte ermöglicht. Für jeden unscharfen Parameter kann dessen Schwankungsbreite definiert und in der Simulation berücksichtigt werden. Die Ausweisung der Unschärfen im Ergebnis erfolgt unter Angabe des sich ergebenden Minimums und Maximums der berechneten Größe (Druckhöhe, Konzentration). Anhand verschiedener Beispielanwendungen werden die in der Arbeit vorgestellten Problemstellungen durch Einsatz von PCSiWaPro behandelt. Die Arbeit gibt ebenso einen Ausblick auf weiterführenden Forschungs- und Entwicklungsbedarf, der sich aus den in der Arbeit erzielten Ergebnissen und Betrachtungen ableiten lässt.:Abbildungsverzeichnis Tabellenverzeichnis Abkürzungsverzeichnis Symbolverzeichnis 1 Einleitung 2 Wasserhaushaltsberechnung in variabel gesättigten porösen Medien 2.1 Zugrundeliegende Gleichung 2.2 Numerische Lösung 3 Transport- und Umsetzungsprozesse 3.1 Erhaltungsgleichung 3.2 Transportprozesse 3.3 Umsetzungsprozesse 3.4 Basisgleichung für den Stofftransport in PCSiWaPro 3.5 Numerische Lösung 4 Entwicklung des Programms PCSiWaPro 4.1 Softwarearchitektur 4.2 Datenbankkonzept 4.3 Benutzeroberfläche für das Preprocessing 4.4 Ergebnisvisualisierung und Postprocessing 4.5 Parallelisierung des Rechenkernels 4.6 Dual-Porosität nach DURNER 4.7 Strömungsrandbedingung als zeitvariable Polygonfunktion 4.8 Berücksichtigung von Unsicherheiten in den Eingangsdaten 5 Anwendungsbeispiele 5.1 Deichdurchströmung 5.2 Modellgestützte Sickerwasserprognose mit unscharfen Eingangsdaten 5.3 Test der Parallelisierung am synthetischen Beispiel 5.4 Zusammenfassung Anwendungsbeispiele 6 Zusammenfassung und Ausblick 7 Literaturverzeichnis 8 Anhang / The soil zone, often referred to as root zone, aeration zone or unsaturated zone in the literature, is characterized by variably saturated conditions and is of particular importance in many disciplines. From the groundwater point of view, it is a zone for protection and buffering of environmental processes at the surface. Penetrating hazardous substances can be retarded or even completely decayed due to the transport, degradation and sorption processes which occur and thus, can prevent a contamination of the groundwater. In order to estimate potential threats to the groundwater based on a contaminated site or a damage, a leachate forecast is required in Germany according to the Federal Soil Protection Act (BBodSchG) and the Federal Soil Protection Ordinance (BBodSchV). The unsaturated zone takes on the function of the source and transport term for the pollutant. The source term function is used to describe the temporal discharge behavior of pollutants from the contaminant source with the leachate, the transport term describes the action path in the soil from the top of the site to the groundwater surface. The requirements and tasks of the BMBF-funded research project 'Prognosis of Pollutant Infiltration into Groundwater with Leachate' (“Prognose des Schadstoffeintrags in das Grundwasser mit dem Sickerwasser”) (SiWaP) motivated the development of the PCSiWaPro program. Within the project, the possibility should be created to be able to carry out a model-based leachate forecast with little effort, taking into account the research results from the SiWaP project. Commercially available software had to be left out, since the implementation of new processes, databases and parameters is not possible. At the same time, a total consideration of the complex processes taking place was necessary, as was the documentation of the input and output data to provide evidence. This led to the development of a graphical user interface (GUI) with an assistant that leads the user in 5 sequential steps to a physically based result including a protocol. All necessary input data are pre-assigned with useful values and checked for plausibility when changed by the user. At the same time, the functionality should not be limited to the possibilities of the assistant and the GUI must provide all available options of a numerical simulation to advanced users. The documentation of the input and output data is ensured by using databases. The GUI provides multiple languages for use in engineering offices, authorities or international projects. These requirements justified the development of a simulation system to be able to calculate the water balance and solute transport in unsaturated soils even under complex conditions. The PCSiWaPro program, emerged from the BMBF joint project SiWaP mentioned above, was an integral part of subsequent research projects, the results of which were incorporated into the further development of the program and expanded its fields of application outside of the leachate forecast. Required input data such as soil hydraulic and solute transport parameters are often subject to uncertainties or can only be captured in value ranges. In order to show such blurring in the calculation results of numerical simulations, the fuzzy set theory was used, which enables the uncertainties to be assigned using so-called α-cuts. The fluctuation range for each uncertain parameter can be defined individually and considered in the simulation. The blurring in the result is indicated by specifying the resulting minimum and maximum of the calculated quantity (pressure level, concentration). Using various sample applications, the problems presented in the thesis are dealt with by using PCSiWaPro. The thesis also gives an outlook on further research and development perspectives, which are derived from the results achieved in this thesis and the demands from the daily practice.:Abbildungsverzeichnis Tabellenverzeichnis Abkürzungsverzeichnis Symbolverzeichnis 1 Einleitung 2 Wasserhaushaltsberechnung in variabel gesättigten porösen Medien 2.1 Zugrundeliegende Gleichung 2.2 Numerische Lösung 3 Transport- und Umsetzungsprozesse 3.1 Erhaltungsgleichung 3.2 Transportprozesse 3.3 Umsetzungsprozesse 3.4 Basisgleichung für den Stofftransport in PCSiWaPro 3.5 Numerische Lösung 4 Entwicklung des Programms PCSiWaPro 4.1 Softwarearchitektur 4.2 Datenbankkonzept 4.3 Benutzeroberfläche für das Preprocessing 4.4 Ergebnisvisualisierung und Postprocessing 4.5 Parallelisierung des Rechenkernels 4.6 Dual-Porosität nach DURNER 4.7 Strömungsrandbedingung als zeitvariable Polygonfunktion 4.8 Berücksichtigung von Unsicherheiten in den Eingangsdaten 5 Anwendungsbeispiele 5.1 Deichdurchströmung 5.2 Modellgestützte Sickerwasserprognose mit unscharfen Eingangsdaten 5.3 Test der Parallelisierung am synthetischen Beispiel 5.4 Zusammenfassung Anwendungsbeispiele 6 Zusammenfassung und Ausblick 7 Literaturverzeichnis 8 Anhang