Refine
Year of publication
Document Type
- Doctoral Thesis (407) (remove)
Language
- English (407) (remove)
Keywords
- Fernerkundung (14)
- Erdbeben (12)
- Klimawandel (12)
- climate change (12)
- remote sensing (12)
- Anden (11)
- Andes (11)
- Tektonik (10)
- Erosion (9)
- Geomorphologie (9)
Institute
- Institut für Geowissenschaften (407) (remove)
Unlike today’s prevailing terrestrial features, the geologic past of Central Asia witnessed marine environments and conditions as well. A vast, shallow sea, known as proto-Paratethys, extended across Eurasia from the Mediterranean Tethys to the Tarim Basin in western China during Cretaceous to Paleogene times. This sea formed about 160 million years ago (during Jurassic times) when the waters of the Tethys Ocean flooded into Eurasia. It drastically retreated to the west and became isolated as the Paratethys during the Late Eocene-Oligocene (ca. 34 Ma).
Having well-constrained timing and paleogeography for the Cretaceous-Paleogene proto-Paratethys sea incursions in Central Asia is essential to properly understand and distinguish the controlling mechanisms and their link to Asian paleoenvironmental and paleoclimatic change. The Cretaceous-Paleogene tectonic evolution of the Pamir and Tibet and their far-field effects play a significant role on the sedimentological and structural evolution of the Central Asian basins and on the evolution of the proto-Paratethys sea fluctuations as well. Comparing the records of the sea incursions to the tectonic and eustatic events has paramount importance to reveal the controlling mechanisms behind the sea incursions. However, due to inaccuracies in the dating of rocks (mostly continental rocks and marine rocks with benthic microfossils providing low-resolution biostratigraphic constraints) and conflicting results, there has been no consensus on the timing of the sea incursions and interpretation of their records has been in question. Here, we present a new chronostratigraphic framework based on biostratigraphy and magnetostratigraphy as well as a detailed paleoenvironmental analysis for the Cretaceous and Paleogene proto-Paratethys Sea incursions in the Tajik and Tarim basins, in Central Asia. This enables us to identify the major drivers of marine fluctuations and their potential consequences on regional and global climate, particularly Asian aridification and the global carbon cycle perturbations such as the Paleocene-Eocene Thermal Maximum (PETM). To estimate the paleogeographic evolution of the proto-Paratethys Sea, the refined age constraints and detailed paleoenvironmental interpretations are combined with successive paleogeographic maps. Regional coastlines and depositional environments during the Cretaceous-Paleogene sea advances and retreats were drawn based on the results of this thesis and integrated with existing literature to generate new paleogeographic maps.
Before its final westward retreat in the Eocene, a total of six Cretaceous and Paleogene major sea incursions have been distinguished from the sedimentary records of the Tajik and Tarim basins in Central Asia. All have been studied and documented here.
We identify the presence of marine conditions already in the Early Cretaceous in the western Tajik Basin, followed by the Cenomanian (ca. 100 Ma) and Santonian (ca. 86 Ma) major marine incursions far into the eastern Tajik and Tarim basins separated by a Turonian-Coniacian (ca. 92-86 Ma) regression. Basin-wide tectonic subsidence analyses imply that the Early Cretaceous invasion of the sea into the Tajik Basin is related to increased Pamir tectonism (at ca. 130 – 90 Ma) in a retro-arc basin setting inferred to be linked to collision and subduction. This tectonic event mainly governed the Cenomanian (ca. 100 Ma) sea incursion in conjunction with a coeval global eustatic high resulting in the maximum geographic extent of the sea. The following Turonian-Coniacian (ca. 92-86 Ma) major regression, driven by eustasy, coincides with a sharp slowdown in tectonic subsidence related to a regime change in Pamir tectonism from compression to extension. The Santonian (ca. 86 Ma) major sea incursion was more likely controlled dominantly by eustasy as also evidenced by the coeval fluctuations in the west Siberian Basin. During the early Maastrichtian, the global Late Cretaceous cooling is inferred from the disappearance of mollusk-rich limestones and the dominance of bryozoan-rich and echinoderm-rich limestones in the Tajik Basin documenting the first evidence for the Late Cretaceous cooling event in Central Asia.
Following the last Cretaceous sea incursion, a major regional restriction event, marked by the exceptionally thick (≤ 400 m) shelf evaporites is assigned a Danian-Selandian age (ca. 63-59 Ma). This is followed by the largest recorded proto-Paratethys sea incursion with a transgression estimated as early Thanetian (ca. 59-57 Ma) and a regression within the Ypresian (ca. 53-52 Ma). The transgression of the next incursion is now constrained as early Lutetian (ca. 47-46 Ma), whereas its regression is constrained as late Lutetian (ca. 41 Ma) and is associated with a drastic increase in both tectonic subsidence and basin infilling. The age of the final and least pronounced sea incursion restricted to the westernmost margin of the Tarim Basin is assigned as Bartonian–Priabonian (ca. 39.7-36.7 Ma). We interpret the long-term westward retreat of the proto-Paratethys Sea starting at ca. 41 Ma to be associated with far-field tectonic effects of the Indo-Asia collision and Pamir/Tibetan plateau uplift. Short-term eustatic sea level transgressions are superimposed on this long-term regression and seem coeval with the transgression events in the other northern Peri-Tethyan sedimentary provinces for the 1st and 2nd Paleogene sea incursions. However, the last Paleogene sea incursion is interpreted as related to tectonism. The transgressive and regressive intervals of the proto-Paratethys Sea correlate well with the reported humid and arid phases, respectively in the Qaidam and Xining basins, thus demonstrating the role of the proto-Paratethys Sea as an important moisture source for the Asian interior and its regression as a contributor to Asian aridification.
We lastly study the mechanics, relative contribution and preservation efficiency of ancient epicontinental seas as carbon sinks with new and existing data, using organic rich (sapropel) deposits dated to the PETM from the extensive epicontinental proto-Paratethys and West Siberian seas. We estimate ca. 1390±230 Gt organic C burial, a substantial amount compared to previously estimated global total excess organic C burial (ca. 1700-2900 Gt) is focused in the proto-Paratethys and West Siberian seas alone. We also speculate that enhanced organic carbon burial later over much of the proto-Paratethys (and later Paratethys) basin (during the deposition of the Kuma Formation and Maikop series, repectively) may have majorly contributed to drawdown of atmospheric carbon dioxide before and during the EOT cooling and glaciation of Antarctica. For past periods with smaller epicontinental seas, the effectiveness of this negative carbon cycle feedback was arguably diminished, and the same likely applies to the present-day.
Carbonates play a key role in the chemistry and dynamics of our planet. They are directly connected to the CO2 budget of our atmosphere and have a great impact on the deep carbon cycle. Moreover, recent studies have shown that carbonates are stable along the geothermal gradient down to Earth's lower mantle conditions, changing their crystal structure and related properties. Subducted carbonates may also react with silicates to form new phases. These reactions will redistribute elements, such as calcium (Ca), magnesium (Mg), iron (Fe) and carbon in the form of carbon dioxide (CO2), but also trace elements, that are carried by the carbonates. The trace elements of most interest are strontium (Sr) and rare earth elements (REE) which have been found to be important constituents in the composition of the primitive lower mantle and in mineral inclusions found in super-deep diamonds. However, the stability of carbonates in presence of mantle silicates at relevant temperatures is far from being well understood. Related to this, very little is known about distribution processes of trace elements between carbonates and mantle silicates. To shed light on these processes, we studied reactions between Sr- and REE-containing CaCO3 and Mg/Fe-bearing silicates of the system (Mg,Fe)2SiO4 - (Mg,Fe)SiO3 at high pressure and high temperature using synchrotron radiation based μ-X-ray diffraction (μ-XRD) and μ-X-ray fluorescence (μ-XRF) with μm-resolution in a laser-heated diamond anvil cell. X-ray diffraction is used to derive the structural changes of the phase reactions whereas X-ray fluorescence gives information on the chemical changes in the sample. In-situ experiments at high pressure and high temperature were performed at beamline P02.2 at PETRA III (Hamburg, Germany) and at beamline ID27 at ESRF (Grenoble, France). In addition to μ-XRD and μ-XRF, ex-situ measurements were made on the recovered sample material using transmission electron microscopy (TEM) and provided further insights into the reaction kinetics of carbonate-silicate reactions.
Our investigations show that CaCO3 is unstable in presence of mantle silicates above 1700 K and a reaction takes place in which magnesite plus CaSiO3-perovskite are formed. In addition, we observed that a high content of iron in the carbonate-silicate system favours dolomite formation during the reaction. The subduction of natural carbonates with significant amounts of Sr leads to a comprehensive investigation of the stability not only of CaCO3 phases in contact with mantle silicates but also of SrCO3 (and of Sr-bearing CaCO3). We found that SrCO3 reacts with (Mg,Fe)SiO3-perovskite to form magnesite and gained evidence for the formation of SrSiO3-perovskite.
To complement our study on the stability of SrCO3 at conditions of the Earth's lower mantle, we performed powder X-ray diffraction and single crystal X-ray diffraction experiments at ambient temperature and up to 49 GPa. We observed a transformation from SrCO3-I into a new high-pressure phase SrCO3-II at around 26 GPa with Pmmn crystal structure and a bulk modulus of 103(10) GPa. This information is essential to fully understand the phase behaviour and stability of carbonates in the Earth's lower mantle and to elucidate the possibility of introducing Sr into mantle silicates by carbonate-silicate reactions.
Simultaneous recording of μ-XRD and μ-XRF in the μm-range over the heated areas provides spatial information not only about phase reactions but also on the elemental redistribution during the reactions. A comparison of the spatial intensity distribution of the XRF signal before and after heating indicates a change in the elemental distribution of Sr and an increase in Sr-concentration was found around the newly formed SrSiO3-perovskite. With the help of additional TEM analyses on the quenched sample material the elemental redistribution was studied at a sub-micrometer scale. Contrary to expectations from combined μ-XRD and μ-XRF measurements, we found that La and Eu were not incorporated into the silicate phases, instead they tend to form either isolated oxide phases (e.g. Eu2O3, La2O3) or hydroxyl-bastnäsite (La(CO3)(OH)). In addition, we observed the transformation from (Mg,Fe)SiO3-perovskite to low-pressure clinoenstatite during pressure release. The monoclinic structure (P21/c) of this phase allows the incorporation of Ca as shown by additional EDX analyses and, to a minor extent, Sr too.
Based on our experiments, we can conclude that a detection of the trace elements in-situ at high pressure and high temperature remains challenging. However, our first findings imply that silicates may incorporate the trace elements provided by the carbonates and indicate that carbonates may have a major effect on the trace element contents of mantle phases.
The transfer of particulate organic carbon from continents to the ocean is an important component of the global carbon cycle. Transfer to and burial of photosynthetically fixed biospheric organic carbon in marine sediments can effectively sequester atmospheric carbon dioxide over geological timescales. The exhumation and erosion of fossil organic carbon contained in sedimentary rocks, i.e. petrogenic carbon, can result in remineralization, releasing carbon to the atmosphere. In contrast, eroded petrogenic organic carbon that gets transferred back to the ocean and reburied does not affect atmospheric carbon content.
Mountain ranges play a key role in this transfer since they can source vast amounts of sediment including particulate organic carbon. Globally, the export of both, biospheric and petrogenic organic carbon has been linked to sediment export. Additionally, short transfer times from mountains to the ocean and high sediment concentrations have been shown to increase the likelihood of organic carbon burial. While the importance of mountain ranges in the organic carbon cycle is now widely recognized, the processes acting within mountain ranges to influence the storage, cycling and mobilization of organic carbon, as well as carbon fluxes from mountain ranges remain poorly constrained.
In this thesis, I employ different methods to assess the nature and fate of particulate organic carbon in mountain belts, ranging from the molecular to regional landscape scale. These studies are located along the Trans-Himalayan Kali Gandaki River in Central Nepal. This river traverses all major geological and climatic zones of the Himalaya, from the dry northern Tibetan plateau to the high-relief, monsoon dominated steep High Himalaya and the lower relief and abundant vegetation of the Lesser Himalayan region.
First, I document how biospheric organic matter has accumulated during the Holocene in the headwaters of the Kali Gandaki River valley, by combining compound specific isotope measurements with different dating methods and grain size data, and investigate the stability of this organic carbon reservoir on millennial timescales. I show, that around 1.6 ka an eco-geomorphic tipping point occurred leading to a destabilization of the landscape resulting in today’s high erosion rates and the excavation of the aged organic carbon reservoir. This study highlights the climatic and geomorphic controls on biospheric organic carbon storage and release from mountain ranges.
Second, I systematically investigate the spatial variation of particulate organic carbon fluxes across the Himalaya along the Kali Gandaki River, using bulk stable and radioactive isotopes combined with a new Bayesian modeling approach. The detailed dataset allows the distinction of aged and modern biospheric organic carbon as well as petrogenic organic carbon across the Himalayan mountain range and the investigation of the role of climatic and geomorphic factors in their riverine export. The data suggest a decoupling of the particulate organic carbon from the sediment yield along the Kali Gandaki River, partially driven by climatic and geomorphic processes. In contrast to the suspended sediment, a large part of the particulate organic carbon exported by the river originates from the Tibetan part of the catchment and is dominated by petrogenic organic carbon derived from Jurassic shales with only minor contributions of modern and aged biospheric organic carbon. These findings emphasize the importance of organic carbon source distribution and erosion mechanisms in determining the organic carbon export from mountain ranges.
In a third step, I explore the potential of ultra-high resolution mass spectrometry for particulate organic carbon transport studies. I have generated a novel and unprecedented high-resolution molecular dataset, which contains up to 103 molecular formulas of the lipid fraction of particulate organic matter for modern and aged biospheric carbon, petrogenic organic carbon and river sediments. First, I test if this dataset can be used to better resolve different organic carbon sources and to identify new geochemical tracers. Using multivariate statistics, I identify up to 10² characteristic molecular formulas for the major organic carbon sources in the upper part of the Kali Gandaki catchment, and trace their transfer from the surrounding landscape into the river sediment. Second, I test the potential of the molecular dataset to trace molecular transformations along source-to-sink pathways. I identify changes in molecular metrics derived from the dataset, which are characteristic of transformation processes during incorporation of litter into soil, the aging of soil material, and the mobilization of the organic carbon into the river. These two studies demonstrate that high-resolution molecular datasets open a promising analytical window on particulate organic carbon and can provide novel insights into the composition, sourcing and transformation of riverine particulate organic carbon.
Collectively, these studies advance our understanding of the processes contributing to the storage and mobilization of organic carbon in the Central Himalaya, the mountain belt that dominates global erosional fluxes. They do so by identifying the major sources of particulate organic carbon to the Trans-Himalayan Kali Gandaki River, by elucidating their sensitivity to climate and geomorphic processes, and by identifying some of the transformations of this material on the molecular scale. As a result, the thesis demonstrates that the amount and composition of organic carbon routed from mountain belts is a function of the dynamic interactions of geologic, biologic, geomorphic and climatic processes within the mountain belt. This understanding will ultimately help in answering whether the build-up and erosion of mountain ranges over geological time represents a net carbon source or sink to the atmosphere. Beyond this, the thesis contributes to our technical ability to characterize organic matter and attribute it to sources by scoping the potential of high-end molecular analysis.
Ferruginous conditions were a prominent feature of the oceans throughout the Precambrian Eons and thus throughout much of Earth’s history. Organic matter mineralization and diagenesis within the ferruginous sediments that deposited from Earth’s early oceans likely played a key role in global biogeochemical cycling. Knowledge of organic matter mineralization in ferruginous sediments, however, remains almost entirely conceptual, as modern analogue environments are extremely rare and largely unstudied, to date. Lake Towuti on the island of Sulawesi, Indonesia is such an analogue environment and the purpose of this PhD project was to investigate the rates and pathways of organic matter mineralization in its ferruginous sediments.
Lake Towuti is the largest tectonic lake in Southeast Asia and is hosted in the mafic and ultramafic rocks of the East Sulawesi Ophiolite. It has a maximum water depth of 203 m and is weakly thermally stratified. A well-oygenated surface layer extends to 70 m depth, while waters below 130 m are persistently anoxic. Intensive weathering of the ultramafic catchment feeds the lake with large amounts of iron(oxy)hydroxides while the runoff contains only little sulfate, leading to sulfate-poor (< 20 µM) lake water and anoxic ferruginous conditions below 130 m. Such conditions are analogous to the ferruginous water columns that persisted throughout much of the Archean and Proterozoic eons. Short (< 35 cm) sediment cores were collected from different water depths corresponding to different bottom water redox conditions. Also, a drilling campaign of the International Continental Scientific Drilling Program (ICDP) retrieved a 114 m long sediment core dedicated for geomicrobiological investigations from a water depth of 153 m, well below the depth of oxygen penetration at the time of sampling. Samples collected from these sediment cores form the fundament of this thesis and were used to perform a suite of biogeochemical and microbiological analyses.
Geomirobiological investigations depend on uncontaminated samples. However, exploration of subsurface environments relies on drilling, which requires the use of a drilling fluid. Drilling fluid infiltration during drilling can not be avoided. Thus, in order to trace contamination of the sediment core and to identify uncontaminated samples for further analyses a simple and inexpensive technique for assessing contamination during drilling operations was developed and applied during the ICDP drilling campaign. This approach uses an aqeous fluorescent pigment dispersion commonly used in the paint industry as a particulate tracer. It has the same physical properties as conventionally used particulate tracers. However, the price is nearly four orders of magnitude lower solving the main problem of particulate tracer approaches. The approach requires only a minimum of equipment and allows for a rapid contamination assessment potentially even directly on site, while the senstitivity is in the range of already established approaches. Contaminated samples in the drill core were identified and not included for further geomicrobiological investigations.
Biogeochemical analyses of short sediment cores showed that Lake Towutis sediments are strongly depleted in electron acceptors commonly used in microbial organic matter mineralization (i.e. oxygen, nitrate, sulfate). Still, the sediments harbor high microbial cell densities, which are a function of redox conditions of Lake Towuti’s bottom water. In shallow water depths bottom water oxygenation leads to a higher input of labile organic matter and electron acceptors like sulfate and iron, which promotes a higher microbial abundance. Microbial analyses showed that a versatile microbial community with a potential to perform metabolisms related to iron and sulfate reduction, fermentation as well as methanogenesis inhabits Lake Towuti’s surface sediments.
Biogeochemical investigations of the upper 12 m of the 114 m sediment core showed that Lake Towuti’s sediment is extremely rich in iron with total concentrations up to 2500 µmol cm-3 (20 wt. %), which makes it the natural sedimentary environment with the highest total iron concentrations studied to date. In the complete or near absence of oxygen, nitrate and sulfate, organic matter mineralization in ferruginous sediments would be expected to proceed anaerobically via the energetically most favorable terminal electron acceptors available - in this case ferric iron. Astonishingly, however, methanogenesis is the dominant (>85 %) organic matter mineralization process in Lake Towuti’s sediment. Reactive ferric iron known to be available for microbial iron reduction is highly abundant throughout the upper 12 m and thus remained stable for at least 60.000 years. The produced methane is not oxidized anaerobically and diffuses out of the sediment into the water column. The proclivity towards methanogenesis, in these very iron-rich modern sediments, implies that methanogenesis may have played a more important role in organic matter mineralization thoughout the Precambrian than previously thought and thus could have been a key contributor to Earth’s early climate dynamics.
Over the whole sequence of the 114 m long sediment core siderites were identified and characterized using high-resolution microscopic and spectroscopic imaging together with microchemical and geochemical analyses. The data show early diagenetic growth of siderite crystals as a response to sedimentary organic matter mineralization. Microchemical zoning was identified in all siderite crystals. Siderite thus likely forms during diagenesis through growth on primary existing phases and the mineralogical and chemical features of these siderites are a function of changes in redox conditions of the pore water and sediment over time. Identification of microchemical zoning in ancient siderites deposited in the Precambrian may thus also be used to infer siderite growth histories in ancient sedimentary rocks including sedimentary iron formations.
Salt pans also termed playas are common landscape features of hydrologically closed basins in arid and semiarid zones, where evaporation significantly exceeds the local precipitation. The analysis and monitoring of salt pan environments is important for the evaluation of current and future impact of these landscape features. Locally, salt pans have importance for the ecosystem, wildlife and human health, and through dust emissions they influence the climate on regional and global scales. Increasing economic exploitation of these environments in the last years, e.g. by brine extraction for raw materials, as well as climate change severely affect the water, material and energy balance of these systems. Optical remote sensing has the potential to characterise salt pan environments and to increase the understanding of processes in playa basins, as well as to assess wider impacts and feedbacks that exist between climate forcing and human intervention in their regions. Remote sensing techniques can provide information for extensive regions on a high temporal basis compared to traditional field samples and ground observations. Specifically, for salt pans that are often challenging to study because of their large size, remote location, and limited accessibility due to missing infrastructure and ephemeral flooding. Furthermore, the availability of current and upcoming hyperspectral remote sensing data opened the opportunity for the analyses of the complex reflectance signatures that relate to the mineralogical mixtures found in the salt pan sediments. However, these new advances in sensor technology, as well as increased data availability currently have not been fully explored for the study of salt pan environments. The potential of new sensors needs to be assessed and state of the art methods need to be adapted and improved to provide reliable information for in depth analysis of processes and characterisation of the recent condition, as well as to support long-term monitoring and to evaluate environmental impacts of changing climate and anthropogenic activity.
This thesis provides an assessment of the capabilities of optical remote sensing for the study of salt pan environments that combines the information of hyperspectral data with the increased temporal coverage of multispectral observations for a more complete understanding of spatial and temporal complexity of salt pan environments using the Omongwa salt pan located in the south-west Kalahari as a test site. In particular, hyperspectral data are used for unmixing of the mineralogical surface composition, spectral feature-based modelling for quantification of main crust components, as well as time-series based classification of multispectral data for the assessment of the long-term dynamic and the analysis of the seasonal process regime. The results show that the surface of the Omongwa pan can be categorized into three major crust types based on diagnostic absorption features and mineralogical ground truth data. The mineralogical crust types can be related to different zones of surface dynamic as well as pan morphology that influences brine flow during the pan inundation and desiccation cycles. Using current hyperspectral imagery, as well as simulated data of upcoming sensors, robust quantification of the gypsum component could be derived. For the test site the results further indicate that the crust dynamic is mainly driven by flooding events in the wet season, but it is also influenced by temperature and aeolian activity in the dry season. Overall, the scientific outcomes show that optical remote sensing can provide a wide range of information helpful for the study of salt pan environments. The thesis also highlights that remote sensing approaches are most relevant, when they are adapted to the specific site conditions and research scenario and that upcoming sensors will increase the potential for mineralogical, sedimentological and geomorphological analysis, and will improve the monitoring capabilities with increased data availability.
To investigate the reliability and stability of spherical harmonic models based on archeo/-paleomagnetic data, 2000 Geomagnetic models were calculated. All models are based on the same data set but with randomized uncertainties. Comparison of these models to the geomagnetic field model gufm1 showed that large scale magnetic field structures up to spherical harmonic degree 4 are stable throughout all models. Through a ranking of all models by comparing the dipole coefficients to gufm1 more realistic uncertainty estimates were derived than the authors of the data provide.
The derived uncertainty estimates were used in further modelling, which combines archeo/-paleomagnetic and historical data. The huge difference in data count, accuracy and coverage of these two very different data sources made it necessary to introduce a time dependent spatial damping, which was constructed to constrain the spatial complexity of the model. Finally 501 models were calculated by considering that each data point is a Gaussian random variable, whose mean is the original value and whose standard deviation is its uncertainty. The final model arhimag1k is calculated by taking the mean of the 501 sets of Gauss coefficients. arhimag1k fits different dependent and independent data sets well. It shows an early reverse flux patch at the core-mantle boundary between 1000 AD and 1200 AD at the location of the South Atlantic Anomaly today. Another interesting feature is a high latitude flux patch over Greenland between 1200 and 1400 AD. The dipole moment shows a constant behaviour between 1600 and 1840 AD.
In the second part of the thesis 4 new paleointensities from 4 different flows of the island Fogo, which is part of Cape Verde, are presented. The data is fitted well by arhimag1k with the exception of the value at 1663 of 28.3 microtesla, which is approximately 10 microtesla lower than the model suggest.
Orogenic peridotites represent portions of upper subcontinental mantle now incorporated in mountain belts. They often contain layers, lenses and irregular bodies of pyroxenite and eclogite. The origin of this heterogeneity and the nature of these layers is still debated but it is likely to involve processes such as transient melts coming from the crust or the mantle and segregating in magma conduits, crust-mantle interaction, upwelling of the asthenosphere and metasomatism. All these processes occur in the lithospheric mantle and are often related with the subduction of crustal rocks to mantle depths. In fact, during subduction, fluids and melts are released from the slab and can interact with the overlying mantle, making the study of deep melts in this environment crucial to understand mantle heterogeneity and crust-mantle interaction. The aim of this thesis is precisely to better constrain how such processes take place studying directly the melt trapped as primary inclusions in pyroxenites and eclogites. The Bohemian Massif, crystalline core of the Variscan belt, is targeted for these purposes because it contains orogenic peridotites with layers of pyroxenite and eclogite and other mafic rocks enclosed in felsic high pressure and ultra-high pressure crustal rocks. Within this Massif mafic rocks from two areas have been selected: the garnet clinopyroxenite in orogenic peridotite of the Granulitgebirge and the ultra-high pressure eclogite in the diamond-bearing gneisses of the Erzgebirge. In both areas primary melt inclusions were recognized in the garnet, ranging in size between 2-25 µm and with different degrees of crystallization, from glassy to polycrystalline. They have been investigated with Micro Raman spectroscopy and EDS mapping and the mineral assemblage is kumdykolite, phlogopite, quartz, kokchetavite, phase with a main Raman peak at 430 cm-1, phase with a main Raman peak at 412 cm-1, white mica and calcite with some variability in relative abundance depending on the case study. In the Granulitgebirge osumilite and pyroxene are also present, whereas calcite is one of the main phases in the Erzgebirge. The presence of glass and the mineral assemblage in the nanogranitoids suggest that they were former droplets of melt trapped in the garnet while it was growing. Glassy inclusions and re-homogenized nanogranitoids show a silicate melt that is granitic, hydrous, high in alkalis and weakly peraluminous. The melt is also enriched in both case studies in Cs, Pb, Rb, U, Th, Li and B suggesting the involvement of crustal component, i.e. white mica (main carrier of Cs, Pb, Rb, Li and B), and a fluid (Cs, Th and U) in the melt producing reaction. The whole rock in both cases mainly consists of garnet and clinopyroxene with, in Erzgebirge samples, the additional presence of quartz both in the matrix and as a polycrystalline inclusion in the garnet. The latter is interpreted as a quartz pseudomorph after coesite and occurs in the same microstructural position as the melt inclusions. Both rock types show a crustal and subduction zone signature with garnet and clinopyroxene in equilibrium. Melt was likely present during the metamorphic peak of the rock, as it occurs in garnet.
Our data suggest that the processes most likely responsible for the formation of the investigated rocks in both areas is a metasomatic reaction between a melt produced in the crust and mafic layers formerly located in the mantle wedge for the Granulitgebirge and in the subducted continental crust itself in the Erzgebirge. Thus metasomatism in the first case took place in the mantle overlying the slab, whereas in the second case metasomatism took place in the continental crust that already contained, before subduction, mafic layers. Moreover, the presence of former coesite in the same microstructural position of the melt inclusions in the Erzgebirge garnets suggest that metasomatism took place at ultra-high pressure conditions.
Summarizing, in this thesis we provide new insights into the geodynamic evolution of the Bohemian Massif based on the study of melt inclusions in garnet in two different mafic rock types, combining the direct microstructural and geochemical investigation of the inclusions with the whole-rock and mineral geochemistry. We report for the first time data, directly extracted from natural rocks, on the metasomatic melt responsible for the metasomatism of several areas of the Bohemian Massif. Besides the two locations here investigated, belonging to the Saxothuringian Zone, a signature similar to the investigated melt is clearly visible in pyroxenite and peridotite of the T-7 borehole (again Saxothuringian Zone) and the durbachite suite located in the Moldanubian Zone.
Seismological and seismotectonic analysis of the northwestern Argentine Central Andean foreland
(2020)
After a severe M W 5.7 earthquake on October 17, 2015 in El Galpón in the province of Salta NW Argentina, I installed a local seismological network around the estimated epicenter. The network covered an area characterized by inherited Cretaceous normal faults and neotectonic faults with unknown recurrence intervals, some of which may have been reactivated normal faults. The 13 three-component seismic stations recorded data continuously for 15 months.
The 2015 earthquake took place in the Santa Bárbara System of the Andean foreland, at about 17km depth. This region is the easternmost morphostructural region of the central Andes. As a part of the broken foreland, it is bounded to the north by the Subandes fold-and-thrust belt and the Sierras Pampeanas to the south; to the east lies the Chaco-Paraná basin.
A multi-stage morphotectonic evolution with thick-skinned basement uplift and coeval thin-skinned deformation in the intermontane basins is suggested for the study area. The release of stresses associated with the foreland deformation can result in strong earthquakes, as the study area is known for recurrent and historical, destructive earthquakes. The available continuous record reaches back in time, when the strongest event in 1692 (magnitude 7 or intensity IX) destroyed the city of Esteco. Destructive earthquakes and surface deformation are thus a hallmark of this part of the Andean foreland.
With state-of-the-art Python packages (e.g. pyrocko, ObsPy), a semi-automatic approach is followed to analyze the collected continuous data of the seismological network. The resulting 1435 hypocenter locations consist of three different groups: 1.) local crustal earthquakes (nearly half of the events belong to this group), 2.) interplate activity, of regional distance in the slab of the Nazca-plate, and 3.) very deep earthquakes at about 600km depth. My major interest focused on the first event class. Those crustal events are partly aftershock events of the El Galpón earthquake and a second earthquake, in the south of the same fault. Further events can be considered as background seismicity of other faults within the study area. Strikingly, the seismogenic zone encompass the whole crust and propagates brittle deformation down, close to the Moho.
From the collected seismological data, a local seismic velocity model is estimated, using VELEST. After the execution of various stability tests, the robust minimum 1D-velocity model implies guiding values for the composition of the local, subsurface structure of the crust. Afterwards, performing a hypocenter relocation enables the assignment of individual earthquakes to aftershock clusters or extended seismotectonic structures. This allows the mapping of previously unknown seismogenic faults.
Finally, focal mechanisms are modeled for events with acurately located hypocenters, using the newly derived local velocity model. A compressive regime is attested by the majority of focal mechanisms, while the strike direction of the individual seismogenic structures is in agreement with the overall north – south orientation of the Central Andes, its mountain front, and individual mountain ranges in the southern Santa-Bárbara-System.
The current thesis contains the results from two experimental and one modelling study focused on the topic of ductile strain localization in the presence of material heterogeneities. Localization of strain in the high temperature regime is a well known feature of rock deformation occurring in nature at different scales and in a variety of lithologies. Large scale shear zones at the roots of major crustal fault zones are considered responsible for the activity of plate tectonics on our planet. A large number of mechanisms are suggested to be associated with strain softening and nucleation of localization. Among these, the presence of material heterogeneities within homogeneous host rocks is frequently observed in field examples to trigger shear zone development. Despite a number of studies conducted on the topic, the mechanisms controlling initiation and evolution of localization are not fully understood yet. We investigated, experimentally and by means of numerical modelling, phenomenological and microphysical aspects of high temperature strain localization in a homogeneous body containing single and paired inclusions of weaker material. A monomineralic carbonate system composed of Carrara marble (homogeneous, strong matrix) and Solnhofen limestone (weak planar inclusions) is selected for our studies based on its versatility as an experimental material and on the frequent occurrence of carbonate rocks at the core of natural shear zones.
To explore the influence of different loading conditions on heterogeneity-induced high temperature shear zones we conducted torsion experiments under constant twist (deformation) rate and constant torque (stress) conditions in a Paterson-type deformation apparatus on hollow cylinders of marble containing single planar inclusions of limestone. At the imposed experimental conditions (900 ◦C temperature and 400 MPa confining pressure) both materials deform plastically and the marble is ≈ 9 times stronger than the limestone. The viscosity contrast between the two materials induces a perturbation of the stress field within the marble matrix at the tip of the planar inclusion. Early on along the deformation path (at bulk shear strains ≈ 0.3), heterogeneous distribution of strain can be observed under both loading conditions and a small area of incipient strain localization is formed at the tip of the weak limestone inclusion. Strongly deformed grains, incipient dynamic recrystallization and a weak crystallographic preferred orientation characterize the marble within an area a few mm in front of the inclusion. As the bulk strain is increased (up to γ ≈ 1), the area of microstructural modification is expanded along the inclusion plane, the texture strengthens and grain size refinement by dynamic recrystallization becomes pervasive. Locally, evidences for coexisting brittle deformation are also observed regardless of the imposed loading conditions. A shear zone is effectively formed within the deforming Carrara marble, its geometry controlled by the plane containing the thin plate of limestone. Thorough microstructural and textural analysis, however, do not reveal substantial differences in the mechanisms or magnitude of strain localization at the different loading conditions. We conclude that, in the presence of material heterogeneities capable of inducing strain softening, the imposed loading conditions do not affect ductile localization in its nucleating and transient stages.
As the ultimate goal of experimental rock deformation is the extrapolation of results to geologically relevant time and space scales, we developed 2D numerical models reproducing (and benchmarked to) our experimental results. Our cm-scaled models have been implemented with a first-order strain-dependent softening law to reproduce the effect of rheological weakening in the deforming material. We successfully reproduced the local stress concentration at the inclusion tips and the strain localization initiated in the marble matrix. The heterogeneous distribution of strain and its evolution with imposed bulk deformation (i.e. the shape and extent of the nucleating shear zone) are observed to depend on the degree of softening imposed to the deforming matrix. When a second (artificial) softening step is introduced at elevated bulk strains in the model, the formation of a secondary high strain layer is observed at the core of the initial shear zone, analogous to the development of ultramylonite bands in high strain natural shear zones. Our results do not only reproduce the nucleation and transient evolution of a heterogeneity-induced high temperature shear zone with high accuracy, but also confirm the importance of introducing reliable softening laws capable of mimicking strain weakening to numerical models of crustal scale ductile processes.
Material heterogeneities inducing strain localization in the field are often consisting of brittle precursors (joints and fractures). More generally, the interaction of brittle and ductile deformation mechanisms and its effect on the localization of strain have been a key topic in the structural geology community for a long time. The positive feedback between (micro)fracturing and ductile strain localization is a well recognized effect in a number of field examples. We experimentally investigated the influence of brittle deformation on the initiation and evolution of high temperature shear zones in a strong matrix containing pairs of weak material heterogeneities. Our Carrara marble-Solnhofen limestone inclusions system was tested in triaxial compression under constant strain rate and high temperature (900 ◦C) conditions in a Paterson deformation apparatus. The inclusion pairs were arranged in non-overlapping step-over geometries of either compressional or extensional nature. Experimental runs were conducted at different confining pressures (30, 50, 100 and 300 MPa) to induce various amounts of brittle deformation within the marble matrix. At low confinement (30 and 50 MPa) abundant brittle deformation is observed in all configurations, but the spatial distribution of cracks is dependent on the kinematics of the step-over region: concentrated along the shearing plane between the inclusions in the extensional samples, or broadly distributed around the inclusions but outside the step-over region in the compressional configuration. Accordingly, brittle-assisted ductile processes tend to localize deformation along the inclusions plane in the extensional geometry or to distribute widely across large areas of the matrix in the compressional step-over. At pressures of 100 and 300 MPa fracturing is mostly suppressed in both configurations and strain is accommodated almost entirely by viscous creep. In extensional samples this leads to progressive de-localization with increasing confinement. Our results show that, while ductile localization of strain is indeed more efficient where assisted by brittle processes, these latter are only effective if themselves heterogeneously distributed, ultimately a function of the local stress perturbations.
The significant environmental and socioeconomic consequences of hydrometeorological extreme events, such as extreme rainfall, are constituted as a most important motivation for analyzing these events in the south-central Andes of NW Argentina. The steep topographic and climatic gradients and their interactions frequently lead to the formation of deep convective storms and consequently trigger extreme rainfall generation.
In this dissertation, I focus on identifying the dominant climatic variables and atmospheric conditions and their spatiotemporal variability leading to deep convection and extreme rainfall in the south-central Andes.
This dissertation first examines the significant contribution of temperature on atmospheric humidity (dew-point temperature, Td) and on convection (convective available potential energy, CAPE) for deep convective storms and hence, extreme rainfall along the topographic and climatic gradients. It was found that both climatic variables play an important role in extreme rainfall generation. However, their contributions differ depending on topographic and climatic sub-regions, as well as rainfall percentiles.
Second, this dissertation explores if (near real-time) the measurements conducted by the Global Navigation Satellite System (GNSS) on integrated water vapor (IWV) provide reliable data for explaining atmospheric humidity. I argue that GNSS-IWV, in conjunction with other atmospheric stability parameters such as CAPE, is able to decipher the extreme rainfall in the eastern central Andes. In my work, I rely on a multivariable regression analysis described by a theoretical relationship and fitting function analysis.
Third, this dissertation identifies the local impact of convection on extreme rainfall in the eastern Andes. Relying on a Principal Component Analysis (PCA) it was found that during the existence of moist and warm air, extreme rainfall is observed more often during local night hours. The analysis includes the mechanisms for this observation.
Exploring the atmospheric conditions and climatic variables leading to extreme rainfall is one of the main findings of this dissertation. The conditions and variables are a prerequisite for understanding the dynamics of extreme rainfall and predicting these events in the eastern Andes.