Processes of the Earth’s surface occur at different scales of time and intensity. Climate in particular determines the activity and seasonal development of vegetation. These dynamics are predominantly driven by temperature in the humid mid-latitudes and by the availability of water in semi-arid regions. Human activities are a modifying parameter for many ecosystems and can become the prime force in well-developed regions with an intensively managed environment. Accounting for these dynamics, i.e. seasonal dynamics of ecosystems and short- to long-term changes in land-cover composition, requires multiple measurements in time. With respect to the characterization of the Earth surface and its transformation due to global warming and human-induced global change, there is a need for appropriate data and methods to determine the activity of vegetation and the change of land cover. Space-borne remote sensing is capable of monitoring the activity and development of vegetation as well as changes of the land surface. In many instances, satellite images are the only means to comprehensively assess the surface characteristics of large areas. A high temporal frequency of image acquisition, forming a time series of satellite data, can be employed for mapping the development of vegetation in space and time. Time series allow for detecting and assessing changes and multi-year transformation processes of high and low intensity, or even abrupt events such as fire and flooding. The operational processing of satellite data and automated information-extraction techniques are the basis for consistent and continuous long-term product generation. This provides the potential for directly using remote-sensing data and products for analyzing the land surface in relation to global warming and global change, including deforestation and land transformation. This study aims at the development of an advanced approach to time-series generation using data-quality indicators. A second goal focuses on the application of time series for automated land-cover classification and update, using fractional cover estimates to accommodate for the comparatively coarse spatial resolution. Requirements of this study are the robustness and high accuracy of the approaches as well as the full transferability to other regions and datasets. In this respect, the developments of this study form a methodological framework, which can be filled with appropriate modules for a specific sensor and application. In order to attain the first goal, time-series compilation, a stand-alone software application called TiSeG (Time Series Generator) has been developed. TiSeG evaluates the pixel-level quality indicators provided with each MODIS land product. It computes two important data-availability indicators, the number of invalid pixels and the maximum gap length. Both indices are visualized in time and space, indicating the feasibility of temporal interpolation. The level of desired data quality can be modified spatially and temporally to account for distinct environments in a larger study area and for seasonal differences. Pixels regarded as invalid are either masked or interpolated with spatial or temporal techniques. / Prozesse an der Erdoberfläche finden auf verschiedenen Intensitätsskalen und in unterschiedlichen Zeiträumen statt. Dabei steuert das Klima die saisonale Aktivität der Vegetation, welche in den humiden Mittelbreiten hauptsächlich durch die Temperatur bestimmt wird. In semi-ariden Gebieten hingegen ist die Verfügbarkeit von Wasser als Haupteinflussfaktor für das Pflanzenwachstum zu betrachten. Andererseits greift auch der Mensch modifizierend in das Ökosystem ein. Dies gilt insbesondere für die stark besiedelten und intensiver genutzten Räume der Erde, in denen die Umwelt nahezu ausschließlich durch den Menschen gesteuert wird. Zur Beurteilung dieser Dynamiken, sowohl der natürlichen saisonalen Muster als auch der kurz- bis langfristigen Änderungen der Landschaft, ist die Aufnahme einer Vielzahl von Messungen über eine längere Periode erforderlich. Insbesondere im Zusammenhang mit der Charakterisierung der Landoberfläche und deren Veränderung im Rahmen der Erderwärmung aber auch des wachsenden Einflusses des Menschen auf die Umwelt besteht somit ein Bedarf an geeigneten Daten und Methoden zur Bestimmung der jährlichen Aktivität von Vegetationseinheiten und der wiederholbaren Kartierung der Landoberfläche. Die Satellitenfernerkundung ist in der Lage, durch Messung von Strahlung die Aktivität der Vegetation zu bestimmen sowie die Klassifikation der Landoberfläche abzuleiten. In vielen Fällen sind Satellitenaufnahmen die einzige Möglichkeit, große Flächen der Erde umfassend und einheitlich zu beurteilen. Dabei kann durch eine Vielzahl aufeinander folgender Aufnahmen, d.h. eine Zeitreihe aus Satellitendaten, die Entwicklung der Vegetation in Raum und Zeit beobachtet werden. Zeitreihen bieten das Potential, Veränderungen der Landoberfläche über mehrere Jahre zu dokumentieren und somit Prozesse sowohl hoher als auch niedriger Intensität abzuleiten. Neben diesen gerichteten Veränderungen können auch plötzliche Ereignisse, wie z.B. Hochwasser oder Brände, mit Zeitreihen erfasst und in Bezug auf normale Bedingungen ausgewertet werden. Insbesondere die operationelle Prozessierung der Satellitendaten und die automatisierte Ableitung von Informationen bilden die Basis für konsistente und kontinuierliche Produkte über längere Zeiträume. Somit besteht das Potential, die Ergebnisse direkt in die Erforschung der Landoberfläche und deren Veränderung, z.B. durch die Erderwärmung, Walddegradation, oder die Nutzung vormals natürlicher Flächen, einzubinden. Diese Dissertation befasst sich mit der Entwicklung von Methoden zur Zeitreihengenerierung unter Verwendung der Qualitätsindikatoren einzelner Aufnahmen. Ein zweites Ziel der Arbeit ist die Anwendung der optimierten Zeitreihen zur automatisierten und reproduzierbaren Kartierung der Landoberfläche, wobei unscharfe Klassifikationsverfahren zur genaueren Charakterisierung der räumlich nur grob aufgelösten Daten eingesetzt werden. Damit erfordert diese Arbeit sowohl Robustheit der eingesetzten Methoden als auch eine hohe Genauigkeit der Ergebnisse. Ebenso maßgeblich ist die Übertragbarkeit der Verfahren, einerseits auf verschiedene Regionen als auch auf verschiedene Datensätze. Daher sind die hier vorgenommenen Entwicklungen als ein Rahmen zu verstehen, der je nach Sensor oder Anwendung mit verschiedenen Modulen besetzt werden kann. Zum Erreichen des ersten Zieles, der Zeitreihengenerierung, wurde das eigenständige Softwareprodukt TiSeG (Time Series Generator) entwickelt. TiSeG dient der Auswertung der Qualitätsindikatoren die mit jedem MODIS-Produkt für terrestrische Applikationen zur Verfügung gestellt werden. Dabei werden in Bezug auf die Generierung von Zeitreihen zwei Indizes der Datenverfügbarkeit ermittelt: erstens die Anzahl der ungültigen Pixel und zweitens die längste zeitliche Datenlücke. Beide Indizes werden räumlich und zeitlich dargestellt und geben so dem Bearbeiter die Information, ob mit den aktuellen Qualitätsangaben die Generierung einer sinnvollen Zeitreihe durch zeitliche Dateninterpolation möglich ist. Die Qualitätseinstellungen können sowohl zeitlich als auch räumlich angepasst werden. Eine zeitliche Änderung kann beispielsweise für bestimmte Jahreszeiten sinnvoll sein. Räumlich unterschiedliche Qualitätseinstellungen eignen sich für größere Untersuchungsgebiete mit differenzierten physisch-geographischen Charakteristika. Als ungültig betrachtete Pixel können durch einen Fehlwert maskiert oder durch zeitliche und räumliche Interpolation neu errechnet werden. Das zweite Ziel dieser Arbeit ist die automatische Klassifikation von Zeitreihen. Hierzu wurde ein modulares Verfahren der überwachten Klassifikation entwickelt. Aufgrund der groben räumlichen Auflösung von MODIS-Daten erschien es besonders wichtig, ein unscharfes Verfahren aufzubauen, das die Heterogenität der Klassen in vielen Räumen besser abbilden kann. Dabei wurde der Klassenanteil eines jeden Pixel ermittelt. Die Schlüsselmodule zur erfolgreichen Durchführung der Klassifikation waren eine Multiskalenanalyse und die geeignete Auswahl von Merkmalen und Stichproben zum Trainieren des Klassifikators. Der eigentliche Klassifikationsschritt wurde durch eine Erweiterung des Entscheidungsbaumklassifikators durchgeführt. Diese Erweiterung kann in den bestehenden Rahmen von „random forest“ und „bagging“ (ein Akronym für „bootstrap aggregation“) eingeordnet werden. Jedoch wurden die dort angewendeten Verfahren in dieser Arbeit zu einem deutlich strategischen Vorgehen modifiziert, d.h. es wurde auf das Prinzip der Zufallsauswahl und der Wiederverwendung von Stichproben (Ziehen mit Zurücklegen) verzichtet. In Bezug auf die anfangs geschilderten Anforderungen, Robustheit, Genauigkeit und Übertragbarkeit, kann an dieser Stelle festgestellt werden, dass sich die in dieser Arbeit entwickelten Methodiken als geeignet erwiesen haben. Insbesondere die Übertragbarkeit auf andere Regionen und Daten war eine große Herausforderung, da hierdurch kein zusätzliches a priori Wissen außer den Trainingsdaten zum überwachten Klassifizieren benutzt werden konnte. Die regionale Übertragbarkeit ist für mehrere Untersuchungsräume mit sehr unterschiedlichen physisch-geographischen Eigenschaften demonstriert worden. Obwohl TiSeG derzeit nicht direkt auf andere Datensätze außer MODIS angewendet wird, bildet die zugrunde liegende Idee der Auswertung von Qualitätsdaten zur Zeitreihengenerierung sowie der entwickelte Rahmen der Software die Möglichkeit der Erweiterung, z.B. auf Daten von MERIS und auf das zukünftige VIIRS-Instrument. Des Weiteren ist das grundlegende Konzept der Qualitätsauswertung in einem separaten Prozessor für AVHRR NDVI Daten zu einem vollautomatischen, schrittweise interpolierenden Verfahren erweitert worden. Das in dieser Arbeit vorgestellte modulare Klassifikationsverfahren erfordert keine besonderen Eingangsdaten, wie z.B. bestimmte MODIS Zeitreihen. Dies wurde durch den Gebrauch unterschiedlicher Eingangsdaten zur Generierung der Maßzahlen bei der Sensitivitätsanalyse bestätigt. Damit ist der Prozess weder auf MODIS-Zeitserien noch auf Zeitreihen generell beschränkt. Der gesamte automatische Klassifikationsprozess ist datengesteuert, sowohl was die zu klassifizierenden Daten als auch die Trainingsdaten angeht. Die Unschärfe in den Ergebnissen ermöglicht die detaillierte Auswertung der Klassenzusammensetzung, was ein besonders wichtiger Aspekt bei grob aufgelösten Datenprodukten und deren Anpassungsfähigkeit auf andere Anwendungen ist.

