The 11 July 1889 Chilik earthquake (M-w 8.0-8.3) forms part of a remarkable sequence of large earthquakes in the late nineteenth and early twentieth centuries in the northern Tien Shan. Despite its importance, the source of the 1889 earthquake remains unknown, though the macroseismic epicenter is sited in the Chilik valley, similar to 100 km southeast of Almaty, Kazakhstan (similar to 2 million population). Several short fault segments that have been inferred to have ruptured in 1889 are too short on their own to account for the estimated magnitude. In this paper we perform detailed surveying and trenching of the similar to 30 km long Saty fault, one of the previously inferred sources, and find that it was formed in a single earthquake within the last 700 years, involving surface slip of up to 10 m. The scarp-forming event, likely to be the 1889 earthquake, was the only surface-rupturing event for at least 5000 years and potentially for much longer. From satellite imagery we extend the mapped length of fresh scarps within the 1889 epicentral zone to a total of similar to 175 km, which we also suggest as candidate ruptures from the 1889 earthquake. The 175 km of rupture involves conjugate oblique left-lateral and right-lateral slip on three separate faults, with step overs of several kilometers between them. All three faults were essentially invisible in the Holocene geomorphology prior to the last slip. The recurrence interval between large earthquakes on any of these faults, and presumably on other faults of the Tien Shan, may be longer than the timescale over which the landscape is reset, providing a challenge for delineating sources of future hazard.

This case study evaluates the suitability of radar-based quantitative precipitation estimates (QPEs) for the simulation of streamflow in the Marikina River Basin (MRB), the Philippines. Hourly radar-based QPEs were produced from reflectivity that had been observed by an S-band radar located about 90 km from the MRB. Radar data processing and precipitation estimation were carried out using the open source library wradlib. To assess the added value of the radar-based QPE, we used spatially interpolated rain gauge observations (gauge-only (GO) product) as a benchmark. Rain gauge observations were also used to quantify rainfall estimation errors at the point scale. At the point scale, the radar-based QPE outperformed the GO product in 2012, while for 2013, the performance was similar. For both periods, estimation errors substantially increased from daily to the hourly accumulation intervals. Despite this fact, both rainfall estimation methods allowed for a good representation of observed streamflow when used to force a hydrological simulation model of the MRB. Furthermore, the results of the hydrological simulation were consistent with rainfall verification at the point scale: the radar-based QPE performed better than the GO product in 2012, and equivalently in 2013. Altogether, we could demonstrate that, in terms of streamflow simulation, the radar-based QPE can perform as good as or even better than the GO product - even for a basin such as the MRB which has a comparatively dense rain gauge network. This suggests good prospects for using radar-based QPE to simulate and forecast streamflow in other parts of the Philippines where rain gauge networks are not as dense.

Six N-alkylpyridinium salts [CnPy](2)[MCl4] (n = 4 or 12 and M = Co, Cu, Zn) were synthesized, and their structure and thermal properties were studied. The [C4Py](2)[MCl4] compounds are monoclinic and crystallize in the space group P2(1)/n. The crystals of the longer chain analogues [C12Py](2)[MCl4] are triclinic and crystallize in the space group P (1) over bar. Above the melting temperature, all compounds are ionic liquids (ILs). The derivatives with the longer C12 chain exhibit liquid crystallinity and the shorter chain compounds only show a melting transition. Consistent with single-crystal analysis, electron paramagnetic resonance spectroscopy suggests that the [CuCl4](2-) ions in the Cu-based ILs have a distorted tetrahedral geometry.

Metabolically active microbial communities are present in a wide range of subsurface environments. Techniques like enumeration of microbial cells, activity measurements with radiotracer assays and the analysis of porewater constituents are currently being used to explore the subsurface biosphere, alongside with molecular biological analyses. However, many of these techniques reach their detection limits due to low microbial activity and abundance. Direct measurements of microbial turnover not just face issues of insufficient sensitivity, they only provide information about a single specific process but in sediments many different process can occur simultaneously. Therefore, the development of a new technique to measure total microbial activity would be a major improvement. A new tritium-based hydrogenase-enzyme assay appeared to be a promising tool to quantify total living biomass, even in low activity subsurface environments. In this PhD project total microbial biomass and microbial activity was quantified in different subsurface sediments using established techniques (cell enumeration and pore water geochemistry) as well as a new tritium-based hydrogenase enzyme assay. By using a large database of our own cell enumeration data from equatorial Pacific and north Pacific sediments and published data it was shown that the global geographic distribution of subseafloor sedimentary microbes varies between sites by 5 to 6 orders of magnitude and correlates with the sedimentation rate and distance from land. Based on these correlations, global subseafloor biomass was estimated to be 4.1 petagram-C and ~0.6 % of Earth's total living biomass, which is significantly lower than previous estimates. Despite the massive reduction in biomass the subseafloor biosphere is still an important player in global biogeochemical cycles. To understand the relationship between microbial activity, abundance and organic matter flux into the sediment an expedition to the equatorial Pacific upwelling area and the north Pacific Gyre was carried out. Oxygen respiration rates in subseafloor sediments from the north Pacific Gyre, which are deposited at sedimentation rates of 1 mm per 1000 years, showed that microbial communities could survive for millions of years without fresh supply of organic carbon. Contrary to the north Pacific Gyre oxygen was completely depleted within the upper few millimeters to centimeters in sediments of the equatorial upwelling region due to a higher supply of organic matter and higher metabolic activity. So occurrence and variability of electron acceptors over depth and sites make the subsurface a complex environment for the quantification of total microbial activity. Recent studies showed that electron acceptor processes, which were previously thought to thermodynamically exclude each other can occur simultaneously. So in many cases a simple measure of the total microbial activity would be a better and more robust solution than assays for several specific processes, for example sulfate reduction rates or methanogenesis. Enzyme or molecular assays provide a more general approach as they target key metabolic compounds. Since hydrogenase enzymes are ubiquitous in microbes, the recently developed tritium-based hydrogenase radiotracer assay is applied to quantify hydrogenase enzyme activity as a parameter of total living cell activity. Hydrogenase enzyme activity was measured in sediments from different locations (Lake Van, Barents Sea, Equatorial Pacific and Gulf of Mexico). In sediment samples that contained nitrate, we found the lowest cell specific enzyme activity around 10^(-5) nmol H_(2) cell^(-1) d^(-1). With decreasing energy yield of the electron acceptor used, cell-specific hydrogenase activity increased and maximum values of up to 1 nmol H_(2) cell^(-1) d^(-1) were found in samples with methane concentrations of >10 ppm. Although hydrogenase activity cannot be converted directly into a turnover rate of a specific process, cell-specific activity factors can be used to identify specific metabolism and to quantify the metabolically active microbial population. In another study on sediments from the Nankai Trough microbial abundance and hydrogenase activity data show that both the habitat and the activity of subseafloor sedimentary microbial communities have been impacted by seismic activities. An increase in hydrogenase activity near the fault zone revealed that the microbial community was supplied with hydrogen as an energy source and that the microbes were specialized to hydrogen metabolism.

Subsurface microbial communities undertake many terminal electron-accepting processes, often simultaneously. Using a tritium-based assay, we measured the potential hydrogen oxidation catalyzed by hydrogenase enzymes in several subsurface sedimentary environments (Lake Van, Barents Sea, Equatorial Pacific, and Gulf of Mexico) with different predominant electron-acceptors. Hydrogenases constitute a diverse family of enzymes expressed by microorganisms that utilize molecular hydrogen as a metabolic substrate, product, or intermediate. The assay reveals the potential for utilizing molecular hydrogen and allows qualitative detection of microbial activity irrespective of the predominant electron-accepting process. Because the method only requires samples frozen immediately after recovery, the assay can be used for identifying microbial activity in subsurface ecosystems without the need to preserve live material. We measured potential hydrogen oxidation rates in all samples from multiple depths at several sites that collectively span a wide range of environmental conditions and biogeochemical zones. Potential activity normalized to total cell abundance ranges over five orders of magnitude and varies, dependent upon the predominant terminal electron acceptor. Lowest per-cell potential rates characterize the zone of nitrate reduction and highest per-cell potential rates occur in the methanogenic zone. Possible reasons for this relationship to predominant electron acceptor include (i) increasing importance of fermentation in successively deeper biogeochemical zones and (ii) adaptation of H(2)ases to successively higher concentrations of H-2 in successively deeper zones.

