Refine
Year of publication
Document Type
- Doctoral Thesis (510) (remove)
Keywords
- Fernerkundung (15)
- climate change (14)
- Klimawandel (13)
- remote sensing (13)
- Anden (12)
- Andes (12)
- Erdbeben (12)
- Seismologie (10)
- Tektonik (10)
- Erosion (9)
Institute
- Institut für Geowissenschaften (510) (remove)
Recent large earthquakes put in evidence the need of improving and developing robust and rapid procedures to properly calculate the magnitude of an earthquake in a short time after its occurrence. The most famous example is the 26 December 2004 Sumatra earthquake, when the limitations of the standard procedures adopted at that time by many agencies failed to provide accurate magnitude estimates of this exceptional event in time to launch early enough warnings and appropriate response. Being related to the radiated seismic energy ES, the energy magnitude ME is a good estimator of the high frequency content radiated by the source which goes into the seismic waves. However, a procedure to rapidly determine ME (that is to say, within 15 minutes after the earthquake occurrence) was required. Here it is presented a procedure able to provide in a rapid way the energy magnitude ME for shallow earthquakes by analyzing teleseismic P‑waves in the distance range 20-98. To account for the energy loss experienced by the seismic waves from the source to the receivers, spectral amplitude decay functions obtained from numerical simulations of Greens functions based on the average global model AK135Q are used. The proposed method has been tested using a large global dataset (~1000 earthquakes) and the obtained rapid ME estimations have been compared to other magnitude scales from different agencies. Special emphasis is given to the comparison with the moment magnitude MW, since the latter is very popular and extensively used in common seismological practice. However, it is shown that MW alone provide only limited information about the seismic source properties, and that disaster management organizations would benefit from a combined use of MW and ME in the prompt evaluation of an earthquake’s tsunami and shaking potential. In addition, since the proposed approach for ME is intended to work without knowledge of the fault plane geometry (often available only hours after an earthquake occurrence), the suitability of this method is discussed by grouping the analyzed earthquakes according to their type of mechanism (strike-slip, normal faulting, thrust faulting, etc.). No clear trend is found from the rapid ME estimates with the different fault plane solution groups. This is not the case for the ME routinely determined by the U.S. Geological Survey, which uses specific radiation pattern corrections. Further studies are needed to verify the effect of such corrections on ME estimates. Finally, exploiting the redundancy of the information provided by the analyzed dataset, the components of variance on the single station ME estimates are investigated. The largest component of variance is due to the intra-station (record-to-record) error, although the inter-station (station-to-station) error is not negligible and is of several magnitude units for some stations. Moreover, it is shown that the intra-station component of error is not random but depends on the travel path from a source area to a given station. Consequently, empirical corrections may be used to account for the heterogeneities of the real Earth not considered in the theoretical calculations of the spectral amplitude decay functions used to correct the recorded data for the propagation effects.
The seismically active Alborz mountains of northern Iran are an integral part of the Arabia-Eurasia collision. Linked strike-slip and thrust/reverse-fault systems in this mountain belt are characterized by slow loading rates, and large earthquakes are highly disparate in space and time. Similar to other intracontinental deformation zones such a pattern of tectonic activity is still insufficiently understood, because recurrence intervals between seismic events may be on the order of thousands of years, and are thus beyond the resolution of short term measurements based on GPS or instrumentally recorded seismicity. This study bridges the gap of deformation processes on different time scales. In particular, my investigation focuses on deformation on the Quaternary time scale, beyond present-day deformation rates, and it uses present-day and paleotectonic characteristics to model fault behavior. The study includes data based on structural and geomorphic mapping, faultkinematic analysis, DEM-based morphometry, and numerical fault-interaction modeling. In order to better understand the long- to short term behavior of such complex fault systems, I used geomorphic surfaces as strain markers and dated fluvial and alluvial surfaces using terrestrial cosmogenic nuclides (TCN, 10Be, 26Al, 36Cl) and optically stimulated luminescence (OSL). My investigation focuses on the seismically active Mosha-Fasham fault (MFF) and the seismically virtually inactive North Tehran Thrust (NTT), adjacent to the Tehran metropolitan area. Fault-kinematic data reveal an early mechanical linkage of the NTT and MFF during an earlier dextral transpressional stage, when the shortening direction was oriented northwest. This regime was superseded by Pliocene to Recent NE-oriented shortening, which caused thrusting and sinistral strike-slip faulting. In the course of this kinematic changeover, the NTT and MFF were reactivated and incorporated into a nascent transpressional duplex, which has significantly affected landscape evolution in this part of the range. Two of three distinctive features which characterize topography and relief in the study area can be directly related to their location inside the duplex array and are thus linked to interaction between eastern MFF and NTT, and between western MFF and Taleghan fault, respectively. To account for inferred inherited topography from the previous dextral-transpression regime, a new concept of tectonic landscape characterization has been used. Accordingly, I define simple landscapes as those environments, which have developed during the influence of a sustained tectonic regime. In contrast, composite landscapes contain topographic elements inherited from previous tectonic conditions that are inconsistent with the regional present-day stress field and kinematic style. Using numerical fault-interaction modeling with different tectonic boundary conditions, I calculated synoptic snapshots of artificial topography to compare it with the real topographic metrics. However, in the Alborz mountains, E-W faults are favorably oriented to accommodate the entire range of NW- to NE-directed compression. These faults show the highest total displacement which might indicate sustained faulting under changing boundary conditions. In contrast to the fault system within and at the flanks of the Alborz mountains, Quaternary deformation in the adjacent Tehran plain is characterized by oblique motion and thrust and strike-slip fault systems. In this morphotectonic province fault-propagation folding along major faults, limited strike-slip motion, and en-échelon arrays of second-order upper plate thrusts are typical. While the Tehran plain is characterized by young deformation phenomena, the majority of faulting took place in the early stages of the Quaternary and during late Pliocene time. TCN-dating, which was performed for the first time on geomorphic surfaces in the Tehran plain, revealed that the oldest two phases of alluviation (units A and B) must be older than late Pleistocene. While urban development in Tehran increasingly covers and obliterates the active fault traces, the present-day kinematic style, the vestiges of formerly undeformed Quaternary landforms, and paleo earthquake indicators from the last millennia attest to the threat that these faults and their related structures pose for the megacity.
In den letzten drei Jahrzehnten wurden in einigen Seen und Feuchtgebieten in bewaldeten Einzugsgebieten Nordost-Brandenburgs sinkende Wasserstände beobachtet. In diesen Gebieten bestimmt die Grundwasserneubildung im Einzugsgebiet maßgeblich das Wasserdargebot der Seen und Feuchtgebiete, die deshalb hier als grundwasserabhängige Landschaftselemente bezeichnet werden. Somit weisen die sinkenden Wasserstände auf einen Rückgang der wegen des geringen Niederschlagsdargebotes ohnehin schon geringen Grundwasserneubildung hin. Die Höhe der Grundwasserneubildung ist neben den hydroklimatischen Randbedingungen auch von der Landnutzung abhängig. Veränderungen in der Waldvegetation und der hydroklimatischen Randbedingungen bewirken Änderungen der Grundwasserneubildung und beeinflussen somit auch den Wasserhaushalt der Seen und Feuchtgebiete. Aktuell wird die Waldvegetation durch Kiefernmonokulturen dominiert, mit im Vergleich zu anderen Baumarten höherer Evapotranspiration. Entwicklungen in der Forstwirtschaft streben die Verringerung von Kiefernmonokulturen an. Diese sollen langfristig auf geeigneten Standorten durch Laubmischwälder ersetzt werden. Dadurch lassen sich eine geringere Evapotranspiration und damit eine höhere Grundwasserneubildung erreichen. In der vorliegenden Arbeit werden am Beispiel des Redernswalder Sees und des Briesensees die Ursachen der beobachteten sinkenden Wasserstände analysiert. Ihre Wasserstände nahmen in den letzten 25 Jahren um mehr als 3 Meter ab. Weiterhin wird untersucht, wie die erwarteten Klimaänderungen und Veränderungen in der Waldbewirtschaftung die zukünftige Grundwasserneubildung und den Wasserhaushalt von Seen beeinflussen können. Die Entwicklung der Grundwasserneubildung im Untersuchungsgebiet wurde mit dem Wasserhaushaltsmodell WaSiM-ETH simuliert. Die Analyse der Wechselwirkungen der Seen mit dem regionalen quartären Grundwasserleitersystem erfolgte mit dem 3D-Grundwassermodell FEFLOW. Mögliche zukünftige Veränderungen der Grundwasserneubildung und der Seewasserstände durch Klimaänderungen und Waldumbau wurden mit Szenarienrechnungen bis zum Jahr 2100 analysiert. Die modellgestützte Analyse zeigte, dass die beobachteten abnehmenden Wasserstände zu etwa gleichen Anteilen durch Veränderungen der hydroklimatischen Randbedingungen sowie durch Veränderungen in der Waldvegetation und damit abnehmenden Grundwasserneubildungsraten zu erklären sind. Die zukünftigen Entwicklungen der Grundwasserneubildung und der Wasserstände sind geprägt von sich ändernden hydroklimatischen Randbedingungen und einem sukzessiven Wandel der Kiefernbestände zu Laubwäldern. Der Waldumbau hat positive Wirkungen auf die Grundwasserneubildung und damit auf die Wasserstände. Damit können die Einflüsse des eingesetzten REMO-A1B-Klimaszenarios zum Ende des Modellzeitraumes durch den Waldumbau nicht kompensiert werden, das Sinken des Wasserstandes wird jedoch wesentlich reduziert. Bei dem moderateren REMO-B1-Klimaszenario werden die Wasserstände des Jahres 2008 durch den Waldumbau bis zum Jahr 2100 überschritten.
Based on technological advances made within the past decades, ground-penetrating radar (GPR) has become a well-established, non-destructive subsurface imaging technique. Catalyzed by recent demands for high-resolution, near-surface imaging (e.g., the detection of unexploded ordnances and subsurface utilities, or hydrological investigations), the quality of today's GPR-based, near-surface images has significantly matured. At the same time, the analysis of oil and gas related reflection seismic data sets has experienced significant advances. Considering the sensitivity of attribute analysis with respect to data positioning in general, and multi-trace attributes in particular, trace positioning accuracy is of major importance for the success of attribute-based analysis flows. Therefore, to study the feasibility of GPR-based attribute analyses, I first developed and evaluated a real-time GPR surveying setup based on a modern tracking total station (TTS). The combination of current GPR systems capability of fusing global positioning system (GPS) and geophysical data in real-time, the ability of modern TTS systems to generate a GPS-like positional output and wireless data transmission using radio modems results in a flexible and robust surveying setup. To elaborate the feasibility of this setup, I studied the major limitations of such an approach: system cross-talk and data delays known as latencies. Experimental studies have shown that when a minimal distance of ~5 m between the GPR and the TTS system is considered, the signal-to-noise ratio of the acquired GPR data using radio communication equals the one without radio communication. To address the limitations imposed by system latencies, inherent to all real-time data fusion approaches, I developed a novel correction (calibration) strategy to assess the gross system latency and to correct for it. This resulted in the centimeter trace accuracy required by high-frequency and/or three-dimensional (3D) GPR surveys. Having introduced this flexible high-precision surveying setup, I successfully demonstrated the application of attribute-based processing to GPR specific problems, which may differ significantly from the geological ones typically addressed by the oil and gas industry using seismic data. In this thesis, I concentrated on archaeological and subsurface utility problems, as they represent typical near-surface geophysical targets. Enhancing 3D archaeological GPR data sets using a dip-steered filtering approach, followed by calculation of coherency and similarity, allowed me to conduct subsurface interpretations far beyond those obtained by classical time-slice analyses. I could show that the incorporation of additional data sets (magnetic and topographic) and attributes derived from these data sets can further improve the interpretation. In a case study, such an approach revealed the complementary nature of the individual data sets and, for example, allowed conclusions about the source location of magnetic anomalies by concurrently analyzing GPR time/depth slices to be made. In addition to archaeological targets, subsurface utility detection and characterization is a steadily growing field of application for GPR. I developed a novel attribute called depolarization. Incorporation of geometrical and physical feature characteristics into the depolarization attribute allowed me to display the observed polarization phenomena efficiently. Geometrical enhancement makes use of an improved symmetry extraction algorithm based on Laplacian high-boosting, followed by a phase-based symmetry calculation using a two-dimensional (2D) log-Gabor filterbank decomposition of the data volume. To extract the physical information from the dual-component data set, I employed a sliding-window principle component analysis. The combination of the geometrically derived feature angle and the physically derived polarization angle allowed me to enhance the polarization characteristics of subsurface features. Ground-truth information obtained by excavations confirmed this interpretation. In the future, inclusion of cross-polarized antennae configurations into the processing scheme may further improve the quality of the depolarization attribute. In addition to polarization phenomena, the time-dependent frequency evolution of GPR signals might hold further information on the subsurface architecture and/or material properties. High-resolution, sparsity promoting decomposition approaches have recently had a significant impact on the image and signal processing community. In this thesis, I introduced a modified tree-based matching pursuit approach. Based on different synthetic examples, I showed that the modified tree-based pursuit approach clearly outperforms other commonly used time-frequency decomposition approaches with respect to both time and frequency resolutions. Apart from the investigation of tuning effects in GPR data, I also demonstrated the potential of high-resolution sparse decompositions for advanced data processing. Frequency modulation of individual atoms themselves allows to efficiently correct frequency attenuation effects and improve resolution based on shifting the average frequency level. GPR-based attribute analysis is still in its infancy. Considering the growing widespread realization of 3D GPR studies there will certainly be an increasing demand towards improved subsurface interpretations in the future. Similar to the assessment of quantitative reservoir properties through the combination of 3D seismic attribute volumes with sparse well-log information, parameter estimation in a combined manner represents another step in emphasizing the potential of attribute-driven GPR data analyses.
In the high mountains of Asia, glaciers cover an area of approximately 115,000 km² and constitute one of the largest continental ice accumulations outside Greenland and Antarctica. Their sensitivity to climate change makes them valuable palaeoclimate archives, but also vulnerable to current and predicted Global Warming. This is a pressing problem as snow and glacial melt waters are important sources for agriculture and power supply of densely populated regions in south, east, and central Asia. Successful prediction of the glacial response to climate change in Asia and mitigation of the socioeconomic impacts requires profound knowledge of the climatic controls and the dynamics of Asian glaciers. However, due to their remoteness and difficult accessibility, ground-based studies are rare, as well as temporally and spatially limited. We therefore lack basic information on the vast majority of these glaciers. In this thesis, I employ different methods to assess the dynamics of Asian glaciers on multiple time scales. First, I tested a method for precise satellite-based measurement of glacier-surface velocities and conducted a comprehensive and regional survey of glacial flow and terminus dynamics of Asian glaciers between 2000 and 2008. This novel and unprecedented dataset provides unique insights into the contrasting topographic and climatic controls of glacial flow velocities across the Asian highlands. The data document disparate recent glacial behavior between the Karakoram and the Himalaya, which I attribute to the competing influence of the mid-latitude westerlies during winter and the Indian monsoon during summer. Second, I tested whether such climate-related longitudinal differences in glacial behavior also prevail on longer time scales, and potentially account for observed regionally asynchronous glacial advances. I used cosmogenic nuclide surface exposure dating of erratic boulders on moraines to obtain a glacial chronology for the upper Tons Valley, situated in the headwaters of the Ganges River. This area is located in the transition zone from monsoonal to westerly moisture supply and therefore ideal to examine the influence of these two atmospheric circulation regimes on glacial advances. The new glacial chronology documents multiple glacial oscillations during the last glacial termination and during the Holocene, suggesting largely synchronous glacial changes in the western Himalayan region that are related to gradual glacial-interglacial temperature oscillations with superimposed monsoonal precipitation changes of higher frequency. In a third step, I combine results from short-term satellite-based climate records and surface velocity-derived ice-flux estimates, with topographic analyses to deduce the erosional impact of glaciations on long-term landscape evolution in the Himalayan-Tibetan realm. The results provide evidence for the long-term effects of pronounced east-west differences in glaciation and glacial erosion, depending on climatic and topographic factors. Contrary to common belief the data suggest that monsoonal climate in the central Himalaya weakens glacial erosion at high elevations, helping to maintain a steep southern orographic barrier that protects the Tibetan Plateau from lateral destruction. The results of this thesis highlight how climatic and topographic gradients across the high mountains of Asia affect glacier dynamics on time scales ranging from 10^0 to 10^6 years. Glacial response times to climate changes are tightly linked to properties such as debris cover and surface slope, which are controlled by the topographic setting, and which need to be taken into account when reconstructing mountainous palaeoclimate from glacial histories or assessing the future evolution of Asian glaciers. Conversely, the regional topographic differences of glacial landscapes in Asia are partly controlled by climatic gradients and the long-term influence of glaciers on the topographic evolution of the orogenic system.
Flood design necessitates discharge estimates for large recurrence intervals. However, in a flood frequency analysis, the uncertainty of discharge estimates increases with higher recurrence intervals, particularly due to the small number of available flood data. Furthermore, traditional distribution functions increase unlimitedly without consideration of an upper bound discharge. Hence, additional information needs to be considered which is representative for high recurrence intervals. Envelope curves which bound the maximum observed discharges of a region are an adequate regionalisation method to provide additional spatial information for the upper tail of a distribution function. Probabilistic regional envelope curves (PRECs) are an extension of the traditional empirical envelope curve approach, in which a recurrence interval is estimated for a regional envelope curve (REC). The REC is constructed for a homogeneous pooling group of sites. The estimation of this recurrence interval is based on the effective sample years of data considering the intersite dependence among all sites of the pooling group. The core idea of this thesis was an improvement of discharge estimates for high recurrence intervals by integrating empirical and probabilistic regional envelope curves into the flood frequency analysis. Therefore, the method of probabilistic regional envelope curves was investigated in detail. Several pooling groups were derived by modifying candidate sets of catchment descriptors and settings of two different pooling methods. These were used to construct PRECs. A sensitivity analysis shows the variability of discharges and the recurrence intervals for a given site due to the different assumptions. The unit flood of record which governs the intercept of PREC was determined as the most influential aspect. By separating the catchments into nested and unnested pairs, the calculation algorithm for the effective sample years of data was refined. In this way, the estimation of the recurrence intervals was improved, and therefore the use of different parameter sets for nested and unnested pairs of catchments is recommended. In the second part of this thesis, PRECs were introduced into a distribution function. Whereas in the traditional approach only discharge values are used, PRECs provide a discharge and its corresponding recurrence interval. Hence, a novel approach was developed, which allows a combination of the PREC results with the traditional systematic flood series while taking the PREC recurrence interval into consideration. An adequate mixed bounded distribution function was presented, which in addition to the PREC results also uses an upper bound discharge derived by an empirical envelope curve. By doing so, two types of additional information which are representative for the upper tail of a distribution function were included in the flood frequency analysis. The integration of both types of additional information leads to an improved discharge estimation for recurrence intervals between 100 and 1000 years.
Large-scale volcanic deformation recently detected by radar interferometry (InSAR) provides new information and thus new scientific challenges for understanding volcano-tectonic activity and magmatic systems. The destabilization of such a system at depth noticeably affects the surrounding environment through magma injection, ground displacement and volcanic eruptions. To determine the spatiotemporal evolution of the Lazufre volcanic area located in the central Andes, we combined short-term ground displacement acquired by InSAR with long-term geological observations. Ground displacement was first detected using InSAR in 1997. By 2008, this displacement affected 1800 km2 of the surface, an area comparable in size to the deformation observed at caldera systems. The original displacement was followed in 2000 by a second, small-scale, neighbouring deformation located on the Lastarria volcano. We performed a detailed analysis of the volcanic structures at Lazufre and found relationships with the volcano deformations observed with InSAR. We infer that these observations are both likely to be the surface expression of a long-lived magmatic system evolving at depth. It is not yet clear whether Lazufre may trigger larger unrest or volcanic eruptions; however, the second deformation detected at Lastarria and the clear increase of the large-scale deformation rate make this an area of particular interest for closer continuous monitoring.
Das Parallel-Seismik-Verfahren dient vor allem der nachträglichen Längenmessung von Fundamentpfählen oder ähnlichen Elementen zur Gründung von Bauwerken. Eine solche Messung wird beispielsweise notwendig, wenn ein Gebäude verstärkt, erhöht oder anders als bisher genutzt werden soll, aber keine Unterlagen mehr über die Fundamente vorhanden sind. Das Messprinzip des schon seit einigen Jahrzehnten bekannten Verfahrens ist relativ einfach: Auf dem Pfahlkopf wird meist durch Hammerschlag eine Stoßwelle erzeugt, die durch den Pfahl nach unten läuft. Dabei wird Energie in den Boden abgegeben. Die abgestrahlten Wellen werden von Sensoren in einem parallel zum Pfahl hergestellten Bohrloch registriert. Aus den Laufzeiten lassen sich die materialspezifischen Wellengeschwindigkeiten im Pfahl und im Boden sowie die Pfahllänge ermitteln. Bisher wurde meist ein sehr einfaches Verfahren zur Datenauswertung verwendet, das die Länge der Pfähle systematisch überschätzt. In der vorliegenden Dissertation wurden die mathematisch-physikalischen Grundlagen beleuchtet und durch Computersimulation die Wellenausbreitung in Pfahl und Boden genau untersucht. Weitere Simulationen klärten den Einfluss verschiedener Mess- und Strukturparameter, beispielsweise den Einfluss von Bodenschichtung oder Fehlstellen im Pfahl. So konnte geklärt werden, in welchen Fällen mit dem Parallel-Seismik-Verfahren gute Ergebnisse erzielt werden können (z. B. bei Fundamenten in Sand oder Ton) und wo es an seine Grenzen stößt (z. B. bei Gründung im Fels). Auf Basis dieser Ergebnisse entstand ein neuer mathematischer Formalismus zur Auswertung der Laufzeiten. In Verbindung mit einem Verfahren zur Dateninversion, d. h. der automatischen Anpassung der Unbekannten in den Gleichungen an die Messergebnisse, lassen sich sehr viel genauere Werte für die Pfahllänge ermitteln als mit allen bisher publizierten Verfahren. Zudem kann man nun auch mit relativ großen Abständen zwischen Bohrloch und Pfahl (2 - 3 m) arbeiten. Die Methode wurde an simulierten Daten ausführlich getestet. Die Messmethode und das neue Auswerteverfahren wurden in einer Reihe praktischer Anwendungen getestet – und dies fast immer erfolgreich. Nur in einem Fall komplizierter Fundamentgeometrie bei gleichzeitig sehr hoher Anforderung an die Genauigkeit war schon nach Simulationen klar, dass hier ein Einsatz nicht sinnvoll ist. Dafür zeigte es sich, dass auch die Länge von Pfahlwänden und Spundwänden ermittelt werden kann. Die Parallel-Seismik-Methode funktioniert als einziges verfügbares Verfahren zur Fundamentlängenermittlung zugleich in den meisten Bodenarten sowie an metallischen und nichtmetallischen Fundamenten und kommt ohne Kalibrierung aus. Sie ist nun sehr viel breiter einsetzbar und liefert sehr viel genauere Ergebnisse. Die Simulationen zeigten noch Potential für Erweiterungen, zum Beispiel durch den Einsatz spezieller Sensoren, die zusätzliche Wellentypen empfangen und unterscheiden können.
Temporal gravimeter observations, used in geodesy and geophysics to study variation of the Earth’s gravity field, are influenced by local water storage changes (WSC) and – from this perspective – add noise to the gravimeter signal records. At the same time, the part of the gravity signal caused by WSC may provide substantial information for hydrologists. Water storages are the fundamental state variable of hydrological systems, but comprehensive data on total WSC are practically inaccessible and their quantification is associated with a high level of uncertainty at the field scale. This study investigates the relationship between temporal gravity measurements and WSC in order to reduce the hydrological interfering signal from temporal gravity measurements and to explore the value of temporal gravity measurements for hydrology for the superconducting gravimeter (SG) of the Geodetic Observatory Wettzell, Germany. A 4D forward model with a spatially nested discretization domain was developed to simulate and calculate the local hydrological effect on the temporal gravity observations. An intensive measurement system was installed at the Geodetic Observatory Wettzell and WSC were measured in all relevant storage components, namely groundwater, saprolite, soil, top soil and snow storage. The monitoring system comprised also a suction-controlled, weighable, monolith-filled lysimeter, allowing an all time first comparison of a lysimeter and a gravimeter. Lysimeter data were used to estimate WSC at the field scale in combination with complementary observations and a hydrological 1D model. Total local WSC were derived, uncertainties were assessed and the hydrological gravity response was calculated from the WSC. A simple conceptual hydrological model was calibrated and evaluated against records of a superconducting gravimeter, soil moisture and groundwater time series. The model was evaluated by a split sample test and validated against independently estimated WSC from the lysimeter-based approach. A simulation of the hydrological gravity effect showed that WSC of one meter height along the topography caused a gravity response of 52 µGal, whereas, generally in geodesy, on flat terrain, the same water mass variation causes a gravity change of only 42 µGal (Bouguer approximation). The radius of influence of local water storage variations can be limited to 1000 m and 50 % to 80 % of the local hydro¬logical gravity signal is generated within a radius of 50 m around the gravimeter. At the Geodetic Observatory Wettzell, WSC in the snow pack, top soil, unsaturated saprolite and fractured aquifer are all important terms of the local water budget. With the exception of snow, all storage components have gravity responses of the same order of magnitude and are therefore relevant for gravity observations. The comparison of the total hydrological gravity response to the gravity residuals obtained from the SG, showed similarities in both short-term and seasonal dynamics. However, the results demonstrated the limitations of estimating total local WSC using hydrological point measurements. The results of the lysimeter-based approach showed that gravity residuals are caused to a larger extent by local WSC than previously estimated. A comparison of the results with other methods used in the past to correct temporal gravity observations for the local hydrological influence showed that the lysimeter measurements improved the independent estimation of WSC significantly and thus provided a better way of estimating the local hydrological gravity effect. In the context of hydrological noise reduction, at sites where temporal gravity observations are used for geophysical studies beyond local hydrology, the installation of a lysimeter in combination with complementary hydrological measurements is recommended. From the hydrological view point, using gravimeter data as a calibration constraint improved the model results in comparison to hydrological point measurements. Thanks to their capacity to integrate over different storage components and a larger area, gravimeters provide generalized information on total WSC at the field scale. Due to their integrative nature, gravity data must be interpreted with great care in hydrological studies. However, gravimeters can serve as a novel measurement instrument for hydrology and the application of gravimeters especially designed to study open research questions in hydrology is recommended.
In 1915, Alfred Wegener published his hypotheses of plate tectonics that revolutionised the world for geologists. Since then, many scientists have studied the evolution of continents and especially the geologic structure of orogens: the most visible consequence of tectonic processes. Although the morphology and landscape evolution of mountain belts can be observed due to surface processes, the driving force and dynamics at lithosphere scale are less well understood despite the fact that rocks from deeper levels of orogenic belts are in places exposed at the surface. In this thesis, such formerly deeply-buried (ultra-) high-pressure rocks, in particular eclogite facies series, have been studied in order to reveal details about the formation and exhumation conditions and rates and thus provide insights into the geodynamics of the most spectacular orogenic belt in the world: the Himalaya. The specific area investigated was the Kaghan Valley in Pakistan (NW Himalaya). Following closure of the Tethyan Ocean by ca. 55-50 Ma, the northward subduction of the leading edge of India beneath the Eurasian Plate and subsequent collision initiated a long-lived process of intracrustal thrusting that continues today. The continental crust of India – granitic basement, Paleozoic and Mesozoic cover series and Permo-Triassic dykes, sills and lavas – has been buried partly to mantle depths. Today, these rocks crop out as eclogites, amphibolites and gneisses within the Higher Himalayan Crystalline between low-grade metamorphosed rocks (600-640°C/ ca. 5 kbar) of the Lesser Himalaya and Tethyan sediments. Beside tectonically driven exhumation mechanisms the channel flow model, that describes a denudation focused ductile extrusion of low viscosity material developed in the middle to lower crust beneath the Tibetan Plateau, has been postulated. To get insights into the lithospheric and crustal processes that have initiated and driven the exhumation of this (ultra-) high-pressure rocks, mineralogical, petrological and isotope-geochemical investigations have been performed. They provide insights into 1) the depths and temperatures to which these rocks were buried, 2) the pressures and temperatures the rocks have experienced during their exhumation, 3) the timing of these processes 4) and the velocity with which these rocks have been brought back to the surface. In detail, through microscopical studies, the identification of key minerals, microprobe analyses, standard geothermobarometry and modelling using an effective bulk rock composition it has been shown that published exhumation paths are incomplete. In particular, the eclogites of the northern Kaghan Valley were buried to depths of 140-100 km (36-30 kbar) at 790-640°C. Subsequently, cooling during decompression (exhumation) towards 40-35 km (17-10 kbar) and 630-580°C has been superseded by a phase of reheating to about 720-650°C at roughly the same depth before final exhumation has taken place. In the southern-most part of the study area, amphibolite facies assemblages with formation conditions similar to the deduced reheating phase indicate a juxtaposition of both areas after the eclogite facies stage and thus a stacking of Indian Plate units. Radiometric dating of zircon, titanite and rutile by U-Pb and amphibole and micas by Ar-Ar reveal peak pressure conditions at 47-48 Ma. With a maximum exhumation rate of 14 cm/a these rocks reached the crust-mantle boundary at 40-35 km within 1 Ma. Subsequent exhumation (46-41 Ma, 40-35 km) decelerated to ca. 1 mm/a at the base of the continental crust but rose again to about 2 mm/a in the period of 41-31 Ma, equivalent to 35-20 km. Apatite fission track (AFT) and (U-Th)/He ages from eclogites, amphibolites, micaschists and gneisses yielded moderate Oligocene to Miocene cooling rates of about 10°C/Ma in the high altitude northern parts of the Kaghan Valley using the mineral-pair method. AFT ages are of 24.5±3.8 to 15.6±2.1 Ma whereas apatite (U-Th)/He analyses yielded ages between 21.0±0.6 and 5.3±0.2 Ma. The southern-most part of the Valley is dominated by younger late Miocene to Pliocene apatite fission track ages of 7.6±2.1 and 4.0±0.5 Ma that support earlier tectonically and petrologically findings of a juxtaposition and stack of Indian Plate units. As this nappe is tectonically lowermost, a later distinct exhumation and uplift driven by thrusting along the Main Boundary Thrust is inferred. A multi-stage exhumation path is evident from petrological, isotope-geochemical and low temperature thermochronology investigations. Buoyancy driven exhumation caused an initial rapid exhumation: exhumation as fast as recent normal plate movements (ca. 10 cm/a). As the exhuming units reached the crust-mantle boundary the process slowed down due to changes in buoyancy. Most likely, this exhumation pause has initiated the reheating event that is petrologically evident (e.g. glaucophane rimmed by hornblende, ilmenite overgrowth of rutile). Late stage processes involved widespread thrusting and folding with accompanied regional greenschist facies metamorphism, whereby contemporaneous thrusting on the Batal Thrust (seen by some authors equivalent to the MCT) and back sliding of the Kohistan Arc along the inverse reactivated Main Mantle Thrust caused final exposure of these rocks. Similar circumstances have been seen at Tso Morari, Ladakh, India, 200 km further east where comparable rock assemblages occur. In conclusion, as exhumation was already done well before the initiation of the monsoonal system, climate dependent effects (erosion) appear negligible in comparison to far-field tectonic effects.
The Tibetan Plateau is the largest elevated landmass in the world and profoundly influences atmospheric circulation patterns such as the Asian monsoon system. Therefore this area has been increasingly in focus of palaeoenvironmental studies. This thesis evaluates the applicability of organic biomarkers for palaeolimnological purposes on the Tibetan Plateau with a focus on aquatic macrophyte-derived biomarkers. Submerged aquatic macrophytes have to be considered to significantly influence the sediment organic matter due to their high abundance in many Tibetan lakes. They can show highly 13C-enriched biomass because of their carbon metabolism and it is therefore crucial for the interpretation of δ13C values in sediment cores to understand to which extent aquatic macrophytes contribute to the isotopic signal of the sediments in Tibetan lakes and in which way variations can be explained in a palaeolimnological context. Additionally, the high abundance of macrophytes makes them interesting as potential recorders of lake water δD. Hydrogen isotope analysis of biomarkers is a rapidly evolving field to reconstruct past hydrological conditions and therefore of special relevance on the Tibetan Plateau due to the direct linkage between variations of monsoon intensity and changes in regional precipitation / evaporation balances. A set of surface sediment and aquatic macrophyte samples from the central and eastern Tibetan Plateau was analysed for composition as well as carbon and hydrogen isotopes of n-alkanes. It was shown how variable δ13C values of bulk organic matter and leaf lipids can be in submerged macrophytes even of a single species and how strongly these parameters are affected by them in corresponding sediments. The estimated contribution of the macrophytes by means of a binary isotopic model was calculated to be up to 60% (mean: 40%) to total organic carbon and up to 100% (mean: 66%) to mid-chain n-alkanes. Hydrogen isotopes of n-alkanes turned out to record δD of meteoric water of the summer precipitation. The apparent enrichment factor between water and n-alkanes was in range of previously reported ones (≈-130‰) at the most humid sites, but smaller (average: -86‰) at sites with a negative moisture budget. This indicates an influence of evaporation and evapotranspiration on δD of source water for aquatic and terrestrial plants. The offset between δD of mid- and long-chain n-alkanes was close to zero in most of the samples, suggesting that lake water as well as soil and leaf water are affected to a similar extent by those effects. To apply biomarkers in a palaeolimnological context, the aliphatic biomarker fraction of a sediment core from Lake Koucha (34.0° N; 97.2° E; eastern Tibetan Plateau) was analysed for concentrations, δ13C and δD values of compounds. Before ca. 8 cal ka BP, the lake was dominated by aquatic macrophyte-derived mid-chain n-alkanes, while after 6 cal ka BP high concentrations of a C20 highly branched isoprenoid compound indicate a predominance of phytoplankton. Those two principally different states of the lake were linked by a transition period with high abundances of microbial biomarkers. δ13C values were relatively constant for long-chain n-alkanes, while mid-chain n-alkanes showed variations between -23.5 to -12.6‰. Highest values were observed for the assumed period of maximum macrophyte growth during the late glacial and for the phytoplankton maximum during the middle and late Holocene. Therefore, the enriched values were interpreted to be caused by carbon limitation which in turn was induced by high macrophyte and primary productivity, respectively. Hydrogen isotope signatures of mid-chain n-alkanes have been shown to be able to track a previously deduced episode of reduced moisture availability between ca. 10 and 7 cal ka BP, indicated by a 20‰ shift towards higher δD values. Indications for cooler episodes at 6.0, 3.1 and 1.8 cal ka BP were gained from drops of biomarker concentrations, especially microbial-derived hopanoids, and from coincidental shifts towards lower δ13C values. Those episodes correspond well with cool events reported from other locations on the Tibetan Plateau as well as in the Northern Hemisphere. To conclude, the study of recent sediments and plants improved the understanding of factors affecting the composition and isotopic signatures of aliphatic biomarkers in sediments. Concentrations and isotopic signatures of the biomarkers in Lake Koucha could be interpreted in a palaeolimnological context and contribute to the knowledge about the history of the lake. Aquatic macrophyte-derived mid-chain n-alkanes were especially useful, due to their high abundance in many Tibetan Lakes and their ability to record major changes of lake productivity and palaeo-hydrological conditions. Therefore, they have the potential to contribute to a fuller understanding of past climate variability in this key region for atmospheric circulation systems.
The origin and evolution of granites has been widely studied because granitoid rocks constitute a major portion of the Earth ́s crust. The formation of granitic magma is, besides temperature mainly triggered by the water content of these rocks. The presence of water in magmas plays an important role due to the ability of aqueous fluids to change the chemical composition of the magma. The exsolution of aqueous fluids from melts is closely linked to a fractionation of elements between the two phases. Then, aqueous fluids migrate to shallower parts of the Earth ́s crust because of it ́s lower density compared to that of melts and adjacent rocks. This process separates fluids and melts, and furthermore, during the ascent, aqueous fluids can react with the adjacent rocks and alter their chemical signature. This is particularly impor- tant during the formation of magmatic-hydrothermal ore deposits or in the late stages of the evolution of magmatic complexes. For a deeper insight to these processes, it is essential to improve our knowledge on element behavior in such systems. In particular, trace elements are used for these studies and petrogenetic interpretations because, unlike major elements, they are not essential for the stability of the phases involved and often reflect magmatic processes with less ambiguity. However, for the majority of important trace elements, the dependence of the geochemical behavior on temperature, pressure, and in particular on the composition of the system are only incompletely or not at all experimentally studied. Former studies often fo- cus on the determination of fluid−melt partition coefficients (Df/m=cfluid/cmelt) of economically interesting elements, e.g., Mo, Sn, Cu, and there are some partitioning data available for ele- ments that are also commonly used for petrological interpretations. At present, no systematic experimental data on trace element behavior in fluid−melt systems as function of pressure, temperature, and chemical composition are available. Additionally, almost all existing data are based on the analysis of quenched phases. This results in substantial uncertainties, particularly for the quenched aqueous fluid because trace element concentrations may change upon cooling. The objective of this PhD thesis consisted in the study of fluid−melt partition coefficients between aqueous solutions and granitic melts for different trace elements (Rb, Sr, Ba, La, Y, and Yb) as a function of temperature, pressure, salinity of the fluid, composition of the melt, and experimental and analytical approach. The latter included the refinement of an existing method to measure trace element concentrations in fluids equilibrated with silicate melts di- rectly at elevated pressures and temperatures using a hydrothermal diamond-anvil cell and synchrotron radiation X-ray fluorescence microanalysis. The application of this in-situ method enables to avoid the main source of error in data from quench experiments, i.e., trace element concentration in the fluid. A comparison of the in-situ results to data of conventional quench experiments allows a critical evaluation of quench data from this study and literature data. In detail, starting materials consisted of a suite of trace element doped haplogranitic glasses with ASI varying between 0.8 and 1.4 and H2O or a chloridic solution with m NaCl/KCl=1 and different salinities (1.16 to 3.56 m (NaCl+KCl)). Experiments were performed at 750 to 950◦C and 0.2 or 0.5 GPa using conventional quench devices (externally and internally heated pressure vessels) with different quench rates, and at 750◦C and 0.2 to 1.4 GPa with in-situ analysis of the trace element concentration in the fluids. The fluid−melt partitioning data of all studied trace elements show 1. a preference for the melt (Df/m < 1) at all studied conditions, 2. one to two orders of magnitude higher Df/m using chloridic solutions compared to experiments with H2O, 3. a clear dependence on the melt composition for fluid−melt partitioning of Sr, Ba, La, Y, and Yb in experiments using chloridic solutions, 4. quench rate−related differences of fluid−melt partition coefficients of Rb and Sr, and 5. distinctly higher fluid−melt partitioning data obtained from in-situ experiments than from comparable quench runs, particularly in the case of H2O as starting solution. The data point to a preference of all studied trace elements for the melt even at fairly high salinities, which contrasts with other experimental studies, but is supported by data from studies of natural co-genetically trapped fluid and melt inclusions. The in-situ measurements of trace element concentrations in the fluid verify that aqueous fluids will change their composition upon cooling, which is in particular important for Cl free systems. The distinct differences of the in-situ results to quench data of this study as well as to data from the literature signify the im- portance of a careful fluid sampling and analysis. Therefore, the direct measurement of trace element contents in fluids equilibrated with silicate melts at elevated PT conditions represents an important development to obtain more reliable fluid−melt partition coefficients. For further improvement, both the aqueous fluid and the silicate melt need to be analyzed in-situ because partitioning data that are based on the direct measurement of the trace element content in the fluid and analysis of a quenched melt are still not completely free of quench effects. At present, all available data on element complexation in aqueous fluids in equilibrium with silicate melts at high PT are indirectly derived from partitioning data, which involves in these experiments assumptions on the species present in the fluid. However, the activities of chemical components in these partitioning experiments are not well constrained, which is required for the definition of exchange equilibria between melt and fluid species. For example, the melt-dependent variation of partition coefficient observed for Sr imply that this element can not only be complexed by Cl− as suggested previously. The data indicate a more complicated complexation of Sr in the aqueous fluid. To verify this hypothesis, the in-situ setup was also used to determine strontium complexation in fluids equilibrated with silicate melts at desired PT conditions by the application of X-ray absorption near edge structure (XANES) spectroscopy. First results show a strong effect of both fluid and melt composition on the resulting XANES spectra, which indicates different complexation environments for Sr.
Foreland-basin systems are excellent archives to decipher the feedbacks between surface and tectonic processes in orogens. The sedimentary architecture of a foreland-basin system reflects the balance between tectonic subsidence causing long-term accommodation space and sediment influx corresponding to efficiency of erosion and mass-redistribution processes. In order to explore the effects of climatic and tectonic forcing in such a system, I investigated the Oligo-Miocene foreland-basin sediments of the southern Alborz mountains, an intracontinental orogen in northern Iran, related to the Arabia-Eurasia continental collision. This work includes absolute dating methods such as 40Ar/39Ar and zircon (U-Th)/He thermochronology, magnetostratigraphy, sedimentological analysis, sandstone and conglomerate provenance study, carbon and oxygen isotope analysis, and clay mineralogy study. Results show a systematic correlation between coarsening-upward cycles and sediment accumulation rates in the basin on 105 to 106yr time scales. During thrust loading phases, the coarse-grained fraction supplied by the uplifting range is stored in the proximal part of the basin (sedimentary facies retrogradation), while fine-grained sediments are deposited in distal sectors. Variations in sediment provenance during these phases of enhanced tectonic activity give evidence for erosional unroofing phases and/or drainage-reorganization events. In addition, enhanced tectonic activity promoted the growth of topography and associated orographic barrier effects, as demonstrated by sedimentologic indicators and the analysis of stable C and O isotopes from calcareous paleosols and lacustrine/palustrine samples. Extensive progradation of coarse-grained deposits occurs during phases of decreased subsidence, when the coarse-grained fraction supplied by the uplifting range cannot be completely stored in the proximal part of the basin. In this environment, a reduction in basin subsidence is associated with laterally stacked fluvial channel deposits, and is related to intra-foreland uplift, as documented by growth strata, tectonic tilting, and sediment reworking. Increase in sediment accumulation rate associated with progradation of vertically-stacked coarse-grained fluvial channels also occurs. Paleosol O-isotope data shows that this increase is related to wetter climatic phases, suggesting that surface processes are more efficient and exhumation rates increase, giving rise to a positive feedback. Furthermore, isotopic and sedimentologic data show that starting from 10-9 Ma, climate became less arid with an increase in seasonality of precipitation. Because important changes were also recorded in the Mediterranean Sea and Asia at that time, the evidence for climatic variability observed in the Alborz mountains most likely reflects changes in Northern Hemisphere atmospheric circulation patterns. This study has additional implications for the evolution of the Alborz mountains and the Arabia-Eurasia continental collision zone. At the orogenic scale, the locus of deformation did not move steadily southward, but stepped forward and backward since Oligocene time. In particular, from ~ 17.5 to 6.2 Ma the orogen grew by a combination of frontal accretion and wedge-internal deformation on time scales of ca. 0.7 to 2 m.y. Moreover, the provenance data suggest that prior to 10-9 Ma the shortening direction changed from NW-SE to NNE-SSW, in agreement with structural data. On the scale of the entire collision zone, the evolution of the studied basins and adjacent mountain ranges suggests a new geodynamic model for the evolution of the Arabia-Eurasia continental collision zone. Numerous sedimentary basins in the Alborz mountains and in other locations of the Arabia-Eurasia collision zone record a change from a tensional (transtensional) to a compressional (transpressional) tectonic setting by ~ 36 Ma. I interpret this to reflect the onset of subduction of the stretched Arabian continental lithosphere beneath central Iran, leading to moderate plate coupling and lower- and upper-plate deformation (soft continental collision). The increase in deformation rates in the southern Alborz mountains from ~ 17.5 Ma suggests that significant upper-plate deformation must have started by the early Miocene most likely in response to an increase in degree of plate coupling. I suggest that this was related to the subduction of thicker Arabian continental lithosphere and the consequent onset of hard continental collision. This model reconciles the apparent lag time of 15-20 m.y between the late Eocene to early Oligocene age for the initial Arabia-Eurasia continental collision and the onset of widespread deformation across the collision zone to the north in early to late Miocene time.
This study presents noble gas compositions (He, Ne, Ar, Kr, and Xe) of lavas from several Hawaiian volcanoes. Lavas from the Hawaii Scientific Drilling Project (HSDP) core, surface samples from Mauna Kea, Mauna Loa, Kilauea, Hualalai, Kohala and Haleakala as well as lavas from a deep well on the summit of Kilauea were investigated. Noble gases, especially helium, are used as tracers for mantle reservoirs, based on the assumption that high 3He/4He ratios (>8 RA) represent material from the deep and supposedly less degassed mantle, whereas lower ratios (~ 8 RA) are thought to represent the upper mantle. Shield stage Mauna Kea, Kohala and Kilauea lavas yielded MORB-like to moderately high 3He/4He ratios, while 3He/4He ratios in post-shield stage Haleakala lavas are MORB-like. Few samples show 20Ne/22Ne and 21Ne/22Ne ratios different from the atmospheric values, however, Mauna Kea and Kilauea lavas with excess in mantle Ne agree well with the Loihi-Kilauea line in a neon three-isotope plot, whereas one Kohala sample plots on the MORB correlation line. The values in the 4He/40Ar* (40Ar* denotes radiogenic Ar) versus 4He diagram imply open system fractionation of He from Ar, with a deficiency in 4He. Calculated 4He/40Ar*, 3He/22Nes (22NeS denotes solar Ne) and 4He/21Ne ratios for the sample suite are lower than the respective production and primordial ratios, supporting the observation of a fractionation of He from the heavier noble gases, with a depletion of He with respect to Ne and Ar. The depletion of He is interpreted to be partly due to solubility controlled gas loss during magma ascent. However, the preferential He loss suggests that He is more incompatible than Ne and Ar during magmatic processes. In a binary mixing model, the isotopic He and Ne pattern are best explained by a mixture of a MORB-like end-member with a plume like or primordial end-member with a fractionation in 3He/22Ne, represented by a curve parameter r of 15 (r=(³He/²²Ne)MORB/(³He/²²Ne)PLUME or PRIMORDIAL). Whether the high 3He/4He ratios in Hawaiian lavas are indicative of a primitive component within the Hawaiian plume or are rather a product of the crystal-melt- partitioning behavior during partial melting remains to be resolved.
Recent years witnessed a vast advent of stalagmites as palaeoclimate archives. The multitude of geochemical and physical proxies and a promise of a precise and accurate age model greatly appeal to palaeoclimatologists. Although substantial progress was made in speleothem-based palaeoclimate research and despite high-resolution records from low-latitudinal regions, proving that palaeo-environmental changes can be archived on sub-annual to millennial time scales our comprehension of climate dynamics is still fragmentary. This is in particular true for the summer monsoon system on the Indian subcontinent. The Indian summer monsoon (ISM) is an integral part of the intertropical convergence zone (ITCZ). As this rainfall belt migrates northward during boreal summer, it brings monsoonal rainfall. ISM strength depends however on a variety of factors, including snow cover in Central Asia and oceanic conditions in the Indic and Pacific. Presently, many of the factors influencing the ISM are known, though their exact forcing mechanism and mutual relations remain ambiguous. Attempts to make an accurate prediction of rainfall intensity and frequency and drought recurrence, which is extremely important for South Asian countries, resemble a puzzle game; all interaction need to fall into the right place to obtain a complete picture. My thesis aims to create a faithful picture of climate change in India, covering the last 11,000 ka. NE India represents a key region for the Bay of Bengal (BoB) branch of the ISM, as it is here where the monsoon splits into a northwestward and a northeastward directed arm. The Meghalaya Plateau is the first barrier for northward moving air masses and receives excessive summer rainfall, while the winter season is very dry. The proximity of Meghalaya to the Tibetan Plateau on the one hand and the BoB on the other hand make the study area a key location for investigating the interaction between different forcings that governs the ISM. A basis for the interpretation of palaeoclimate records, and a first important outcome of my thesis is a conceptual model which explains the observed pattern of seasonal changes in stable isotopes (d18O and d2H) in rainfall. I show that although in tropical and subtropical regions the amount effect is commonly called to explain strongly depleted isotope values during enhanced rainfall, alone it cannot account for observed rainwater isotope variability in Meghalaya. Monitoring of rainwater isotopes shows no expected negative correlation between precipitation amount and d18O of rainfall. In turn I find evidence that the runoff from high elevations carries an inherited isotopic signature into the BoB, where during the ISM season the freshwater builds a strongly depleted plume on top of the marine water. The vapor originating from this plume is likely to memorize' and transmit further very negative d18O values. The lack of data does not allow for quantication of this plume effect' on isotopes in rainfall over Meghalaya but I suggest that it varies on seasonal to millennial timescales, depending on the runoff amount and source characteristics. The focal point of my thesis is the extraction of climatic signals archived in stalagmites from NE India. High uranium concentration in the stalagmites ensured excellent age control required for successful high-resolution climate reconstructions. Stable isotope (d18O and d13C) and grey-scale data allow unprecedented insights into millennial to seasonal dynamics of the summer and winter monsoon in NE India. ISM strength (i. e. rainfall amount) is recorded in changes in d18Ostalagmites. The d13C signal, reflecting drip rate changes, renders a powerful proxy for dry season conditions, and shows similarities to temperature-related changes on the Tibetan Plateau. A sub-annual grey-scale profile supports a concept of lower drip rate and slower stalagmite growth during dry conditions. During the Holocene, ISM followed a millennial-scale decrease of insolation, with decadal to centennial failures resulting from atmospheric changes. The period of maximum rainfall and enhanced seasonality corresponds to the Holocene Thermal Optimum observed in Europe. After a phase of rather stable conditions, 4.5 kyr ago, the strengthening ENSO system dominated the ISM. Strong El Nino events weakened the ISM, especially when in concert with positive Indian Ocean dipole events. The strongest droughts of the last 11 kyr are recorded during the past 2 kyr. Using the advantage of a well-dated stalagmite record at hand I tested the application of laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) to detect sub-annual to sub-decadal changes in element concentrations in stalagmites. The development of a large ablation cell allows for ablating sample slabs of up to 22 cm total length. Each analyzed element is a potential proxy for different climatic parameters. Combining my previous results with the LAICP- MS-generated data shows that element concentration depends not only on rainfall amount and associated leaching from the soil. Additional factors, like biological activity and hydrogeochemical conditions in the soil and vadose zone can eventually affect the element content in drip water and in stalagmites. I present a theoretical conceptual model for my study site to explain how climatic signals can be transmitted and archived in stalagmite carbonate. Further, I establish a first 1500 year long element record, reconstructing rainfall variability. Additionally, I hypothesize that volcanic eruptions, producing large amounts of sulfuric acid, can influence soil acidity and hence element mobilization.
Active continental margins are affected by complex feedbacks between tectonic, climate and surface processes, the intricate relations of which are still a matter of discussion. The Chilean convergent margin, forming the outstanding Andean subduction orogen, constitutes an ideal natural laboratory for the investigation of climate, tectonics and their interactions. In order to study both processes, I examined marine and lacustrine sediments from different depositional environments on- and offshore the south-central Chilean coast (38-40°S). I combined sedimentological, geochemical and isotopical analyses to identify climatic and tectonic signals within the sedimentary records. The investigation of marine trench sediments (ODP Site 1232, SONNE core 50SL) focused on frequency changes of turbiditic event layers since the late Pleistocene. In the active margin setting of south-central Chile, these layers were considered to reflect periodically occurring earthquakes and to constitute an archive of the regional paleoseismicity. The new results indicate glacial-interglacial changes in turbidite frequencies during the last 140 kyr, with short recurrence times (~200 years) during glacial and long recurrence times (~1000 years) during interglacial periods. Hence, the generation of turbidites appears to be strongly influenced by climate and sea level changes, which control on the amount of sediment delivered to the shelf edge and therewith the stability of the continental slope: more stable slope conditions during interglacial periods entail lower turbidite frequencies than in glacial periods. Since glacial turbidite recurrence times are congruent with earthquake recurrence times derived from the historical record and other paleoseismic archives of the region, I concluded that only during cold stages the sediment availability and slope instability enabled the complete series of large earthquakes to be recorded. The sediment transport to the shelf region is not only driven by climate conditions but also influenced by local forearc tectonics. Accelerating uplift rates along major tectonic structures involved drainage anomalies and river flow inversions, which seriously altered the sediment supply to the Pacific Ocean. Two examples for the tectonic hindrance of fluvial systems are the coastal lakes Lago Lanalhue and Lago Lleu Lleu. Both lakes developed within former river valleys, which once discharged towards the Pacific and were dammed by tectonically uplifted sills at ~8000 yr BP. Analyses of sediment cores from the lakes showed similar successions of marine/brackish deposits at the bottom, covered by lacustrine sediments on top. Dating of the transitions between these different units and the comparison with global sea level curves allowed me to calculate local Holocene uplift rates, which are distinctly higher for the upraised sills (Lanalhue: 8.83 ± 2.7 mm/yr, Lleu Lleu: 11.36 ± 1.77 mm/yr) than for the lake basins (Lanalhue: 0.42 ± 0.71 mm/yr, Lleu Lleu: 0.49 ± 0.44 mm/yr). I hence considered the sills to be the surface expression of a blind thrust associated with a prominent inverse fault that is controlling regional uplift and folding. After the final separation of Lago Lanalhue and Lago Lleu Lleu from the Pacific, a constant deposition of lacustrine sediments preserved continuous records of local environmental changes. Sequences from both lakes indicate a long-term climate trend with a significant shift from more arid conditions during the Mid-Holocene (8000 – 4200 cal yr BP) to more humid conditions during the Late Holocene (4200 cal yr BP – present). This trend is consistent with other regional paleoclimatic data and interpreted to reflect changes in the strength/position of the Southern Westerly Winds. Since ~5000 years, sediments of Lago Lleu Lleu are marked by numerous intercalated detrital layers that recur with a mean frequency of ~210 years. Deposition of these layers may be triggered by local tectonics (i.e. earthquakes), but may also originate from changes in the local climate (e.g. onset of modern ENSO conditions). During the last 2000 years, pronounced variations in the terrigenous sediment supply to both lakes suggest important hydrological changes on the centennial time-scale as well. A lower input of terrigenous matter points to less humid phases between 200 cal yr B.C. - 150 cal yr A.D., 900 - 1350 cal yr A.D. and 1850 cal yr A.D. to present (broadly corresponding to the Roman, Medieval, and Modern Warm Periods). More humid periods persisted from 150 - 900 cal yr A.D. and 1350 - 1850 cal yr A.D. (broadly corresponding to the Dark Ages and the Little Ice Age). In conclusion, the combined investigation of marine and lacustrine sediments is a feasible method for the reconstruction of climatic and tectonic processes on different time scales. My approach allows exploring both climate and tectonics in one and the same archive, and is largely transferable to other active margins worldwide.
Submarine landslides can generate local tsunamis posing a hazard to human lives and coastal facilities. Two major related problems are: (i) quantitative estimation of tsunami hazard and (ii) early detection of the most dangerous landslides. This thesis focuses on both those issues by providing numerical modeling of landslide-induced tsunamis and by suggesting and justifying a new method for fast detection of tsunamigenic landslides by means of tiltmeters. Due to the proximity to the Sunda subduction zone, Indonesian coasts are prone to earthquake, but also landslide tsunamis. The aim of the GITEWS-project (German-Indonesian Tsunami Early Warning System) is to provide fast and reliable tsunami warnings, but also to deepen the knowledge about tsunami hazards. New bathymetric data at the Sunda Arc provide the opportunity to evaluate the hazard potential of landslide tsunamis for the adjacent Indonesian islands. I present nine large mass movements in proximity to Sumatra, Java, Sumbawa and Sumba, whereof the largest event displaced 20 km³ of sediments. Using numerical modeling, I compute the generated tsunami of each event, its propagation and runup at the coast. Moreover, I investigate the age of the largest slope failures by relating them to the Great 1977 Sumba earthquake. Continental slopes off northwest Europe are well known for their history of huge underwater landslides. The current geological situation west of Spitsbergen is comparable to the continental margin off Norway after the last glaciation, when the large tsunamigenic Storegga slide took place. The influence of Arctic warming on the stability of the Svalbard glacial margin is discussed. Based on new geophysical data, I present four possible landslide scenarios and compute the generated tsunamis. Waves of 6 m height would be capable of reaching northwest Europe threatening coastal areas. I present a novel technique to detect large submarine landslides using an array of tiltmeters, as a possible tool in future tsunami early warning systems. The dislocation of a large amount of sediment during a landslide produces a permanent elastic response of the earth. I analyze this response with a mathematical model and calculate the theoretical tilt signal. Applications to the hypothetical Spitsbergen event and the historical Storegga slide show tilt signals exceeding 1000 nrad. The amplitude of landslide tsunamis is controlled by the product of slide volume and maximal velocity (slide tsunamigenic potential). I introduce an inversion routine that provides slide location and tsunamigenic potential, based on tiltmeter measurements. The accuracy of the inversion and of the estimated tsunami height near the coast depends on the noise level of tiltmeter measurements, the distance of tiltmeters from the slide, and the slide tsunamigenic potential. Finally, I estimate the applicability scope of this method by employing it to known landslide events worldwide.
It has always been enigmatic which processes control the accretion of the North American terranes towards the Pacific plate and the landward migration of the San Andreas plate boundary. One of the theories suggests that the Pacific plate first cools and captures the uprising mantle in the slab window, and then it causes the accretion of the continental crustal blocks. The alternative theory attributes the accretion to the capture of Farallon plate fragments (microplates) stalled in the ceased Farallon-North America subduction zone. Quantitative judgement between these two end-member concepts requires a 3D thermomechanical numerical modeling. However, the software tool required for such modeling is not available at present in the geodynamic modeling community. The major aim of the presented work is comprised basically of two interconnected tasks. The first task is the development and testing of the research Finite Element code with sufficiently advanced facilities to perform the three-dimensional geological time scale simulations of lithospheric deformation. The second task consists in the application of the developed tool to the Neogene deformations of the crust and the mantle along the San Andreas Fault System in Central and northern California. The geological time scale modeling of lithospheric deformation poses numerous conceptual and implementation challenges for the software tools. Among them is the necessity to handle the brittle-ductile transition within the single computational domain, adequately represent the rock rheology in a broad range of temperatures and stresses, and resolve the extreme deformations of the free surface and internal boundaries. In the framework of this thesis the new Finite Element code (SLIM3D) has been successfully developed and tested. This code includes a coupled thermo-mechanical treatment of deformation processes and allows for an elasto-visco-plastic rheology with diffusion, dislocation and Peierls creep mechanisms and Mohr-Coulomb plasticity. The code incorporates an Arbitrary Lagrangian Eulerian formulation with free surface and Winkler boundary conditions. The modeling technique developed is used to study the aspects influencing the Neogene lithospheric deformation in central and northern California. The model setup is focused on the interaction between three major tectonic elements in the region: the North America plate, the Pacific plate and the Gorda plate, which join together near the Mendocino Triple Junction. Among the modeled effects is the influence of asthenosphere upwelling in the opening slab window on the overlying North American plate. The models also incorporate the captured microplate remnants in the fossil Farallon subduction zone, simplified subducting Gorda slab, and prominent crustal heterogeneity such as the Salinian block. The results show that heating of the mantle roots beneath the older fault zones and the transpression related to fault stepping, altogether, render cooling in the slab window alone incapable to explain eastward migration of the plate boundary. From the viewpoint of the thermomechanical modeling, the results confirm the geological concept, which assumes that a series of microplate capture events has been the primary reason of the inland migration of the San Andreas plate boundary over the recent 20 Ma. The remnants of the Farallon slab, stalled in the fossil subduction zone, create much stronger heterogeneity in the mantle than the cooling of the uprising asthenosphere, providing the more efficient and direct way for transferring the North American terranes to Pacific plate. The models demonstrate that a high effective friction coefficient on major faults fails to predict the distinct zones of strain localization in the brittle crust. The magnitude of friction coefficient inferred from the modeling is about 0.075, which is far less than typical values 0.6 – 0.8 obtained by variety of borehole stress measurements and laboratory data. Therefore, the model results presented in this thesis provide additional independent constrain which supports the “weak-fault” hypothesis in the long-term ongoing debate over the strength of major faults in the SAFS.
This thesis presents investigations on sediments from two African lakes which have been recording changes in their surrounding environmental and climate conditions since more than 200,000 years. Focus of this work is the time of the last Glacial and the Holocene (the last ~100,000 years before present [in the following 100 kyr BP]). One important precondition for this kind of research is a good understanding of the present ecosystems in and around the lakes and of the sediment formation under modern climate conditions. Both studies therefore include investigations on the modern environment (including organisms, soils, rocks, lake water and sediments). A 90 m long sediment sequence was investigated from Lake Tswaing (north-eastern South Africa) using geochemical analyses. These investigations document alternating periods of high detrital input and low (especially autochthonous) organic matter content and periods of low detrital input, carbonatic or evaporitic sedimentation and high autochthonous organic matter content. These alternations are interpreted as changes between relatively humid and arid conditions, respectively. Before c. 75 kyr BP, they seem to follow changes in local insolation whereas afterwards they appear to be acyclic and are probably caused by changes in ocean circulation and/or in the mean position of the Inter-Tropical Convergence Zone (ITCZ). Today, these factors have main influence on precipitation in this area where rainfall occurs almost exclusively during austral summer. All modern organisms were analysed for their biomarker and bulk organic and compound-specific stable carbon isotope composition. The same investigations on sediments from the modern lake floor document the mixed input of the investigated individual organisms and reveal additional influences by methanotrophic bacteria. A comparison of modern sediment characteristics with those of sediments covering the time 14 to 2 kyr BP shows changes in the productivity of the lake and the surrounding vegetation which are best explained by changes in hydrology. More humid conditions are indicated for times older than 10 kyr BP and younger than 7.5 kyr BP, whereas arid conditions prevailed in between. These observations agree with the results from sediment composition and indications from other climate archives nearby. The second lake study deals with Lake Challa, a small, deep crater lake on the foot of Mount Kilimanjaro. In this lake form mm-scale laminated sediments which were analyses with micro-XRF scanning for changes in the element composition. By comparing these results with investigations on thin sections, results from ongoing sediment trap studies, meteorological data, and investigations on the surrounding rocks and soils, I develop a model for seasonal variability in the limnology and sedimentation of Lake Challa. The lake appears to be stratified during the warm rain seasons (October – December and March – May) during which detrital material is delivered to the lake and carbonates precipitate. On the lake floor forms a dark lamina with high contents of Fe and Ti and high Ca/Al and low Mn/Fe ratios. Diatoms bloom during the cool and windy season (June – September) when mixing down to c. 60 m depth provides easily bio-available nutrients. Contemporaneously, Fe and Mn-oxides are precipitating which cause high Mn/Fe ratios in the light diatom-rich laminae of the sediments. Trends in the Mn/Fe ratio of the sediments are interpreted to reflect changes in the intensity or duration of seasonal mixing in Lake Challa. This interpretation is supported by parallel changes in the organic matter and biogenic silica content observed in the 22 m long profile recovered from Lake Challa. This covers the time of the last 25 kyr BP. It documents a transition around 16 kyr BP from relatively well-mixed conditions with high detrital input during glacial times to stronger stratified conditions which are probably related to increasing lake levels in Challa and generally more humid conditions in East Africa. Intensified mixing is recorded for the time of the Younger Dryas and the period between 11.4 and 10.7 kyr BP. For these periods, reduced intensity of the SW monsoon and intensified NE monsoon are reported from archives of the Indian-Asian Monsoon region, arguing for the latter as a probable source for wind mixing in Lake Challa. This connection is probably also responsible for contemporaneous events in the Mn/Fe ratios of the Lake Challa sediments and in other records of northern hemisphere monsoon intensity during the Holocene and underlines the close interaction of global low latitude atmospheric circulation.
The Earth’s magnetic field (EMF) is generated by convections in the electrically conducting liquid iron-rich outer core, modified by the Earth’s rotation. A drastic manifestation of the dynamics of this fluid body is the occurrence of geomagnetic field reversals in the Earth’s history but also geomagnetic excursions, which are more frequent features of otherwise stable polarity chrons, but often poorly constrained in the geological record. To better understand the origin of the field, we need to know how the field has varied on different geological timescales. This includes not only information about changes in the ancient field’s direction but also about the absolute intensity (palaeointensity) and the age. This palaeointensity record is needed for compiling a full-vector description of the field. A palaeomagnetic and palaeointensity study on lava flows allows gaining insights about the evolution of the EMF through time and space. However, constraining the EMF evolution over different geological timescales remains a difficult objective due to the paucity of available palaeointensity data. One new alternative approach in palaeointensity studies is the recently proposed multispecimen parallel differential pTRM (MS) method, which has potentially several advantages over the commonly used Thellier method, because it is in theory independent of magnetic domain state, less prone to biasing effects, such as thermal alteration and significantly faster to perform in the laboratory. A study of highly active volcanic regions, such as the Trans-Mexican Volcanic Belt, seems promising when attempting a full-vector reconstruction or when looking for field excursions. One aim of this thesis was to gain new information about the occurrence and global validity of geomagnetic excursions from the Brunhes- or Matuyama Chron. For this purpose some 75 lava flows from within the Trans-Mexican Volcanic Belt were sampled for palaeomagnetic analyses. The scatter of virtual geomagnetic poles from lavas younger than 1.7 Ma was used for estimating palaeosecular variation and was found to be consistent with latitude dependent Model G and other high quality palaeomagnetic data from Mexico. The palaeomagnetic mean-vectors of 56 lavas were correlated to the Geomagnetic Polarity Timescale supplemented with information on geomagnetic excursions. On the grounds of their associated radioisotopic ages, four lavas were tentatively correlated with known excursions from marine records. Two lava flows dating of Brunhes Chron were associated with the Big Lost and Delts/Stage 17 excursions, respectively. From further two flows dating of Matuyama Chron, one flow was associated with either the Santa Rosa- or Kamikatsura excursions, while the other could have been emplaced during the Gilsa excursion. The most significant outcome was the finding that both Brunhes excursional flows display nearly fully reversed directions that deviate almost 180°C from the expected normal polarity direction. This observation could indicate that in particular the Big Lost and Delta/Stage17 excursions may represent other short periods during which the field completed a full reversal for a short time, such as was previously found for other older cryptochrons or tiny wiggles. Another focus of this thesis was set on estimating the feasibility of the new MS method for routine palaeointensity determination. This was accomplished by applying the MS method to samples from 11 historical lava flows from Mexico and Iceland from which the actual field intensity was either known from contemporary observatory data, or deduced from magnetic field models. Comparing observed with expected intensity values allowed to test the accuracy of the MS method. It a was found that the majority of palaeointensity estimates after the MS method yielded results that were very close or indistinguishable within the range of uncertainty from the expected values. However, a general trend towards an overestimate in the palaeointensity was also observed, which, on the grounds of corroborating rock magnetic analyses, was associated with multidomain material. This observation was taken as first evidence that the MS method is not entirely independent of magnetic domain state, as was originally claimed. However, a second experiment in which a modification of the most widely used Thellier method was applied to sister samples from 5 Icelandic flows revealed that, in comparison to the MS method, the latter produced more accurate and statistically better defined palaeointensities. Thus, from these first results, the MS method appeared as a viable alternative for future palaeointensity studies. Subsequently it was attempted to corroborate the directional record from Mexican lavas with palaeointensity data. It was possible to acquire palaeointensity estimates for 32 out of 51 investigated lava flows. These new results revealed that the new MS palaeointensities for Mexico are, with a high degree of statistical significance, around 30% higher than expected. The generally high palaeointensities seem to corroborate the results obtained from historical lava flows in this study and other previous studies on synthetic samples where domain state effects were found to cause overestimates in the palaeointensity of up to 30 per cent in the MS method. The primary process that leads to this overestimate is assigned to an asymmetry in the demagnetisation and remagnetisation process. Yet, this overestimate is expected to be no larger than what might be expected from Thellier experiments performed on samples with a given degree of multidomain behaviour.
Modern anthropogenic forcing of atmospheric chemistry poses the question of how the Earth System will respond as thousands of gigatons of greenhouse gas are rapidly added to the atmosphere. A similar, albeit nonanthropogenic, situation occurred during the early Paleogene, when catastrophic release of carbon to the atmosphere triggered abrupt increase in global temperatures. The best documented of these events is the Paleocene-Eocene Thermal Maximum (PETM, ~55 Ma) when the magnitude of carbon addition to the oceans and atmosphere was similar to those expected for the future. This event initiated global warming, changes in hydrological cycles, biotic extinction and migrations. A recently proposed hypothesis concerning changes in marine ecosystems suggests that this global warming strongly influenced the shallow-water biosphere, triggering extinctions and turnover in the Larger Foraminifera (LF) community and the demise of corals. The successions from the Adriatic Carbonate Platform (SW Slovenia) represent an ideal location to test the hypothesis of a possible causal link between the PETM and evolution of shallow-water organisms because they record continuous sedimentation from the Late Paleocene to the Early Eocene and are characterized by a rich biota, especially LF, fundamental for detailed biostratigraphic studies. In order to reconstruct paleoenvironmental conditions during deposition, I focused on sedimentological analysis and paleoecological study of benthic assemblages. During the Late Paleocene-earliest Eocene, sedimentation occurred on a shallow-water carbonate ramp system characterized by enhanced nutrient levels. LF represent the common constituent of the benthic assemblages that thrived in this setting throughout the Late Paleocene to the Early Eocene. With detailed biostratigraphic and chemostratigraphic analyses documenting the most complete record to date available for the PETM event in a shallow-water marine environment, I correlated chemostratigraphically for the first time the evolution of LF with the δ¹³C curves. This correlation demonstrated that no major turnover in the LF communities occurred synchronous with the PETM; thus the evolution of LF was mainly controlled by endogenous biotic forces. The study of Late Thanetian metric-sized microbialite-coral mounds which developed in the middle part of the ramp, documented the first Cenozoic occurrence of microbially-cemented mounds. The development of these mounds, with temporary dominance of microbial communities over corals, suggest environmentally-triggered “phase shifts” related to frequent fluctuations of nutrient/turbidity levels during recurrent wet phases which preceding the extreme greenhouse conditions of the PETM. The paleoecological study of the coral community in the microbialites-coral mounds, the study of corals from Early Eocene platform from SW France, and a critical, extensive literature research of Late Paleocene – Early Eocene coral occurrences from the Tethys, the Atlantic, the Caribbean realms suggested that these corals types, even if not forming extensive reefs, are common in the biofacies as small isolated colonies, piles of rubble or small patch-reefs. These corals might have developed ‘alternative’ life strategies to cope with harsh conditions (high/fluctuating nutrients/turbidity, extreme temperatures, perturbation of aragonite saturation state) during the greenhouse times of the early Paleogene, representing a good fossil analogue to modern corals thriving close to their thresholds for survival. These results demonstrate the complexity of the biological responses to extreme conditions, not only in terms of temperature but also nutrient supply, physical disturbance and their temporal variability and oscillating character.
This study investigated the warm-temperate to tropical shallow-water Miocene carbonates of the Perfugas basin (Anglona area), northern Sardinia, Italy (Central Mediterranean). The aim of this study was to identify and document the existence and significance of early diagenesis in this carbonate system, especially the diagenetic history, which reflects the diagenetic potential in terms of skeletal mineralogy. The motivation behind the present study was to investigate the role that early cementation has over facies stabilization linked to differences in biotic associations in shallow-water settings. Principal to this was to unravel the amount, kind and distribution of early cements in this type of carbonates, in order to complement previous studies, and hence acquire a more global perspective on non-tropical carbonate settings. The shallow-buried Sedini Limestone Unit was investigated for variations on early diagenetic features, as well as for the type of biotic association, and oxygen and carbon stable isotope stratigraphy. Results showed, that particularly at the Perfugas basin (< 15 km2), which evolves in time from a ramp into a steep-flanked platform, shallow-water facies are characterized by a “transitional” type of biotic association. The biotic assemblages change gradually over time from a heterozoan-rich into a photozoan-rich depositional system. This transition implies a change in the depositional environmental control factors such as temperature. It is considered that sedimentation took place under warm-temperate waters, which shifted to more warmer or tropical waters through time. Moreover, it was noticed that along with these changes, marine early syn-depositional cements (high-Mg calcite), with particular fabrics (e.g. fibrous), gradually contributed to the early lithification of rocks, favoring a steepening of the platform relief. The major controls for the shift of the depositional geometry was triggered by the change of the type of biotic associations (carbonate factory), related with the shift towards warmer conditions, and the development of early marine cementation. The identification of the amount and distribution of different cement phases, porosities and early diagenetic features, within facies and stratigraphy, showed that diagenesis is differential along depth, and within the depositional setting. High-Mg calcite cements (micrite, fibrous and syntaxial inclusion-rich) are early syn-depositional, facies-related (shallow-water), predominant at the platform phase, and marine in origin. Low-Mg calcite cements (bladed, syntaxial inclusionpoor and blocky) are early to late post-depositional, non-facies related (shallow- to deep-water) and shallow-burial marine in origin. However, a particular difference exists when looking at the amount and distribution of low-Mg calcite bladed cements. They become richer in shallow-water facies at the platform phase, suggesting that the enrichment of bladed cementation is linked to the appearance of metastable grains (e.g. aragonite). In both depositional profiles, the development of secondary porosity is the product of fabric-selective dissolution of grains (aragonite, high-Mg calcite) and/or cements (syntaxial inclusion-rich). However, stratigraphy and stable isotopes (oxygen and carbon), indicate that the molds found at shallower facies located beneath, and close to stratigraphic boundaries, have been produced by the infiltration of meteoric-derived water, which caused recrystallization without calcite cementation. Away from these stratigraphic locations, shallow- and deep-water facies show molds, and recrystallization, as well as low-Mg calcite cementation, interpreted as occurring during burial of these sediments by marine waters. The main cement source is suggested to be aragonite. Our results indicate that the Sedini Limestone Unit was transformed in three different diagenetic environments (marine, meteoric and shallow-burial marine); however, the degree of transformation in each diagenetic environment differs in the heterozoan-dominated ramp from the photozoan-dominated platform. It is suggested that the sediments from the ramp follow a diagenetic pathway similar to their heterozoan counterparts (i.e. lack of marine cementation, and loss of primary porosity by compaction), and the sediments from the platform follow a diagenetic pathway similar to their photozoan counterparts (i.e. marine cementation occluding primary porosity). However, in this carbonate setting, cements are Mg-calcite, no meteoric cementation was produced, and secondary porosity at shallow-water facies of the platform phase is mostly open and preserved. Despite the temporal and transitional change in biotic associations, ramp and platform facies (shallow- to deep-water facies) showed an oxygen isotope record overprinted by diagenesis. Oxygen primary marine signatures were not found. It is believed that burial diagenesis (recrystallization and low-Mg calcite cementation) was the main reason. This was unexpected at the ramp, since heterozoan-rich carbonates can hold isotope values close to primary marine signals due to their low-Mg calcite original composition. Ramp and platform facies (shallow- to deep-water facies) showed a carbon isotope record that was less affected by diagenesis. However, only at deep-water facies, did the carbon record show positive values comparable with carbon primary marine signals. The positive carbon values were noticed with major frequency at the platform deep-water facies. Moreover, these values usually showed a covariant trend with the oxygen isotope record; even that the latter did not hold positive values. The main conclusion of this work is that carbonates, deposited under warm-temperate to tropical conditions, have a unique facies, diagenesis and chemostratigraphic expression, which is different from their cool-water heterozoan or warm-water photozoan counterparts, reflecting the “transitional” nature of biotic association.
Modern acquisition of seismic data on receiver networks worldwide produces an increasing amount of continuous wavefield recordings. Hence, in addition to manual data inspection, seismogram interpretation requires new processing utilities for event detection, signal classification and data visualization. Various machine learning algorithms, which can be adapted to seismological problems, have been suggested in the field of pattern recognition. This can be done either by means of supervised learning using manually defined training data or by unsupervised clustering and visualization. The latter allows the recognition of wavefield patterns, such as short-term transients and long-term variations, with a minimum of domain knowledge. Besides classical earthquake seismology, investigations of temporal patterns in seismic data also concern novel approaches such as noise cross-correlation or ambient seismic vibration analysis in general, which have moved into focus within the last decade. In order to find records suitable for the respective approach or simply for quality control, unsupervised preprocessing becomes important and valuable for large data sets. Machine learning techniques require the parametrization of the data using feature vectors. Applied to seismic recordings, wavefield properties have to be computed from the raw seismograms. For an unsupervised approach, all potential wavefield features have to be considered to reduce subjectivity to a minimum. Furthermore, automatic dimensionality reduction, i.e. feature selection, is required in order to decrease computational cost, enhance interpretability and improve discriminative power. This study presents an unsupervised feature selection and learning approach for the discovery, imaging and interpretation of significant temporal patterns in seismic single-station or network recordings. In particular, techniques permitting an intuitive, quickly interpretable and concise overview of available records are suggested. For this purpose, the data is parametrized by real-valued feature vectors for short time windows using standard seismic analysis tools as feature generation methods, such as frequency-wavenumber, polarization, and spectral analysis. The choice of the time window length is dependent on the expected durations of patterns to be recognized or discriminated. We use Self-Organizing Maps (SOMs) for a data-driven feature selection, visualization and clustering procedure, which is particularly suitable for high-dimensional data sets. Using synthetics composed of Rayleigh and Love waves and three different types of real-world data sets, we show the robustness and reliability of our unsupervised learning approach with respect to the effect of algorithm parameters and data set properties. Furthermore, we approve the capability of the clustering and imaging techniques. For all data, we find improved discriminative power of our feature selection procedure compared to feature subsets manually selected from individual wavefield parametrization methods. In particular, enhanced performance is observed compared to the most favorable individual feature generation method, which is found to be the frequency spectrum. The method is applied to regional earthquake records at the European Broadband Network with the aim to define suitable features for earthquake detection and seismic phase classification. For the latter, we find that a combination of spectral and polarization features favor S wave detection at a single receiver. However, SOM-based visualization of phase discrimination shows that clustering applied to the records of two stations only allows onset or P wave detection, respectively. In order to improve the discrimination of S waves on receiver networks, we recommend to consider additionally the temporal context of feature vectors. The application to continuous recordings of seismicity close to an active volcano (Mount Merapi, Java, Indonesia) shows that two typical volcano-seismic events (VTB and Guguran) can be detected and distinguished by clustering. In contrast, so-called MP events cannot be discriminated. Comparable results are obtained for selected features and recognition rates regarding a previously implemented supervised classification system. Finally, we test the reliability of wavefield clustering to improve common ambient vibration analysis methods such as estimation of dispersion curves and horizontal to vertical spectral ratios. It is found, that in general, the identified short- and long-term patterns have no significant impact on those estimates. However, for individual sites, effects of local sources can be identified. Leaving out the corresponding clusters, yields reduced uncertainties or allows for improving estimation of dispersion curves.
The modern foreland basin straddling the eastern margin of the Andean orogen is the prime example of a retro-arc foreland basin system adjacent to a subduction orogen. While widely studied in the central and southern Andes, the spatial and temporal evolution of the Cenozoic foreland basin system in the northern Andes has received considerably less attention. This is in part due to the complex geodynamic boundary conditions, such as the oblique subduction and accretion of the Caribbean plates to the already complex interaction between the Nazca and the South American plates. In the Colombian Andes, for example, a foreland basin system has been forming since ~80 Ma over an area previously affected by rift tectonics during the Mesozoic. This setting of Cenozoic contractile deformation superposed on continental crust pre-strained by extensional processes thus represents a natural, yet poorly studied experimental set-up, where the role of tectonic inheritance on the development of foreland basin systems can be evaluated. However, a detailed documentation of the early foreland basin evolution in this part of the Andes has thus far only been accomplished in the more internal sectors of the orogen. In this study, I integrate new structural, sedimentological and biostratigraphic data with low-temperature thermochronology from the eastern sector of the Colombian Andes, in order to provide the first comprehensive account of mountain building and related foreland basin sedimentation in this part of the orogen, and to assess as to what extent pre-existent basement anisotropies have conditioned the locus of foreland deformation in space and time. In the Medina Basin, along the eastern flank of the Eastern Cordillera, I integrated detailed structural mapping and new sedimentological data with a new chronostratigraphic framework based on detailed palynology that links an eastward-thinning early Oligocene to early Miocene syntectonic wedge containing rapid facies changes with an episode of fast tectonic subsidence starting at ~30 Ma. This record represents the first evidence of topographic loading generated by slip along the principal basement-bounding thrusts in the Eastern Cordillera to the west of the basin and thus constrains the onset of mountain building in this area. A comprehensive assessment of exhumation patterns based on zircon fission-track (ZFT), apatite fission-track (AFT) analysis and thermal modelling reveals the location of these thrust loads to have been located along the contractionally reactivated Soapaga Fault in the axial sector of the Eastern Cordillera. Farther to the east, AFT and ZFT data also document the onset of thrust-induced exhumation associated with contractional reactivation of the main range-bounding Servita Fault at ~20 Ma. Associated with this episode of orogenic growth, peak burial temperature estimates based on vitrinite reflectance data in the Cenozoic sedimentary record of the adjacent Medina Basin documents earlier incorporation of the western sector of the basin into the advancing fold and thrust belt. I combined these new thermochronological data with published AFT analyses and known chronologic indicators of brittle deformation in order to evaluate the patterns of orogenic-front migration in the Andes of central Colombia. This spatiotemporal analysis of deformation reveals an episodic pattern of eastward migration of the orogenic front at an average rate of 2.5-2.7 mm/yr during the Late Cretaceous-Cenozoic. I identified three major stages of orogen propagation. First, following initiation of mountain building in the Central Cordillera during the Late Cretaceous, the orogenic front propagate eastward at slow rates (0.5-3.1 mm/yr) until early Eocene times. Such slow orogenic advance would have resulted from limited accretionary flux related to slow and oblique (SW-NE-oriented) convergence of the Farallon and South American plates during that time. A second stage of rapid orogenic advance (4.0-18.0 mm/yr) during the middle-late Eocene, and locally of at least 100 mm/yr in the middle Eocene, resulted from initial tectonic inversion of the Eastern Cordillera. I correlate this episode of rapid orogen-front migration with an increase in the accretionary flux triggered by acceleration in convergence and a rotation of the convergence vector to a more orogen-perpendicular direction. Finally, stagnation of the Miocene deformation front along former rift-bounding reactivated faults in the eastern flank of the Eastern Cordillera led to a decrease in the rates of orogenic advance. Post-late Miocene-Pliocene thrusting along the actively deforming front of the Eastern Cordillera at this latitude suggests averaged Miocene-Holocene orogen propagation rates of 1.2-2.1 mm/yr. In addition, ZFT data suggest that exhumation along the eastern flank of the orogen occurred at moderate rates of ~0.3 mm/yr during the Miocene, prior to an acceleration of exhumation since the Pliocene, as suggested by recently published AFT data. In order to evaluate the relations between thrust loading and sedimentary facies evolution in the foreland, I analyzed gravel progradation in the foreland basin system. In particular, I compared one-dimensional Eocene to Pliocene sediment accumulation rates in the Medina basin with a three-dimensional sedimentary budget based on the interpretation of ~1800 km of industry-style seismic reflection profiles and borehole data tied to the new chronostratigraphic framework. The sedimentological data from the Medina Basin reveal rapid accumulation of fluvial and lacustrine sediments at rates of up to ~ 0.5 mm/yr during the Miocene. Provenance data based on gravel petrography and paleocurrents reveal that these Miocene fluvial systems were sourced by Upper Cretaceous and Paleocene sedimentary units exposed to the west, in the Eastern Cordillera. Peak sediment-accumulation rates in the upper Carbonera Formation and the Guayabo Group occur during episodes of gravel progradation in the proximal foredeep in the Early and Late Miocene. I interpreted this positive correlation between sediment accumulation and gravel deposition as the direct consequence of thrust activity in the Servita-Lengupá Fault. This contrasts with current models relating gravel progradation to episodes of tectonic quiescence in more distal portions of foreland basin systems and calls for a re-evaluation of tectonic history interpretations inferred from sedimentary units in other mountain belts. In summary, my results document a late Eocene-early Miocene eastward advance of the topographic loads associated with the leading edge of deformation in the northern Andes of Colombia. Crustal thickening of the Eastern Cordillera associated with initiation of thrusting along the Servitá Fault illustrates that this sector of the Andean orogen acquired ~90% of its present width already by the early Miocene (~20 Ma). My data thus demonstrate that inherited crustal anisotropies, such as the former rift-bounding faults of the Eastern Cordillera, favour a non-systematic progression of foreland basin deformation through time by preferentially concentrating accommodation of slip and thrust-loading. These new chronology of exhumation and deformation associated with specific structures in the Colombian Andes also constitutes an important advance towards the understanding of models for hydrocarbon maturation, migration and trap formation along the prolific petroleum province of the Llanos Basin in the modern foredeep area.
Im Rahmen der Dissertation wurden an Wässern und freien Gasen aus Thermalquellen sowie an weniger als 5 Millionen Jahre alten basischen Vulkaniten des zentralandinen Puna-Hochplateaus (NE-Argentinien) umfangreiche element- und isotopengeochemische Untersuchungen durchgeführt und die Edelgasgehalte und -isotopensignaturen in diesen Medien bestimmt. Damit soll ein Beitrag zum besseren Verständnis der jüngeren Subduktionsgeschichte im Bereich der südlichen Zentralanden geleistet, die Wechselwirkungen zwischen ozeanischer Unter- und kontinentaler Oberplatte sichtbar gemacht und die Edelgassystematik verbessert werden. Wie die Ergebnisse der Untersuchungen an Gasen aus den Thermalquellen der Puna-Region zeigen, ist der Anteil an Mantel-Helium in den Thermalquellen dieser Region mit bis zu 67 % wesentlich höher als in der westlich gelegenen vulkanisch aktiven Westkordillere und den anderen angrenzenden Gebieten. In einigen Quellen konnten sogar Anteile an Mantel-Neon nachgewiesen werden, was aufgrund von Überlagerungen mit Neon atmosphärischen und krustalen Ursprungs weltweit bisher nur vereinzelt gelungen ist. Für kontinentale Bereiche mit großer Krustendicke ist ein solch starker Mantelgasfluss äußerst ungewöhnlich und bedeutet, dass Mantelschmelzen bis in die Kruste aufgedrungen sind und tief reichende Wegsamkeiten existieren, so dass die Mantelgase aufsteigen können, ohne stark krustal beeinflusst zu werden. Dass im Bereich der Puna rezent Mantelmaterial in die Kruste aufsteigt, zu diesem Ergebnis kommen auch aktuelle seismologische Untersuchungen. Zudem wurden junge, vorwiegend monogenetische Basalte bis basaltische Andesite geochemisch auf ihre Haupt-, Neben- und Spurenbestandteile sowie ihre Gehalte an Seltenenerdenelementen hin untersucht. Auch wurden die Isotopenverhältnisse von Sr, Nd und Pb in den Gesteinen bestimmt und petrographisch-mineralogische Analysen der darin enthaltenen Olivine und Pyroxene durchgeführt. Wie die Resultate belegen, haben die Magmen bei ihrem Aufstieg durch die Erdkruste insbesondere Material aus der Oberkruste assimiliert und sind zudem durch Fluide aus der abtauchenden Platte beeinflusst worden. Damit konnte gezeigt werden, dass einfache geochemische Methoden allein nicht ausreichen, um die Mantelquelle der Magmen ermitteln oder Aussagen über die Asthenosphärendynamik in der Region machen zu können. Im Gegensatz dazu zeigen die Messungen der Edelgasisotopenverhältnisse in den Fluideinschlüssen der Olivine und Pyroxene, dass deren Edelgaszusammensetzung nicht durch Krustenkontamination beeinflusst wurde, weil die Magmen erst nach der Olivin- bzw. Pyroxen-Kristallisation Schmelzen aus der Oberkruste assimiliert haben. Darüber hinaus konnten durch die Edelgasisotopenmessungen die bisher höchsten magmatischen He- und Ne-Isotopenverhältnisse von ganz Südamerika nachgewiesen werden. Aus der unterschiedlichen Höhe der Messwerte ist zu schließen, dass die im Osten der Puna vorkommenden älteren Laven aus einem nichtkonvektiven (lithosphärischen) Mantel stammen, während die am vulkanischen Bogen und Westrand der Puna gelegenen jüngeren Laven, ihren Ursprung in einer konvektiven (asthenosphärischen) Mantelquelle haben. Zudem konnte gezeigt werden, dass der Mantelgasfluss in der Region in den letzten 5 Millionen Jahren stark zunahm und sich die Eruption von mantelstämmigen basischen Laven in dieser Zeit kontinuierlich in westliche Richtung zum aktiven Vulkanbogen hin verlagerte. Im daraus abgeleiteten Modell beruht dieser Prozess (1) auf einer an die Kontinentalverschiebung gekoppelten W-Drift des Kontinents und (2) auf einem mit der Versteilung der Unterplatte verbundenen Vordringen des subkontinentalen asthenosphärischen Mantels nach W, nach dem Ende der Subduktion des unterseeischen aseismischen Juan Fernández-Rückens in der Region. Zudem gibt es starke Argumente dafür, dass die asthenosphärischen Magmen aus einer fluidreichen Zone in 500 – 600 km Tiefe parallel zur subduzierten Platte aufsteigen und nicht, wie bisher angenommen, durch Schmelzbildung in Bereichen unter 200 km Tiefe, allein durch Entwässerung der abtauchenden Platte erzeugt werden. Zu diesem Resultat führt vor allem die Kombination der He-Isotopenverhältnisse mit Ergebnissen seismologischer Untersuchungen.
Landscapes evolve in a complex interplay between climate and tectonics. Thus, the geomorphic characteristics of a landscape can only be understood if both, climatic and tectonic signals of past and ongoing processes can be identified. In order to evaluate the impact of both forcing factors it is crucial to quantify the evolution of geomorphic markers in natural environments. The Cenozoic Andes are an ideal setting to evaluate tectonic and climatic aspects of landscape evolution at different time and length scales in different natural compartments. The Andean Cordillera constitutes the type subduction orogen and is associated with the subduction of the oceanic Nazca Plate beneath the South American continent since at least 200 million years. In Chile and the adjacent regions this convergent margin is characterized by active tectonics, volcanism, and mountain building. Importantly, along the coast of Chile megathrust earthquakes occur frequently and influence landscape evolution. In fact, the largest earthquake ever recorded occurred in south-central Chile in 1960 and comprised a rupture zone of ~ 1000 km length. However, on longer time scales beyond historic documentation of seismicity it is not well known, how such seismotectonic segments have behaved and how they influence the geomorphic evolution of the coastal realms. With several semi-independent morphotectonic segments, recurrent megathrust earthquakes, and a plethora of geomorphic features indicating sustained tectonism, the margin of Chile is thus a key area to study relationships between surface processes and tectonics. In this study, I combined geomorphology, geochronology, sedimentology, and morphometry to quantify the Pliocene-Pleistocene landscape evolution of the tectonically active south-central Chile forearc. Thereby, I provide (1) new results about the influence of seismotectonic forearc segmentation on the geomorphic evolution and (2) new insights in the interaction between climate and tectonics with respect to the morphology of the Chilean forearc region. In particular, I show that the forearc is characterized by three long-term segments that are not correlated with short-lived earthquake-rupture zones that may. These segments are the Nahuelbuta, Toltén, and Bueno segments, each recording a distinct geomorphic and tectonic evolution. The Nahuelbuta and Bueno segments are undergoing active tectonic uplift. The long-term behavior of these two segments is manifested in form of two doubly plunging, growing antiforms that constitute an integral part of the Coastal Cordillera and record the uplift of marine and river terraces. In addition, these uplifting areas have caused major changes in flow directions or rivers. In contrast, the Toltén segment, situated between the two other segments, appears to be quasi-stable. In order to further quantify uplift and incision in the actively deforming Nahuelbuta segment, I dated an erosion surface and fluvial terraces in the Coastal Cordillera with cosmogenic 10Be and 26Al and optically stimulated luminescence, respectively. According to my results, late Pleistocene uplift rates corresponding to 0.88 mm a-1 are faster than surface-uplift rates averaging over the last 5 Ma, which are in the range of 0.21 mm a-1. This discrepancy suggests that surface uplift is highly variable in time and space and might preferably concentrate along reverse faults as indicated by a late Pleistocene flow reversal. In addition, the results of exposure dating with cosmogenic 10Be and 26Al indicate that the morphotectonic segmentation of this region of the forearc has been established in Pliocene time, coeval with the initiation of uplift of the Coastal Cordillera about 5 Ma ago, inferred to be related to a shift in subduction mode from erosion to accretion. Finally, I dated volcanic clasts obtained from alluvial surfaces in the Central Depression, a low-relief sector separating the Coastal from the Main Cordillera, with stable cosmogenic 3He and 21Ne, in order to reveal the controls of sediment accumulation in the forearc. My results document that these gently sloping surfaces have been deposited 150 to 300 ka ago. This deposition may be related to changes in the erosional regime during glacial episodes. Taken together, the data indicates that the overall geomorphic expression of the forearc is of post-Miocene age and may be intimately related to a climatic overprint of the tectonic system. This climatic forcing is also reflected in the topography and local relief of the Central and Southern Andes that vary considerably along the margin, determined by the dominant surface process that in turn is eventually controlled by climate. However, relief also partly reflects surface processes that have taken place under past climatic conditions. This emphasizes that due care has to be exercised when interpreting landscapes as mirrors of modern climates.
The arctic region is undergoing the most rapid environmental change experienced on Earth, and the rate of change is expected to increase over the coming decades. Arctic coasts are particularly vulnerable because they lie at the interface between terrestrial systems dominated by permafrost and marine systems dominated by sea ice. An increased rise in sea level and degradation of sea-ice as predicted by the Intergovernmental Panel on Climate Change in its most recent report and as observed recently in the Arctic will likely result in greater rates of coastal retreat. An increase in coastal erosion would result in dramatic increases in the volume of sediment, organic carbon and contaminants to the Arctic Ocean. These in turn have the potential to create dramatic changes in the geochemistry and biodiversity of the nearshore zone and affect the Arctic Ocean carbon cycle. To calculate estimates of organic carbon input from coastal erosion to the Arctic Ocean, current methods rely on the length of the coastline in the form of non self-similar line datasets. This thesis however emphasizes that using shorelines drawn at different scales can induce changes in the amount of sediment released by 30% in some cases. It proposes a substitute method of computations of erosion based on areas instead of lengths (i.e. buffers instead of shoreline lengths) which can be easily implemented at the circum-Arctic scale. Using this method, variations in quantities of eroded sediment are, on average, 70% less affected by scale changes and are therefore a more reliable method of calculation. Current estimates of coastal erosion rates in the Arctic are scarce and long-term datasets are a handful, which complicates assessment and prognosis of coastal processes, in particular the occurrence of coastal hazards. This thesis aims at filling the gap by providing the first long-term dataset (1951-2006) of coastal erosion on the Bykovsky Peninsula, North-East Siberia. This study shows that the coastline, which is made of ice-rich permafrost, retreated at a mean annual rate of 0.59 m/yr between 1951and 2006. Rates were highly variable: 97.0 % of the rates observed were less than 2 m/yr and 81.6% were less than 1m/yr. However, no significant trend in erosion could be recorded despite the study of five temporal subperiods within 1951-2006. The juxtaposition of wind records could not help to explain erosion records either and this thesis emphasizes the local controls on erosion, in particular the cryostratigraphy, the proximity of the Peninsula to the Lena River Delta freshwater plume and the local topographical constraints on swell development. On ice-rich coastal stretches of the Artic, the interaction of coastal dynamics and permafrost leads to the occurrence of spectacular “C-shaped” depressions termed retrogressive thaw slumps which can reach lengths of up to 650 m. On Herschel Island and at King Point (Yukon Coastal Plain, northern Canada), topographical, sedimentological and biogeochemical surveys were conducted to investigate the present and past activity of these landforms. In particular, undisturbed tundra areas were compared with zones of former slump activity, now stabilized and re-vegetated. This thesis shows that stabilized areas are drier and less prone to plant growth than undisturbed areas and feature fundamentally different geotechnical properties. Radiocarbon dating and topographical surveys indicated until up to 300 BP a likely period of dramatic slump activity on Herschel Island, similar to the one currently observed, which led to the creation of these surfaces. This thesis hypothesizes the occurrence of a ~250 years cycle of slump activity on the Herschel Island shoreline based on the surveyed topography and cryostratigraphy and anticipates higher frequency of slump activity in the future. The variety of processes described in this thesis highlights the changing nature of the intensity and frequency of physical processes acting upon the arctic coast. It also challenges current perceptions of the threats to existing industry and community infrastructure in the Arctic. The increasing presence of humans on Artic coasts coupled with the expected development of shipping will drive an increase in economical and industrial activity on these coasts which remains to be addressed scientifically.
The South Chilean subduction zone between 41° and 43.5°S : seismicity, structure and state of stress
(2008)
While the northern and central part of the South American subduction zone has been intensively studied, the southern part has attracted less attention, which may be due to its difficult accessibility and lower seismic activity. However, the southern part exhibits strong seismic and tsunamogenic potential with the prominent example of the Mw=9.5 May 22, 1960 Valdivia earthquake. In this study data from an amphibious seismic array (Project TIPTEQ) is presented. The network reached from the trench to the active magmatic arc incorporating the Island of Chiloé and the north-south trending Liquiñe-Ofqui fault zone (LOFZ). 364 local events were observed in an 11-month period from November 2004 until October 2005. The observed seismicity allows to constrain for the first time the current state of stress of the subducting plate and magmatic arc, as well as the local seismic velocity structure. The downgoing Benioff zone is readily identifiable as an eastward dipping plane with an inclination of ~30°. Main seismic activity occurred predominantly in a belt parallel to the coast of Chiloé Island in a depth range of 12-30 km, which is presumably related to the plate interface. The down-dip termination of abundant intermediate depth seismicity at approximately 70 km depth seems to be related to the young age (and high temperature) of the oceanic plate. A high-quality subset of events was inverted for a 2-D velocity model. The vp model resolves the sedimentary basins and the downgoing slab. Increased velocities below the longitudinal valley and the eastern part of Chiloé Island suggest the existence of a mantle bulge. Apart from the events in the Benioff Zone, shallow crustal events were observed mainly in different clusters along the magmatic arc. These crustal clusters of seismicity are related to the LOFZ, as well as to the volcanoes Chaitén, Michinmahuida and Corcovado. Seismic activity up to a magnitude of 3.8 Mw reveals the recent activity of the fault zone. Focal mechanisms for the events along the LOFZ were calculated using a moment tensor inversion of amplitude spectra for body waves which mostly yield strike-slip mechanisms indicating a SW-NE striking of sigma_1 for the LOFZ. Focal mechanism stress inversion indicates a strike-slip regime along the arc and a thrust regime in the Benioff zone. The observed deformation - which is also revealed by teleseismic observations - suggests a confirmation for the proposed northward movement of a forearc sliver acting as a detached continental micro-plate.
Paleoenvironmental records provide ample information on the Late Quaternary climatic evolution. Due to the great diversity of continental mid-latitude environments the synthetic picture of the past mid-latitudinal climate changes is, however, far from being complete. Owing to its significant size and landlocked setting the Black Sea constitutes a perfect location to study patterns and mechanisms of climate change along the continental interior of Central and Eastern Europe and Asia Minor. Presently, the southern drainage area of the Black Sea is characterized by a Mediterranean-type climate while the northern drainage is under the influence of Central and Northern European climate. During the Last Glacial a decrease in the global sea level disconnected the Black Sea from the Mediterranean Sea transforming it into a giant closed lake. At that time atmospheric precipitation and related with it river run-off were the most important factors driving sediment supply and water chemistry of the Black ‘Lake’. Therefore studying properties of the Black Sea sediments provides important information on the interactions and development of the Mediterranean and Central and North European climate in the past. One significant outcome of my thesis is an improved chronostraphigraphical framework for the glacial lacustrine unit of the Black Sea sediment cores, which allowed to refine the environmental history of the Black Sea region and enabled a reliable correlation with data from other marine and terrestrial archives. Data gathered along a N-S transect presented on a common time scale revealed coherent changes in the basin and its surrounding. During the glacial, the southward-shifted Polar Front reduced moisture transport to the northern drainage of the Black Sea and let the southern drainage become dominant in freshwater and sediment supply into the basin. Changes in NW Anatolian precipitation reconstructed from the variability of the terrigenous input imply that during the glacial the regional rainfall variability was strongly influenced by Mediterranean sea surface temperatures and decreased in response to the cooling associated with the North Atlantic Heinrich Events H1 and H2. In contrast to regional precipitation changes, the hydrological properties of the Black Sea remained relatively stable under full glacial conditions. First significant modification in the freshwater/sediment sources reconstructed from changes in the sediment composition, lithology, and 18O of ostracods took place at around 16.4 cal ka BP, simultaneous to the early deglacial northward retreat of the oceanic and atmospheric polar fronts. Meltwater pulses, most probably derived from the disintegrating European ice sheets, changed the isotopic composition of the Black Sea and increased the supply from northern sediment sources. While these changes signalized a mitigation of the Northern European and Mediterranean climate, a decisive increase in local temperature was indicated only later at the transition from the Oldest Dryas to the Bølling around 14.6 cal ka BP. At that time the warming of the Black Sea surface initiated massive phytoplankton blooms, which in turn, induced the precipitation of inorganic carbonates. This biologically triggered process significantly changed the water chemistry and was recorded by simultaneous shifts in the elemental composition of ostracod shells and in the isotopic composition of the inorganically-precipitated carbonates. Starting with the B/A warming and continuing through the YD cold interval and the Early Holocene warming, the Black Sea temperature signal corresponds to the precipitation and temperature changes recorded in the wider Mediterranean region. Early Holocene conditions, similar to those of the Bølling/Allerød, were punctured by the marine inflow from the Mediterranean at ~ 9.3 cal ka BP, which terminated the lacustrine phase of the Black Sea and had a substantial impact on the chemical and physical properties of its water.
The styles of deformation of the fore-arc wedges along the Chilean convergent margin are observed to vary significantly, despite similar plate kinematic conditions. Here, I focus on the analysis of fore-arc deformation on two regions along the Chilean convergent margin at 20°-24°S and 37°-42°S. Although both regions are subjected to the oblique subduction of the oceanic Nazca plate and backstopped by the Andes mountain chain; they display different patterns of deformation. The northern Chilean study area (20° - 24°S) is characterized by an exceptionally thick crust of about 60 km beneath the Altiplano – Puna plateau, lack of an accretionary wedge in the fore-arc due to hyperarid climate, and consequently a sediment starved trench. Two major margin parallel strike slip faults are observed in this area, the Atacama Fault Zone (AFZ) and the Precordilleran Fault System (PFS). Both strike-slip faults do not exhibit significant recent displacement. The southern study area (37° - 42°S), compared to the northern study area, is characterized by lower topography, high precipitation rates (~2000 mm/yr), and a younger subducted oceanic plate. An active strike-slip fault, the Liquiñe-Ofqui-Fault-Zone (LOFZ), shows ~1 cm/yr recent dextral movement and shapes the surface of this area. Thus, the southern Chilean study area exhibits localized strike-slip motion. Within this area the largest earthquake ever recorded, the 1960 Valdivia earthquake, occurred with a moment magnitude of MW=9.5. I have constructed 2D thermal models and 3D mechanical models for both Chilean study areas to study processes related to active subduction. The applied numerical method is the finite element technique by means of the commercial software package ABAQUS. The thermal models are focused on the thermal conditions along the plate interface. The thermal structure along the plate interface reveals the limits of coupling but also the type of transition from coupled to uncoupled and vice versa. The model results show that shear heating at the plate interface is an important mechanism that should be taken into account. The models also show that the thermal condition at the downdip limit of the coupling zone leads to a sharp decrease of friction along the interface. Due to the different geometries of the two Chilean study areas, such as the slab dip and the thickness of the continental crust, the downdip limit of the southern study area is slightly shallower than that of the northern study area. The results of the 2D thermal models are used to constrain the spatial extent of the coupling zone in the 3D mechanical models. 3D numerical simulations are used to investigate how geometry, rheology and mechanical parameters influence strain partitioning and styles of deformation in the Chilean fore-arc. The general outline of the models is based on the fore-arc geometry and boundary conditions as derived from geophysical and geological field data. I examined the influence of different rheological approaches and varying physical properties of the fore-arc to identify and constrain the parameters controlling the difference in surface deformation between the northern and southern study area. The results of numerical studies demonstrate that a small slab dip, a high coefficient of basal friction, a high obliquity of convergence, and a high Young’s modulus favour localisation of deformation in the fore-arc wedge. This parameter study helped me to constrain preferred models for the two Chilean study areas that fit to first order observations. These preferred models explain the difference in styles of deformation as controlled by the angle of obliquity, the dip of subducting slab, and the strength of wedge material. The difference in styles can be even larger if I apply stronger coupling between plates within the southern area; however, several independent observations indicate opposite tendency showing southward decrease of intensity of coupling. The weaker wedge material of the preferred model for the northern study area is associated with advanced development of the adjacent orogen, the Central Andes. Analysis of world-wide examples of oblique subduction zones supports the conclusion that more mature subduction zones demonstrate less pronounced localization of strike-slip motion.
Mafic magmatism in the Eastern Cordillera and Putumayo Basin, Colombia : causes and consequences
(2007)
The Eastern Cordillera of Colombia is mainly composed of sedimentary rocks deposited since early Mesozoic times. Magmatic rocks are scarce. They are represented only by a few locally restricted occurrences of dykes and sills of mafic composition presumably emplaced in the Cretaceous and of volcanic rocks of Neogene age. This work is focused on the study of the Cretaceous magmatism with the intention to understand the processes causing the genesis of these rocks and their significance in the regional tectonic setting of the Northern Andes. The magmatic rocks cut the Cretaceous sedimentary succession of black shales and marlstones that crop out in both flanks of the Eastern Cordillera. The studied rocks were classified as gabbros (Cáceres, Pacho, Rodrigoque), tonalites (Cáceres, La Corona), diorites and syenodiorites (La Corona), pyroxene-hornblende gabbros (Pacho), and pyroxene-hornblendites (Pajarito). The gabbroic samples are mainly composed of plagioclase, clinopyroxene, and/or green to brown hornblende, whereas the tonalitic rocks are mainly composed of plagioclase and quartz. The samples are highly variable in crystal sizes from fine- to coarse-grained. Accessory minerals such as biotite, titanite and zircon are present. Some samples are characterized by moderate to strong alteration, and show the presence of epidote, actinolite and chlorite. Major and trace element compositions of the rocks as well as the rock-forming minerals show significant differences in the geochemical and petrological characteristics for the different localities, suggesting that this magmatism does not result from a single melting process. The wide compositional spectrum of trace elements in the intrusions is characteristic for different degrees of mantle melting and enrichment of incompatible elements. MORB- and OIB-like compositions suggest at least two different sources of magma with tholeiitic and alkaline affinity, respectively. Evidence of slab-derived fluids can be recognized in the western part of the basin reflected in higher Ba/Nb and Sr/P ratios and also in the Sr radiogenic isotope ratios, which is possible a consequence of metasomatism in the mantle due to processes related to the presence of a previously subducted slab. The trace element patterns evidence an extensional setting in the Cretaceous basin producing a continental rift, with continental crust being stretched until oceanic crust was generated in the last stages of this extension. Electron microprobe analyses (EMPA) of the major elements and synchrotron radiation micro-X-ray fluorescence (μ-SRXRF) analyses of the trace element composition of the early crystallized minerals of the intrusions (clinopyroxenes and amphiboles) reflect the same dual character that has been found in the bulk-rock analyses. Despite the observed alteration of the rocks, the mineral composition shows evidences for an enriched and a relative depleted magma source. Even the normalization of the trace element concentrations of clinopyroxenes and amphiboles to the whole rock nearly follows the pattern predicted by published partition coefficients, suggesting that the alteration did not change the original trace element compositions of the investigated minerals. Sr-Nd-Pb isotope data reveal a large isotopic variation but still suggest an initial origin of the magmas in the mantle. Samples have moderate to highly radiogenic compositions of 143Nd/144Nd and high 87Sr/86Sr ratios and follow a trend towards enriched mantle compositions, like the local South American Paleozoic crust. The melts experienced variable degrees of contamination by sediments, crust, and seawater. The age corrected Pb isotope ratios show two separated groups of samples. This suggests that the chemical composition of the mantle below the Northern Andes has been modified by the interaction with other components resulting in a heterogeneous combination of materials of diverse origins. Although previous K/Ar age dating have shown that the magmatism took place in the Cretaceous, the high error of the analyses and the altered nature of the investigated minerals did preclude reliable interpretations. In the present work 40Ar/39Ar dating was carried out. The results show a prolonged history of magmatism during the Cretaceous over more than 60 Ma, from ~136 to ~74 Ma (Hauterivian to Campanian). Pre-Cretaceous rifting phases occurred in the Triassic-Jurassic for the western part of the basin and in the Paleozoic for the eastern part. Those previous rifting phases are decisive mechanisms controlling the localization and composition of the Cretaceous magmatism. Therefore, it is the structural position and not the age of the intrusions which preconditions the kind of magmatism and the degree of melting. The divergences on ages are the consequence of the segmentation of the basin in several sub-basins which stretching, thermal evolution and subsidence rate evolved independently. The first hypothesis formulated at the beginning of this investigation was that the Cretaceous gabbroic intrusions identified in northern Ecuador could be correlated with the intrusions described in the Eastern Cordillera. The mafic occurrences should mark the location of the most subsiding places of the large Cretaceous basin in northern South America. For this reason, the gabbroic intrusions cutting the Cretaceous succession in the Putumayo Basin, southern Colombia, were investigated. The results of the studies were quite unexpected. The petrologic and geochemical character of the magmatic rocks indicates subduction-related magmatism. K/Ar dating of amphibole yields a Late Miocene to Pliocene age (6.1 ± 0.7 Ma) for the igneous event in the basin. Although there is no correlation between this magmatic event and the Cretaceous magmatic event, the data obtained has significant tectonic and economic implications. The emplacement of the Neogene gabbroic rocks coincides with the late Miocene/Pliocene Andean orogenic uplift as well as with a significant pulse of hydrocarbon generation and expulsion.
The Andean orogen is the most outstanding example of mountain building caused by the subduction of oceanic below continental lithosphere. The Andes formed by the subduction of the Nazca and Antarctic oceanic plates under the South American continent over at least ~200 million years. Tectonic and climatic conditions vary markedly along this north-south–oriented plate boundary, which thus represents an ideal natural laboratory to study tectonic and climatic segmentation processes and their possible feedbacks. Most of the seismic energy on Earth is released by earthquakes in subduction zones, like the giant 1960, Mw 9.5 event in south-central Chile. However, the segmentation mechanisms of surface deformation during and between these giant events have remained poorly understood. The Andean margin is a key area to study seismotectonic processes because of its along-strike variability under similar plate kinematic boundary conditions. Active deformation has been widely studied in the central part of the Andes, but the south-central sector of the orogen has gathered less research efforts. This study focuses on tectonics at the Neogene and late Quaternary time scales in the Main Cordillera and coastal forearc of the south-central Andes. For both domains I document the existence of previously unrecognized active faults and present estimates of deformation rates and fault kinematics. Furthermore these data are correlated to address fundamental mountain building processes like strain partitioning and large-scale segmentation. In the Main Cordillera domain and at the Neogene timescale, I integrate structural and stratigraphic field observations with published isotopic ages to propose four main phases of coupled styles of tectonics and distribution of volcanism and magmatism. These phases can be related to the geometry and kinematics of plate convergence. At the late Pleistocene timescale, I integrate field observations with lake seismic and bathymetric profiles from the Lago Laja region, located near the Andean drainage divide. These data reveal Holocene extensional faults, which define the Lago Laja fault system. This fault system has no significant strike-slip component, contrasting with the Liquiñe-Ofqui dextral intra-arc system to the south, where Holocene strike-slip markers are ubiquitous. This contrast in structural style along the arc is coincident with a marked change in along-strike fault geometries in the forearc, across the Arauco Peninsula. Thereon I propose that a net gradient in the degree of partitioning of oblique subduction occurs across the Arauco transition zone. To the north, the margin parallel component of oblique convergence is distributed in a wide zone of diffuse deformation, while to the south it is partitioned along an intra-arc, margin-parallel strike-slip fault zone. In the coastal forearc domain and at the Neogene timescale, I integrate structural and stratigraphic data from field observations, industry reflection-seismic profiles and boreholes to emphasize the influence of climate-driven filling of the trench on the mechanics and kinematics of the margin. I show that forearc basins in the 34-45°S segment record Eocene to early Pliocene extension and subsidence followed by ongoing uplift and contraction since the late Pliocene. I interpret the first stage as caused by tectonic erosion due to high plate convergence rates and reduced trench fill. The subsequent stage, in turn, is related to accretion caused by low convergence rates and the rapid increase in trench fill after the onset of Patagonian glaciations and climate-driven exhumation at ~6-5 Ma. On the late Quaternary timescale, I integrate off-shore seismic profiles with the distribution of deformed marine terraces from Isla Santa María, dated by the radiocarbon method, to show that inverted reverse faulting controls the coastal geomorphology and segmentation of surface deformation. There, a cluster of microearthquakes illuminates one of these reverse faults, which presumingly reaches the plate interface. Furthermore, I use accounts of coseismic uplift during the 1835 M>8 earthquake made by Charles Darwin, to propose that this active reverse fault has been mechanically coupled to the megathrust. This has important implications on the assessment of seismic hazards in this, and other similar regions. These results underscore the need to study plate-boundary deformation processes at various temporal and spatial scales and to integrate geomorphologic, structural, stratigraphic, and geophysical data sets in order to understand the present distribution and causes of tectonic segmentation.
Earthquakes form by sudden brittle failure of rock mostly as shear ruptures along a rupture plane. Beside this, mechanisms other than pure shearing have been observed for some earthquakes mainly in volcanic areas. Possible explanations include complex rupture geometries and tensile earthquakes. Tensile earthquakes occur by opening or closure of cracks during rupturing. They are likely to be often connected with fluids that cause pressure changes in the pore space of rocks leading to earthquake triggering. Tensile components have been reported for swarm earthquakes in West Bohemia in 2000. The aim and subject of this work is an assessment and the accurate determination of such tensile components for earthquakes in anisotropic media. Currently used standard techniques for the retrieval of earthquake source mechanisms assume isotropic rock properties. By means of moment tensors, equivalent forces acting at the source are used to explain the radiated wavefield. Conversely, seismic anisotropy, i.e. directional dependence of elastic properties, has been observed in the earth's crust and mantle such as in West Bohemia. In comparison to isotropy, anisotropy causes modifications in wave amplitudes and shear-wave splitting. In this work, effects of seismic anisotropy on true or apparent tensile source components of earthquakes are investigated. In addition, earthquake source parameters are determined considering anisotropy. It is shown that moment tensors and radiation patterns due to shear sources in anisotropic media may be similar to those of tensile sources in isotropic media. In contrast, similarities between tensile earthquakes in anisotropic rocks and shear sources in isotropic media may exist. As a consequence, the interpretation of tensile source components is ambiguous. The effects that are due to anisotropy depend on the orientation of the earthquake source and the degree of anisotropy. The moment of an earthquake is also influenced by anisotropy. The orientation of fault planes can be reliably determined even if isotropy instead of anisotropy is assumed and if the spectra of the compressional waves are used. Greater difficulties may arise when the spectra of split shear waves are additionally included. Retrieved moment tensors show systematic artefacts. Observed tensile source components determined for events in West Bohemia in 1997 can only partly be attributed to the effects of moderate anisotropy. Furthermore, moment tensors determined earlier for earthquakes induced at the German Continental Deep Drilling Program (KTB), Bavaria, were reinterpreted under assumptions of anisotropic rock properties near the borehole. The events can be consistently identified as shear sources, although their moment tensors comprise tensile components that are considered to be apparent. These results emphasise the necessity to consider anisotropy to uniquely determine tensile source parameters. Therefore, a new inversion algorithm has been developed, tested, and successfully applied to 112 earthquakes that occurred during the most recent intense swarm episode in West Bohemia in 2000 at the German-Czech border. Their source mechanisms have been retrieved using isotropic and anisotropic velocity models. Determined local magnitudes are in the range between 1.6 and 3.2. Fault-plane solutions are similar to each other and characterised by left-lateral faulting on steeply dipping, roughly North-South oriented rupture planes. Their dip angles decrease above a depth of about 8.4km. Tensile source components indicating positive volume changes are found for more than 60% of the considered earthquakes. Their size depends on source time and location. They are significant at the beginning of the swarm and at depths below 8.4km but they decrease in importance later in the course of the swarm. Determined principle stress axes include P axes striking Northeast and Taxes striking Southeast. They resemble those found earlier in Central Europe. However, depth-dependence in plunge is observed. Plunge angles of the P axes decrease gradually from 50° towards shallow angles with increasing depth. In contrast, the plunge angles of the T axes change rapidly from about 8° above a depth of 8.4km to 21° below this depth. By this thesis, spatial and temporal variations in tensile source components and stress conditions have been reported for the first time for swarm earthquakes in West Bohemia in 2000. They also persist, when anisotropy is assumed and can be explained by intrusion of fluids into the opened cracks during tensile faulting.
Lakustrine Sedimente als Archive des spätquartären Umweltwandels in der Amery-Oase, Ostantarktis
(2006)
Im Rahmen einer deutsch-australischen Forschungskooperation erfolgte im Südsommer 2001/2002 eine Expedition in die Amery-Oase (70°50’S, 68°00’E), die im Einzugsgebiet des Lambert-Gletscher/Amery-Schelfeis-Systems, dem größten ostantarktischen Eis-Drainagesystem, liegt. Von deutscher Seite wurden im Zuge der Geländekampagne erstmals lakustrine Sedimentsequenzen gewonnen, um die bislang wenig erforschte spätquartäre Klima- und Umweltgeschichte dieser rund 1800 km<sup>2 großen eisfreien Region zu rekonstruieren. Die drei untersuchten Glazialseen Beaver, Radok und Terrasovoje unterscheiden sich sowohl deutlich in ihrer Größe, Bathymetrie und den hydrologischen Merkmalen sowie in ihren Sedimentabfolgen. Einen Schwerpunkte dieser Doktorarbeit bildet die Rekonstruktion der Sedimentationsprozesse und des Ablagerungsmilieus sowie Untersuchungen zur Herkunft des detritischen Sedimentmaterials in den Seebecken. Der methodische Ansatz verfolgt die Charakterisierung der klastischen Sedimentfazies an Hand lithologisch-granulometrischer Merkmale sowie mineralogisch-geochemischer Analysen der Sedimentherkunft. Ein weiterer Schwerpunkt ist die Rekonstruktion der holozänen biogen gesteuerten Ablagerungsbedingungen im Terrasovoje-See, die Rückschlüsse auf den kurzfristigen postglazialen Klima- und Umweltwandel in der Amery-Oase gestattet. Dabei wurden mikrofazielle Untersuchungsmethoden und hochauflösende Elementscannermessungen angewandt. Die klastische Sedimentherkunft in den drei Seen unterscheidet sich räumlich deutlich voneinander und spiegelt den komplexen geologischen Aufbau der Amery-Oase wider. Als Sedimentquellen konnten präkambrische Metamorphite, permotriassische Sedimentgesteine und tertiäre Lockersedimente identifiziert werden. Die Varibilität der Herkunftssignale ist zeitlich weniger deutlich als räumlich ausgeprägt und deutet auf relativ konstante Liefergebiete in den einzelnen Seen hin. Das glaziolakustrine Ablagerungsmilieu der drei untersuchten Seen zeigt klare räumliche und zeitliche Unterschiede. In allen drei Seen setzen sich die älteren Sedimente aus grobkörnigem, häufig diamiktischem Material zusammen, während die jüngeren Sedimente aus feinkörnigen Laminiten bestehen. Die lithofazielle Zweiteilung in den Sedimentabfolgen deutet auf einen Rückzug der Gletscher und/oder einen Anstieg der Wassertiefen im Übergang von den grobkörnigen zu den feinkörnigen Ablagerungseinheiten hin. Die oberen feinkörnigen Kernabschnitte spiegeln in allen drei Seen die postglaziale lakustrine Sedimentation wider. Im Beaver-See wird die postglaziale Fazies durch laminierte klastische Stillwassersedimente repräsentiert, im Radok-See durch Turbiditsequenzen und im Terrasovoje-See durch Algenlaminite. Abgesehen vom Terrasovoje-See ist die zeitliche Einordnung der Fazieswechsel auf Grund mangelnder Altersinformationen schwer erfassbar. Im Terrasovoje-See setzte die postglaziale Sedimentation um rund 12,4 cal. ka ein. Somit weisen die darunterliegenden glazigenen Klastika mindestens ein spätpleistozänes Alter auf. Die sedimentologischen Eigenschaften, Änderungen der Sedimentationsraten und organogene Zusammensetzung der postglazialen Biogenlaminite des Terrasovoje-Sees deuten auf Variationen der paläolimnologischen Bedingungen hinsichtlich Eisbedeckung, biologischer Produktivität, Wasserstand, Redoxbedingungen und Salinität hin, die mit regionalen holozänen Klimaänderungen in Verbindung gebracht werden können. Weitere Anhaltspunkte ergeben sich aus der Zusammensetzung und den Mächtigkeitsvariationen der Laminae, die generell aus Wechsellagerungen von Cyanobakterienmatten mit feinklastischen Lagen bestehen. Lagenzählungen der Laminae belegen Änderungen des Ablagerungsmilieus auf subdekadischen Zeitskalen, wobei zeitweilige jährliche Signale nicht ausgeschlossen werden können. Unter Berücksichtigung aller faziellen Indikatoren lässt sich aus der Sedimentabfolge des Terrasovoje-Sees ein frühholozänes Klimaoptimum zwischen 9 und 7 cal. ka sowie weitere Wärmephasen zwischen 3,2 und 2,3 cal. ka sowie 1,5 und 1,0 cal. ka ableiten. Im Vergleich mit Eiskernarchiven und anderen Seesedimentabfolgen aus ostantarktischen Oasen zeigt sich, dass das Auftreten postglazialer Warmphasen nicht allenorts einem allgemein gültigen räumlich-zeitlichen Muster folgt. Die Ursachen hierfür liegen vermutlich in den lokalen geographischen Gegebenheiten. Es lässt sich daraus schliessen, dass die bisher vorliegenden Klimarekonstruktionen eher das Lokalklima an einem Untersuchungsstandort als das Großklima der Ostantarktis reflektieren. Daraus ergibt sich die Notwendigkeit weiterer Untersuchungen von antarktischen Klimaarchiven und Untersuchungsstandorten, um örtliche von überregionalen Klimasignalen besser unterscheiden zu können.
The intracontinental endorheic Aral Sea, remote from oceanic influences, represents an excellent sedimentary archive in Central Asia that can be used for high-resolution palaeoclimate studies. We performed palynological, microfacies and geochemical analyses on sediment cores retrieved from Chernyshov Bay, in the NW part of the modern Large Aral Sea. The most complete sedimentary sequence, whose total length is 11 m, covers approximately the past 2000 years of the late Holocene. High-resolution palynological analyses, conducted on both dinoflagellate cysts assemblages and pollen grains, evidenced prominent environmental change in the Aral Sea and in the catchment area. The diversity and the distribution of dinoflagellate cysts within the assemblages characterized the sequence of salinity and lake-level changes during the past 2000 years. Due to the strong dependence of the Aral Sea hydrology to inputs from its tributaries, the lake levels are ultimately linked to fluctuations in meltwater discharges during spring. As the amplitude of glacial meltwater inputs is largely controlled by temperature variations in the Tien Shan and Pamir Mountains during the melting season, salinity and lake-level changes of the Aral Sea reflect temperature fluctuations in the high catchment area during the past 2000 years. Dinoflagellate cyst assemblages document lake lowstands and hypersaline conditions during ca. 0–425 AD, 920–1230 AD, 1500 AD, 1600–1650 AD, 1800 AD and since the 1960s, whereas oligosaline conditions and higher lake levels prevailed during the intervening periods. Besides, reworked dinoflagellate cysts from Palaeogene and Neogene deposits happened to be a valuable proxy for extreme sheet-wash events, when precipitation is enhanced over the Aral Sea Basin as during 1230–1450 AD. We propose that the recorded environmental changes are related primarily to climate, but may have been possibly amplified during extreme conditions by human-controlled irrigation activities or military conflicts. Additionally, salinity levels and variations in solar activity show striking similarities over the past millennium, as during 1000–1300 AD, 1450–1550 and 1600–1700 AD when low lake levels match well with an increase in solar activity thus suggesting that an increase in the net radiative forcing reinforced past Aral Sea’s regressions. On the other hand, we used pollen analyses to quantify changes in moisture conditions in the Aral Sea Basin. High-resolution reconstruction of precipitation (mean annual) and temperature (mean annual, coldest versus warmest month) parameters are performed using the “probability mutual climatic spheres” method, providing the sequence of climate change for the past 2000 years in western Central Asia. Cold and arid conditions prevailed during ca. 0–400 AD, 900–1150 AD and 1500–1650 AD with the extension of xeric vegetation dominated by steppe elements. Conversely, warmer and less arid conditions occurred during ca. 400–900 AD and 1150–1450 AD, where steppe vegetation was enriched in plants requiring moister conditions. Change in the precipitation pattern over the Aral Sea Basin is shown to be predominantly controlled by the Eastern Mediterranean (EM) cyclonic system, which provides humidity to the Middle East and western Central Asia during winter and early spring. As the EM is significantly regulated by pressure modulations of the North Atlantic Oscillation (NAO) when the system is in a negative phase, a relationship between humidity over western Central Asia and the NAO is proposed. Besides, laminated sediments record shifts in sedimentary processes during the late Holocene that reflect pronounced changes in taphonomic dynamics. In Central Asia, the frequency of dust storms occurring during spring when the continent is heating up is mostly controlled by the intensity and the position of the Siberian High (SH) Pressure System. Using titanium (Ti) content in laminated sediments as a proxy for aeolian detrital inputs, changes in wind dynamics over Central Asia is documented for the past 1500 years, offering the longest reconstruction of SH variability to date. Based on high Ti content, stronger wind dynamics are reported from 450–700 AD, 1210–1265 AD, 1350–1750 AD and 1800–1975 AD, reporting a stronger SH during spring. In contrast, lower Ti content from 1750–1800 AD and 1980–1985 AD reflect a diminished influence of the SH and a reduced atmospheric circulation. During 1180–1210 AD and 1265–1310 AD, considerably weakened atmospheric circulation is evidenced. As a whole, though climate dynamics controlled environmental changes and ultimately modulated changes in the western Central Asia’s climate system, it is likely that changes in solar activity also had an impact by influencing to some extent the Aral Sea’s hydrology balance and also regional temperature patterns in the past. <hr> The appendix of the thesis is provided via the HTML document as ZIP download.
One of the most difficult issues when dealing with optical water remote-sensing is its acceptance as a useful application for environmental research. This problem is, on the one hand, concerned with the optical complexity and variability of the investigated natural media, and therefore the question arises as to the plausibility of the parameters derived from remote-sensing techniques. Detailed knowledge about the regional bio- and chemico-optical properties is required for such studies, however such information is seldom available for the sites of interest. On the other hand, the primary advantage of remote-sensing information, which is the provision of a spatial overview, may not be exploited fully by the disciplines that would benefit most from such information. It is often seen in a variety of disciplines that scientists have been primarily trained to look at discrete data sets, and therefore have no experience of incorporating information dealing with spatial heterogeneity. In this thesis, the opportunity was made available to assess the potential of Ocean Colour data to provide spatial and seasonal information about the surface waters of Lake Baikal (Siberia). While discrete limnological field data is available, the spatial extension of Lake Baikal is enormous (ca. 600 km), while the field data are limited to selected sites and expedition time windows. Therefore, this remote-sensing investigation aimed to support a multi-disciplinary limnological investigation within the framework of the paleoclimate EU-project ‘High Resolution CONTINENTal Paleoclimate Record in Lake Baikal, Siberia (CONTINENT)’ using spatial and seasonal information from the SeaWiFS satellite (NASA). From this, the SeaWiFS study evolved to become the first efficient bio-optical satellite study of Lake Baikal. During the course of three years, field work including spectral field measurements and water sampling, was carried out at Lake Baikal in Southern Siberia, and at the Mecklenburg and Brandenburg lake districts in Germany. The first step in processing the SeaWiFS satellite data involved adapting the SeaDAS (NASA) atmospheric-correction processing to match as close as possible the specific conditions of Lake Baikal. Next, various Chl-a algorithms were tested on the atmospherically-corrected optimized SeaWiFS data set (years 2001 to 2002), comparing the CONTINENT pigment ground-truth data with the Chl-a concentrations derived from the satellite data. This showed the high performance of the global Chl-a products OC2 and OC4 for the oligotrophic, transparent waters (bio-optical Case 1) of Lake Baikal. However, considerable Chl-a overestimation prevailed in bio-optical Case 2 areas for the case of discharge events. High-organic terrigenous input into Lake Baikal could be traced and information extracted using the SeaWiFS spectral data. Suspended Particulate Matter (SPM) was quantified by the regression of the SeaDAS attenuation coefficient as the optical parameter with SPM field data. Finally, the Chl-a and terrigenous input maps derived from the remote sensing data were used to assist with analyzing the relationships between the various discrete data obtained during the CONTINENT field work. Hence, plausible spatial and seasonal information describing autochthonous and allochthonous material in Lake Baikal could be provided by satellite data.Lake Baikal, with its bio-optical complexity and its different areas of Case 1 and Case 2 waters, is a very interesting case study for Ocean Colour analyses. Proposals for future Ocean Colour studies of Lake Baikal are discussed, including which bio-optical parameters for analytical models still need to be clarified by field investigations.
Subduction zones are regions of intense earthquake activity up to great depth. Sources are located inside the subducting lithosphere and, as a consequence, seismic radiation from subduction zone earthquakes is strongly affected by the interior slab structure. The wave field of these intraslab events observed in the forearc region is profoundly influenced by a seismically slow layer atop the slab surface. This several kilometer thick low-velocity channel (wave guide) causes the entrapment of seismic energy producing strong guided wave phases that appear in P onsets in certain regions of the forearc. Observations at the Chile-Peru subduction zone presented here, as well as observations at several other circum-pacific subduction zones show such signals. Guided wave analysis contributes details of immense value regarding the processes near the slab surface, such as layering of subducted lithosphere, source locations of intraslab seismicity and most of all, range and manner of mineralogical phase transitions. Seismological data stem from intermediate depth events (depth range 70 km - 300 km) recorded in northern Chile near 21 Grad S during the collaborative research initiative " Deformation Processes in the Andes" (SFB 267). A subset of stations - all located within a slab-parallel transect close to 69 Grad W - show low-frequency first arrivals (2 Hz), sometimes followed by a second high-frequency phase. We employ 2-dimensional finite-difference simulations of complete P-SV wave propagation to explore the parameter space of subduction zone wave guides and explain the observations. Key processes underlying the guided wave propagation are studied: Two distinct mechanisms of decoupling of trapped energy from the wave guide are analyzed - a prerequisite to observe the phases at stations located at large distances from the wave guide (up to 100 km). Variations of guided wave effects perpendicular to the strike of the subduction zone are investigated, such as the influence of phases traveling in the fast slab. Further, the merits and limits of guided wave analysis are assessed. Frequency spectra of the guided wave onsets prove to be a robust quantity that captures guided wave characteristics at subduction zones including higher mode excitation. They facilitate the inference of wave guide structure and source positioning: The peak frequency of the guided wave fundamental mode is associated with a certain combination of layer width and velocity contrast. The excitation strength of the guided wave fundamental mode and higher modes is associated with source position and orientation relative to the low-velocity layer. The guided wave signals at the Chile-Peru subduction zone are caused by energy that leaks from the subduction zone wave guide. On the one hand, the bend shape of the slab allows for leakage at a depth of 100 km. On the other, equalization of velocities between the wave guide and the host rocks causes further energy leakage at the contact zone between continental and oceanic crust (70 km depth). Guided waves bearing information on deep slab structure can therefore be recorded at specific regions in the forearc. These regions are determined based on slab geometry, and their locations coincide with the observations. A number of strong constraints on the structure of the Chile-Peru slab are inferred: The deep wave guide for intraslab events is formed by a layer of 2 km average width that remains seismically slow (7 percent velocity reduction compared to surrounding mantle). This low-velocity layer at the top of the Chile-Peru slab is imaged from a depth of 100 km down to at least 160 km. Intermediate depth events causing the observed phases are located inside the layer or directly beneath it in the slab mantle. The layer is interpreted as partially eclogized lower oceanic crust persisting to depth beyond the volcanic arc.
About 24 % of the land surface in the northern hemisphere are underlayed by permafrost in various states. Permafrost aggradation occurs under special environmental conditions with overall low annual precipitation rates and very low mean annual temperatures. Because the general permafrost occurrence is mainly driven by large-scale climatic conditions, the distribution of permafrost deposits can be considered as an important climate indicator. The region with the most extensive continuous permafrost is Siberia. In northeast Siberia, the ice- and organic-rich permafrost deposits of the Ice Complex are widely distributed. These deposits consist mostly of silty to fine-grained sandy sediments that were accumulated during the Late Pleistocene in an extensive plain on the then subaerial Laptev Sea shelf. One important precondition for the Ice Complex sedimentation was, that the Laptev Sea shelf was not glaciated during the Late Pleistocene, resulting in a mostly continuous accumulation of permafrost sediments for at least this period. This shelf landscape became inundated and eroded in large parts by the Holocene marine transgression after the Last Glacial Maximum. Remnants of this landscape are preserved only in the present day coastal areas. Because the Ice Complex deposits contain a wide variety of palaeo-environmental proxies, it is an excellent palaeo-climate archive for the Late Quaternary in the region. Furthermore, the ice-rich Ice Complex deposits are sensible to climatic change, i.e. climate warming. Because of the large-scale climatic changes at the transition from the Pleistocene to the Holocene, the Ice Complex was subject to extensive thermokarst processes since the Early Holocene. Permafrost deposits are not only an environmental indicator, but also an important climate factor. Tundra wetlands, which have developed in environments with aggrading permafrost, are considered a net sink for carbon, as organic matter is stored in peat or is syn-sedimentary frozen with permafrost aggradation. Contrary, the Holocene thermokarst development resulted in permafrost degradation and thus the release of formerly stored organic carbon. Modern tundra wetlands are also considered an important source for the climate-driving gas methane, originating mainly from microbial activity in the seasonal active layer. Most scenarios for future global climate development predict a strong warming trend especially in the Arctic. Consequently, for the understanding of how permafrost deposits will react and contribute to such scenarios, it is necessary to investigate and evaluate ice-rich permafrost deposits like the widespread Ice Complex as climate indicator and climate factor during the Late Quaternary. Such investigations are a pre-condition for the precise modelling of future developments in permafrost distribution and the influence of permafrost degradation on global climate. The focus of this work, which was conducted within the frame of the multi-disciplinary joint German-Russian research projects "Laptev Sea 2000" (1998-2002) and "Dynamics of Permafrost" (2003-2005), was twofold. First, the possibilities of using remote sensing and terrain modelling techniques for the observation of periglacial landscapes in Northeast Siberia in their present state was evaluated and applied to key sites in the Laptev Sea coastal lowlands. The key sites were situated in the eastern Laptev Sea (Bykovsky Peninsula and Khorogor Valley) and the western Laptev Sea (Cape Mamontovy Klyk region). For this task, techniques using CORONA satellite imagery, Landsat-7 satellite imagery, and digital elevation models were developed for the mapping of periglacial structures, which are especially indicative of permafrost degradation. The major goals were to quantify the extent of permafrost degradation structures and their distribution in the investigated key areas, and to establish techniques, which can be used also for the investigation of other regions with thermokarst occurrence. Geographical information systems were employed for the mapping, the spatial analysis, and the enhancement of classification results by rule-based stratification. The results from the key sites show, that thermokarst, and related processes and structures, completely re-shaped the former accumulation plain to a strongly degraded landscape, which is characterised by extensive deep depressions and erosional remnants of the Late Pleistocene surface. As a results of this rapid process, which in large parts happened within a short period during the Early Holocene, the hydrological and sedimentological regime was completely changed on a large scale. These events resulted also in a release of large amounts of organic carbon. Thermokarst is now the major component in the modern periglacial landscapes in terms of spatial extent, but also in its influence on hydrology, sedimentation and the development of vegetation assemblages. Second, the possibilities of using remote sensing and terrain modelling as a supplementary tool for palaeo-environmental reconstructions in the investigated regions were explored. For this task additionally a comprehensive cryolithological field database was developed for the Bykovsky Peninsula and the Khorogor Valley, which contains previously published data from boreholes, outcrops sections, subsurface samples, and subsurface samples, as well as additional own field data. The period covered by this database is mainly the Late Pleistocene and the Holocene, but also the basal deposits of the sedimentary sequence, interpreted as Pliocene to Early Pleistocene, are contained. Remote sensing was applied for the observation of periglacial strucures, which then were successfully related to distinct landscape development stages or time intervals in the investigation area. Terrain modelling was used for providing a general context of the landscape development. Finally, a scheme was developed describing mainly the Late Quaternary landscape evolution in this area. A major finding was the possibility of connecting periglacial surface structures to distinct landscape development stages, and thus use them as additional palaeo-environmental indicator together with other proxies for area-related palaeo-environmental reconstructions. In the landscape evolution scheme, i.e. of the genesis of the Late Pleistocene Ice Complex and the Holocene thermokarst development, some new aspects are presented in terms of sediment source and general sedimentation conditions. This findings apply also for other sites in the Laptev Sea region.
Diagenetic studies of carbonate rocks focused for a long time on photozoan carbonate assemblages deposited in tropical climates. The results of these investigations were taken as models for the diagenetic evolution of many fossil carbonates. Only in recent years the importance of heterozoan carbonates, generally formed out of the tropics or in deeper waters, was realized. Diagenetic studies focusing on this kind of rocks are still scarce, but indicate that the diagenetic evolution of these rocks might be a better model for many fossil carbonate settings ("calcite-sea" carbonates) than the photozoan model used before. This study deals with the determination of the diagenetic pathways and environments in such shallow-water heterozoan carbonate assemblages. Special emphasis is put on the identification of early, near-seafloor diagenetic processes and on the evaluation of the amount of constructive diagenesis in form of cementation in this diagenetic environment. As study area the Central Mediterranean, the Maltese Islands and Sicily, was chosen. Here two sections were logged in Olio-Miocene shallow-water carbonates consisting of different kinds of heterozoan assemblages. The study area is very suitable for the investigation of constructive early diagenetic processes, as the rocks were never deeply buried and burial diagenetic pressure solution and cementation as cause of lithification could be ruled out. Nevertheless, the carbonate rocks are well lithified and form steep cliffs, implying cementation/lithification in another, shallower diagenetic environment. To determine the diagenetic pathways and environments, detailed transmitted light and cathodoluminescence petrography was carried out on thin sections. Furthermore the stable isotope (δ<sup>18O and δ<sup>13C) composition of the bulk rock, single biota and single cement phases was determined, as well as the major and trace element composition of the single cement phases. Petrographically three (Sicily) to four (Maltese Islands) cementation phases, two phases of fabric selective and one of non-fabric selective dissolution, one phase of neomorphism and one of chemical compaction could be distinguished. The stable isotope measurements of the single cement phases pointed to cement precipitation from marine, marine-derived and meteoric waters. The trace element analysis indicated precipitation under reducing conditions, (A) in an open system with low rock-water interaction on the Maltese Islands and (B) in a closed system with high rock-water interaction on Sicily. For the closed systems case, aragonite as cement source could be concluded because its chemical composition was preserved in the newly formed cements. By integrating these results, diagenetic pathways and environments for the investigated locations were established, and the cement source(s) in the different environments were determined. The diagenetic evolution started in the marine environment with the precipitation of fibrous/fibrous-bladed and epitaxial cement I. These cements formed as High Mg Calcite (HMC) directly out of marine waters. The paleoenvironmentally shallowest part of the section on the Maltese Islands was also exposed to meteoric diagenetic fluids. This meteoric influence lead to the dissolution of aragonitic and HMC skeletons, which sourced the cementation by Low Mg Calcitic (LMC) epitaxial cement II in this part of the Maltese section. Entering the burial-marine environment the main part of dissolution, cementation and neomorphism started to take place. The elevated CO2 content in this environment, caused by the decay of organic matter, lead to the dissolution of aragonitic skeletons, which sourced the cementation by LMC epitaxial cement II, bladed and blocky cements. The earlier precipitated HMC cement phases were either partly dissolved (epitaxial cement I) or neomorphosed to LMC (fibrous/fibrous-bladed and epitaxial cement I). In the burial environment weak chemical compaction took place without sourcing significant amounts of cementation. In a last phase the rocks entered the meteoric realm by uplift, which caused non-fabric selective dissolution. This study shows that early diagenetic processes, taking place at or just below the sediment-water-interface, are very important for the mineralogical stabilization of heterozoan carbonate strata. The main amount of constructive diagenesis in form of cementation takes place in this environment, sourced by dissolution of aragonitic and, to a lesser degree, of HMC skeletons. The results of this study imply that the primary amount of aragonitic skeletons in heterozoan carbonate sediments must be carefully assessed, as they are the main early diagenetic cement source. In fossil heterozoan carbonate rocks, aragonitic skeletons might be the cement source even when no relict structures like micritic envelops or biomolds are preserved. In general, the diagenetic evolution of heterozoan carbonate rocks is a good model for the diagenesis of "calcite-sea" time carbonate rocks.
Fault planes of large earthquakes incorporate inhomogeneous structures. This can be observed in teleseismic studies through the spatial distribution of slip and seismic moment release caused by the mainshock. Both parameters are often concentrated on patches on the fault plane with much higher values for slip and moment release than their adjacent areas. These patches are called asperities which obviously have a strong influence on the mainshock rupture propagation. Condition and properties of structures in the fault plane area, which are responsible for the evolution of such asperities or their significance on damage distributions of future earthquakes, are still not well understood and subject to recent geo-scientific studies. In the presented thesis asperity structures are identified on the fault plane of the Mw=8.0 Antofagasta earthquake in northern Chile which occurred on 30th of July, 1995. It was a thrust-type event in the seismogenic zone between the subducting pacific Nazca plate and the overriding South American plate. In cooperation of the German Task Force for Earthquakes and the CINCA'95 project a network of up to 44 seismic stations was set up to record the aftershock sequence. The seaward extension of the network with 9 OBH stations increased significantly the precision of hypocenter determinations. They were distributed mainly on the fault plane itself around the city of Antofagasta and Mejillones Peninsula. The asperity structures were recognized here by the spatial variations of local seismological parameters; at first by the spatial distribution of the seismic b-value on the fault plane, derived from the magnitude-frequency relation of Gutenberg-Richter. The correlation of this b-value map with other parameters like the mainshock source time function, the gravity isostatic residual anomalies, the aftershock radiated seismic energy distribution and the vp/vs ratios from a local earthquake tomograhpy study revealed some ideas about the composition and asperity generating processes. The investigation of 295 aftershock focal mechanism solutions supported the resulting fault plane structure and proposed a 3D similar stress state in the area of the Antofagasta fault plane.
This thesis work describes a new experimental method for the determination of Mode II (shear) fracture toughness, KIIC of rock and compares the outcome to results from Mode I (tensile) fracture toughness, KIC, testing using the International Society of Rock Mechanics Chevron-Bend method.Critical Mode I fracture growth at ambient conditions was studied by carrying out a series of experiments on a sandstone at different loading rates. The mechanical and microstructural data show that time- and loading rate dependent crack growth occurs in the test material at constant energy requirement.The newly developed set-up for determination of the Mode II fracture toughness is called the Punch-Through Shear test. Notches were drilled to the end surfaces of core samples. An axial load punches down the central cylinder introducing a shear load in the remaining rock bridge. To the mantle of the cores a confining pressure may be applied. The application of confining pressure favours the growth of Mode II fractures as large pressures suppress the growth of tensile cracks.Variation of geometrical parameters leads to an optimisation of the PTS- geometry. Increase of normal load on the shear zone increases KIIC bi-linear. High slope is observed at low confining pressures; at pressures above 30 MPa low slope increase is evident. The maximum confining pressure applied is 70 MPa. The evolution of fracturing and its change with confining pressure is described.The existence of Mode II fracture in rock is a matter of debate in the literature. Comparison of the results from Mode I and Mode II testing, mainly regarding the resulting fracture pattern, and correlation analysis of KIC and KIIC to physico-mechanical parameters emphasised the differences between the response of rock to Mode I and Mode II loading. On the microscale, neither the fractures resulting from Mode I the Mode II loading are pure mode fractures. On macroscopic scale, Mode I and Mode II do exist.
The role of feedback between erosional unloading and tectonics controlling the development of the Himalaya is a matter of current debate. The distribution of precipitation is thought to control surface erosion, which in turn results in tectonic exhumation as an isostatic compensation process. Alternatively, subsurface structures can have significant influence in the evolution of this actively growing orogen. Along the southern Himalayan front new 40Ar/39Ar white mica and apatite fission track (AFT) thermochronologic data provide the opportunity to determine the history of rock-uplift and exhumation paths along an approximately 120-km-wide NE-SW transect spanning the greater Sutlej region of the northwest Himalaya, India. 40Ar/39Ar data indicate, consistent with earlier studies that first the High Himalayan Crystalline, and subsequently the Lesser Himalayan Crystalline nappes were exhumed rapidly during Miocene time, while the deformation front propagated to the south. In contrast, new AFT data delineate synchronous exhumation of an elliptically shaped, NE-SW-oriented ~80 x 40 km region spanning both crystalline nappes during Pliocene-Quaternary time. The AFT ages correlate with elevation, but show within the resolution of the method no spatial relationship to preexisting major tectonic structures, such as the Main Central Thrust or the Southern Tibetan Fault System. Assuming constant exhumation rates and geothermal gradient, the rocks of two age vs. elevation transects were exhumed at ~1.4 ±0.2 and ~1.1 ±0.4 mm/a with an average cooling rate of ~50-60 °C/Ma during Pliocene-Quaternary time. The locus of pronounced exhumation defined by the AFT data coincides with a region of enhanced precipitation, high discharge, and sediment flux rates under present conditions. We therefore hypothesize that the distribution of AFT cooling ages might reflect the efficiency of surface processes and fluvial erosion, and thus demonstrate the influence of erosion in localizing rock-uplift and exhumation along southern Himalayan front, rather than encompassing the entire orogen.Despite a possible feedback between erosion and exhumation along the southern Himalayan front, we observe tectonically driven, crustal exhumation within the arid region behind the orographic barrier of the High Himalaya, which might be related to and driven by internal plateau forces. Several metamorphic-igneous gneiss dome complexes have been exhumed between the High Himalaya to the south and Indus-Tsangpo suture zone to the north since the onset of Indian-Eurasian collision ~50 Ma ago. Although the overall tectonic setting is characterized by convergence the exhumation of these domes is accommodated by extensional fault systems.Along the Indian-Tibetan border the poorly described Leo Pargil metamorphic-igneous gneiss dome (31-34°N/77-78°E) is located within the Tethyan Himalaya. New field mapping, structural, and geochronologic data document that the western flank of the Leo Pargil dome was formed by extension along temporally linked normal fault systems. Motion on a major detachment system, referred to as the Leo Pargil detachment zone (LPDZ) has led to the juxtaposition of low-grade metamorphic, sedimentary rocks in the hanging wall and high-grade metamorphic gneisses in the footwall. However, the distribution of new 40Ar/39Ar white mica data indicate a regional cooling event during middle Miocene time. New apatite fission track (AFT) data demonstrate that subsequently more of the footwall was extruded along the LPDZ in a brittle stage between 10 and 2 Ma with a minimum displacement of ~9 km. Additionally, AFT-data indicate a regional accelerated cooling and exhumation episode starting at ~4 Ma. Thus, tectonic processes can affect the entire orogenic system, while potential feedbacks between erosion and tectonics appear to be limited to the windward sides of an orogenic systems.
The India-Eurasia continental collision zone provides a spectacular example of active mountain building and climatic forcing. In order to quantify the critically important process of mass removal, I analyzed spatial and temporal precipitation patterns of the oscillating monsoon system and their geomorphic imprints. I processed passive microwave satellite data to derive high-resolution rainfall estimates for the last decade and identified an abnormal monsoon year in 2002. During this year, precipitation migrated far into the Sutlej Valley in the northwestern part of the Himalaya and reached regions behind orographic barriers that are normally arid. There, sediment flux, mean basin denudation rates, and channel-forming processes such as erosion by debris-flows increased significantly. Similarly, during the late Pleistocene and early Holocene, solar forcing increased the strength of the Indian summer monsoon for several millennia and presumably lead to analogous precipitation distribution as were observed during 2002. However, the persistent humid conditions in the steep, high-elevation parts of the Sutlej River resulted in deep-seated landsliding. Landslides were exceptionally large, mainly due to two processes that I infer for this time: At the onset of the intensified monsoon at 9.7 ka BP heavy rainfall and high river discharge removed material stored along the river, and lowered the baselevel. Second, enhanced discharge, sediment flux, and increased pore-water pressures along the hillslopes eventually lead to exceptionally large landslides that have not been observed in other periods. The excess sediments that were removed from the upstream parts of the Sutlej Valley were rapidly deposited in the low-gradient sectors of the lower Sutlej River. Timing of downcutting correlates with centennial-long weaker monsoon periods that were characterized by lower rainfall. I explain this relationship by taking sediment flux and rainfall dynamics into account: High sediment flux derived from the upstream parts of the Sutlej River during strong monsoon phases prevents fluvial incision due to oversaturation the fluvial sediment-transport capacity. In contrast, weaker monsoons result in a lower sediment flux that allows incision in the low-elevation parts of the Sutlej River.
Following work is embedded in the multidisciplinary study DESERT (DEad SEa Rift Transect) that has been carried out in the Middle East since the beginning of the year 2000. It focuses on the structure of the southern Dead Sea Transform (DST), the transform plate boundary between Africa (Sinai) and the Arabian microplate. The left-lateral displacement along this major active strike-slip fault amounts to more than 100 km since Miocene times. The DESERT near-vertical seismic reflection (NVR) experiment crossed the DST in the Arava Valley between Red Sea and Dead Sea, where its main fault is called Arava Fault. The 100 km long profile extends in a NW—SE direction from Sede Boqer/Israel to Ma'an/Jordan and coincides with the central part of a wide-angle seismic refraction/reflection line. Near-vertical seismic reflection studies are powerful tools to study the crustal architecture down to the crust/mantle boundary. Although they cannot directly image steeply dipping fault zones, they can give indirect evidence for transform motion by offset reflectors or an abrupt change in reflectivity pattern. Since no seismic reflection profile had crossed the DST before DESERT, important aspects of this transform plate boundary and related crustal structures were not known. Thus this study aimed to resolve the DST's manifestation in both the upper and the lower crust. It was to show, whether the DST penetrates into the mantle and whether it is associated with an offset of the crust/mantle boundary, which is observed at other large strike-slip zones. In this work a short description of the seismic reflection method and the various processing steps is followed by a geological interpretation of the seismic data, taking into account relevant information from other studies. Geological investigations in the area of the NVR profile showed, that the Arava Fault can partly be recognized in the field by small scarps in the Neogene sediments, small pressure ridges or rhomb-shaped grabens. A typical fault zone architecture with a fault gauge, fault-related damage zone, and undeformed host rock, that has been reported from other large fault zones, could not be found. Therefore, as a complementary part to the NVR experiment, which was designed to resolve deeper crustal structures, ASTER (Advanced Spacebourne Thermal Emission and Reflection Radiometer) satellite images were used to analyze surface deformation and determine neotectonic activity.
Paleomagnetic dating of climatic events in Late Quaternary sediments of Lake Baikal (Siberia)
(2004)
Lake Baikal provides an excellent climatic archive for Central Eurasia as global climatic variations are continuously depicted in its sediments. We performed continuous rock magnetic and paleomagnetic analyses on hemipelagic sequences retrieved from 4 underwater highs reaching back 300 ka. The rock magnetic study combined with TEM, XRD, XRF and geochemical analyses evidenced that a magnetite of detrital origin dominates the magnetic signal in glacial sediments whereas interglacial sediments are affected by early diagenesis. HIRM roughly quantifies the hematite and goethite contributions and remains the best proxy for estimating the detrital input in Lake Baikal. Relative paleointensity records of the earth′s magnetic field show a reproducible pattern, which allows for correlation with well-dated reference curves and thus provides an alternative age model for Lake Baikal sediments. Using the paleomagnetic age model we observed that cooling in the Lake Baikal region and cooling of the sea surface water in the North Atlantic, as recorded in planktonic foraminifera δ18 O, are coeval. On the other hand, benthic δ18 O curves record mainly the global ice volume change, which occurs later than the sea surface temperature change. This proves that a dating bias results from an age model based on the correlation of Lake Baikal sedimentary records with benthic δ18 O curves. The compilation of paleomagnetic curves provides a new relative paleointensity curve, “Baikal 200”. With a laser-assisted grain size analysis of the detrital input, three facies types, reflecting different sedimentary dynamics can be distinguished. (1) Glacial periods are characterised by a high clay content mostly due to wind activity and by occurrence of a coarse fraction (sand) transported over the ice by local winds. This fraction gives evidence for aridity in the hinterland. (2) At glacial/interglacial transitions, the quantity of silt increases as the moisture increases, reflecting increased sedimentary dynamics. Wind transport and snow trapping are the dominant process bringing silt to a hemipelagic site (3) During the climatic optimum of the Eemian, the silt size and quantity are minimal due to blanketing of the detrital sources by the vegetal cover.
Gangschwärme nehmen eine bedeutende Stellung im Verständnis zur kontinentalen Fragmentierung ein. Einerseits markieren sie das Paläo-Spannungsfeld und helfen bei der Rekonstruktion der strukturellen Entwicklung der gedehnten Lithosphäre, andererseits gibt ihre petrologische Beschaffenheit Aufschluß über die Entstehung des Magmas, Aufstieg und Platznahme und schließlich erlaubt ihre Altersbestimmung die Rekonstruktion einer chronologischen Reihenfolge magmatischer und struktureller Ereignisse. Das Arbeitsgebiet im namibianischen Henties Bay-Outjo Dike swarm (HOD) war zur Zeit der Unterkreide einem Rifting mit intensiver Platznahme von überwiegend mafischen Gängen unterworfen. Geochemische Signaturen weisen die Gänge als erodierte Förderkanäle der Etendeka Plateaubasalte aus. Durch den Einsatz von hochauflösenden Aeromagnetik- und Satellitendaten war es möglich, die Geometrie des Gangschwarmes erstmals detailliert synoptisch zu erfassen. Viele zu den Schichten des Grundgebirges foliationsparallel verlaufende magnetische Anomalien können unaufgeschlossenen kretazischen Intrusionen zugeordnet werden. Bei der nach Norden propagierenden Südatlantiköffnung spielte die unterschiedliche strukturelle Vorzeichnung durch die neoproterozoischen Faltengürtel sowie Lithologie und Spannungsfeld des Angola Kratons eine bedeutende Rolle. Im küstennahen zentralen Bereich war dank der Vorzeichnung des Nordost streichenden Damara-Faltengürtels ein Rifting in Nordwest-Südost-Richtung dominierend, bis das Angola Kraton ein weiteres Fortscheiten nach Nordosten hemmte und die Vorzeichnung des Nordwest streichenden Kaoko-Faltengürtels an der Westgrenze den weiteren Riftverlauf und die letztendlich erfolgreiche Öffnung des Südatlantiks bestimmte. Aus diesem Grund kann das Gebiet des HOD als ein failed rift betrachtet werden. Die Entwicklung des Spannungsfeldes im HOD kann folgendermaßen skizziert werden: 1. Platznahme von Gängen bei gleichzeitig hoher Dehnungsrate und hohem Magmenfluß. 2. Platznahme von Zentralvulkanen entlang reaktivierter paläozoischer Lineamente bei Abnahme der Dehnungsrate und fortbestehendem hohen Magmenfluß. 3. Abnahme/Versiegen des Magmenflusses und neotektonische Bewegungen führen zur Bildung von Halbgräben.
Western Anatolia that represents the eastward lateral continuation of the Aegean domain is composed of several tectono-metamorphic units showing occurrences of high-pressure/low-temperature (HP-LT) rocks. While some of these metamorphic rocks are vestiges of the Pan-African or Cimmerian orogenies, others are the result of the more recent Alpine orogenesis. In southwest Turkey, the Menderes Massif occupies an extensive area tectonically overlain by nappe units of the Izmir-Ankara Suture Zone in the north, the Afyon Zone in the east, and the Lycian Nappes in the south. In the present study, investigations in the metasediments of the Lycian Nappes and underlying southern Menderes Massif revealed widespread occurrences of Fe-Mg-carpholite-bearing rocks. This discovery leads to the very first consideration that both nappe complexes recorded HP-LT metamorphic conditions during the Alpine orogenesis. P-T conditions for the HP metamorphic peak are about 10-12 kbar/400°C in the Lycian Nappes, and 12-14 kbar/470-500°C in the southern Menderes Massif, documenting a burial of at least 30 km during subduction and nappe stacking. Ductile deformation analysis in concert with multi-equilibrium thermobarometric calculations reveals that metasediments from the Lycian Nappes recorded distinct exhumation patterns after a common HP metamorphic peak. The rocks located far from the contact separating the Lycian Nappes and the Menderes Massif, where HP parageneses are well preserved, retained a single HP cooling path associated with top-to-the-NNE shearing related to the Akçakaya shear zone. This zone of strain localization is an intra-nappe contact that was active in the early stages of exhumation of HP rocks, within the stability field of Fe-Mg-carpholite. The rocks located close to the contact with the Menderes Massif, where HP parageneses are completely retrogressed into chlorite and mica, recorded warmer exhumation paths associated with top-to-the-E intense shearing. This deformation occurred after the southward emplacement of Lycian Nappes, and is contemporaneous with the reactivation of the ’Lycian Nappes-Menderes Massif′ contact as a major shear zone (the Gerit shear zone) that allowed late exhumation of HP parageneses under warmer conditions. The HP rocks from the southern Menderes Massif recorded a simple isothermal decompression at about 450°C during exhumation, and deformation during HP event and its exhumation is characterized by a severe N-S to NE-SW stretching. The age of the HP metamorphism recorded in the Lycian Nappes is assumed to range between the Latest Cretaceous (age of the youngest sediments in the Lycian allochthonous unit) and the Eocene (age of the Cycladic Blueschists). A probable Palaeocene age is suggested. The age of the HP metamorphism that affected the cover series of the Menderes Massif is constrained between the Middle Palaeocene (age of the uppermost metaolistostrome of the Menderes ’cover′) and the Middle Eocene (age of the HP metamorphism in the Dilek-Selçuk region that belongs to the Cycladic Complex). Apatite fission track data for the rocks on both sides of the ’Lycian Nappes/Menderes Massif’ contact suggest that these rocks were very close to the paleo-Earth surface in the Late Oligocene-Early Miocene time. This study in the Lycian Nappes and in the Menderes Massif establishes the existence of an extensive Alpine HP metamorphic belt in southwest Turkey. HP rocks were involved in the accretionary complex related to northward-verging subduction of the Neo-Tethys Ocean, Late Cretaceous obduction and subsequent Early Tertiary continental collision of the passive margin (Anatolide-Tauride block) beneath the active margin of the northern plate (Sakarya micro-continent). During the Eocene, the accretionary complex was made of three stacked HP units. The lowermost corresponds to the imbricated ’core′ and HP ’cover′ of the Menderes Massif, the intermediate one consists of the Cycladic Blueschist Complex (Dilek-Selçuk unit), and the uppermost unit is made of the HP Lycian Nappes. Whereas the basement units of both Aegean and Anatolian regions underwent a different pre-Mesozoic tectonic history, they were probably juxtaposed by the end of the Paleozoic and underwent a common Mesozoic history. Then, the basements and their cover, as well as the Cycladic Blueschists and the Lycian Nappes were involved in similar evolutional accretionary complexes during the Eocene and Oligocene times.
In this work, an approach of paleoclimate reconstruction for tropical East Africa is presented. After giving a short summary of modern climate conditions in the tropics and the East African climate peculiarity, the potential of reconstructing climate from paleolake sediments is discussed. As demonstrated, the hydrologic sensitivity of high-elevated closed-basin lakes in the Central Kenya Rift yields valuable guaranties for the establishment of long-term climate records. Temporal fluctuations of the limnological characteristics saved in the lake sediments are used to define variations in the Quaternary climate history. Based on diatom analyses in radiocarbon- and 40Ar/39Ar-dated sediments, a chronology of paleoecologic fluctuations is developed for the Central Kenya Rift -lakes Nakuru, Elmenteita and Naivasha. At least during the penultimate interglacial (around 140 to 60 kyr BP) and during the last interglacial (around 12 to 4 kyr BP), these lakes experienced several transgression-regression cycles on time intervals of about 11,000 years. Additionally, a long-term trend of lake evolution is found suggesting the general succession from deep freshwater lakes towards more saline waters during the last million years. Using ecologic transfer functions and a simple lake-balance model, the observed paleohydrologic fluctuations are linked to potential precipitation-evaporation changes in the lake basins. Though also tectonic influences on the drainage pattern and the effect of varied seepage are investigated, it can be shown that already a small increase in precipitation of about 30±10 % may have affected the hydrologic budget of the intra-rift lakes within the reconstructed range. The findings of this study help to assess the natural climate variability of East Africa. They furthermore reflect the sensitivity of the Central Kenya Rift -lakes to fluctuations of large-scale climate parameters, such as solar radiation and sea-surface temperatures of the Indian Ocean.
The Dead Sea Transform (DST) is a prominent shear zone in the Middle East. It separates the Arabian plate from the Sinai microplate and stretches from the Red Sea rift in the south via the Dead Sea to the Taurus-Zagros collision zone in the north. Formed in the Miocene about 17 Ma ago and related to the breakup of the Afro-Arabian continent, the DST accommodates the left-lateral movement between the two plates. The study area is located in the Arava Valley between the Dead Sea and the Red Sea, centered across the Arava Fault (AF), which constitutes the major branch of the transform in this region. A set of seismic experiments comprising controlled sources, linear profiles across the fault, and specifically designed receiver arrays reveals the subsurface structure in the vicinity of the AF and of the fault zone itself down to about 3-4 km depth. A tomographically determined seismic P velocity model shows a pronounced velocity contrast near the fault with lower velocities on the western side than east of it. Additionally, S waves from local earthquakes provide an average P-to-S velocity ratio in the study area, and there are indications for a variations across the fault. High-resolution tomographic velocity sections and seismic reflection profiles confirm the surface trace of the AF, and observed features correlate well with fault-related geological observations. Coincident electrical resistivity sections from magnetotelluric measurements across the AF show a conductive layer west of the fault, resistive regions east of it, and a marked contrast near the trace of the AF, which seems to act as an impermeable barrier for fluid flow. The correlation of seismic velocities and electrical resistivities lead to a characterisation of subsurface lithologies from their physical properties. Whereas the western side of the fault is characterised by a layered structure, the eastern side is rather uniform. The vertical boundary between the western and the eastern units seems to be offset to the east of the AF surface trace. A modelling of fault-zone reflected waves indicates that the boundary between low and high velocities is possibly rather sharp but exhibits a rough surface on the length scale a few hundreds of metres. This gives rise to scattering of seismic waves at this boundary. The imaging (migration) method used is based on array beamforming and coherency analysis of P-to-P scattered seismic phases. Careful assessment of the resolution ensures reliable imaging results. The western low velocities correspond to the young sedimentary fill in the Arava Valley, and the high velocities in the east reflect mainly Precambrian igneous rocks. A 7 km long subvertical scattering zone reflector is offset about 1 km east of the AF surface trace and can be imaged from 1 km to about 4 km depth. The reflector marks the boundary between two lithological blocks juxtaposed most probably by displacement along the DST. This interpretation as a lithological boundary is supported by the combined seismic and magnetotelluric analysis. The boundary may be a strand of the AF, which is offset from the current, recently active surface trace. The total slip of the DST may be distributed spatially and in time over these two strands and possibly other faults in the area.
Late Miocene to Quaternary volcanic rocks from the frontal arc to the back-arc region of the Central Volcanic Zone in the Andes show a wide range of delta 11B values (+4 to -7 ‰) and boron concentrations (6 to 60 ppm). Positive delta 11B values of samples from the volcanic front indicate involvement of a 11B-enriched slab component, most likely derived from altered oceanic crust, despite the thick Andean continental lithosphere, and rule out a pure crust-mantle origin for these lavas. The delta 11B values and B concentrations in the lavas decrease systematically with increasing depth of the Wadati-Benioff Zone. This across-arc variation in delta 11B values and decreasing B/Nb ratios from the arc to the back-arc samples are attributed to the combined effects of B-isotope fractionation during progressive dehydration in the slab and a steady decrease in slab-fluid flux towards the back arc, coupled with a relatively constant degree of crustal contamination as indicated by similar Sr, Nd and Pb isotope ratios in all samples. Modelling of fluid-mineral B-isotope fractionation as a function of temperature fits the across-arc variation in delta 11B and we conclude that the B-isotope composition of arc volcanics is dominated by changing delta 11B composition of B-rich slab-fluids during progressive dehydration. Crustal contamination becomes more important towards the back-arc due to the decrease in slab-derived fluid flux. Because of this isotope fractionation effect, high delta 11B signatures in volcanic arcs need not necessarily reflect differences in the initial composition of the subducting slab. Three-component mixing calculations for slab-derived fluid, the mantle wedge and the continental crust based on B, Sr and Nd isotope data indicate that the slab-fluid component dominates the B composition of the fertile mantle and that the primary arc magmas were contaminated by an average addition of 15 to 30 % crustal material.