Weltweit sind Trockengebiete in ständiger Veränderung, verursacht durch natürliche klimatische Schwankungen und oftmals durch Prozesse der Landdegradation. Auch weisen die meisten semi-ariden Naturräume eine große räumliche Heterogenität auf, hervorgerufen durch ein kleinräumiges Mosaik aus Gräsern, kleineren Sträuchern und Bereichen offenliegenden Bodens. Die Dichte der Vegetation wird primär vom pflanzenverfügbaren Wasser bestimmt, aber auch der Entwicklungs- und Degradationszustand der Böden sowie anthropogen bedingte Faktoren spielen hierbei eine Rolle. Zur Charakterisierung und Kartierung der Vegetation sowie zur Bewertung des Bodenerosionsrisikos und des Degradationszustands hat sich die Erhebung der Bedeckungsgrade von vitaler, photosynthetisch aktiver Vegetation (PV), von abgestorbener oder zeitweise vertrockneter und somit nicht photosynthetisch aktiver Vegetation (NPV) sowie von offenliegendem Boden als zweckmäßig herausgestellt. Die Nutzung der Fernerkundung für diese Aufgabe erfolgt zumeist nur für kleinmaßstäbige Kartierungen und – im Falle von Multispektralsensoren – unter Vernachlässigung nicht-photosynthetisch aktiver Vegetation. Die räumliche Variabilität der Vegetation-Boden-Mosaike liegt oftmals in der Größenordnung von wenigen Metern und somit unterhalb des räumlichen Auflösungsvermögens von Fernerkundungssystemen. Um dennoch die verschiedenen Anteile innerhalb eines Pixels identifizieren und quantifizieren zu können, sind Methoden der Subpixel-Klassifikation notwendig. In dieser Arbeit wird eine Methodik zur verbesserten und automatisierbaren Ableitung von Bodenbedeckungsgraden in semi-ariden Naturräumen vorgestellt. Hierzu wurde ein Verfahren zur linearen spektralen Entmischung in Form einer Multiple Endmember Spectral Mixture Analysis (MESMA) entwickelt und umgesetzt. Durch diese Methodik kann explizit die spektrale Variabilität von Vegetation und Boden in das Mischungsmodellmiteinbezogen werden, und quantitative Anteile für die funktionalen Klassen PV, NPV und Boden innerhalb eines Pixels erfasst werden. Durch die räumliche Kartierung der verwendeten EM wird weiterhin eine thematische Klassifikation erreicht. Die hierfür benötigten Informationen können – wie im Falle der Spektren reiner Materialien (EM-Spektren) – aus den Bilddaten selbst abgeleitet werden, oder können – wie ein Geländemodell und die Information über den Scanwinkel – im Zuge der Vorprozessierung aus weiteren Datenquellen erzeugt werden. Hinsichtlich der automatisierten EM-Ableitung wird eine zweistufige Methodik eingesetzt, welche auf einer angepassten Version des Sequential Maximum Angle Convex Cone (SMACC)-Verfahrens sowie der Analyse einer ersten Entmischungsiteration basiert. Die Klassifikation der gefundenen potentiellen EM erfolgt durch ein merkmalsbasiertes Verfahren. Weiterhin weisen nicht-photosynthetisch aktive Vegetation und Boden eine hohe spektrale Ähnlichkeit auf. Zur sicheren Trennung kann die Identifikation schmaler Absorptionsbanden dienen. Zu diesen zählen beispielsweise die Absorptionsbanden von Holozellulose und – je nach Bodentyp – Absorptionsbanden von Bodenmineralen. Auch die spektrale Variabilität der Klassen PV und NPV erfordert zur sicheren Unterscheidung die Verwendung biophysikalisch erklärbarer Merkmale im Spektrum. Hierzu zählen unter anderem die Stärke der Chlorophyll-Absorption, die Form und Lage der ’RedEdge’ und das Auftreten von Holozellulosebanden. Da diese spektrale Information bei herkömmlichen Entmischungsansätzen nicht berücksichtigt wird, erfolgt überwiegend eine Optimierung der Gesamtalbedo, was zu einer schlechten Trennung der Klassen führen kann. Aus diesem Grund wird in der vorliegenden Arbeit der MESMA-Ansatz dahingehend erweitert, dass spektrale Information in Form von identifizierten und parametrisierten Absorptionsbanden in den Entmischungsprozess mit einfließt und hierdurch das Potential hyperspektraler Datensätze besser genutzt werden kann. Auch wird in einer zusätzlichen Entmischungsiteration die räumliche Nachbarschaft betrachtet, um insbesondere die Verwendung des sinnvollsten Boden-EMs zu gewährleisten. Ein zusätzliches Problemfeld stellt die numerische Lösung des überbestimmten und oftmals schlecht konditionierten linearen Mischungsmodells dar. Hierzu kann durch die Verwendung des BVLS-Algorithmus und des Ausschlusses kritischer EM-Kombinationen eine numerisch stabile Lösung gefunden werden. Um die oftmals immense Rechenzeit von MESMA-Verfahren zu verkürzen, besteht die Möglichkeit einer iterativen EM-Auswahl und somit die Vermeidung einer Lösung des Mischungssystems durch Berechnung aller EM-Kombinationen (’Brute-Force’-Ansatz). Ein weiterer wichtiger Punkt ist die explizite pixelweise Angabe zur Zuverlässigkeit der Entmischungsergebnisse. Dies erfolgt auf Basis des Mischungsmodells selbst, durch den Vergleich zu empirischen Regressionsmodellen, durch die Berücksichtigung des lokalen Einfallswinkels sowie durch die Integration von Qualitätsangaben der Ausgangsdaten. Um das Verfahren systematisch und unter kontrollierten Bedingungen zu verifizieren und um den Einfluss verschiedener externer Parameter sowie die typischen Genauigkeiten auf einer breiten Datenbasis zu ermitteln, wird eine Simulationskette zur Erzeugung synthetischer Mischungen erstellt. In diese Simulationen fließen Feldspektren von Böden und Pflanzen verschiedener semi-arider Gebiete mit ein, um möglichst viele Fälle abdecken zu können. Die eigentliche Validierung erfolgt auf HyMap-Datensätzen des Naturparks ’Cabo de Gata’ in der andalusischen Provinz Almería sowie auf Messungen, die begleitend im Feld durchgeführt wurden. Hiermit konnte die Methodik auf ihre Genauigkeit unter den konkreten Anforderungen des Anwendungsbeispiels überprüft werden. Die erzielbare Genauigkeit dieser automatisierten Methodik liegt mit einem mittleren Fehler um rund 10% Abundanz absolut im selben Wertebereich oder nur geringfügig höher als die Ergebnisse publizierter manueller MESMA-Ansätze. Weiterhin konnten die typischen Genauigkeiten der Verifikation im Zuge der Validierung bestätigt werden. Den limitierenden Faktor des Ansatzes stellen in der Praxis fehlerhafte oder unvollständige EM-Modelle dar. Mit der vorgestellten Methodik ist somit die Möglichkeit gegeben, die Bedeckungsgrade quantitativ und automatisiert im Subpixelbereich zu erfassen. / Hyperspectral remote sensing has shown the potential to directly identify and parameterize absorption features and form parameters of spectral signatures. Using this information bio- and geo-physical properties of the sensed surfaces can be retrieved. Conventional spectral unmixing approaches disregard this information since all data values are equally incorporated despite their spectral position. The approach introduced in this thesis for the first time combines an automated spectral unmixing methodology with methods of spectral identification in order to use the full potential of spectroscopic data. An intrinsic property of semi-arid environments all over the world is large spatial heterogeneity caused by a small-scaled mosaic of grasses, small bushes and patches of bare soil. Vegetation cover is primarily limited by water, but also the status of the soil as well as anthropogenetic factors are important. The cover percentage is a standard parameter in vegetation mapping, for soil risk assessment, and for the assessment of land degradation. Therefore, it was found appropriate in a large number of studies to retrieve cover percentage for the functional classes of bare soil, vital photosynthetic active vegetation (PV) as well as dead or senescent non-photosynthetic Vegetation (NPV). For these tasks remote sensing is often applied only for small-scaled applications, and, when using multispectral sensors, commonly neglects non-photosynthetic vegetation. Also, the spatial extent of vegetation and soil patches is in many cases smaller than the spatial resolution of the sensor. Thus subpixel classifiers are required in order to identify and quantify the fractions within a single pixel. The aim of this study is to develop a methodology to provide accurate cover fractions for the mentioned functional classes at the patch scale. Therefore an automated methodology for quantitative subpixel classification based on linear spectral unmixing of hyperspectral data is introduced. This automated thematic processor, denoted muMESMA in the following, can be linked to automated pre-processing chains for hyperspectral data. All parts of muMESMA are verified using a series of spectral simulations in order to systematically test the consistency, general applicability and dependency on input variables of the approach. Also, typical error margins of the approach are estimated based on this large database of test scenarios. Validation of the thematic processor is then carried out for an application example based on multiple HyMap scenes of the ’Cabo de Gata’ natural park, province Almería, southeast Spain. The area is dominated by tussock grasses, small bushes and abandoned agricultural fields, and can be considered typical for semi-arid grasslands. To ensure an automated and widely applicable methodology, all required information should be retrieved from the imagery itself or from standard pre-processing products. In the case of spectral unmixing, the detection of all required spectral signatures of pure scene materials (denoted as Endmembers, EM) from the image itself is critical. For this purpose the Sequential Maximum Angle Convex Cone (SMACC) method is adjusted for automated usage. In order to increase the EM detection rate an additional approach is included in muMESMA. Potential new EM can be identified during a first unmixing iteration based on the assumption that pixels which can not be accurately modelled have a higher chance of being EM themselves. Using this two-step approach for endmember detection, about 70% of all EM spectra could be automatically retrieved in a series of simulated test cases. As for most unsupervised classification methods these two approaches do not include an identification of spectra. Thus muMESMA includes a decision-based spectral classifier for the classes of interest of this study. Based on diagnostic features in the spectra, the use of threshold values can be mostly avoided. The classification accuracy of this approach was 93% when using pure field spectra of the relevant classes. But this chance of ~ 1 in 10 for an inappropriate EM is critical since the accurate unmixing of most scenes requires more than 10 EM. Thus in order to guarantee suitable EM spectra and therefore accurate unmixing results a manual control of retrieved EM is highly recommended. After that the final EM set is selected based on different criteria and data tests, such as Endmember Average RMSE (EAR). To include the spectral variability of all classes in the unmixing process a large number of EM is required. But linear dependencies between spectrally similar classes significantly reduce the mathematical accuracy of the solution when using an increasing number of EMs. To avoid an ’ill-condition’ of the unmixing a Multiple Endmember Spectral Mixture Analysis (MESMA) approach is used. In this approach each image spectrum is linearly unmixed using all combinations of one spectrum for each of the classes PV, NPV and Soil. Afterwards, the model with the smallest modelling RMS-error is selected and the corresponding fractions are assigned to the pixel. For simulation scenarios the mean unmixing error can be reduced from 11.3% when using conventional manual unmixing to 9.6% using an automated MESMA. For the test area of ’Cabo de Gata’ the MESMA approach was able to accurately model the whole image with one comprehensive EM set despite changing spectral characteristics of surfaces (e. g., changes in vegetation status or composition, variation in soil type). In addition, a good condition of the mixing model without linear dependencies was achieved. The mean error over all test plots could be reduced by 0.6% abundance absolute. Using only the modelling RMS error as model selection criterion is not optimal since this measure is dominated by the overall albedo and small but important spectral features are often neglected. This might cause an increased inaccuracy of estimated ground cover fractions for the spectrally similar classes of NPV and bare soil, since the overall shape of the spectra is highly similar. Thus the residual spectrum (i. e., the part of the signal which cannot be modelled with the given EMs) is analysed using a knowledge-based approach in order to identify a meaningful EM model for each pixel. Using a combined model selection criterion based on this residual analysis, on modelling RMS and on physical constraints, the unmixing error can be further reduced to an average of 8.6% abundance absolute in simulation scenarios. The robustness of the approach is also increased, especially when introducing additional sources of error like partial shadowing, lichens, increased sensor noise or inexact EM. This increased stability is highly important for the HyMap imagery of Cabo de Gata since in all real-world applications the whole spectral variability cannot be completely represented by EMs. Thus the mean unmixing error can be significantly lowered from 15.1% abundance absolute when only optimizing the model RMS error to 11.6% abundance absolute when using the combined optimization criterion including residual analysis. The correlation between field measurements and unmixing results was significant (p < 0.005) with correlation coefficients between 0.76 and 0.86 (R2 values). This average accuracy of the automated muMESMA approach is in the same range as published values of manual field cover measurements, and is equal or only slightly higher than published results using manual MESMA approaches with similar thematic. Also, an additional unmixing iteration incorporating spatial neighbourhood information is included to ensure that a reasonable EM model is selected for each pixel. To reduce the calculation time of the three MESMA unmixing iterations muMESMA can iteratively propose EM models to avoid a ’brute force’ calculation. This methodology based on the residual analysis can significantly reduce calculation time by a factor of 7, but also reduces unmixing accuracy. Another important part of this automated thematic processor is a reliability measure of unmixing quality for every pixel. Based on the residual analysis, modelling error, comparison with empirical regression models and local incidence angle, pixels which are likely error-prone can be identified with the muMESMA approach. For areas with rough terrain the unmixing accuracy can further be increased. Local incidence angle effects significantly affect the accuracy of ground cover estimates for all measurements off-nadir where terrain is rough and vegetation patchy. When the line of sight from bare soil patches to the sensor is (partially) blocked by vegetation, the cover percentage of bare soil is underestimated, while fractional cover by vegetation is overestimated. Using the local incidence angle and empirical vegetation parameters (average height and basal area) a pixel-based empirical correction can be carried out. For the ’Cabo de Gata’ test site the average accuracy of ground increased by an additional 2% abundance absolute. Within this thesis it could be shown that the estimation of ground cover fractions for the functional classes PV, NPV and bare soil can be carried out with increased accuracy using hyperspectral data. With the exception of a recommended manual check the muMESMA approach is fully automated and can retrieve all required information from standardized input data generated by most preprocessing chains. The results of the verification carried out on a large database of field spectra from different semi-arid natural environments and different phenological conditions could be confirmed in the validation for the specific application case ’Cabo de Gata’ based on three multitemporal HyMap data sets. Thus the wide applicability of the approach can be assumed.