The behaviour of individuals, businesses, and government entities before, during, and immediately after a disaster can dramatically affect the impact and recovery time. However, existing risk-assessment methods rarely include this critical factor. In this Perspective, we show why this is a concern, and demonstrate that although initial efforts have inevitably represented human behaviour in limited terms, innovations in flood-risk assessment that integrate societal behaviour and behavioural adaptation dynamics into such quantifications may lead to more accurate characterization of risks and improved assessment of the effectiveness of risk-management strategies and investments. Such multidisciplinary approaches can inform flood-risk management policy development.

The climate is a complex dynamical system involving interactions and feedbacks among different processes at multiple temporal and spatial scales. Although numerous studies have attempted to understand the climate system, nonetheless, the studies investigating the multiscale characteristics of the climate are scarce. Further, the present set of techniques are limited in their ability to unravel the multi-scale variability of the climate system. It is completely plausible that extreme events and abrupt transitions, which are of great interest to climate community, are resultant of interactions among processes operating at multi-scale. For instance, storms, weather patterns, seasonal irregularities such as El Niño, floods and droughts, and decades-long climate variations can be better understood and even predicted by quantifying their multi-scale dynamics. This makes a strong argument to unravel the interaction and patterns of climatic processes at different scales. With this background, the thesis aims at developing measures to understand and quantify multi-scale interactions within the climate system. In the first part of the thesis, I proposed two new methods, viz, multi-scale event synchronization (MSES) and wavelet multi-scale correlation (WMC) to capture the scale-specific features present in the climatic processes. The proposed methods were tested on various synthetic and real-world time series in order to check their applicability and replicability. The results indicate that both methods (WMC and MSES) are able to capture scale-specific associations that exist between processes at different time scales in a more detailed manner as compared to the traditional single scale counterparts. In the second part of the thesis, the proposed multi-scale similarity measures were used in constructing climate networks to investigate the evolution of spatial connections within climatic processes at multiple timescales. The proposed methods WMC and MSES, together with complex network were applied to two different datasets. In the first application, climate networks based on WMC were constructed for the univariate global sea surface temperature (SST) data to identify and visualize the SSTs patterns that develop very similarly over time and distinguish them from those that have long-range teleconnections to other ocean regions. Further investigations of climate networks on different timescales revealed (i) various high variability and co-variability regions, and (ii) short and long-range teleconnection regions with varying spatial distance. The outcomes of the study not only re-confirmed the existing knowledge on the link between SST patterns like El Niño Southern Oscillation and the Pacific Decadal Oscillation, but also suggested new insights into the characteristics and origins of long-range teleconnections. In the second application, I used the developed non-linear MSES similarity measure to quantify the multivariate teleconnections between extreme Indian precipitation and climatic patterns with the highest relevance for Indian sub-continent. The results confirmed significant non-linear influences that were not well captured by the traditional methods. Further, there was a substantial variation in the strength and nature of teleconnection across India, and across time scales. Overall, the results from investigations conducted in the thesis strongly highlight the need for considering the multi-scale aspects in climatic processes, and the proposed methods provide robust framework for quantifying the multi-scale characteristics.

Sea surface temperature (SST) patterns can – as surface climate forcing – affect weather and climate at large distances. One example is El Niño-Southern Oscillation (ENSO) that causes climate anomalies around the globe via teleconnections. Although several studies identified and characterized these teleconnections, our understanding of climate processes remains incomplete, since interactions and feedbacks are typically exhibited at unique or multiple temporal and spatial scales. This study characterizes the interactions between the cells of a global SST data set at different temporal and spatial scales using climate networks. These networks are constructed using wavelet multi-scale correlation that investigate the correlation between the SST time series at a range of scales allowing instantaneously deeper insights into the correlation patterns compared to traditional methods like empirical orthogonal functions or classical correlation analysis. This allows us to identify and visualise regions of – at a certain timescale – similarly evolving SSTs and distinguish them from those with long-range teleconnections to other ocean regions. Our findings re-confirm accepted knowledge about known highly linked SST patterns like ENSO and the Pacific Decadal Oscillation, but also suggest new insights into the characteristics and origins of long-range teleconnections like the connection between ENSO and Indian Ocean Dipole.

Hydrologic regionalization deals with the investigation of homogeneity in watersheds and provides a classification of watersheds for regional analysis. The classification thus obtained can be used as a basis for mapping data from gauged to ungauged sites and can improve extreme event prediction. This paper proposes a wavelet power spectrum (WPS) coupled with the self-organizing map method for clustering hydrologic catchments. The application of this technique is implemented for gauged catchments. As a test case study, monthly streamflow records observed at 117 selected catchments throughout the western United States from 1951 through 2002. Further, based on WPS of each station, catchments are classified into homogeneous clusters, which provides a representative WPS pattern for the streamflow stations in each cluster.

Regionalization and pooling stations to form homogeneous regions or communities are essential for reliable parameter transfer, prediction in ungauged basins, and estimation of missing information. Over the years, several clustering methods have been proposed for regional analysis. Most of these methods are able to quantify the study region in terms of homogeneity but fail to provide microscopic information about the interaction between communities, as well as about each station within the communities. We propose a complex network-based approach to extract this valuable information and demonstrate the potential of our approach using a rainfall network constructed from the Indian gridded daily precipitation data. The communities were identified using the network-theoretical community detection algorithm for maximizing the modularity. Further, the grid points (nodes) were classified into universal roles according to their pattern of within- and between-community connections. The method thus yields zoomed-in details of individual rainfall grids within each community.

Hydrometric networks play a vital role in providing information for decision-making in water resource management. They should be set up optimally to provide as much information as possible that is as accurate as possible and, at the same time, be cost-effective. Although the design of hydrometric networks is a well-identified problem in hydrometeorology and has received considerable attention, there is still scope for further advancement. In this study, we use complex network analysis, defined as a collection of nodes interconnected by links, to propose a new measure that identifies critical nodes of station networks. The approach can support the design and redesign of hydrometric station networks. The science of complex networks is a relatively young field and has gained significant momentum over the last few years in different areas such as brain networks, social networks, technological networks, or climate networks. The identification of influential nodes in complex networks is an important field of research. We propose a new node-ranking measure – the weighted degree–betweenness (WDB) measure – to evaluate the importance of nodes in a network. It is compared to previously proposed measures used on synthetic sample networks and then applied to a real-world rain gauge network comprising 1229 stations across Germany to demonstrate its applicability. The proposed measure is evaluated using the decline rate of the network efficiency and the kriging error. The results suggest that WDB effectively quantifies the importance of rain gauges, although the benefits of the method need to be investigated in more detail.

This study aims to compare impacts of climate change on streamflow in four large representative African river basins: the Niger, the Upper Blue Nile, the Oubangui and the Limpopo. We set up the eco-hydrological model SWIM (Soil and Water Integrated Model) for all four basins individually. The validation of the models for four basins shows results from adequate to very good, depending on the quality and availability of input and calibration data. For the climate impact assessment, we drive the model with outputs of five bias corrected Earth system models of Coupled Model Intercomparison Project Phase 5 (CMIP5) for the representative concentration pathways (RCPs) 2.6 and 8.5. This climate input is put into the context of climate trends of the whole African continent and compared to a CMIP5 ensemble of 19 models in order to test their representativeness. Subsequently, we compare the trends in mean discharges, seasonality and hydrological extremes in the 21st century. The uncertainty of results for all basins is high. Still, climate change impact is clearly visible for mean discharges but also for extremes in high and low flows. The uncertainty of the projections is the lowest in the Upper Blue Nile, where an increase in streamflow is most likely. In the Niger and the Limpopo basins, the magnitude of trends in both directions is high and has a wide range of uncertainty. In the Oubangui, impacts are the least significant. Our results confirm partly the findings of previous continental impact analyses for Africa. However, contradictory to these studies we find a tendency for increased streamflows in three of the four basins (not for the Oubangui). Guided by these results, we argue for attention to the possible risks of increasing high flows in the face of the dominant water scarcity in Africa. In conclusion, the study shows that impact intercomparisons have added value to the adaptation discussion and may be used for setting up adaptation plans in the context of a holistic approach.

