Exploring physical processes related to past climate proxies: lake sediments and stable water isotopes

Proxy data are essential for the investigation of climate variability on time scales larger than the historical meteorological observation period. The potential value of a proxy depends on our ability to understand and quantify the physical processes that relate the corresponding climate parameter and the signal in the proxy archive. These processes can be explored under present-day conditions. In this thesis, both statistical and physical models are applied for their analysis, focusing on two specific types of proxies, lake sediment data and stable water isotopes.rnIn the first part of this work, the basis is established for statistically calibrating new proxies from lake sediments in western Germany. A comprehensive meteorological and hydrological data set is compiled and statistically analyzed. In this way, meteorological times series are identified that can be applied for the calibration of various climate proxies. A particular focus is laid on the investigation of extreme weather events, which have rarely been the objective of paleoclimate reconstructions so far. Subsequently, a concrete example of a proxy calibration is presented. Maxima in the quartz grain concentration from a lake sediment core are compared to recent windstorms. The latter are identified from the meteorological data with the help of a newly developed windstorm index, combining local measurements and reanalysis data. The statistical significance of the correlation between extreme windstorms and signals in the sediment is verified with the help of a Monte Carlo method. This correlation is fundamental for employing lake sediment data as a new proxy to reconstruct windstorm records of the geological past.rnThe second part of this thesis deals with the analysis and simulation of stable water isotopes in atmospheric vapor on daily time scales. In this way, a better understanding of the physical processes determining these isotope ratios can be obtained, which is an important prerequisite for the interpretation of isotope data from ice cores and the reconstruction of past temperature. In particular, the focus here is on the deuterium excess and its relation to the environmental conditions during evaporation of water from the ocean. As a basis for the diagnostic analysis and for evaluating the simulations, isotope measurements from Rehovot (Israel) are used, provided by the Weizmann Institute of Science. First, a Lagrangian moisture source diagnostic is employed in order to establish quantitative linkages between the measurements and the evaporation conditions of the vapor (and thus to calibrate the isotope signal). A strong negative correlation between relative humidity in the source regions and measured deuterium excess is found. On the contrary, sea surface temperature in the evaporation regions does not correlate well with deuterium excess. Although requiring confirmation by isotope data from different regions and longer time scales, this weak correlation might be of major importance for the reconstruction of moisture source temperatures from ice core data. Second, the Lagrangian source diagnostic is combined with a Craig-Gordon fractionation parameterization for the identified evaporation events in order to simulate the isotope ratios at Rehovot. In this way, the Craig-Gordon model can be directly evaluated with atmospheric isotope data, and better constraints for uncertain model parameters can be obtained. A comparison of the simulated deuterium excess with the measurements reveals that a much better agreement can be achieved using a wind speed independent formulation of the non-equilibrium fractionation factor instead of the classical parameterization introduced by Merlivat and Jouzel, which is widely applied in isotope GCMs. Finally, the first steps of the implementation of water isotope physics in the limited-area COSMO model are described, and an approach is outlined that allows to compare simulated isotope ratios to measurements in an event-based manner by using a water tagging technique. The good agreement between model results from several case studies and measurements at Rehovot demonstrates the applicability of the approach. Because the model can be run with high, potentially cloud-resolving spatial resolution, and because it contains sophisticated parameterizations of many atmospheric processes, a complete implementation of isotope physics will allow detailed, process-oriented studies of the complex variability of stable isotopes in atmospheric waters in future research.rn

Einfluss der Landnutzungsänderung auf die Variabilität von Grundwasserneubildungs- und Sickerwasserrate - Qualitative Grundwassergefährdungsabschätzung mit einem GIS -