Verbleibende Unsicherheiten im Kohlenstoffhaushalt in Ökosystemen der hohen nördlichen Breiten können teilweise auf die Schwierigkeiten bei der Erfassung der räumlich und zeitlich hoch variablen Methanemissionsraten von Permafrostböden zurückgeführt werden. Methan ist ein global abundantes atmosphärisches Spurengas, welches signifikant zur Erwärmung der Atmosphäre beiträgt. Aufgrund der hohen Sensibilität des arktischen Bodenkohlenstoffreservoirs sowie der großen von Permafrost unterlagerten Landflächen sind arktische Gebiete am kritischsten von einem globalen Klimawandel betroffen. Diese Dissertation adressiert den Bedarf an Modellierungsansätzen für die Bestimmung der Quellstärke nordsibirischer permafrostbeeinflusster Ökosysteme der nassen polygonalen Tundra mit Hinblick auf die Methanemissionen auf regionalem Maßstab. Die Arbeit präsentiert eine methodische Struktur in welcher zwei prozessbasierte Modelle herangezogen werden, um die komplexen Wechselwirkungen zwischen den Kompartimenten Pedosphäre, Biosphäre und Atmosphäre, welche zu Methanemissionen aus Permafrostböden führen, zu erfassen. Es wird ein Upscaling der Gesamtmethanflüsse auf ein größeres, von Permafrost unterlagertes Untersuchungsgebiet auf Basis eines prozessbasierten Modells durchgeführt. Das prozessbasierte Vegetationsmodell Biosphere Energy Hydrology Transfer Model (BETHY/DLR) wird für die Berechnung der Nettoprimärproduktion (NPP) arktischer Tundravegetation herangezogen. Die NPP ist ein Maß für die Substratverfügbarkeit der Methanproduktion und daher ein wichtiger Eingangsparameter für das zweite Modell: Das prozessbasierte Methanemissionsmodell wird anschließend verwendet, um die Methanflüsse einer gegebenen Bodensäule explizit zu berechnen. Dabei werden die Prozesse der Methanogenese, Methanotrophie sowie drei verschiedene Transportmechanismen – molekulare Diffusion, Gasblasenbildung und pflanzengebundener Transport durch vaskuläre Pflanzen – berücksichtigt. Das Methanemissionsmodell ist für Permafrostbedingungen modifiziert, indem das tägliche Auftauen des Permafrostbodens in der kurzen arktischen Vegetationsperiode berücksichtigt wird. Der Modellantrieb besteht aus meteorologischen Datensätzen des European Center for Medium-Range Weather Forecasts (ECMWF). Die Eingangsdatensätze werden mit Hilfe von in situ Messdaten validiert. Zusätzliche Eingangsdaten für beide Modelle werden aus Fernerkundungsdaten abgeleitet, welche mit Feldspektralmessungen validiert werden. Eine modifizierte Landklassifikation auf der Basis von Landsat-7 Enhanced Thematic Mapper Plus (ETM+) Daten wird für die Ableitung von Informationen zu Feuchtgebietsverteilung und Vegetationsbedeckung herangezogen. Zeitserien der Auftautiefe werden zur Beschreibung des Auftauens bzw. Rückfrierens des Bodens verwendet. Diese Faktoren sind die Haupteinflussgrößen für die Modellierung von Methanemissionen aus permafrostbeeinflussten Tundraökosystemen. Die vorgestellten Modellergebnisse werden mittels Eddy-Kovarianz-Messungen der Methanflüsse validiert, welche während der Vegetationsperioden der Jahre 2003-2006 im südlichen Teil des Lena Deltas (72°N, 126°E) vom Alfred Wegener Institut für Polar- und Meeresforschung (AWI) durchgeführt wurden. Das Untersuchungsgebiet Lena Delta liegt an der Laptewsee in Nordostsibirien und ist durch Ökosysteme der arktischen nassen polygonalen Tundra sowie kalten kontinuierlichen Permafrost charakterisiert. Zeitlich integrierte Werte der modellierten Methanflüsse sowie der in situ Messungen zeigen gute Übereinstimmungen und weisen auf eine leichte Modellunterschätzung von etwa 10%. / Remaining uncertainties in the carbon budget of high latitude ecosystems are partly due to difficulties in assessing methane emission rates from permafrost soils the source strengths of which are highly variable in space and time. Methane is a globally abundant atmospheric trace gas that contributes significantly to the warming of the atmosphere. Due to the high sensitivity of the arctic soil carbon reservoir and the large surface area underlain by permafrost, arctic regions are most critically influenced by a changing climate. This dissertation addresses the need for modelling approaches to determine the source strength of northern Siberian permafrost affected wet polygonal tundra ecosystems with regard to methane emission on the regional scale. It presents a methodical structure wherein two process-based models are used to capture the complex interrelated processes between pedosphere, biosphere and atmosphere that lead to methane emission from permafrost soils on the regional scale. Upscaling of methane fluxes for a larger permafrost site is performed using results of a process-based model. The process-based vegetation model Biosphere Energy Transfer Hydrology Model (BETHY/DLR) is applied to estimate net primary productivity (NPP) of arctic tundra vegetation. NPP is parameterized as a measure for substrate availability for methane production and thus an important input parameter for the second model: the process-based wetland methane emission model is subsequently used to explicitly model methane fluxes for a given soil column, taking into account methanogenesis, methane oxidation and three different transport mechanisms, namely molecular diffusion, ebullition and plant-mediated transport through vascular plants. The methane emission model is modified for permafrost conditions by explicitly considering daily thawing of permafrost during the short arctic growing season. Model forcing consists of meteorological data sets obtained from the European Center for Medium-Range Weather Forecasts (ECMWF). Input data sets are validated against field measurements. Auxiliary input data for both models are derived from satellite imagery and validated by field spectral measurements. A modified land use/land classification (LULC) scheme based on Landsat-7 Enhanced Thematic Mapper Plus (ETM+) data is used to derive information on wetland distribution and vegetation cover. Time series of active layer thickness are used to describe thawing/freezing of soils. These parameters are key factors in modelling methane emissions from permafrost influenced tundra ecosystems. Validation of presented model results is performed using eddy covariance measurements of methane flux on the landscape scale carried out during the growing seasons 2003-2006 in the southern part of the Lena Delta (72°N, 126°E) by Alfred Wegener Institute for Polar and Marine Research (AWI). The Lena Delta study site is located at the Laptev Sea in northeast Siberia and is characterized by arctic wet polygonal tundra ecosystems and cold continuous permafrost. Timeintegrated values for modelled methane fluxes and in situ measurements compare reasonably well and indicate a moderate model underestimation of about 10%.