This study intends to contribute to the ongoing discussion on whether land use and land cover changes (LULC) or climate trends have the major influence on the observed increase of flood magnitudes in the Sahel. A simulation-based approach is used for attributing the observed trends to the postulated drivers. For this purpose, the ecohydrological model SWIM (Soil and Water Integrated Model) with a new, dynamic LULC module was set up for the Sahelian part of the Niger River until Niamey, including the main tributaries Sirba and Goroul. The model was driven with observed, reanalyzed climate and LULC data for the years 1950-2009. In order to quantify the shares of influence, one simulation was carried out with constant land cover as of 1950, and one including LULC. As quantitative measure, the gradients of the simulated trends were compared to the observed trend. The modeling studies showed that for the Sirba River only the simulation which included LULC was able to reproduce the observed trend. The simulation without LULC showed a positive trend for flood magnitudes, but underestimated the trend significantly. For the Goroul River and the local flood of the Niger River at Niamey, the simulations were only partly able to reproduce the observed trend. In conclusion, the new LULC module enabled some first quantitative insights into the relative influence of LULC and climatic changes. For the Sirba catchment, the results imply that LULC and climatic changes contribute in roughly equal shares to the observed increase in flooding. For the other parts of the subcatchment, the results are less clear but show, that climatic changes and LULC are drivers for the flood increase; however their shares cannot be quantified. Based on these modeling results, we argue for a two-pillar adaptation strategy to reduce current and future flood risk: Flood mitigation for reducing LULC-induced flood increase, and flood adaptation for a general reduction of flood vulnerability.

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.

Central Asia is located at the confluence of large-scale atmospheric circulation systems. It is thus likely to be highly susceptible to changes in the dynamics of those systems; however, little is still known about the regional paleoclimate history. Here we present carbon and hydrogen isotopic compositions of n-alkanoic acids from a late Holocene sediment core from Lake Karakuli (eastern Pamir, Xinjiang Province, China). Instrumental evidence and isotopeenabled climate model experiments with the Laboratoire de Meteorologie Dynamique Zoom model version 4 (LMDZ4) demonstrate that delta D values of precipitation in the region are influenced by both temperature and precipitation amount. We find that these parameters are inversely correlated on an annual scale, i.e., the climate has varied between relatively cool and wet and more warm and dry over the last 50 years. Since the isotopic signals of these changes are in the same direction and therefore additive, isotopes in precipitation are sensitive recorders of climatic changes in the region. Additionally, we infer that plants use year-round precipitation (including snowmelt), and thus leaf wax delta D values must also respond to shifts in the proportion of moisture derived from westerly storms during late winter and early spring. Downcore results give evidence for a gradual shift to cooler and wetter climates between 3.5 and 2.5 cal kyr BP, interrupted by a warm and dry episode between 3.0 and 2.7 kyr BP. Further cool and wet episodes occur between 1.9 and 1.5 and between 0.6 and 0.1 kyr BP, the latter coeval with the Little Ice Age. Warm and dry episodes from 2.5 to 1.9 and 1.5 to 0.6 kyr BP coincide with the Roman Warm Period and Medieval Climate Anomaly, respectively. Finally, we find a drying tend in recent decades. Regional comparisons lead us to infer that the strength and position of the westerlies, and wider northern hemispheric climate dynamics, control climatic shifts in arid Central Asia, leading to complex local responses. Our new archive from Lake Karakuli provides a detailed record of the local signatures of these climate transitions in the eastern Pamir.

Sedimentary lipid biomarkers have become widely used tools for reconstructing past climatic and ecological changes due to their ubiquitous occurrence in lake sediments. In particular, the hydrogen isotopic composition (expressed as delta D values) of leaf wax lipids derived from terrestrial plants has been a focus of research during the last two decades and the understanding of competing environmental and plant physiological factors influencing the delta D values has greatly improved. Comparatively less attention has been paid to lipid biomarkers derived from aquatic plants, although these compounds are abundant in many lacustrine sediments. We therefore conducted a field and laboratory experiment to study the effect of salinity and groundwater discharge on the isotopic composition of aquatic plant biomarkers. We analyzed samples of the common submerged plant species, Potamogeton pectinatus (sago pondweed), which has a wide geographic distribution and can tolerate high salinity. We tested the effect of groundwater discharge (characterized by more negative delta D values relative to lake water) and salinity on the delta D values of n-alkanes from P. pectinatus by comparing plants (i) collected from the oligotrophic freshwater Lake Stechlin (Germany) at shallow littoral depth from locations with and without groundwater discharge, and (ii) plants grown from tubers collected from the eutrophic Lake Muggelsee in nutrient solution at four salinity levels. Isotopically depleted groundwater did not have a significant influence on the delta D values of n-alkanes in Lake Stechlin P. pectinatus and calculated isotopic fractionation factors epsilon(l/w) between lake water and n-alkanes averaged -137 +/- 9%(n-C-23), -136 +/- 7%(n-C-25) and -131 +/- 6%(n-C-27), respectively. Similar epsilon values were calculated for plants from Lake Muggelsee grown in freshwater nutrient solution (-134 +/- 11% for n-C-23), while greater fractionation was observed at increased salinity values of 10 (163 +/- 12%) and 15(-172 +/- 15%). We therefore suggest an average e value of -136 +/- 9% between source water and the major n-alkanes in P. pectinatus grown under freshwater conditions. Our results demonstrate that isotopic fractionation can increase by 30-40% at salinity values 10 and 15. These results could be explained either by inhibited plant growth at higher salinity, or by metabolic adaptation to salt stress that remain to be elucidated. A potential salinity effect on dD values of aquatic lipids requires further examination, since this would impact on the interpretation of downcore isotopic data in paleohydrologic studies. (C) 2017 Elsevier Ltd. All rights reserved.

Python is at the forefront of scientific computation for seismologists and therefore should be introduced to students interested in becoming seismologists. On its own, Python is open source and well designed with extensive libraries. However, Python code can also be executed, visualized, and communicated to others with "Jupyter Notebooks". Thus, Jupyter Notebooks are ideal for teaching students Python and scientific computation. In this article, we designed an openly available Python library and collection of Jupyter Notebooks based on defined scientific computation learning goals for seismology students. The Notebooks cover topics from an introduction to Python to organizing data, earthquake catalog statistics, linear regression, and making maps. Our Python library and collection of Jupyter Notebooks are meant to be used as course materials for an upper-division data analysis course in an Earth Science Department, and the materials were tested in a Probabilistic Seismic Hazard course. However, seismologists or anyone else who is interested in Python for data analysis and map making can use these materials.

Afyon Zone, which was derived from the Anatolide-Tauride platform during closure of the Neo-Tethys, is made up of pre-Mesozoic basement and unconformably overlying Triassic-Early Tertiary cover series. The Afyon Zone contains widespread metavolcanic rocks, which are dominated by rhyolite, dacite, and trachyandesite. They form a distinct volcanic succession, which is separated from the underlying Silurian-Lower Carboniferous metacarbonates and meta-siliciclastics by a regional unconformity. Trachyandesitic metavolcanics are made up of massive lava flows, pyroclastics and epiclastics, less frequently, domes and dikes, which were developed on a deeply eroded subaerial landmass. U/Pb and Pb/Pb zircon geochronology yielded Lower Triassic (similar to 250 Ma) ages, which are interpreted as extrusion age of trachyandesitic volcanics. Based on the stratigraphic, geochronological, and geochemical data, we suggest that these Lower Triassic magmatic rocks represent an extensional tectonic setting on the northern active margin of the Gondwana, which led to the development of the northern branch of the Neo-Tethys.