Ziel der Arbeit ist die Entwicklung eines allgemeingültigen, im mittleren Maßstabsbereich durchgeführten, GIS-Bewertungsverfahrens zur Ermittlung der Grundwasserneubildungsrate unter Berücksichtigung einer langzeitlichen Landnutzungsänderung sowie die Darstellung eines Szenarios zur Grundwassergefährdungsabschätzung für ausgewählte rüstungsspezifi-sche organische Verbindungen. Das Untersuchungsgebiet befindet sich in einem während der beiden Weltkriege stark militär- und rüstungstechnisch genutzten Raum. Die Beurtei-lungsszenarien der Grundwassergefährdung erfolgen rein qualitativ.Die Abhängigkeit der Grundwasserneubildungsrate von der Landnutzungsänderung ist deut-lich zu erkennen. Mit zunehmender Flächenversiegelung, einher gehend mit veränderten Sickerwasserraten und des oberflächennahen Abflusses, reduziert sich die Grundwasser-neubildungsrate deutlich. Die Abschätzung der Grundwassergefährdung für die sprengstofftypischen Verbindungen erfolgt im ersten Schritt emissionsbezogen für die ungesättigte Zone unter Berücksichtigung der Standort- und Bodeneigenschaften und der physikalisch-chemischen Stoffeigenschaften. Im zweiten Schritt erfolgt eine immissionsbezogene Betrachtung der potenziellen Belastun-gen von Trinkwasserbrunnen.Trotz starker Generalisierungseffekte bietet dieser stoffspezifische Bewertungsansatz in Verbindung mit physikalischen Methoden für organische Verbindungen einen guten ersten Überblick über eine potenzielle Gefährdungen des Grundwassers. Die Datenbestände der entwickelten GIS-Anwendung sind leicht erweiterbar, so dass eine Abschätzung der potenzi-ellen Gefährdung auch für andere organische Stoffgruppen sowie auch eine Übertragung auf andere Untersuchungsgebiete möglich ist.

Polycyclic aromatic hydrocarbons (PAHs) in river bank soils

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous and found in the atmosphere, aquatic environment, sediments and soils. For environmental risk assessments and the allocation of the polluter it is important to know the PAH sources. PAH contamination sites are usually the result of anthropogenic processes. Three major sources are known: i) petroleum, including crude oil and its refined products and coals (petrogenic PAHs), ii) burning of organic matter (pyrogenic PAHs) and iii) transformation products of natural organic precursors present in the environment (diagenetic processes). In one case elevated PAH concentrations were found in river bank soils when building a retention area along the Mosel River. The source of the PAHs in this area was unclear and required the investigation of possible sources. To evaluate the PAH distribution along the Mosel River, a section of ~ 160 km along the river and a short section along the Saar River were investigated within this study. Concentrations of the Σ16 EPA PAHs were as high as 81 mg kg-1 dry weight (dw). Additionally, coal particles were identified in some soils, which originated from mining activities in the Saarland region. PAH distribution patterns of the 16 EPA PAHs suggest a mainly pyrogenic origin and in some cases a mixture of pyrogenic and petrogenic origin. For a comprehensive investigation five sampling sites were selected. Two sites were located before the confluence of the Mosel and Saar River, one site at the confluence and two sites after the confluence. The examination included typical forensic methods such as PAH distribution patterns of 45 PAHs (including alkylated PAHs), calculation of PAH ratios, determination of PAH alkyl homologues, n-alkanes, principal component analysis (PCA) and coal petrography. The results revealed a mainly pyrogenic source at sampling sites before the confluence of the two rivers. At and after the confluence, a mixture of pyrogenic and petrogenic inputs were present. With the help of coal petrography, coal derived particles could be identified in these soils. Therefore, coal was suggested to be the petrogenic source. It could be shown that sites with diffuse sources of contaminants, like the bank soils of the Mosel River, are difficult to characterize. As previously mentioned for detailed source identifications, the use of various forensic methods is essential. Determination of PAH alkyl homologue series, biomarkers and isotopes are often recommended. Source identification was evaluated using three different methods (i.e. PAH distribution patterns of an extended PAH spectrum, PAH ratios and analyses of n-alkanes). It was assessed if these methods were sufficient for the initial steps in identifying sources of PAHs in selected samples, and if they could be used for decision-making purposes. Point- and non-point sources were identified by applying the three methods and it could be shown that these relatively simple methods are sufficient in determining the primary source. In a last step of this study two soils (one before the confluence of the Mosel and Saar rivers and one after the confluence), and one sediment of the Mosel River were evaluated by investigating the mutagenic potential of the soils and the sediment with a fluctuation version of the Ames-test. The study showed that coal bearing soils at the Mosel River do not exhibit a greater mutagenic potential than other soils or sediments without coal particles.