Im Rahmen der vorliegenden Arbeit wurden die Druck-Temperatur-Bildungsbedingungen metapelitischer Gesteine aus verschiedenen lithostratigraphischen Formationen des Spessart-Kristallins untersucht. Geologisch stellt das Spessart-Kristallin einen Teil der Mitteldeutschen Kristallinzone dar, die sich innerhalb des Orogens der mitteleuropäischen Varisziden am nördlichen Rand des Saxothuringikums erstreckt. Das wesentliche Ziel der Untersuchungen bestand darin, mittels phasenpetrologischer und geothermobarometrischer Methoden die Metamorphose-Entwicklung des Spessart-Kristallins zu rekonstruieren und in Beziehung zur geodynamischen Geschichte der Varisziden zu setzen. Der Schwerpunkt der Arbeiten lag auf den Staurolith-Glimmerschiefern der Mömbris-Formation. Darüber hinaus wurden Gesteine der Geiselbach-, Alzenau- und Elterhof-Formation einbezogen. Als Grundlage für die Phasenpetrologie wurden petrographische, mineralchemische und geochemische Untersuchungen durchgeführt. Prograde Metamophose-Zonen können im Spessart-Kristallin nicht kartiert werden. Die Protolithe der untersuchten Metasedimente stellten vermutlich häufig saure bis intermediäre Magmatite dar, für die Geiselbach- und Elterhof-Formation wohl auch quarzreiche Sedimente. Die geo¬chemischen Daten lassen für die Mömbris-, Alzenau- und Elterhof-Formation Grauwacken bis Pelite als sedimentäre Edukte der Metamorphite annehmen, die Gesteine der Geiselbach-Formation könnten auf Arkosen zurückgehen. Eine Ablagerung der sedimentären Edukte im Bereich eines Kontinen¬talen Inselbogens bis Aktiven Kontinentalrandes ist für die Mömbris- und Alzenau-Formation wahrscheinlich, für die Geiselbach- und Elterhof-Formation liegt kein eindeutiges Bild vor. Zur Abschätzung der Metamorphosebedingungen wurden verschiedene Phasendiagramme verwendet, die auf den metapelitischen Modellsystemen KMnFMASH (K2O-MnO-FeO-MgO-Al2O3-SiO2-H2O) und KFMASH (K2O-FeO-MgO-Al2O3-SiO2-H2O) basieren, insbesondere P-T-Pseudoschnitte und T-X-Schnitte. Weiterhin wurden konventionelle Geothermobarometer berechnet und Abschätzungen mittels intern-konsistenter thermodynamischer Datensätze vorgenommen. Die theoretischen Grundlagen dieses phasenpetrologischen Ansatzes werden kurz erläutert. Für den Metamorphose-Höhepunkt der Gesteine ergaben sich Temperaturen im Bereich von ca. 600 - 615 °C und Drucke um 6.5 - 8 kbar. Diese Daten weisen eine recht gute Übereinstimmung zu den bisher in der Literatur bekannten Werten für das Spessart-Kristallin auf. Im Anschluß an die amphibolitfazielle Metamorphose wurden die Gesteine mehr oder minder stark retrograd überprägt. Anzeichen für eine polymetamorphe Entwicklung dieses Teils der Mitteldeutschen Kristallinzone liegen nicht vor. Die rekonstruierten P-T-Pfade bzw. P-T-Pfad-Segmente dokumentieren eine recht einheitliche metamorphe Entwicklung im Uhrzeigersinn („clockwise“) und weisen auf eine Barrow-type Metamorphose hin. Die P-T-Pfade der meisten Proben zeigen einen charakteristischen Verlauf mit einer Phase nahezu isothermaler Dekompression. Demgegenüber konnte für einige Disthen-führende Proben ein etwas flacherer P-T-Pfad mit einer offenbar geringfügig stärker Temperatur-betonten Entwicklung differenziert werden. Das Metamorphose-Maximum ist für diese Gesteine durch Temperaturen von ca. 620 - 630 °C und Drucke von etwa 6 - 8 kbar gekennzeichnet. Damit wird eine leichte Zunahme des Metamorphosegrades nach Süden innerhalb der Mömbris-Formation, die verschiedentlich vermutet worden war, nachgewiesen. Die neu erarbeiteten Pfade sind aufgrund des methodischen Ansatzes, der die Zusammensetzung und Mineralparagenese der jeweiligen Probe berücksichtigt, im Vergleich zu früheren Arbeiten deutlich besser abgesichert. Sie dokumentieren erstmals in dieser Form die Druck-Temperatur-Geschichte des Spessart-Kristallins. Die P-T-Pfade lassen auf eine relativ schnelle Versenkung der Gesteine bei einem recht niedrigen geothermischen Gradienten und eine anschließende rasche Heraushebung aus einer Tiefe von etwa 25 - 28 km auf etwa 15 - 18 km bei einer eher geringen Temperaturabnahme schließen. Die damit für das Spessart-Kristallin dokumentierte Entwicklung fügt sich gut in das aktuelle geotektonische Modell einer Kollision eines passiven Kontinentalrandes mit einem kontinentalen Bogen ein und steht in Analogie zur derzeit gängigen Vorstellung, die Mitteldeutsche Kristallinschwelle repräsentiere einen variszischen aktiven Plattenrand. / Metamorphic Evolution of the Spessart Crystalline Complex, Central European Variscides: studies in phase-petrology, mineral-chemistry, and geochemistry of metapelitic rocks. – 211 pp. + appendix, Würzburg, 2008. Within the Central European Variscides the Spessart Crystalline Complex (SCC) is part of the Mid German Crystalline Rise which nowadays is interpreted as an active Variscan margin. This thesis aimed at the metamorphic evolution of the SCC and looked for possible contributions to the geodynamic evolution of the area in question. Studies, therefore, concentrated on metapelitic rocks from different lithostratigraphic formations of the SCC, especially on the staurolite-micaschists of the Mömbris-formation, but samples from the (very badly exposed) Geiselbach-, Alzenau-, and Elterhof-formations were analyzed as well. Emphasis was mainly given to phase-petrology and to geothermobarometry, preceded by studies in petrology, mineral chemistry and geochemistry the results of which are documented in the appendices. The evaluation of the metamorphic conditions used different phase-diagrams, i.e. P-T-pseudosections and T-X-sections, based on the metapelitic model-systems KMnFMASH (K2O-MnO-FeO-MgO-Al2O3-SiO2-H2O) and KFMASH (K2O-FeO-MgO-Al2O3-SiO2-H2O). Further estimations base on geothermobarometers and internally-consistent thermodynamic datasets. The theoretical background of these methods is shortly explained. The metamorphic maximum is given at temperatures of about 600 - 615 °C and pressures at some 6.5 - 8 kbar. These data are in good accordance with published evidence for the Spessart Crystalline Complex. Metamorphism under conditions of amphibolite facies was followed by a retrograde overprint. There is no indication of polymetamorphic development in this part of the Mid German Crystalline Rise. The P-T-paths reconstructed in this thesis document a more or less uniform clockwise metamorphic development and, thus, indicate a Barrow-type metamorphism. The P-T-paths of most of the samples show a characteristic course with a phase of nearly isothermal decompression. Some kyanite-bearing samples with a metamorphic peak at temperatures of 620 - 630 °C and pressures of some 6 - 8 kbar are characterized by a somewhat less steep P-T-path. In summary, these new P-T-paths document for the first time the pressure-temperature-history of the Spessart Crystalline Complex this way. The pressure-temperature-paths indicate a relatively quick burial of the rocks at a comparatively low geothermal gradient and a following rapid exhumation from a depth of some 25 - 28 km to 15 - 18 km with a rather low decrease of temperature. This development of the Spessart-Crystalline-Complex fits well into the actual geodynamic model of a collision of a passive margin with a continental arc and agrees with the current conception of the Mid German Crystalline Rise.