Sinkholes and depressions are typical landforms of karst regions. They pose a considerable natural hazard to infrastructure, agriculture, economy and human life in affected areas worldwide. The physio-chemical processes of sinkholes and depression formation are manifold, ranging from dissolution and material erosion in the subsurface to mechanical subsidence/failure of the overburden. This thesis addresses the mechanisms leading to the development of sinkholes and depressions by using complementary methods: remote sensing, distinct element modelling and near-surface geophysics. In the first part, detailed information about the (hydro)-geological background, ground structures, morphologies and spatio-temporal development of sinkholes and depressions at a very active karst area at the Dead Sea are derived from satellite image analysis, photogrammetry and geologic field surveys. There, clusters of an increasing number of sinkholes have been developing since the 1980s within large-scale depressions and are distributed over different kinds of surface materials: clayey mud, sandy-gravel alluvium and lacustrine evaporites (salt). The morphology of sinkholes differs depending in which material they form: Sinkholes in sandy-gravel alluvium and salt are generally deeper and narrower than sinkholes in the interbedded evaporite and mud deposits. From repeated aerial surveys, collapse precursory features like small-scale subsidence, individual holes and cracks are identified in all materials. The analysis sheds light on the ongoing hazardous subsidence process, which is driven by the base-level fall of the Dead Sea and by the dynamic formation of subsurface water channels. In the second part of this thesis, a novel, 2D distinct element geomechanical modelling approach with the software PFC2D-V5 to simulating individual and multiple cavity growth and sinkhole and large-scale depression development is presented. The approach involves a stepwise material removal technique in void spaces of arbitrarily shaped geometries and is benchmarked by analytical and boundary element method solutions for circular cavities. Simulated compression and tension tests are used to calibrate model parameters with bulk rock properties for the materials of the field site. The simulations show that cavity and sinkhole evolution is controlled by material strength of both overburden and cavity host material, the depth and relative speed of the cavity growth and the developed stress pattern in the subsurface. Major findings are: (1) A progressively deepening differential subrosion with variable growth speed yields a more fragmented stress pattern with stress interaction between the cavities. It favours multiple sinkhole collapses and nesting within large-scale depressions. (2) Low-strength materials do not support large cavities in the material removal zone, and subsidence is mainly characterised by gradual sagging into the material removal zone with synclinal bending. (3) High-strength materials support large cavity formation, leading to sinkhole formation by sudden collapse of the overburden. (4) Large-scale depression formation happens either by coalescence of collapsing holes, block-wise brittle failure, or gradual sagging and lateral widening. The distinct element based approach is compared to results from remote sensing and geophysics at the field site. The numerical simulation outcomes are generally in good agreement with derived morphometrics, documented surface and subsurface structures as well as seismic velocities. Complementary findings on the subrosion process are provided from electric and seismic measurements in the area. Based on the novel combination of methods in this thesis, a generic model of karst landform evolution with focus on sinkhole and depression formation is developed. A deepening subrosion system related to preferential flow paths evolves and creates void spaces and subsurface conduits. This subsequently leads to hazardous subsidence, and the formation of sinkholes within large-scale depressions. Finally, a monitoring system for shallow natural hazard phenomena consisting of geodetic and geophysical observations is proposed for similarly affected areas.

Mechanical and/or chemical removal of material from the subsurface may generate large subsurface cavities, the destabilisation of which can lead to ground collapse and the formation of sinkholes. Numerical simulation of the interaction of cavity growth, host material deformation and overburden collapse is desirable to better understand the sinkhole hazard but is a challenging task due to the involved high strains and material discontinuities. Here, we present 2-D distinct element method numerical simulations of cavity growth and sinkhole development. Firstly, we simulate cavity formation by quasi-static, stepwise removal of material in a single growing zone of an arbitrary geometry and depth. We benchmark this approach against analytical and boundary element method models of a deep void space in a linear elastic material. Secondly, we explore the effects of properties of different uniform materials on cavity stability and sinkhole development. We perform simulated biaxial tests to calibrate macroscopic geotechnical parameters of three model materials representative of those in which sinkholes develop at the Dead Sea shoreline: mud, alluvium and salt. We show that weak materials do not support large cavities, leading to gradual sagging or suffusion-style subsidence. Strong materials support quasi-stable to stable cavities, the overburdens of which may fail suddenly in a caprock or bedrock collapse style. Thirdly, we examine the consequences of layered arrangements of weak and strong materials. We find that these are more susceptible to sinkhole collapse than uniform materials not only due to a lower integrated strength of the overburden but also due to an inhibition of stabilising stress arching. Finally, we compare our model sinkhole geometries to observations at the Ghor Al-Haditha sinkhole site in Jordan. Sinkhole depth ∕ diameter ratios of 0.15 in mud, 0.37 in alluvium and 0.33 in salt are reproduced successfully in the calibrated model materials. The model results suggest that the observed distribution of sinkhole depth ∕ diameter values in each material type may partly reflect sinkhole growth trends.

The Najd, Oman, is located in one of the most arid environments in the world. The groundwater in this region is occurring in four different aquifers A to D of the Hadhramaut Group consisting mainly of different types of limestone and dolomite. The quality of the groundwater is dominated by the major ions sodium, calcium, magnesium, sulphate, and chloride, but the hydrochemical character is varying among the four aquifers. Mineralization within the separate aquifers increases along the groundwater flow direction from south to north-northeast up to high saline sodium-chloride water in aquifer D in the northeast area of the Najd. Environmental isotope analyses of hydrogen and oxygen were conducted to monitor the groundwater dynamics and to evaluate the recharge conditions of groundwater into the Najd aquifers. Results suggest an earlier recharge into these aquifers as well as ongoing recharge takes place in the region down to present day. Mixing of modern and submodern waters was detected by water isotopes in aquifer D in the mountain chain (Jabal) area and along the northern side of the mountain range. In addition, delta H-2 and delta O-18 variations suggest that aquifers A, B, and C are assumed to be connected by faults and fractures, and interaction between the aquifers may occur. Low tritium concentrations support the mixing assumption in the recharge area. The knowledge about the groundwater development is an important factor for the sustainable use of water resources in the Dhofar region.

A parameterization for the motion of ice-shelf fronts on a Cartesian grid in finite-difference land-ice models is presented. The scheme prevents artificial thinning of the ice shelf at its edge, which occurs due to the finite resolution of the model. The intuitive numerical implementation diminishes numerical dispersion at the ice front and enables the application of physical boundary conditions to improve the calculation of stress and velocity fields throughout the ice-sheet-shelf system. Numerical properties of this subgrid modification are assessed in the Potsdam Parallel Ice Sheet Model (PISM-PIK) for different geometries in one and two horizontal dimensions and are verified against an analytical solution in a flow-line setup.

We have redetermined focal depths of moderate and major earthquakes with reported lower-crust and upper-mantle depths that have occurred in Tien-Shan, since the availability of broad-band array data. Records of earthquakes at global arrays have been used for identification and modelling of depth phases in order to make accurate estimation of focal depths. Our results show that half of the purportedly deep earthquakes are indeed originating from depths attributable to middle-crust and lower-crust regions. Also one exceptional event in the northern foreland of Tien-Shan in Junggar Basin is located in the upper mantle at the depth of 64 km. Such unusually deep earthquakes for intraplate continental tectonic domain are all located at the margin of Tien-Shan with its adjacent stable blocks and at least some of them have occurred where the brittle behaviour of continental rocks is not highly expected. The reverse mechanisms of all these earthquakes and their proximity to formerly subducting and later colliding and underplating stable blocks and their interactions with overlying Tien-Shan are clues to explain this extremity.

Analysis of time-lapse ground-penetrating radar (GPR) data can provide information regarding subsurface hydrological processes, such as preferential flow. However, the analysis of time-lapse data is often limited by data quality; for example, for noisy input data, the interpretation of difference images is often difficult. Motivated by modern image-processing tools, we have developed two robust GPR attributes, which allow us to distinguish amplitude (contrast similarity) and time-shift (structural similarity) variations related to differences between individual time-lapse GPR data sets. We tested and evaluated our attributes using synthetic data of different complexity. Afterward, we applied them to a field data example, in which subsurface flow was induced by an artificial rainfall event. For all examples, we identified our structural similarity attribute to be a robust measure for highlighting time-lapse changes also in data with low signal-to-noise ratios. We determined that our new attribute-based workflow is a promising tool to analyze time-lapse GPR data, especially for imaging subsurface hydrological processes.

Flow phenomena in the unsaturated zone are highly variable in time and space. Thus, it is challenging to measure and monitor such processes under field conditions. Here, we present a new setup and interpretation approach for combining a dye tracer experiment with a 4D ground-penetrating radar (GPR) survey. Therefore, we designed a rainfall experiment during which we measured three surface-based 3D GPR surveys using a pair of 500 MHz antennas. Such a survey setup requires accurate acquisition and processing techniquesto extract time-lapse information supporting the interpretation of selected cross-sections photographed after excavating the site. Our results reveal patterns of traveltime changes in the measured GPR data, which are associated with soil moisture changes. As distinct horizons are present at our site, such changes can be quantified and transferred into changes in total soil moisture content. Our soil moisture estimates are similar to the amount of infiltrated water, which confirms our experimental approach and makes us confident for further developing this strategy, especially, with respect to improving the temporal and spatial resolution. (C) 2015 Elsevier B.V. All rights reserved.