Sustainable improvement of water networks : to achieve best decisions for rehabilitation plans using advanced GIS and statistical techniques

Many of developing countries are facing crisis in water management due to increasing of population, water scarcity, water contaminations and effects of world economic crisis. Water distribution systems in developing countries are facing many challenges of efficient repair and rehabilitation since the information of water network is very limited, which makes the rehabilitation assessment plans very difficult. Sufficient information with high technology in developed countries makes the assessment for rehabilitation easy. Developing countries have many difficulties to assess the water network causing system failure, deterioration of mains and bad water quality in the network due to pipe corrosion and deterioration. The limited information brought into focus the urgent need to develop economical assessment for rehabilitation of water distribution systems adapted to water utilities. Gaza Strip is subject to a first case study, suffering from severe shortage in the water supply and environmental problems and contamination of underground water resources. This research focuses on improvement of water supply network to reduce the water losses in water network based on limited database using techniques of ArcGIS and commercial water network software (WaterCAD). A new approach for rehabilitation water pipes has been presented in Gaza city case study. Integrated rehabilitation assessment model has been developed for rehabilitation water pipes including three components; hydraulic assessment model, Physical assessment model and Structural assessment model. WaterCAD model has been developed with integrated in ArcGIS to produce the hydraulic assessment model for water network. The model have been designed based on pipe condition assessment with 100 score points as a maximum points for pipe condition. As results from this model, we can indicate that 40% of water pipeline have score points less than 50 points and about 10% of total pipes length have less than 30 score points. By using this model, the rehabilitation plans for each region in Gaza city can be achieved based on available budget and condition of pipes. The second case study is Kuala Lumpur Case from semi-developed countries, which has been used to develop an approach to improve the water network under crucial conditions using, advanced statistical and GIS techniques. Kuala Lumpur (KL) has water losses about 40% and high failure rate, which make severe problem. This case can represent cases in South Asia countries. Kuala Lumpur faced big challenges to reduce the water losses in water network during last 5 years. One of these challenges is high deterioration of asbestos cement (AC) pipes. They need to replace more than 6500 km of AC pipes, which need a huge budget to be achieved. Asbestos cement is subject to deterioration due to various chemical processes that either leach out the cement material or penetrate the concrete to form products that weaken the cement matrix. This case presents an approach for geo-statistical model for modelling pipe failures in a water distribution network. Database of Syabas Company (Kuala Lumpur water company) has been used in developing the model. The statistical models have been calibrated, verified and used to predict failures for both networks and individual pipes. The mathematical formulation developed for failure frequency in Kuala Lumpur was based on different pipeline characteristics, reflecting several factors such as pipe diameter, length, pressure and failure history. Generalized linear model have been applied to predict pipe failures based on District Meter Zone (DMZ) and individual pipe levels. Based on Kuala Lumpur case study, several outputs and implications have been achieved. Correlations between spatial and temporal intervals of pipe failures also have been done using ArcGIS software. Water Pipe Assessment Model (WPAM) has been developed using the analysis of historical pipe failure in Kuala Lumpur which prioritizing the pipe rehabilitation candidates based on ranking system. Frankfurt Water Network in Germany is the third main case study. This case makes an overview for Survival analysis and neural network methods used in water network. Rehabilitation strategies of water pipes have been developed for Frankfurt water network in cooperation with Mainova (Frankfurt Water Company). This thesis also presents a methodology of technical condition assessment of plastic pipes based on simple analysis. This thesis aims to make contribution to improve the prediction of pipe failures in water networks using Geographic Information System (GIS) and Decision Support System (DSS). The output from the technical condition assessment model can be used to estimate future budget needs for rehabilitation and to define pipes with high priority for replacement based on poor condition. rn