This work presents a new method to measure model independent viscosities of inhomogeneous materials at high temperatures. Many mechanisms driving volcanic eruptions are strongly influenced by the viscous properties of the participating materials. Since an eruption takes place at temperatures at which these materials (predominantly silicate melts) are not completely molten, typically inhomogeneities, like e.g. equilibrium and non-equilibrium crystals, are present in the system. In order to incorporate such inhomogeneities into objective material parameters the viscosity measurement is based on a rotational viscometer in a wide gap Couette setup. The gap size between the two concentric cylinders was designed as large as possible in order to account for the inhomogeneities. The emerging difficulties concerning the model independent data reduction from measured values to viscosities are solved using an appropriate interpolation scheme. The method was applied to a material representative for the majority of volcanic eruptions on earth: a typical continental basaltic rock (Billstein/Rhön/Germany). The measured viscosities show a strong shear rate dependency, which surprises, because basaltic melt has been, until now, assumed to behave as a Newtonian fluid. Since a non-Newtonian material shows a very different relaxation behavior in the Couette motion compared to a Newtonian one (which, ultimately, does not show any), and a strong relaxation signal was recorded during viscosity measurements, the equations of Couette motion were investigated. The time dependent stress distribution in a material due to a quasi step-like velocity change at the inner Couette radius (i.e. the spindle) was considered. The results show that a material combining a linear shear modulus and a Newtonian viscosity -- a Maxwell material -- cannot quantify the relaxation behavior. This could be considered as a hint, that the widely used Maxwell relaxation times cannot be applied as a 1:1 mapping from microscopic considerations to macroscopic situations. / Die vorliegende Arbeit beschreibt eine neue Methode zur modellunabhängigen Messung von Viskositäten bei hohen Temperaturen. Viele der Mechanismen, welche vulkanischer Aktivität zugrunde liegen, werden stark durch die viskosen Eigenschaften der beteiligten Materialien beeinflusst. Die eruptierten Materialien -- zum überwiegenden Teil Silikatschmelzen -- sind bei Eruptionstemperatur nicht komplett geschmolzen. Deshalb sind Gleichgewichts- und Nichtgleichgewichtskristalle in den betrachteten Systemen vorhanden. Um diese Inhomogenitäten in objektive Materialparameter einzubeziehen, basiert die vorgestellte Viskositätsmessung auf auf einem Rotationsviskosimeter in einer wide gap''-Anordnung. Die Spaltbreite zwischen den beiden konzentrischen Zylindern wurde so groß wie möglich gemacht um Inhomogenitäten zu berücksichtigen. Die aufkommenden Schwierigkeiten bezüglich der modellunabhängigen Bestimmung der Viscositäten aus den gemessenen Daten wurden mit einer geeigneten Interpolationsmethode gelöst. Mit dieser Methode wurden die Viskositäten eines, für die Mehrheit vulkanischer Eruptionen auf der Erde typisches Material gemessen: eines kontinentalen Basaltes aus Billstein (Rhön, Deutschland). Die gemessenen Viskositäten zeigen bei konstanter Temperatur eine starke Abhängigkeit von der Deformationsrate. Dies überrascht, da basaltische Schmelzen bis heute bei vergleichbaren Temperaturen als Newtonsche Flüssigkeiten betrachtet wurden. Da ein nicht-Newtonsches Material, im Vergleich mit einem Newtonschen, ein deutlich anderes Relaxationsverhalten aufweist (das Newtonsche zeigt ultimativ keine Relaxation), und da ein deutliches Relaxationssignal während der Viskositätsmessung gemessen wurde, wurden die Bewegungsgleichungen der Couette Bewegung untersucht. Die zeitabhängige Spannungeverteilung in einem Material, verursacht durch eine quasi-stufenartige Geschwindigkeitsänderung am inneren Couette-Radius (d. h. am Drehkörper des Viskosimeters) wurde betrachtet. Die Ergebnisse zeigen, dass ein Material, welches ein linear elastisches Schermodul und eine newtonsche Viskosität kombiniert -- ein Maxwell-Material -- das Relaxationverhalten quantitativ nicht beschreiben kann. Dies könnte als Hinweis betrachtet werden, dass die weitverbreiteten Maxwell-Relaxationszeiten nicht 1:1 von mikroskopischen Betrachtungen auf makroskopische Situationen angewendet werden können.