In urban areas, leaf area index (LAI) is a key ecosystem structural attribute with implications for energy and water balance, gas exchange, and anthropogenic energy use. In this study, we estimated LAI spatially using airborne lidar in downtown Santa Barbara, California, USA. We implemented two different modeling approaches. First, we directly estimated effective LAI (LAIe) using scan angle- and clump-corrected lidar laser penetration metrics (LPM). Second, we adapted existing allometric equations to estimate crown structural metrics including tree height and crown base height using lidar. The latter approach allowed for LAI estimates at the individual tree-crown scale. The LPM method, at both high and decimated point densities, resulted in good linear agreement with estimates from ground-based hemispherical photography (r(2) = 0.82, y = 0.99x) using a model that assumed a spherical leaf angle distribution. Within individual tree crown segments, the lidar estimates of crown structure closely paralleled field measurements (e.g., r(2) = 0.87 for crown length). LAI estimates based on the lidar crown measurements corresponded well with estimates from field measurements (r(2) = 0.84, y = 0.97x + 0.10). Consistency of the LPM and allometric lidar methods was also strong at 71 validation plots (r(2) = 0.88) and at 450 additional sample locations across the entire study area (r(2) = 0.72). This level of correspondence exceeded that of the canopy hemispherical photography and allometric, ground-based estimates (r(2) = 0.53). The first-order alignment of these two disparate methods may indicate that the error bounds for mapping LAI in cities are small enough to pursue large scale, spatially explicit estimation. (C) 2015 Elsevier Inc All rights reserved.

Orthopyroxenes of a high temperature protomylonite of the Ivrea Zone, Northern Italy show twin like polysynthetic lamellae parallel to {210} of the hypersthene host. The transformation is caused by plastic deformation under high metamorphic conditions which has resulted in dynamic recrystallization of pyroxene and plagioclase. The lamellae consist of clinohypersthene. The twin plane and the lamellar clino-orthoinversion of hypersthene due to natural deformation have not been described hitherto.

The Garzn Complex of the Garzn Massif in SW Colombia is composed of the Vergel Granulite Unit (VG) and the Las Margaritas Migmatite Unit (LMM). Previous studies reveal peak temperature conditions for the VG of about 740 A degrees C. The present study considers the remarkable exsolution phenomena in feldspars and pyroxenes and titanium-in-quartz thermometry. Recalculated ternary feldspar compositions indicate temperatures around 900-1,000 A degrees C just at or above the ultra-high temperature-metamorphism (UHTM) boundary of granulites. The calculated temperatures range of exsolved ortho- and clinopyroxenes also supports the existence of an UHTM event. In addition, titanium-in-quartz thermometry points towards ultra-high temperatures. It is the first known UHTM crustal segment in the northern part of South America. Although a mean geothermal gradient of ca 38 A degrees C km(-1) could imply additional heat supply in the lower crust controlling this extreme of peak metamorphism, an alternative model is suggested. The formation of the Vergel Granulite Unit is supposed to be formed in a continental back-arc environment with a thinned and weakened crust behind a magmatic arc (Guapotn-Mancagua Gneiss) followed by collision. In contrast, rocks of the adjacent Las Margaritas Migmatite Unit display "normal" granulite facies temperatures and are formed in a colder lower crust outside the arc, preserved by the Guapotn-Mancagu Gneiss. Back-arc formation was followed by inversion and thickening of the basin. The three units that form the modern-day Garzn Massif, were juxtaposed upon each other during collision (at ca. 1,000 Ma) and exhumation. The collision leading to the deformation of the studied area is part of the Grenville orogeny leading to the amalgamation of Rodinia.

Pseudotachylyte veins frequently associated with mylonites and ultramylonites occur within migmatitic paragneisses, metamonzodiorites, as well as felsic and mafic granulites at the base of the section of the Hercynian lower crust exposed in Calabria (Southern Italy). The crustal section is tectonically superposed on lower grade units. Ultramylonites and pseudotachylytes are particularly well developed in migmatitic paragneisses, whereas sparse fault-related pseudotachylytes and thin mylonite/ultramylonite bands occur in granulite-facies rocks. The presence of sillimanite and clinopyroxene in ultramylonites and mylonites indicates that relatively high-temperature conditions preceded the formation of pseudotachylytes. We have analysed pseudotachylytes from different rock types to ascertain their deep crustal origin and to better understand the relationships between brittle and ductile processes during deformation of the deeper crust. Different protoliths were selected to test how lithology controls pseudotachylyte composition and textures. In migmatites and felsic granulites, euhedral or cauliflower-shaped garnets directly crystallized from pseudotachylyte melts of near andesitic composition. This indicates that pseudotachylytes originated at deep crustal conditions (> 0.75 GPa). In mafic protoliths, quenched needle-to-feather-shaped high-alumina orthopyroxene occurs in contact with newly crystallized plagioclase. The pyroxene crystallizes in garnet-free and garnet-bearing veins. The simultaneous growth of orthopyroxene and plagioclase as well as almandine, suggests lower crustal origin, with pressures in excess of 0.85 GPa. The existence of melts of different composition in the same vein indicates the stepwise, non-equilibrium conditions of frictional melting. Melt formed and intruded into pre-existing anisotropies. In mafic granulites, brittle faulting is localized in a previously formed thin high-temperature mylonite bands. migmatitic gneisses are deformed into ultramylonite domains characterized by s-c fabric. Small grain size and fluids lowered the effective stress on the c planes favouring a seismic event and the consequent melt generation. Microstructures and ductile deformation of pseudotachylytes suggest continuous ductile flow punctuated by episodes of high-strain rate, leading to seismic events and melting.

The formation of the Pamir is a key component of the India-Asia collision with major implications for lithospheric processes, plateau formation, land-sea configurations and associated climate changes. Although the formation of the Pamir is traditionally linked to Cenozoic processes associated with the India-Asia collision, the contribution of the Mesozoic tectonic evolution remains poorly understood. The Pamir was formed by the suturing of Gondwanan terranes to the south margin of Eurasia, however, the timing and tectonic mechanisms associated with this Mesozoic accretion remain poorly constrained. These processes are recorded by several igneous belts within these terranes, which are not well studied. Within the Southern Pamir, the Albian-Turonian volcanic rocks and comagmatic plutons of the Kyzylrabat Igneous Complex (KIC) provide an important and still unconstrained record of the Pamir evolution. Here we provide the age, origin and the geodynamic setting of the KIC volcanics by studying their petrology, zircon U-Pb geochronology, geochemistry and isotope composition.17 samples from the KIC volcanics yield U-Pb ages spanning from 92 to 110 Ma. The volcanics are intermediate to acidic in composition (SiO2 = 56-69 wt%) and exhibit high-K calc-alkaline and shoshonitic affinity (K2O/Na2O = 12.2 wt%). They show enrichment in LILE and LREE and depletion in HFSE and HREE with negative Ta, Ti and Nb anomalies, suggesting an arc-related tectonic setting for their formation. Low sNd(t) values (from 9.1 to 4.7), relatively high Sr-87/Sr-86(i) ratios (0.7069-0.7096) and broad range of zircon stif values (from 22.6 to 1.5) suggest a mixture of different magma sources. These features suggest that volcanics were derived by crustal under- or intraplating of an enriched subduction-related mantle shoshonitic magmas, by heating and partial melting of the lower crust, and by mixing of both magma components. Our results further imply that the KIC volcanics represent a shoshonitic suite typical of an evolution from active continental arc to post-collisional setting with a steepening of the Benioff zone and thickening of the crust toward the back-arc. This setting is best explained by the subduction- collision transition along the Shyok suture due to accretion of the Kohistan island arc to the Karakoram. This suggests that a significant part of the crustal shortening and thickening accommodated in the Pamir occurred in the Mesozoic before the India-Asia collision with implications for regional tectonic models. This further suggests the Pamir was already a major topographic feature with potentially important paleoclimate forcing such as the monsoonal circulation. (C) 2017 International Association for Gondwana Research. Published by Elsevier B.V. All rights reserved.

