Refine
Year of publication
Document Type
- Article (2826)
- Doctoral Thesis (508)
- Postprint (139)
- Other (76)
- Review (52)
- Monograph/Edited Volume (34)
- Preprint (17)
- Conference Proceeding (13)
- Habilitation Thesis (12)
- Master's Thesis (6)
Keywords
- climate change (52)
- Holocene (44)
- erosion (28)
- permafrost (27)
- Himalaya (26)
- remote sensing (24)
- Climate change (23)
- Tibetan Plateau (22)
- Earthquake source observations (21)
- Pollen (21)
Institute
- Institut für Geowissenschaften (3691) (remove)
Lamprophyre sind porphyrische, aus Mantelschmelzen gebildete Gesteine, die meist in Form von Gängen auftreten. Sie zeichnen sich durch auffällige und charakteristische texturelle, chemische und mineralogische Eigenschaften aus. Als ehemalige Mantelschmelzen liefern sie Information sowohl über Bedingungen der Schmelzbildung im Mantel als auch über geodynamische Prozesse, die zu metasomatischer Veränderung des Mantels geführt haben. Im Saxothuringikum Mitteleuropas, am Nordrand des Böhmischen Massivs, gibt es zahlreiche Lamprophyrvorkommen, die hier zur Charakterisierung der Mantelentwicklung während der variszischen Orogenese dienen. Die vorliegende Arbeit befaßt sich mit den mineralogischen, geochemischen und isotopischen (Sr-Nd-Pb) Signaturen von spätvariszischen kalkalkalischen Lamprophyren, von postvariszischen ultramafischen Lamprophyren, von Alkalibasalten der Lausitz und, zum Vergleich, von prävariszischen Gabbros. Darüberhinaus nutzt die Arbeit Lithium-Isotopensignaturen kombiniert mit Sr-Nd-Pb–Isotopendaten spätvariszischer kalkalkalischer Lamprophyre aus drei variszischen Domänen (Erzgebirge, Lausitz, Sudeten) zur Erkundung der lokalen Mantelüberprägungen während der variszischen Orogenese.
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.
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.
The mechanisms by which climate and vegetation affect erosion rates over various time scales lie at the heart of understanding landscape response to climate change. Plot-scale field experiments show that increased vegetation cover slows erosion, implying that faster erosion should occur under low to moderate vegetation cover. However, demonstrating this concept over long time scales and across landscapes has proven to be difficult, especially in settings complicated by tectonic forcing and variable slopes. We investigate this problem by measuring cosmogenic Be-10-derived catchment-mean denudation rates across a range of climate zones and hillslope gradients in the Kenya Rift, and by comparing our results with those published from the Rwenzori Mountains of Uganda. We find that denudation rates from sparsely vegetated parts of the Kenya Rift are up to 0.13 mm/yr, while those from humid and more densely vegetated parts of the Kenya Rift flanks and the Rwenzori Mountains reach a maximum of 0.08 mm/yr, despite higher median hillslope gradients. While differences in lithology and recent land-use changes likely affect the denudation rates and vegetation cover values in some of our studied catchments, hillslope gradient and vegetation cover appear to explain most of the variation in denudation rates across the study area. Our results support the idea that changing vegetation cover can contribute to complex erosional responses to climate or land-use change and that vegetation cover can play an important role in determining the steady-state slopes of mountain belts through its stabilizing effects on the land surface.
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.
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 subsurface harbors a large fraction of Earth's living biomass, forming complex microbial ecosystems. Without a profound knowledge of the ongoing biologically mediated processes and their reaction to anthropogenic changes it is difficult to assess the long-term stability and feasibility of any type of geotechnical utilization, as these influence subsurface ecosystems. Despite recent advances in many areas of subsurface microbiology, the direct quantification of turnover processes is still in its infancy, mainly due to the extremely low cell abundances. We provide an overview of the currently available techniques for the quantification of microbial turnover processes and discuss their specific strengths and limitations. Most techniques employed so far have focused on specific processes, e.g. sulfate reduction or methanogenesis. Recent studies show that processes that were previously thought to exclude each other can occur simultaneously, albeit at very low rates. Without the identification of the respective processes it is impossible to quantify total microbial activity. Even in cases where all simultaneously occurring processes can be identified, the typically very low rates prevent quantification. In many cases a simple measure of total microbial activity would be a better and more robust measure than assays for several specific processes. Enzyme or molecular assays provide a more general approach as they target key metabolic compounds. Depending on the compound targeted a broader spectrum of microbial processes can be quantified. The two most promising compounds are ATP and hydrogenase, as both are ubiquitous in microbes. Technical constraints limit the applicability of currently available ATP-assays for subsurface samples. A recently developed hydrogenase radiotracer assay has the potential to become a key tool for the quantification of subsurface microbial activity.
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.
Numerical simulation of fluid-flow processes in a 3D high-resolution carbonate reservoir analogue
(2014)
A high-resolution three-dimensional (3D) outcrop model of a Jurassic carbonate ramp was used in order to perform a series of detailed and systematic flow simulations. The aim of this study was to test the impact of small- and large-scale geological features on reservoir performance and oil recovery. The digital outcrop model contains a wide range of sedimentological, diagenetic and structural features, including discontinuity surfaces, shoal bodies, mud mounds, oyster bioherms and fractures. Flow simulations are performed for numerical well testing and secondary oil recovery. Numerical well testing enables synthetic but systematic pressure responses to be generated for different geological features observed in the outcrops. This allows us to assess and rank the relative impact of specific geological features on reservoir performance. The outcome documents that, owing to the realistic representation of matrix heterogeneity, most diagenetic and structural features cannot be linked to a unique pressure signature. Instead, reservoir performance is controlled by subseismic faults and oyster bioherms acting as thief zones. Numerical simulations of secondary recovery processes reveal strong channelling of fluid flow into high-permeability layers as the primary control for oil recovery. However, appropriate reservoir-engineering solutions, such as optimizing well placement and injection fluid, can reduce channelling and increase oil recovery.
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.
Im Graduiertenkolleg NatRiskChange der Universität Potsdam und anderen Forschungseinrichtungen werden beobachtete sowie zukünftig mögliche Veränderungen von Naturgefahren untersucht. Teil des strukturierten Doktorandenprogramms sind sogenannte Task-Force-Einsätze, bei denen die Promovierende zeitlich begrenzt ein aktuelles Ereignis auswerten. Im Zuge dieser Aktivität wurde die Sturzflut vom 29.05.2016 in Braunsbach (Baden-Württemberg) untersucht.
In diesem Bericht werden erste Auswertungen zur Einordnung der Niederschläge, zu den hydrologischen und geomorphologischen Prozessen im Einzugsgebiet des Orlacher Bachs sowie zu den verursachten Schäden beleuchtet.
Die Region war Zentrum extremer Regenfälle in der Größenordnung von 100 mm innerhalb von 2 Stunden. Das 6 km² kleine Einzugsgebiet hat eine sehr schnelle Reaktionszeit, zumal bei vorgesättigtem Boden. Im steilen Bachtal haben mehrere kleinere und größere Hangrutschungen über 8000 m³ Geröll, Schutt und Schwemmholz in das Gewässer eingetragen und möglicherweise kurzzeitige Aufstauungen und Durchbrüche verursacht. Neben den großen Wassermengen mit einer Abflussspitze in einer Größenordnung von 100 m³/s hat gerade die Geschiebefracht zu großen Schäden an den Gebäuden entlang des Bachlaufs in Braunsbach geführt.
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.
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.
The quantification of spatial propagation of extreme precipitation events is vital in water resources planning and disaster mitigation. However, quantifying these extreme events has always been challenging as many traditional methods are insufficient to capture the nonlinear interrelationships between extreme event time series. Therefore, it is crucial to develop suitable methods for analyzing the dynamics of extreme events over a river basin with a diverse climate and complicated topography. Over the last decade, complex network analysis emerged as a powerful tool to study the intricate spatiotemporal relationship between many variables in a compact way. In this study, we employ two nonlinear concepts of event synchronization and edit distance to investigate the extreme precipitation pattern in the Ganga river basin. We use the network degree to understand the spatial synchronization pattern of extreme rainfall and identify essential sites in the river basin with respect to potential prediction skills. The study also attempts to quantify the influence of precipitation seasonality and topography on extreme events. The findings of the study reveal that (1) the network degree is decreased in the southwest to northwest direction, (2) the timing of 50th percentile precipitation within a year influences the spatial distribution of degree, (3) the timing is inversely related to elevation, and (4) the lower elevation greatly influences connectivity of the sites. The study highlights that edit distance could be a promising alternative to analyze event-like data by incorporating event time and amplitude and constructing complex networks of climate extremes.
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.
The temporal dynamics of climate processes are spread across different timescales and, as such, the study of these processes at only one selected timescale might not reveal the complete mechanisms and interactions within and between the (sub-) processes. To capture the non-linear interactions between climatic events, the method of event synchronization has found increasing attention recently. The main drawback with the present estimation of event synchronization is its restriction to analysing the time series at one reference timescale only. The study of event synchronization at multiple scales would be of great interest to comprehend the dynamics of the investigated climate processes. In this paper, the wavelet-based multi-scale event synchronization (MSES) method is proposed by combining the wavelet transform and event synchronization. Wavelets are used extensively to comprehend multi-scale processes and the dynamics of processes across various timescales. The proposed method allows the study of spatio-temporal patterns across different timescales. The method is tested on synthetic and real-world time series in order to check its replicability and applicability. The results indicate that MSES is able to capture relationships that exist between processes at different timescales.
Quantifying the roles of single stations within homogeneous regions using complex network analysis
(2018)
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.
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.
In the last decade, the number and dimensions of catastrophic flooding events in the Niger River Basin (NRB) have markedly increased. Despite the devastating impact of the floods on the population and the mainly agriculturally based economy of the riverine nations, awareness of the hazards in policy and science is still low. The urgency of this topic and the existing research deficits are the motivation for the present dissertation.
The thesis is an initial detailed assessment of the increasing flood risk in the NRB. The research strategy is based on four questions regarding (1) features of the change in flood risk, (2) reasons for the change in the flood regime, (3) expected changes of the flood regime given climate and land use changes, and (4) recommendations from previous analysis for reducing the flood risk in the NRB.
The question examining the features of change in the flood regime is answered by means of statistical analysis. Trend, correlation, changepoint, and variance analyses show that, in addition to the factors exposure and vulnerability, the hazard itself has also increased significantly in the NRB, in accordance with the decadal climate pattern of West Africa. The northern arid and semi-arid parts of the NRB are those most affected by the changes.
As potential reasons for the increase in flood magnitudes, climate and land use changes are attributed by means of a hypothesis-testing framework. Two different approaches, based on either data analysis or simulation, lead to similar results, showing that the influence of climatic changes is generally larger compared to that of land use changes. Only in the dry areas of the NRB is the influence of land use changes comparable to that of climatic alterations.
Future changes of the flood regime are evaluated using modelling results. First ensembles of statistically and dynamically downscaled climate models based on different emission scenarios are analyzed. The models agree with a distinct increase in temperature. The precipitation signal, however, is not coherent. The climate scenarios are used to drive an eco-hydrological model. The influence of climatic changes on the flood regime is uncertain due to the unclear precipitation signal. Still, in general, higher flood peaks are expected. In a next step, effects of land use changes are integrated into the model. Different scenarios show that regreening might help to reduce flood peaks. In contrast, an expansion of agriculture might enhance the flood peaks in the NRB. Similarly to the analysis of observed changes in the flood regime, the impacts of climate- and land use changes for the future scenarios are also most severe in the dry areas of the NRB.
In order to answer the final research question, the results of the above analysis are integrated into a range of recommendations for science and policy on how to reduce flood risk in the NRB. The main recommendations include a stronger consideration of the enormous natural climate variability in the NRB and a focus on so called “no-regret” adaptation strategies which account for high uncertainty, as well as a stronger consideration of regional differences. Regarding the prevention and mitigation of catastrophic flooding, the most vulnerable and sensitive areas in the basin, the arid and semi-arid Sahelian and Sudano-Sahelian regions, should be prioritized. Eventually, an active, science-based and science-guided flood policy is recommended. The enormous population growth in the NRB in connection with the expected deterioration of environmental and climatic conditions is likely to enhance the region´s vulnerability to flooding. A smart and sustainable flood policy can help mitigate these negative impacts of flooding on the development of riverine societies in West Africa.
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.
Climate or Land Use?
(2015)
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.
Climate or land use?
(2017)
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.
Sediment-discharge hysteresis loops are frequently analyzed to facilitate the understanding of sediment transport processes. Hysteresis patterns, however, are often complex and their interpretation can be complicated. Particularly, quantifying hysteresis patterns remains a problematic issue. Moreover, it is currently unknown how much data is required for analyzing sediment-discharge hysteresis loops in a given area. These open questions and challenges motivated us to develop a new method for quantifying suspended-sediment hysteresis. Subsequently, we applied the new hysteresis index to three suspended-sediment and discharge datasets from a small tropical rainforest catchment. The datasets comprised a different number of events and sampling sites. Our analyses show three main findings: (1) datasets restricted to only few events, which is typical for rapid assessment surveys, were always sufficient to identify the dominating hysteresis pattern in our research area. Furthermore, some of these small datasets contained multiple-peak events that allowed identifying intra-event exhaustion effects and hence, limitations in sediment supply. (2) Datasets comprising complete hydrological years were particularly useful for analyzing seasonal dynamics of hysteresis. These analyses revealed an exhaustion of hysteresis on the inter-event scale which also points to a limited sediment supply. (3) Datasets comprising measurements from two consecutive gauges installed at the catchment outlet and on a slope within that catchment allowed analyzing the change of hysteresis patterns along the flowpath. On the slope, multiple-peak events showed a stronger intra-event exhaustion of hysteresis than at the catchment outlet. Furthermore, exhaustion of hysteresis on the inter-event scale was not evident on the slope but occurred at the catchment outlet. Our results indicate that even small sediment datasets can provide valuable insights into sediment transport processes of small catchments. Furthermore, our results may serve as a first guideline on what to expect from an analysis of hysteresis patterns for datasets of varying quality and quantity. (c) 2014 Elsevier B.V. All rights reserved.
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.
Samples of 474 forest stands in Germany were analysed for concentrations of polycyclic aromatic hydrocarbons (PAHs) in three sampling depths. Enhanced concentrations were mainly found at spots relatively close to densely industrialized and urbanized regions and at some topographically elevated areas. Average enrichment factors between mineral soil and humic layer depend on humus type i.e. decrease from mull via moder to more Based on their compound-patterns, the observed samples could be assigned to three main clusters. For some parts of our study area a uniform assignment of samples to clusters over larger regions could be identified. For instance, samples taken at vicinity to brown-coal strip-mining districts are characterized by high relative abundances of low-molecular-weight PAHs. These results suggest that PAHs are more likely originated from local and regional emitters rather than from long-range transport and that specific source-regions can be identified based on PAH fingerprints. (C) 2015 Elsevier Ltd. All rights reserved.
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.
Organic geochemical proxy data from surface sediment samples and a sediment core from Lake Donggi Cona were used to infer environmental changes on the northeastern Tibetan Plateau spanning the last 18.4 kyr. Long-chain n-alkanes dominate the aliphatic hydrocarbon fraction of the sediment extract from most surface sediment samples and the sediment core. Unsaturated mid-chain n-alkanes (nC(23:1) and nC(25:1)) have high abundances in some samples, especially in core samples from the late glacial and early Holocene. TOC contents, organic biomarker and non-pollen-palynomorph concentrations and results from organic petrologic analysis on selected samples suggest three major episodes in the history of Lake Donggi Cona. Before ca. 12.6 cal ka BP samples contain low amounts of organic matter due to cold and arid conditions during the late glacial. After 12.6 cal ka BP, relatively high contents of TOC and concentrations of Botryococcus fossils, as well as enhanced concentrations of mid-chain n-alkanes and n-alkenes suggest a higher primary and macrophyte productivity than at present This is supported by high contents of palynomorphs derived from higher plants and algae and was possibly triggered by a decrease of salinity and amelioration of climate during the early Holocene. Since 6.8 cal ka BP Lake Donggi Cona has been an oligotrophic freshwater lake. Proxy data suggest that variations in insolation drive ecological changes in the lake, with increased aquatic productivity during the early Holocene summer insolation maximum. Short-term drops of TOC contents or biomarker concentrations (at 9.9 cal ka BP, after 8.0 and between 3.5 and 1.7 cal ka BP) can possibly be related to relatively cool and dry episodes reported from other sites on the north-eastern Tibetan Plateau, which are hypothesized to occur in phase with Northern Hemisphere cooling events.
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.
The Central Asian Pamir Mountains (Pamirs) are a high-altitude region sensitive to climatic change, with only few paleoclimatic records available. To examine the glacial-interglacial hydrological changes in the region, we analyzed the geochemical parameters of a 31-kyr record from Lake Karakul and performed a set of experiments with climate models to interpret the results. delta D values of terrestrial biomarkers showed insolation-driven trends reflecting major shifts of water vapor sources. For aquatic biomarkers, positive delta D shifts driven by changes in precipitation seasonality were observed at ca. 31-30, 28-26, and 17-14 kyr BP. Multiproxy paleoecological data and modelling results suggest that increased water availability, induced by decreased summer evaporation, triggered higher lake levels during those episodes, possibly synchronous to northern hemispheric rapid climate events. We conclude that seasonal changes in precipitation-evaporation balance significantly influenced the hydrological state of a large waterbody such as Lake Karakul, while annual precipitation amount and inflows remained fairly constant.
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.
The Karaburun Peninsula, which is considered part of the Anatolide-Tauride Block of Turkey, contains clastic and carbonate sequences deposited on the northern margin of Gondwana. The Palaeozoic clastic sequence, which is intruded by the Early Triassic granitoid and tectonically overlies a Mesozoic melange sequence, can be divided into three subunits: a lower clastic subunit consisting of a sandstone-shale alternation, an upper clastic subunit consisting of black chert-bearing shales, sandstone and conglomerate, and a Permo-Carboniferous carbonate subunit. The lower Triassic Karaburun I-type granitoid has a high initial Sr-87/Sr-86 ratio (0.709021-0.709168), and low Nd-143/Nd-144 ratio (0.512004-0.512023) and epsilon Nd (-5.34 to -5.70) isotopic values. Geochronological data indicate a crystallization (intrusion) age of 247.1 +/- 2.0 Ma (Scythian). Geochemically, the acidic magmatism reflects a subduction-related continental-arc basin tectonic setting, which can be linked to the opening of the northern branch of Neo-Tethys as a continental back-arc rifting basin on the northern margin of Gondwana. This can be related to the closure through southward subduction of the Palaeotethys Ocean beneath Gondwana.
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.
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 2-D distinct element method (DEM) code (PFC2D_V5) is used here to simulate the evolution of subsidence-related karst landforms, such as single and clustered sinkholes, and associated larger-scale depressions. Subsurface material in the DEM model is removed progressively to produce an array of cavities; this simulates a network of subsurface groundwater conduits growing by chemical/mechanical erosion. The growth of the cavity array is coupled mechanically to the gravitationally loaded surroundings, such that cavities can grow also in part by material failure at their margins, which in the limit can produce individual collapse sinkholes. Two end-member growth scenarios of the cavity array and their impact on surface subsidence were examined in the models: (1) cavity growth at the same depth level and growth rate; (2) cavity growth at progressively deepening levels with varying growth rates. These growth scenarios are characterised by differing stress patterns across the cavity array and its overburden, which are in turn an important factor for the formation of sinkholes and uvalalike depressions. For growth scenario (1), a stable compression arch is established around the entire cavity array, hindering sinkhole collapse into individual cavities and favouring block-wise, relatively even subsidence across the whole cavity array. In contrast, for growth scenario (2), the stress system is more heterogeneous, such that local stress concentrations exist around individual cavities, leading to stress interactions and local wall/overburden fractures. Consequently, sinkhole collapses occur in individual cavities, which results in uneven, differential subsidence within a larger-scale depression. Depending on material properties of the cavity-hosting material and the overburden, the larger-scale depression forms either by sinkhole coalescence or by widespread subsidence linked geometrically to the entire cavity array. The results from models with growth scenario (2) are in close agreement with surface morphological and subsurface geophysical observations from an evaporite karst area on the eastern shore of the Dead Sea.
Evaluation of groundwater dynamics and quality in the Najd aquifers located in the Sultanate of Oman
(2012)
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.
Extracellular DNA (eDNA) is a ubiquitous biological compound in aquatic sediment and soil. Previous studies suggested that eDNA plays an important role in biogeochemical element cycling, horizontal gene transfer and stabilization of biofilm structures. Previous methods for eDNA extraction were either not suitable for oligotrophic sediments or only allowed quantification but no genetic analyses. Our procedure is based on cell detachment and eDNA liberation from sediment particles by sequential washing with an alkaline sodium phosphate buffer followed by a separation of cells and eDNA. The separated eDNA is then bound onto silica particles and purified, whereas the intracellular DNA from the separated cells is extracted using a commercial kit. The method provides extra- and intracellular DNA of high purity that is suitable for downstream applications like PCR. Extracellular DNA was extracted from organic-rich shallow sediment of the Baltic Sea, glacially influenced sediment of the Barents Sea and from the oligotrophic South Pacific Gyre. The eDNA concentration in these samples varied from 23 to 626 ng g(-1) wet weight sediment. A number of experiments were performed to verify each processing step. Although extraction efficiency is higher than other published methods, it is not fully quantitative. (C) 2014 Elsevier B.V. All rights reserved.
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.
Glacial-cycle simulations of the Antarctic Ice Sheet with the Parallel Ice Sheet Model (PISM)
(2020)
Simulations of the glacial-interglacial history of the Antarctic Ice Sheet provide insights into dynamic threshold behavior and estimates of the ice sheet's contributions to global sea-level changes for the past, present and future. However, boundary conditions are weakly constrained, in particular at the interface of the ice sheet and the bedrock. Also climatic forcing covering the last glacial cycles is uncertain, as it is based on sparse proxy data. <br /> We use the Parallel Ice Sheet Model (PISM) to investigate the dynamic effects of different choices of input data, e.g., for modern basal heat flux or reconstructions of past changes of sea level and surface temperature. As computational resources are limited, glacial-cycle simulations are performed using a comparably coarse model grid of 16 km and various parameterizations, e.g., for basal sliding, iceberg calving, or for past variations in precipitation and ocean temperatures. In this study we evaluate the model's transient sensitivity to corresponding parameter choices and to different boundary conditions over the last two glacial cycles and provide estimates of involved uncertainties. We also discuss isolated and combined effects of climate and sea-level forcing. Hence, this study serves as a "cookbook" for the growing community of PISM users and paleo-ice sheet modelers in general. <br /> For each of the different model uncertainties with regard to climatic forcing, ice and Earth dynamics, and basal processes, we select one representative model parameter that captures relevant uncertainties and motivates corresponding parameter ranges that bound the observed ice volume at present. The four selected parameters are systematically varied in a parameter ensemble analysis, which is described in a companion paper.
The biostratigraphy of Campanian-Maastrichtian carbonate platforms is largely based on the larger foraminiferal genus Orbitoides. However, while the taxonomy and the chronostratigraphic age of the younger species of this genus are well established, there are still many controversies on the earliest species. We have restudied their morphological characters using a large collection of samples from the type-localities and from continuous sections in the southern Pyrenees. Based on these new observations, the long forgotten species O. sanctae-pelagiae is reinstated, while O. dordoniensis is considered a junior synonym. Successive populations of O. hottingeri, O. sanctae-pelagiae and O. douvillei show gradual morphological changes in time marked by an increase in the size and complexity of the macrospheric embryonal apparatus, an increase of the size of the adult specimens of both generations and the progressive appearance and development of true lateral chamberlets. The Font de les Bagasses Unit in the southern Pyrenees preserves a high-resolution archive of the evolution of the earliest Orbitoides. Strontium isotope stratigraphy indicates that the oldest species, O. hottingeri, made its first appearance in the earliest Campanian, close to the Santonian-Campanian boundary, and was replaced by O. sanctae-pelagiae at a level closely corresponding to the boundary between the Placenticeras bidorsatum and Menabites delawarensis ammonite zones. (C) 2014 Elsevier Ltd. All rights reserved.
Successful delivery of the United Nations sustainable development goals and implementation of the Paris Agreement requires technologies that utilize a wide range of minerals in vast quantities. Metal recycling and technological change will contribute to sustaining supply, but mining must continue and grow for the foreseeable future to ensure that such minerals remain available to industry. New links are needed between existing institutional frameworks to oversee responsible sourcing of minerals, trajectories for mineral exploration, environmental practices, and consumer awareness of the effects of consumption. Here we present, through analysis of a comprehensive set of data and demand forecasts, an interdisciplinary perspective on how best to ensure ecologically viable continuity of global mineral supply over the coming decades.
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.