The present study concerened with taxonomy and palaeoecology of the Cenomanian-Turonian macrobenthic fauna which includes bivalves, gastropods, echinoids, and coral. In addtion, cephalopods are also taken in consideration. 144 taxa are identified and systematically described. Palaeoecological and taphonomic anylsis of the statistically sampled macrobenthos are also discussed. The biostratigraphic sequences along the Cenomanian-Turonian rocks were carried out on the basis of ammonites and other macrobenthic fauna such as corals and bivalves. In order to reconstruct benthic association, 41 statistically sampled were subjected to cluster ananlysis by using Past Programm (Hammer et al., 2001). 10 association and three assemblages were described in order to reconstruct the different depositional enviroments. / Diese Arbeit behandelt die Taxonomie und Palökologie von Makro-Invertebraten im Bereich Cenoman-Turon im östlichen Sinai, Ägypten

Der Klimawandel und insbesondere die globale Erwärmung gehören aktuell zu den größten Herausforderungen an Politik und Wissenschaft. Steigende CO2-Emissionen sind hierbei maßgeblich für die Klimaerwärmung verantwortlich. Ein regulierender Faktor beim CO2-Austausch mit der Atmosphäre ist die Vegetation, welche als CO2-Senke aber auch als CO2-Quelle fungieren kann. Diese Funktionen können durch Analysen der Landbedeckungsänderung in Kombination mit Modellierungen der Kohlenstoffbilanz quantifiziert werden, was insbesondere von aktuellen und zukünftigen politischen Instrumenten wie CDM (Clean Development Mechanism) oder REDD (Reducing Emissions from Deforestation and Degradation) gefordert wird. Vor allem in Regionen mit starker Landbedeckungsänderung und hoher Bevölkerungsdichte sowie bei geringem Wissen über die Produktivität und CO2-Speicherpotentiale der Vegetation, bedarf es einer Erforschung und Quantifizierung der terrestrischen Kohlenstoffspeicher. Eine Region, für die dies in besonderem Maße zutrifft, ist Westafrika. Jüngste Studien haben gezeigt, dass sich einerseits die Folgen des Klimawandels und Umweltveränderungen sehr stark in Westafrika auswirken werden und andererseits Bevölkerungswachstum eine starke Änderung der Landbedeckung für die Nutzung als agrarische Fläche bewirkt hat. Folglich sind in dieser Region die terrestrischen Kohlenstoffspeicher durch Ausdehnung der Landwirtschaft und Waldrodung besonders gefährdet. Große Flächen agieren anstelle ihrer ursprünglichen Funktion als CO2-Senke bereits als CO2-Quelle. [...] / Global warming associated with climate change is one of the greatest challenges of today's world. One regulating factor of CO2 exchange with the atmosphere is the vegetation cover. Measurements of land cover changes in combination with modeling of the carbon balance can therefore contribute to determining temporal variations of CO2 sources and sinks, which is an essential necessity of existing and prospective political instruments like CDM (Clean Development Mechanism) or REDD (Reducing Emissions from Deforestation and Degradation). The need for quantifiable terrestrial carbon stocks is especially high for regions, where rates of land cover transformation and population density are high and knowledge on vegetation productivity is low. One region which is characterized by these criteria is West Africa. Therefore, carbon stocks in this region are seriously endangered by land cover change like the expansion of agriculture and forest logging. Large areas already act as carbon sources on a yearly basis instead of their previous function as carbon sink. Since only a few studies have analyzed the terrestrial carbon stocks in Africa and especially regional analysis in West Africa are missing, the following study focuses on regional scale modeling of the actual terrestrial carbon stocks. Additionally, the potential carbon stocks of unmanaged land cover and the potential for CO2 payments have been analyzed in this work. To quantify and assess carbon fluxes as well as the loss of carbon, net primary productivity of vegetation has been modeled, based on the plants characteristics to fix carbon from the atmosphere during photosynthesis. Modeling vegetation dynamics and net primary productivity has been realized by using MODIS 250m time series for semi-humid and semi-arid savanna ecosystems in West Africa. This study aimed to quantify CO2 exchanges of the Savanna regions in the Volta basin by applying and adapting the Regional Biomass Model (RBM). The RBM was developed by Jochen Richters (2005) at a resolution of 1000m for the Namibian Kaokoveld. In this study the model was optimized to the scale of 232m to consider the heterogeneous landscape in West Africa (RBM+). New input parameters with higher accuracies and resolution were generated instead of using the global standard products. The most important parameters for the modeling are FPAR and the fractional cover of herbaceous and woody vegetation. To enhance the MODIS FPAR product, linear interpolation and downscaling algorithms were applied. The main objective of the downscaling is a better representation of the finely scattered vegetation by the 232m resolution FPAR. The second optimized parameter, the fractional cover of herbaceous and woody vegetation was represented by the Vegetation Continuous Fields product (VCF) from MODIS in the originally version of the RBM. This global product reflects the vegetation structure of West Africa poorly, since few high resolution training data is available for this region, and the dynamic savanna vegetation can hardly be classified by not regionally adapted methods. Additionally, the data is only available with 500m resolution. Therefore, in this study a new product with 232m resolution was developed which represents the spatial heterogeneity well and, due to the regional adaptation, shows higher accuracies. The percentage cover of woody and herbaceous vegetation and bare soil on 232m MODIS data was calculated in a multi scale approach. Based on very high resolution data, represented by Quickbird and Ikonos with 0.6-4m resolution, and high resolution data from Landsat with 30m resolution, the percentage coverage was estimated for representative focus regions. These classifications were used as a training data set to determine the percentage coverage on the 232m scale with MODIS time series for the whole study region. Based on these optimized and adapted input parameters, the net primary productivity was modeled. Data from a meteorological station and an Eddy-Covariance-Flux allowed a detailed validation of the input parameters and of the model results. The model led to good results as it only overestimated the net primary productivity for the two analyzed years 2005 and 2006 by 8.8 and 2.0 %, respectively. The second aim of the study was an analysis of the potential for long term terrestrial carbon sinks. Classifications of the actual and of the potential land cover were calculated for this analysis. Considering the overall long time CO2 fixation behavior of trees, which depends on their age, longterm carbon stocks for 100 years were simulated. As carbon fixing could be paid by emission trading, which is in future depending on the political Post-Kyoto programs, potential alternative income was calculated with different price scenarios for the three countries. A comparison with the gross domestic products of these countries and with developing aid, showed the significance of CO2 trading in this region.

Die Bevölkerungsexplosion in der Region von Curitiba während der letzten Jahre verursachte eine große Zunahme des Wasserbedarfs; die zunehmend unkontrollierte Besiedlung stellt dabei ein großes Problem für die Umwelt dar. Die größte Bedrohung für die Wasserversorgung dieser Region ist die urbane Invasion in Gebiete, wo die Herkunft der Wasserressourcen liegen. Diese Invasion geht langsam aber stetig voran und bedroht kostbare und nicht ersetzbare Ressourcen. Vor diesem Hintergrund wurde ein Gebiet in der direkten Nähe der Großstadt Curitiba als Studienobjekt ausgewählt. In diesem Gebiet, dem Iraí-Becken, wurde während der Untersuchungszeit ein Trinkwasserspeicher geplant und gebaut. Es besteht die große Gefahr, dass das Iraí-Reservoir kontaminiert wird, obwohl das engere Gebiet rundum den See geschützt werden soll. Die Verschmutzungsgefahr geht hauptsächlich von zwei Nebenflüssen aus, die durch mehr oder weniger besiedeltes Gebiet strömen. Im Arbeitsgebiet befinden sich Brunnen, die der Trinkwasserversorgung dienen. Um die negativen Folgen einer möglichen Verschmutzung des Reservoirs abschätzen zu können, wurde ein Grundwasserfließmodell erstellt. Die erforderliche Wasserbilanz und die räumliche Verteilung der Verschmutzungsempfindlichkeit wurde mit dem hydrologischen Modell "MODBIL" abgeschätzt. Weitere Methoden zur Abschätzung der Verschmutzungs-empfindlichkeit wurden angewandt, um die differierenden Ergebnisse der angewendeten unterschiedlichen Methoden mit einander vergleichen und bewerten zu können. Mit dem kalibrierten Grundwasserfließmodell ist mit der gegebenen hydraulischen Situation vor und nach der Konstruktion des Reservoirs, ein einfaches Particle Tracking Transport Modell eingesetzt worden, um mit unterschiedlichen Szenarien die Beeinflussung vom Reservoirwasser auf das Grundwasser zu simulieren. / The explosive expansion of the population of the Metropolitan Region of Curitiba raised a high increase in the demand for water resources and the uncontrolled settlement poses a large problem for the environment. The greatest menace to the water supply sources of this region is the urban occupation (invasion) into the areas that contain these resources. This occupation continues with its slow, silent, although progressive march, threatening precious and irreplaceable resources. From this background an area in the direct vicinity north-east of Curitiba has been studied. In this area a drinking water reservoir was constructed in the time that the study took place in the Iraí-basin. The Iraí-reservoir even though an area around the lake will be protected may be polluted by two tributaries which flow through more or less densely populated areas. In the study area on the same time wells have been constructed. To estimate what the impact may be from the possibly polluted reservoir on the aquifer a groundwater flow model has been constructed. On the same time to estimate the water balance and the spatial distribution of pollution vulnerability the hydrological model MODBIL has been used. Also other methods have been used to estimate the pollution vulnerability to make a comparison and because none of the methods takes every aspect into account. With the calibrated groundwater flow model for the situation before the construction of the Iraí-reservoir and after its construction, simple particle tracking transport models are constructed as scenarios how the water of the aquifer may be influenced.