The study of outcrop modeling is located at the interface between two fields of expertise, Sedimentology and Computing Geoscience, which respectively investigates and simulates geological heterogeneity observed in the sedimentary record. During the last past years, modeling tools and techniques were constantly improved. In parallel, the study of Phanerozoic carbonate deposits emphasized the common occurrence of a random facies distribution along single depositional domain. Although both fields of expertise are intrinsically linked during outcrop simulation, their respective advances have not been combined in literature to enhance carbonate modeling studies. The present study re-examines the modeling strategy adapted to the simulation of shallow-water carbonate systems, based on a close relationship between field sedimentology and modeling capabilities. In the present study, the evaluation of three commonly used algorithms Truncated Gaussian Simulation (TGSim), Sequential Indicator Simulation (SISim), and Indicator Kriging (IK), were performed for the first time using visual and quantitative comparisons on an ideally suited carbonate outcrop. The results show that the heterogeneity of carbonate rocks cannot be fully simulated using one single algorithm. The operating mode of each algorithm involves capabilities as well as drawbacks that are not capable to match all field observations carried out across the modeling area. Two end members in the spectrum of carbonate depositional settings, a low-angle Jurassic ramp (High Atlas, Morocco) and a Triassic isolated platform (Dolomites, Italy), were investigated to obtain a complete overview of the geological heterogeneity in shallow-water carbonate systems. Field sedimentology and statistical analysis performed on the type, morphology, distribution, and association of carbonate bodies and combined with palaeodepositional reconstructions, emphasize similar results. At the basin scale (x 1 km), facies association, composed of facies recording similar depositional conditions, displays linear and ordered transitions between depositional domains. Contrarily, at the bedding scale (x 0.1 km), individual lithofacies type shows a mosaic-like distribution consisting of an arrangement of spatially independent lithofacies bodies along the depositional profile. The increase of spatial disorder from the basin to bedding scale results from the influence of autocyclic factors on the transport and deposition of carbonate sediments. Scale-dependent types of carbonate heterogeneity are linked with the evaluation of algorithms in order to establish a modeling strategy that considers both the sedimentary characteristics of the outcrop and the modeling capabilities. A surface-based modeling approach was used to model depositional sequences. Facies associations were populated using TGSim to preserve ordered trends between depositional domains. At the lithofacies scale, a fully stochastic approach with SISim was applied to simulate a mosaic-like lithofacies distribution. This new workflow is designed to improve the simulation of carbonate rocks, based on the modeling of each scale of heterogeneity individually. Contrarily to simulation methods applied in literature, the present study considers that the use of one single simulation technique is unlikely to correctly model the natural patterns and variability of carbonate rocks. The implementation of different techniques customized for each level of the stratigraphic hierarchy provides the essential computing flexibility to model carbonate systems. Closer feedback between advances carried out in the field of Sedimentology and Computing Geoscience should be promoted during future outcrop simulations for the enhancement of 3-D geological models.

With accelerating climate cooling in the late Cenozoic, glacial and periglacial erosion became more widespread on the surface of the Earth. The resultant shift in erosion patterns significantly changed the large-scale morphology of many mountain ranges worldwide. Whereas the glacial fingerprint is easily distinguished by its characteristic fjords and U-shaped valleys, the periglacial fingerprint is more subtle but potentially prevails in some mid- to high-latitude landscapes. Previous models have advocated a frost-driven control on debris production at steep headwalls and glacial valley sides. Here we investigate the important role that periglacial processes also play in less steep parts of mountain landscapes. Understanding the influences of frost-driven processes in low-relief areas requires a focus on the consequences of an accreting soil mantle, which characterises such surfaces. We present a new model that quantifies two key physical processes: frost cracking and frost creep, as a function of both temperature and sediment thickness. Our results yield new insights into how climate and sediment transport properties combine to scale the intensity of periglacial processes. The thickness of the soil mantle strongly modulates the relation between climate and the intensity of mechanical weathering and sediment flux. Our results also point to an offset between the conditions that promote frost cracking and those that promote frost creep, indicating that a stable climate can provide optimal conditions for only one of those processes at a time. Finally, quantifying these relations also opens up the possibility of including periglacial processes in large-scale, long-term landscape evolution models, as demonstrated in a companion paper.

Subsea permafrost is perennially cryotic earth material that lies offshore. Most submarine permafrost is relict terrestrial permafrost beneath the Arctic shelf seas, was inundated after the last glaciation, and has been warming and thawing ever since. It is a reservoir and confining layer for gas hydrates and has the potential to release greenhouse gases and affect global climate change. Furthermore, subsea permafrost thaw destabilizes coastal infrastructure. While numerous studies focus on its distribution and rate of thaw over glacial timescales, these studies have not been brought together and examined in their entirety to assess rates of thaw beneath the Arctic Ocean. In addition, there is still a large gap in our understanding of sub-aquatic permafrost processes on finer spatial and temporal scales. The degradation rate of subsea permafrost is influenced by the initial conditions upon submergence. Terrestrial permafrost that has already undergone warming, partial thawing or loss of ground ice may react differently to inundation by seawater compared to previously undisturbed ice-rich permafrost. Heat conduction models are sufficient to model the thaw of thick subsea permafrost from the bottom, but few studies have included salt diffusion for top-down chemical degradation in shallow waters characterized by mean annual cryotic conditions on the seabed. Simulating salt transport is critical for assessing degradation rates for recently inundated permafrost, which may accelerate in response to warming shelf waters, a lengthening open water season, and faster coastal erosion rates. In the nearshore zone, degradation rates are also controlled by seasonal processes like bedfast ice, brine injection, seasonal freezing under floating ice conditions and warm freshwater discharge from large rivers. The interplay of all these variables is complex and needs further research. To fill this knowledge gap, this thesis investigates sub-aquatic permafrost along the southern coast of the Bykovsky Peninsula in eastern Siberia. Sediment cores and ground temperature profiles were collected at a freshwater thermokarst lake and two thermokarst lagoons in 2017. At this site, the coastline is retreating, and seawater is inundating various types of permafrost: sections of ice-rich Pleistocene permafrost (Yedoma) cliffs at the coastline alternate with lagoons and lower elevation previously thawed and refrozen permafrost basins (Alases). Electrical resistivity surveys with floating electrodes were carried out to map ice-bearing permafrost and taliks (unfrozen zones in the permafrost, usually formed beneath lakes) along the diverse coastline and in the lagoons. Combined with the borehole data, the electrical resistivity results permit estimation of contemporary ice-bearing permafrost characteristics, distribution, and occasionally, thickness. To conceptualize possible geomorphological and marine evolutionary pathways to the formation of the observed layering, numerical models were applied. The developed model incorporates salt diffusion and seasonal dynamics at the seabed, including bedfast ice. Even along coastlines with mean annual non-cryotic boundary conditions like the Bykovsky Peninsula, the modelling results show that salt diffusion minimizes seasonal freezing of the seabed, leading to faster degradation rates compared to models without salt diffusion. Seasonal processes are also important for thermokarst lake to lagoon transitions because lagoons can generate cold hypersaline conditions underneath the ice cover. My research suggests that ice-bearing permafrost can form in a coastal lagoon environment, even under floating ice. Alas basins, however, may degrade more than twice as fast as Yedoma permafrost in the first several decades of inundation. In addition to a lower ice content compared to Yedoma permafrost, Alas basins may be pre-conditioned with salt from adjacent lagoons. Considering the widespread distribution of thermokarst in the Arctic, its integration into geophysical models and offshore surveys is important to quantify and understand subsea permafrost degradation and aggradation. Through numerical modelling, fieldwork, and a circum-Arctic review of subsea permafrost literature, this thesis provides new insights into sub-aquatic permafrost evolution in saline coastal environments.