Imaging and modeling of pore scale processes in porous media using X-ray computed tomography and lattice Boltzmann solver

In der Erdöl– und Gasindustrie sind bildgebende Verfahren und Simulationen auf der Porenskala im Begriff Routineanwendungen zu werden. Ihr weiteres Potential lässt sich im Umweltbereich anwenden, wie z.B. für den Transport und Verbleib von Schadstoffen im Untergrund, die Speicherung von Kohlendioxid und dem natürlichen Abbau von Schadstoffen in Böden. Mit der Röntgen-Computertomografie (XCT) steht ein zerstörungsfreies 3D bildgebendes Verfahren zur Verfügung, das auch häufig für die Untersuchung der internen Struktur geologischer Proben herangezogen wird. Das erste Ziel dieser Dissertation war die Implementierung einer Bildverarbeitungstechnik, die die Strahlenaufhärtung der Röntgen-Computertomografie beseitigt und den Segmentierungsprozess dessen Daten vereinfacht. Das zweite Ziel dieser Arbeit untersuchte die kombinierten Effekte von Porenraumcharakteristika, Porentortuosität, sowie die Strömungssimulation und Transportmodellierung in Porenräumen mit der Gitter-Boltzmann-Methode. In einer zylindrischen geologischen Probe war die Position jeder Phase auf Grundlage der Beobachtung durch das Vorhandensein der Strahlenaufhärtung in den rekonstruierten Bildern, das eine radiale Funktion vom Probenrand zum Zentrum darstellt, extrahierbar und die unterschiedlichen Phasen ließen sich automatisch segmentieren. Weiterhin wurden Strahlungsaufhärtungeffekte von beliebig geformten Objekten durch einen Oberflächenanpassungsalgorithmus korrigiert. Die Methode der „least square support vector machine” (LSSVM) ist durch einen modularen Aufbau charakterisiert und ist sehr gut für die Erkennung und Klassifizierung von Mustern geeignet. Aus diesem Grund wurde die Methode der LSSVM als pixelbasierte Klassifikationsmethode implementiert. Dieser Algorithmus ist in der Lage komplexe geologische Proben korrekt zu klassifizieren, benötigt für den Fall aber längere Rechenzeiten, so dass mehrdimensionale Trainingsdatensätze verwendet werden müssen. Die Dynamik von den unmischbaren Phasen Luft und Wasser wird durch eine Kombination von Porenmorphologie und Gitter Boltzmann Methode für Drainage und Imbibition Prozessen in 3D Datensätzen von Böden, die durch synchrotron-basierte XCT gewonnen wurden, untersucht. Obwohl die Porenmorphologie eine einfache Methode ist Kugeln in den verfügbaren Porenraum einzupassen, kann sie dennoch die komplexe kapillare Hysterese als eine Funktion der Wassersättigung erklären. Eine Hysterese ist für den Kapillardruck und die hydraulische Leitfähigkeit beobachtet worden, welche durch die hauptsächlich verbundenen Porennetzwerke und der verfügbaren Porenraumgrößenverteilung verursacht sind. Die hydraulische Konduktivität ist eine Funktion des Wassersättigungslevels und wird mit einer makroskopischen Berechnung empirischer Modelle verglichen. Die Daten stimmen vor allem für hohe Wassersättigungen gut überein. Um die Gegenwart von Krankheitserregern im Grundwasser und Abwässern vorhersagen zu können, wurde in einem Bodenaggregat der Einfluss von Korngröße, Porengeometrie und Fluidflussgeschwindigkeit z.B. mit dem Mikroorganismus Escherichia coli studiert. Die asymmetrischen und langschweifigen Durchbruchskurven, besonders bei höheren Wassersättigungen, wurden durch dispersiven Transport aufgrund des verbundenen Porennetzwerks und durch die Heterogenität des Strömungsfeldes verursacht. Es wurde beobachtet, dass die biokolloidale Verweilzeit eine Funktion des Druckgradienten als auch der Kolloidgröße ist. Unsere Modellierungsergebnisse stimmen sehr gut mit den bereits veröffentlichten Daten überein.