The Skeleton Coast forms part of the Atlantic coastline of NW Namibia comprising several ephemeral rivers, which flow west-southwest towards the Atlantic Ocean. The area is hyper-arid with less than 50 mm average annual rainfall and a rainfall variability of 72%. Therefore, the major catchment areas of the rivers are about 100-200 km further inland in regions with relatively high annual rainfall of 300-600 mm. The coastal plain in the river downstream areas is characterized by a prominent NNW trending, 165 km long belt of 20-50 m high, locally compound, barchanoid and transverse dunes. This dune belt, termed Skeleton Coast Erg, starts abruptly with a series of barchans and large compound dunes 15 km north of the Koigab River and extends from 2-5 km inland sub-parallel to the South Atlantic margin of NW Namibia over a width of 3-20 km. As the SSE-NNW trending dune belt is oriented perpendicular to river flow, the dunefield dams and interacts with the west-southwestward flowing ephemeral river systems. This study focused on three main topics: 1) investigation and classification of the Koigab Fan, 2) the investigation of the Cenozoic succession in the Uniabmond area and 3) comparative studies of fluvio-aeolian interaction between five ephemeral rivers and the Skeleton Coast Erg. Sedimentological and geomorphological investigations show that the Koigab Fan represents a yet undocumented type of a braided fluvial fan system, which operates in an arid climatic, tropical latitude setting, is dominated by ephemeral mixed gravel/sand braided rivers, lacks significant vegetation on the fan surface, has been relatively little affected by human activity, is a perfect study site for recording various types of fluvio-aeolian interaction and thereby acts additionally as a model for certain Precambrian and Early Palaeozoic fan depositional systems deposited prior to the evolution of land plants. The Cenozoic succession in the Uniabmond area consists of three major unconformity-bounded units, which have been subdivided into the Red Canyon, the Whitecliff, and the Uniabmond Formation. The Tertiary Red Canyon Fm. is characterized by continental reddish sediments documenting an alluvial fan and braided river to floodplain depositional environment. The Whitecliff Fm. displays a wide variety of continental and marine facies. This formation provides the possibility to examine fluvio-aeolian interactions and spectacular, steep onlap relationships towards older sediments preserved in ancient seacliffs. The Whitecliff Fm. has been subdivided into four sedimentary cycles, which resulted from sea level changes during the Plio- to Middle Pleistocene. The following Uniabmond Fm. provides a unique insight into the depositional history of the NW Namibian coast during the Last Pleistocene glacial cycle. The formation has been subdivided into four units, which are separated by unconformities controlled by sea level changes. Unit 1 represents deposits of an Eemian palaeo-beach. The overlying Units 2-4 build up the sedimentary body of the Uniab Fan, again a braided river dominated fan, which is nowadays degraded and characterized by deeply incised valleys, deflation surfaces and aeolian landforms. The Uniabmond Fm. is overlain by the dunes of the Skeleton Coast Erg, whose development is related to the Last Glacial Maximum (LGM). The damming of river flow by aeolian landforms has been previously recognized as one of several principal types of fluvio-aeolian interaction. Five ephemeral rivers (from S to N: Koigab, Uniab, Hunkab, Hoanib, Hoarusib), which variously interact with the Skeleton Coast Erg, were chosen for the purpose of this study to consider the variability of parameters within these fluvio-aeolian systems and the resulting differences in the effectiveness of aeolian damming. The fluvio-aeolian interactions between the rivers and the dune field are controlled by the climate characteristics and the geology of the river catchment areas, the sediment load of the rivers, their depositional architecture, the longitudinal river profiles as well as the anatomy of the Skeleton Coast Erg. Resulting processes are 1) aeolian winnowing of fluvially derived sediments and sediment transfer into and deposition in the erg; 2) dune erosion during break-through resulting in hyperconcentrated flow and intra-erg mass flow deposits; 3) the development of extensive flood-reservoir basins caused by dune damming of the rivers during flood; 4) interdune flooding causing stacked mud-pond sequences; and 5) the termination of the erg by more frequent river floods. / Das hyperaride Gebiet der Skelettküste im NW Namibias ist durch eine Vielzahl ephemerer Flußsysteme und einen über 165 km langen, NNW-streichenden Dünengürtel, den Skelettküsten-Erg, charakterisiert. Der Erg verläuft nahezu orthogonal zu den westwärts, zum Atlantik entwässernden ephemeren Flußsystemen, deren Haupteinzugsgebiete etwa 150-300 km östlich der Küste liegen. Bei ihrem unregelmäßigen Abkommen auf dem Weg zum Atlantischen Ozean müssen die meisten dieser westwärts entwässernden Riviere den Dünengürtel passieren. Aus diesem Grund wurde dieses Gebiet ausgewählt, um vergleichende sedimentologische und geomorphologische Studien durchzuführen, die sich auf folgende frei Schwerpunkte konzentrierten: 1) die Untersuchung und Klassifizierung des Koigab Fans und dessen Rolle für den Sedimenteintrag in den Skelettküsten Erg; 2) die detaillierte Aufnahme und Untersuchung der känozoischen Sedimente im Mündungsbereich des Uniab Riviers und deren Rolle bei der Interpretation von Plio- und Pleistozänen Meeresspiegelschwankungen; und 3) die Erlangung eines grundlegenden Verständnis der zugrundeliegenden Kontrollparameter für die Variabilität von fluvio-äolischen Interaktionen und der daraus resultierenden Sedimentkörpergeometrien. Im Mündungsbereich des Koigab Riviers befindet sich der Koigab Fan. Dieser ist mit einer Länge von 15 km und einer Breite von 21 km der größte aktive Fan der Skelettküste, der durch einen der westsüdwestwärts entwässernden ephemeren Flußsystemen gebildet wurde. Sedimentologische und geomorphologische Studien an dessen Ablagerungen ergaben, dass es sich bei ihm um einen bis jetzt noch nicht dokumentierten Fan-Typ innerhalb der bestehenden Fan-Klassifikationen handelt. Der Koigab Fan repräsentiert einen von verflochtenen ephemeren Flußrinnen aufgebauten Fan, der keine Anzeichen von Schuttstömen (debris flows) oder mäandrierenden Flußrinnen aufweist. Der Koigab-Fan stellt daher ein Endglied innerhalb des von Stanistreet & McCarthy (1993) zur Fan-Klassifikation aufgestellten Dreieckdiagramms dar. Dem Koigab-Fan kommt eine besondere Bedeutung zu, da er in Bezug auf Größe und Gradienten eine Brücke zwischen der Gruppe kleiner, steil abfallender Fans und den mit nur geringem Gefälle versehenen Mega-Fans schlägt. Trotz seiner Position in den Tropen liegt der Koigab-Fan in einem ariden Gebiet und ist durch eine ausgeprägte Vegetationsarmut gekennzeichnet. Sedimentologische Untersuchungen an seinen Ablagerungen können daher auch als Analogbeispiele für vegetationsfreie früh-proterozoische und früh-paläozoische Fan-Systeme verwendet werden. Zudem stellt der Koigab Fan eine wesentliche Sedimentquelle für den Skelettküsten Erg dar, was mit Hilfe von Schwermineral- und Korngrößenanalysen nachgewiesen werden konnte. Die im Mündungsbereich des Uniab Riviers anstehenden Sedimente stellen die am vollständigsten erhaltene und am besten aufgeschlossene känozoische Sedimentabfolge entlang der Küste NW-Namibias dar. Im Rahmen des Promotionsprojektes wurde diese Abfolge erstmals vollständig sedimentologisch bearbeitet und anhand markanter Diskordanzen in drei Formationen unterteilt: die Red Canyon, Whitecliff und Uniabmond Formation. Die tertiäre Red Canyon Fm. ist durch Rotsedimente charakterisiert, die in einem alluvial-fluvialen Milieu abgelagert wurden. Die überlagernde Whitecliff Fm. wird aus vier sedimentären Einheiten aufgebaut. Die älteste und die jüngste Einheit enthalten marine und äolische Ablagerungen, in denen fossile Strandniveaus überliefert sind. Diese können herangezogen werden um die Whitekliff Fm. in das späte Pliozän bis ins mittlere Pleistozän einzustufen. Mit den Sedimenten der Uniabmond Fm. ist höchst wahrscheinlich der letzte glaziale Zyklus des Pleistozäns überliefert. Die jüngsten Sedimente der Uniabmond Fm. werden direkt von dem Skelettküsten Erg überlagert, dessen Bildung mit dem letzten glazialen Meeresspiegel-Tiefstand zusammenfällt. Das Aufstauen von Flüssen durch Dünen stellt einen Haupttyp fluvio-äolischer Interaktion dar. Fünf ephemeren Flusssystemen, die unterschiedlich mit dem Skelettküsten Erg in Interaktion treten, wurden innerhalb dieser Studie eingehen untersucht, um die Bandbreite fluvio-äolischer Interaktionen und deren zugrundeliegenden Kontrollparameter zu erfassen. Die fluvio-äolischen Interaktionen zwischen den Flüssen und dem Dünenfeld werden kontrolliert durch Klimacharakteristika, die Geologie des Einzugsgebietes, die Sedimentfracht der Flüsse, ihre Ablagerungsarchitektur und die Anatomie des Skelettküsten Ergs. Die zu beobachtenden Prozesse sind 1) äolische Verfrachtung von fluvial antransportiertem Sediment und dessen Ablagerung in dem Erg; 2) die Erosion und Aufarbeitung von Dünen während Durchbruchsereignissen; 3) die Bildung ausgedehnter Stauseen; 4) die Flutung von Interdünenbereichen mit der Ablagerung von Ton-Silt-Abfolgen; und 5) die Termination des Ergs durch häufiger auftretende Fluten.