The concept of hydrologic connectivity summarizes all flow processes that link separate regions of a landscape. As such, it is a central theme in the field of catchment hydrology, with influence on neighboring disciplines such as ecology and geomorphology. It is widely acknowledged to be an important key in understanding the response behavior of a catchment and has at the same time inspired research on internal processes over a broad range of scales. From this process-hydrological point of view, hydrological connectivity is the conceptual framework to link local observations across space and scales. This is the context in which the four studies this thesis comprises of were conducted. The focus was on structures and their spatial organization as important control on preferential subsurface flow. Each experiment covered a part of the conceptualized flow path from hillslopes to the stream: soil profile, hillslope, riparian zone, and stream. For each study site, the most characteristic structures of the investigated domain and scale, such as slope deposits and peat layers were identified based on preliminary or previous investigations or literature reviews. Additionally, further structural data was collected and topographical analyses were carried out. Flow processes were observed either based on response observations (soil moisture changes or discharge patterns) or direct measurement (advective heat transport). Based on these data, the flow-relevance of the characteristic structures was evaluated, especially with regard to hillslope to stream connectivity. Results of the four studies revealed a clear relationship between characteristic spatial structures and the hydrological behavior of the catchment. Especially the spatial distribution of structures throughout the study domain and their interconnectedness were crucial for the establishment of preferential flow paths and their relevance for large-scale processes. Plot and hillslope-scale irrigation experiments showed that the macropores of a heterogeneous, skeletal soil enabled preferential flow paths at the scale of centimeters through the otherwise unsaturated soil. These flow paths connected throughout the soil column and across the hillslope and facilitated substantial amounts of vertical and lateral flow through periglacial slope deposits. In the riparian zone of the same headwater catchment, the connectivity between hillslopes and stream was controlled by topography and the dualism between characteristic subsurface structures and the geomorphological heterogeneity of the stream channel. At the small scale (1 m to 10 m) highest gains always occurred at steps along the longitudinal streambed profile, which also controlled discharge patterns at the large scale (100 m) during base flow conditions (number of steps per section). During medium and high flow conditions, however, the impact of topography and parafluvial flow through riparian zone structures prevailed and dominated the large-scale response patterns. In the streambed of a lowland river, low permeability peat layers affected the connectivity between surface water and groundwater, but also between surface water and the hyporheic zone. The crucial factor was not the permeability of the streambed itself, but rather the spatial arrangement of flow-impeding peat layers, causing increased vertical flow through narrow “windows” in contrast to predominantly lateral flow in extended areas of high hydraulic conductivity sediments. These results show that the spatial organization of structures was an important control for hydrological processes at all scales and study areas. In a final step, the observations from different scales and catchment elements were put in relation and compared. The main focus was on the theoretical analysis of the scale hierarchies of structures and processes and the direction of causal dependencies in this context. Based on the resulting hierarchical structure, a conceptual framework was developed which is capable of representing the system’s complexity while allowing for adequate simplifications. The resulting concept of the parabolic scale series is based on the insight that flow processes in the terrestrial part of the catchment (soil and hillslopes) converge. This means that small-scale processes assemble and form large-scale processes and responses. Processes in the riparian zone and the streambed, however, are not well represented by the idea of convergence. Here, the large-scale catchment signal arrives and is modified by structures in the riparian zone, stream morphology, and the small-scale interactions between surface water and groundwater. Flow paths diverge and processes can better be represented by proceeding from large scales to smaller ones. The catchment-scale representation of processes and structures is thus the conceptual link between terrestrial hillslope processes and processes in the riparian corridor.

We have undertaken petrographic, mineralogical, geochemical and isotopic investigations on carbonate minerals found within a 10-m-long core from Lonar Lake, central India, with the aim of evaluating their potential as palaeoenvironmental proxies. The core encompasses the entire Holocene and is the first well-dated high-resolution record from central India. While calcite and/or aragonite were found throughout the core, the mineral gaylussite was found only in two specific intervals (46303890 and 2040560 cal a BP). Hydrochemical and isotope data from inflowing streams and lake waters indicate that evaporitic processes play a dominant role in the precipitation of carbonates within this lake. Isotopic (18O and 13C) studies on the evaporative gaylussite crystals and residual bulk carbonates (calcite) from the long core show that evaporation is the major control on 18O enrichment in both the minerals. However, in case of 13C additional mechanisms, for example methanogenesis (gaylussite) and phytoplankton productivity (calcium carbonate), play an additional important role in some intervals. We also discuss the relevance of our investigation for palaeoclimate reconstruction and late Holocene monsoon variability.

Here we present high-resolution delta O-18 records obtained from speleothems collected in the eastern Bolivian Andes. The stable isotope records are related to the regional- to large-scale atmospheric circulation over South America and allow interpreting changes in delta O-18 during the last 1400 yr as a function of changes in precipitation regimes over the southern tropical Andes. Two distinct phases with more negative delta O-18 values, interpreted as periods of increased convective activity over the eastern Andean Cordillera in Bolivia are observed concomitantly with periods of global climate anomalies during the last millennium, such as the Medieval Climate Anomaly (MCA) and the Little Ice Age (LIA) respectively. Changes in the Bolivian delta O-18 record during the LIA are apparently related to a southward displacement of the Intertropical Convergence Zone (ITCZ), which acts as a main moisture driver to intensify convection over the tropical continent. During the MCA, however, the increased convective activity observed in the Bolivian record is likely the result of a different mechanism, which implies moisture sourced mainly from the southern tropical Atlantic. This interpretation is consistent with paleoclimate records further to the north in the tropical Andes that show progressively drier conditions during this time period, indicating a more northerly position of the ITCZ. The transition period between the MCA and the LIA shows a slight tendency toward increased delta O-18 values, indicating weakened convective activity. Our results also reveal a non-stationary anti-phased behavior between the delta O-18 reconstructions from Bolivia and northeastern Brazil that confirms a continental-scale east-west teleconnection across South America during the LIA.

We present a model of the electrical resistivity structure of the lithosphere in the Central Andes between 20 degrees and 24 degrees S from 3-D inversion of 56 long-period magnetotelluric sites. Our model shows a complex resistivity structure with significant variability parallel and perpendicular to the trench direction. The continental forearc is characterized mainly by high electrical resistivity (>1,000m), suggesting overall low volumes of fluids. However, low resistivity zones (LRZs, <5m) were found in the continental forearc below areas where major trench-parallel faults systems intersect NW-SE transverse faults. Forearc LRZs indicate circulation and accumulation of fluids in highly permeable fault zones. The continental crust along the arc shows three distinctive resistivity domains, which coincide with segmentation in the distribution of volcanoes. The northern domain (20 degrees-20.5 degrees S) is characterized by resistivities >1,000m and the absence of active volcanism, suggesting the presence of a low-permeability block in the continental crust. The central domain (20.5 degrees-23 degrees S) exhibits a number of LRZs at varying depths, indicating different levels of a magmatic plumbing system. The southern domain (23 degrees-24 degrees S) is characterized by resistivities >1,000m, suggesting the absence of large magma reservoirs below the volcanic chain at crustal depths. Magma reservoirs located below the base of the crust or in the backarc may fed active volcanism in the southern domain. In the subcontinental mantle, the model exhibits LRZs in the forearc mantle wedge and above clusters of intermediate-depth seismicity, likely related to fluids produced by serpentinization of the mantle and eclogitization of the slab, respectively.

En el presente trabajo se realizó una investigación multidisciplinaria combinando métodos de geomorfología tectónica con estudios geofisicos y estructurales, focalizados principalmente en la caracterización neotectónica de ambos faldeos de la sierra de La Candelaria y del extremo sur de la cuenca de Metán. La zona de estudio se encuentra ubicada en la región limítrofe entre las provincias de Salta y Tucumán y pertenece a la provincia geológica del Sistema Santa Bárbara. El principal objetivo consistió en contextualizar las evidencias de actividad tectónica cuaternaria de la región mediante la propuesta de un modelo estructural novedoso, con el propósito de incrementar la información disponible sobre estructuras neotectónicas y su potencial sismogénico. Con este fin, se aplicaron e integraron diversas técnicas tales como la interpretación de líneas sísmicas de reflexión, construcción de secciones estructurales balanceadas, y métodos geofísicos someros, para constatar el comportamiento en profundidad tanto de las estructuras geológicas identificadas en superficie como de las posibles fallas ciegas corticales involucradas. En primer lugar, se realizó un relevamiento regional del área de estudio empleando imágenes satelitales multiespectrales LANDSAT y SENTINEL 2, que permitieron reconocer diferentes niveles de abanicos aluviales y terrazas fluviales cuaternarios. Mediante la determinación de diferentes indicadores morfométricos en modelos de elevación digital (MED), junto con observaciones de campo, fue posible identificar evidencias de deformación sobre dichos niveles cuaternarios que han sido relacionadas genéticamente con cuatro fallas neotectónicas. Tres de ellas (fallas Arias, El Quemado y Copo Quile) fueron seleccionadas para efectuar estudios de mayor detalle por medio de la aplicación de métodos de geofísica somera (tomografía eléctrica resistiva (ERT) y tomografía sísmica de refracción Sísmica (SRT)), que permitieron corroborar su existencia en profundidad, realizar inferencias geométricas y cinemáticas, y estimar la magnitud de la deformación reciente. Las fallas Arias y El Quemado fueron interpretadas como fallas inversas relacionadas con deslizamiento flexural interstratal, mientras que la falla Copo Quile se interpretó como una falla inversa ciega de bajo ángulo.También se realizó una interpretación conjunta de líneas sísmicas de reflexión y pozos exploratorios pertenecientes a áreas hidrocarburíferas de las cuencas de Choromoro y Metán con el fin de contextualizar las principales estructuras reconocidas en el marco estratigráfico y tectónico regional. Toda la información fue integrada en una sección estructural balanceada mediante técnicas de modelado cinemático. Dicho modelo permite inferir que la deformación cuaternaria reconocida está relacionada al desplazamiento del basamento a lo largo de un corrimiento ciego, responsable del levantamiento de la sierra de La Candelaria y el cerr Cantero. Asimismo, el modelo cinemático permite interpretar la ubicación aproximada de los principales niveles de despegue que controlan el estilo de deformación. El nivel de despegue más somero, que controla la deformación de la cobertura sedimentaria se encuentra a 4 km de profundidad, a 21 km se estima la presencia de otra zona de cizalla subhorizontal dentro del basamento. Finalmente, a partir de la integración de todos los resultados obtenidos, se evaluó el potencial sismogénico de las fallas en la zona de estudio. Las fallas de primer orden que controlan la deformación en la zona son las responsables de los grandes terremotos. Mientras, las fallas Cuaternarias flexodeslizantes e inversas afectan solamente a la cobertura sedimentaria y serían estructuras de segundo orden que acomodan la deformación y fueron activadas durante el cuaternario con movimientos asísmicos y/o sísmicos de muy baja magnitud. Estos resultados permiten inferir que el corrimiento La Candelaria constituye una fuente sismogénica potencial de importancia para la región, donde se ubican numerosas poblaciones y obras civiles de envergadura. Por otra parte, la sección estructural balanceada implica la presencia de otras fallas ciegas de distinto orden de magnitud que podrían ser posibles fuentes sismogénicas profundas adicionales, marcando la necesidad de continuar con el desarrollo de este tipo de estudios en esta región tectónicamente activa.

Uplift in the broken Andean foreland of the Argentine Santa Bárbara System (SBS) is associated with the contractional reactivation of basement anisotropies, similar to those reported from the thick-skinned Cretaceous-Eocene Laramide province of North America. Fault scarps, deformed Quaternary deposits and landforms, disrupted drainage patterns, and medium-sized earthquakes within the SBS suggest that movement along these structures may be a recurring phenomenon, with yet to be defined repeat intervals and rupture lengths. In contrast to the Subandes thrust belt farther north, where eastward-migrating deformation has generated a well-defined thrust front, the SBS records spatiotemporally disparate deformation along structures that are only known to the first order. We present herein the results of geomorphic desktop analyses, structural field observations, and 2D electrical resistivity tomography and seismic-refraction tomography surveys and an interpretation of seismic reflection profiles across suspected fault scarps in the sedimentary basins adjacent to the Candelaria Range (CR) basement uplift, in the south-central part of the SBS. Our analysis in the CR piedmont areas reveals consistency between the results of near-surface electrical resistivity and seismic-refraction tomography surveys, the locations of prominent fault scarps, and structural geometries at greater depth imaged by seismic reflection data. We suggest that this deformation is driven by deep-seated blind thrusting beneath the CR and associated regional warping, while shortening involving Mesozoic and Cenozoic sedimentary strata in the adjacent basins was accommodated by layer-parallel folding and flexural-slip faults that cut through Quaternary landforms and deposits at the surface.

Ring current electrons (1–100 keV) have received significant attention in recent decades, but many questions regarding their major transport and loss mechanisms remain open. In this study, we use the four‐dimensional Versatile Electron Radiation Belt code to model the enhancement of phase space density that occurred during the 17 March 2013 storm. Our model includes global convection, radial diffusion, and scattering into the Earth's atmosphere driven by whistler‐mode hiss and chorus waves. We study the sensitivity of the model to the boundary conditions, global electric field, the electric field associated with subauroral polarization streams, electron loss rates, and radial diffusion coefficients. The results of the code are almost insensitive to the model parameters above 4.5 RERE, which indicates that the general dynamics of the electrons between 4.5 RE and the geostationary orbit can be explained by global convection. We found that the major discrepancies between the model and data can stem from the inaccurate electric field model and uncertainties in lifetimes. We show that additional mechanisms that are responsible for radial transport are required to explain the dynamics of ≥40‐keV electrons, and the inclusion of the radial diffusion rates that are typically assumed in radiation belt studies leads to a better agreement with the data. The overall effect of subauroral polarization streams on the electron phase space density profiles seems to be smaller than the uncertainties in other input parameters. This study is an initial step toward understanding the dynamics of these particles inside the geostationary orbit.

The novel space-borne Global Navigation Satellite System Reflectometry (GNSS-R) technique has recently shown promise in monitoring the ocean state and surface wind speed with high spatial coverage and unprecedented sampling rate. The L-band signals of GNSS are structurally able to provide a higher quality of observations from areas covered by dense clouds and under intense precipitation, compared to those signals at higher frequencies from conventional ocean scatterometers. As a result, studying the inner core of cyclones and improvement of severe weather forecasting and cyclone tracking have turned into the main objectives of GNSS-R satellite missions such as Cyclone Global Navigation Satellite System (CYGNSS). Nevertheless, the rain attenuation impact on GNSS-R wind speed products is not yet well documented. Evaluating the rain attenuation effects on this technique is significant since a small change in the GNSS-R can potentially cause a considerable bias in the resultant wind products at intense wind speeds. Based on both empirical evidence and theory, wind speed is inversely proportional to derived bistatic radar cross section with a natural logarithmic relation, which introduces high condition numbers (similar to ill-posed conditions) at the inversions to high wind speeds. This paper presents an evaluation of the rain signal attenuation impact on the bistatic radar cross section and the derived wind speed. This study is conducted simulating GNSS-R delay-Doppler maps at different rain rates and reflection geometries, considering that an empirical data analysis at extreme wind intensities and rain rates is impossible due to the insufficient number of observations from these severe conditions. Finally, the study demonstrates that at a wind speed of 30 m/s and incidence angle of 30 degrees, rain at rates of 10, 15, and 20 mm/h might cause overestimation as large as approximate to 0.65 m/s (2%), 1.00 m/s (3%), and 1.3 m/s (4%), respectively, which are still smaller than the CYGNSS required uncertainty threshold. The simulations are conducted in a pessimistic condition (severe continuous rainfall below the freezing height and over the entire glistening zone) and the bias is expected to be smaller in size in real environments.

Based on fossil organism remains including plant macrofossils, charcoal, pollen, and invertebrates preserved in syngenetic deposits of the Batagay permafrost sequence in the Siberian Yana Highlands, we reconstructed the environmental history during marine isotope stages (MIS) 6 to 2. Two fossil assemblages, exceptionally rich in plant remains, allowed for a detailed description of the palaeo-vegetation during two climate extremes of the Late Pleistocene, the onset of the last glacial maximum (LGM) and the last interglacial. In addition, altogether 41 assemblages were used to outline the vegetation history since the penultimate cold stage of MIS 6. Accordingly, meadow steppes analogue to modern communities of the phytosociological order Festucetalia lenensis formed the primary vegetation during the Saalian and Weichselian cold stages. Cold-resistant tundra-steppe communities (Carici rupestris-Kobresietea bellardii) as they occur above the treeline today were, in contrast to more northern locations, mostly lacking. During the last interglacial, open coniferous woodland similar to modern larch taiga was the primary vegetation at the site. Abundant charcoal indicates wildfire events during the last interglacial. Zoogenic disturbances of the local vegetation were indicated by the presence of ruderal plants, especially by abundant Urtica dioica, suggesting that the area was an interglacial refugium for large herbivores. Meadow steppes, which formed the primary vegetation during cold stages and provided potentially suitable pastures for herbivores, were a significant constituent of the plant cover in the Yana Highlands also under the full warm stage conditions of the last interglacial. Consequently, meadow steppes occurred in the Yana Highlands during the entire investigated timespan from MIS 6 to MIS 2 documenting a remarkable environmental stability. Thus, the proportion of meadow steppe vegetation merely shifted in response to the respectively prevailing climatic conditions. Their persistence indicates low precipitation and a relatively warm growing season throughout and beyond the late Pleistocene. The studied fossil record also proves that modern steppe occurrences in the Yana Highlands did not establish as late as in the Holocene but instead are relicts of a formerly continuous steppe belt extending from Central Siberia to Northeast Yakutia during the Pleistocene. The persistence of plants and invertebrates characteristic of meadow steppe vegetation in interior Yakutia throughout the late Quaternary indicates climatic continuity and documents the suitability of this region as a refugium also for other organisms of the Pleistocene mammoth steppe including the iconic large herbivores. (C)2018 Elsevier Ltd. All rights reserved.