Sand ramps have been (and still are) neglected in geomorphological research. Only recently any awareness of their potential of being a major source of palaeoenvironmental information, thanks to their multi-process character, has been developed. In Namibia, sand ramps were terra incognita. This study defines, classifies and systematizes sand ramps, investigates the formative processes and examines their palaeoenvironmental significance. The study region is located between the coastal Namib desert and the Great Escarpment, between the Tiras Mountains to the north and the Aus area to the south. Two lines of work were followed: geomorphological and sedimentological investigations in the field, assisted by interpretation of satellite images, aerial photographs and topographic maps, and palaeopedological and sedimentological analytical work in the laboratory. Two generations of sand ramps could be identified. The older generation, represented by a single sand ramp within the study region, is characterized by the presence of old basal sediments. The bulk of the sand ramps is assigned to the young generation, which is divided into three morpho-types: in windward positions voluminous ramps are found, in leeward positions low-volume ramps exist, either of very high or very low slope angle. The most distinct characteristic of sand ramp sediments is their formation by interacting aeolian deposition and fluvial slope wash. The last period of deposition, which shaped all the entire young sand ramps, but also the upper part of the old ramp, is suggested to have occurred after c. 40 ka BP, implying a highly dynamic climatic system during that time, with seasonal aridity and low-frequency, but high-intensity rainfall. A phase of environmental stability followed, most likely around 25 ka BP, supporting growth of vegetation, stabilization and consolidation of the sediments as well as soil formation. Subsequently, the profile was truncated and a desert pavement formed, under climatic conditions comparable to those of the present semi-desert. The ramps were then largely cut off from the bedrock slopes, implying a change towards higher ecosystem variability. As the final major process, recent and modern aeolian sands accumulated on the upper ramp slopes. A luminescence date for the recent sand places their deposition at about 16 ka BP, close to the Last Glacial Maximum. Regarding the source of the sands, a local origin is proposed. For the sand ramp of the old generation the "basic cycle" of initial deposition, stabilization and denudation occurred twelve times, including a phase of calcrete and/or root-cast formation in each of them, adding up to around 60 changes in morphodynamics altogether. At least nine of these cycles took place between 105 ka BP and the LGM, indicating that the general cooling trend during the Late Pleistocene was subject to a high number of oscillations of the environmental conditions not identified before for southern Namibia. Due to the high resolution obtained by the study of sand ramp sediments, but also due to the very special situation of the study area in a desert margin, 100 km from the South Atlantic and in the transition zone between summer and winter rainfall, correlation with stratigraphies (of mostly lower resolution) established for different regions in southern Africa did not appear promising. In conclusion, sand ramps generally serve as a valuable tool for detailed deciphering of past morphodynamics and thereby palaeoenvironmental conditions. For south-west Namibia, sand ramps shed some more light on the Late Quaternary landscape evolution. / In der geomorphologischen Forschung haben Sandrampen bislang wenig Beachtung gefunden, erst in jüngster Zeit ist auf ihr Potential als Speicher von Paläoumweltbedingungen hingewiesen worden. In Namibia waren Sandrampen terra incognita. Die vorliegende Studie definiert, klassifiziert und systematisiert Sandrampen in Südwest-Namibia, sie entschlüsselt ihre Bildungsprozesse und untersucht ihre paläoklimatische Bedeutung. Das Arbeitsgebiet liegt zwischen der küstenparallelen Namib und der Großen Randstufe, zwischen den Tirasbergen im Norden und der Gegend um Aus im Süden. Methodisch standen geomorphologische und sedimentologische Untersuchungen im Gelände im Vordergrund, unterstützt durch die Interpretation von Satellitenbildern, Luftbildern und topographischen Karten. Zahlreiche Substratproben wurden im Labor paläopedologisch und sedimentologisch analysiert. Es konnten zwei Sandrampen-Generationen identifiziert werden. Die ältere Generation, im Arbeitsgebiet durch nur eine Rampe vertreten, zeichnet sich durch das Vorkommen alter basaler Sedimente aus. Alle anderen Rampen sind der jüngeren Generation zuzuordnen und lassen sich in drei Morpho-Typen unterteilen: In Luv-Positionen finden sich voluminöse Rampen, während im Lee geringmächtige Rampen von entweder extrem steiler oder sehr geringer Hangneigung ausgebildet sind. Das auffälligste Merkmal der Sandrampen ist ihre Bildung durch die Interaktion von äolischer Deposition und Hangspülung. Die letzte Depositionsphase, in der die gesamten Körper der jungen Sandrampen abgelagert und die basalen Sedimente der alten Generation überlagert wurden, hat vermutlich nach 40 ka BP stattgefunden. Dies impliziert ein hochdynamisches klimatisches System zu dieser Zeit, mit saisonaler Aridität und seltenen, aber intensiven Regenfällen. Es folgte eine Phase der Ökosystemstabilität, vermutlich um 25 ka BP, in der es zu Vegetationsentwicklung und Bodenbildung sowie zu Stabilisierung und Konsolidierung der Sedimente kam. Eine anschließende Profilkappung mit Wüstenpflasterbildung geschah dann unter ähnlichen klimatischen Bedingungen wie in der heutigen Halbwüste. Danach wurden die Sandrampen fast überall von den Hängen abgeschnitten, was einer Änderung zu höherer Variabilität im Ökosystem zuzusprechen wäre. Die jüngste wesentliche Überprägung bestand in der Ablagerung subrezenter und moderner äolischer Sande in den obersten Bereichen der Sandrampen. Lumineszenzdatierungen stellen die subrezenten Sande ins letzte Hochglazial (~ 16 ka BP). Dem Sandrampensand wird eine lokale Herkunft zugesprochen. Für die Sandrampe der älteren Generation wiederholte sich der Zyklus von Deposition, Stabilisierung und Denudation insgesamt zwölfmal, inklusive je einer Phase von Kalkkrusten- und/oder Wurzelpseudomorphosenbildung. Insgesamt sind rund 60 Prozesswechsel dokumentiert. Mindestens neun dieser Zyklen verliefen zwischen 105 ka BP und dem letzten Hochglazial. Der generelle Abkühlungstrend während des Spätpleistozäns war also einer bedeutenden Anzahl von Schwankungen unterworfen, die bislang für Südnamibia nicht bekannt waren. Aufgrund der hohen zeitlichen Auflösung, die sich aus den Sandrampensedimenten erschließt, aber auch durch die besondere Lage des Arbeitsgebietes in einer Wüstenrandregion, 100 km vom Südatlantik entfernt und in der Übergangszone von Sommer- und Winterniederschlag, wurden Korrelationen mit Stratigraphien (von meist geringerer Auflösung), die für andere Regionen im südlichen Afrika aufgestellt worden sind, als wenig sinnvoll erachtet. Die vorliegende Arbeit zeigt, dass sich Sandrampen generell sehr gut zur Entschlüsselung paläomorphodynamischer Prozesse eignen und damit wesentlich zur Rekonstruktion von Paläoumweltbedingungen beitragen. Für Südwest-Namibia liefern Sandrampen neue, detaillierte Informationen zur spätquartären Landschaftsgeschichte.