Refine
Year of publication
- 2015 (256) (remove)
Document Type
- Article (203)
- Doctoral Thesis (33)
- Postprint (4)
- Review (4)
- Habilitation Thesis (3)
- Other (3)
- Monograph/Edited Volume (2)
- Preprint (2)
- Bachelor Thesis (1)
- Conference Proceeding (1)
Is part of the Bibliography
- yes (256)
Keywords
- Holocene (6)
- climate change (4)
- erosion (4)
- India-Asia collision (3)
- Ostracoda (3)
- Pamir (3)
- Seismicity and tectonics (3)
- ecohydrology (3)
- landscape evolution (3)
- Cerrado (2)
Institute
- Institut für Geowissenschaften (256) (remove)
Several destructive earthquakes have occurred in Tien-Shan region at the beginning of 20th century. However, the detailed seismological characteristics, especially source parameters of those earthquakes are still poorly investigated. The Chon-Kemin earthquake is the strongest instrumentally recorded earthquake in the Tien-Shan region. This earthquake has produced an approximately 200 km long system of surface ruptures along Kemin-Chilik fault zone and killed about similar to 400 people. Several studies presented the different information on the earthquake epicentre location and magnitude, and two different focal mechanisms were also published. The reason for the limited knowledge of the source parameters for the Chon-Kemin earthquake is the complexity of old analogue records processing, digitization and analysis. In this study the data from 23 seismic stations worldwide were collected and digitized. The earthquake epicentre was relocated to 42.996NA degrees and 77.367EA degrees, the hypocentre depth is estimated between 10 and 20 km. The magnitude was recalculated to m(B) 8.05, M-s 7.94 and M-w 8.02. The focal mechanism, determined from amplitude ratios comparison of the observed and synthetic seismograms, was: str = 264A degrees, dip = 52A degrees, rake = 98A degrees. The apparent source time duration was between similar to 45 and similar to 70 s, the maximum slip occurred 25 s after the beginning of the rupture. Two subevents were clearly detected from the waveforms with the scalar moment ratio between them of about 1/3, the third subevent was also detected with less certainty. Taking into account surface rupture information, the fault geometry model with three patches was proposed. Based on scaling relations we conclude that the total rupture length was between similar to 260 and 300 km and a maximum rupture width could reach similar to 70 km.
The Tien-Shan and the neighboring Pamir region are two of the largest mountain belts in the world. Their deformation is dominated by intermontane basins bounded by active thrust and reverse faulting. The Tien-Shan mountain belt is characterized by a very high rate of seismicity along its margins as well as within the Tien-Shan interior. The study area of the here presented thesis, the western part of the Tien-Shan region, is currently seismically active with small and moderate sized earthquakes. However, at the end of the 19th beginning of the 20th century, this region was struck by a remarkable series of large magnitude (M>7) earthquakes, two of them reached magnitude 8.
Those large earthquakes occurred prior to the installation of the global digital seismic network and therefore were recorded only by analog seismic instruments. The processing of the analog data brings several difficulties, for example, not always the true parameters of the recording system are known. Another complicated task is the digitization of those records - a very time-consuming and delicate part. Therefore a special set of techniques is developed and modern methods are adapted for the digitized instrumental data analysis.
The main goal of the presented thesis is to evaluate the impact of large magnitude M≥7.0 earthquakes, which occurred at the turn of 19th to 20th century in the Tien-Shan region, on the overall regional tectonics. A further objective is to investigate the accuracy of previously estimated source parameters for those earthquakes, which were mainly based on macroseismic observations, and re-estimate them based on the instrumental data. An additional aim of this study is to develop the tools and methods for faster and more productive usage of analog seismic data in modern seismology.
In this thesis, the ten strongest and most interesting historical earthquakes in Tien-Shan region are analyzed. The methods and tool for digitizing and processing the analog seismic data are presented. The source parameters of the two major M≥8.0 earthquakes in the Northern Tien-Shan are re-estimated in individual case studies. Those studies are published as peer-reviewed scientific articles in reputed journals. Additionally, the Sarez-Pamir earthquake and its connection with one of the largest landslides in the world, Usoy landslide, is investigated by seismic modeling. These results are also published as a research paper.
With the developed techniques, the source parameters of seven more major earthquakes in the region are determined and their impact on the regional tectonics was investigated. The large magnitudes of those earthquakes are confirmed by instrumental data. The focal mechanism of these earthquakes were determined providing evidence for responsible faults or fault systems.
Shaking water out of soil
(2015)
Moderate to large earthquakes can increase the amount of water flowing in streams. Previous interpretations and models assume that the extra water originates in the saturated zone. Here we show that earthquakes may also release water from the unsaturated zone when the seismic energy is sufficient to overcome the threshold of soil water retention. Soil water may then be released into aquifers, increasing streamflow. After the M8.8 Maule, Chile, earthquake, the discharge in some headwater catchments of the Chilean coastal range increased, and the amount of extra water in the discharge was similar to the total amount of water available for release from the unsaturated zone. Assuming rapid recharge of this water to the water table, a groundwater flow model that accounts for evapotranspiration and water released from soils can reproduce the increase in discharge as well as the enhanced diurnal discharge variations observed after the earthquake. Thus the unsaturated zone may play a previously unappreciated, and potentially significant, role in shallow hydrological responses to earthquakes.
Probabilistic seismic-hazard analysis (PSHA) is the current tool of the trade used to estimate the future seismic demands at a site of interest. A modern PSHA represents a complex framework that combines different models with numerous inputs. It is important to understand and assess the impact of these inputs on the model output in a quantitative way. Sensitivity analysis is a valuable tool for quantifying changes of a model output as inputs are perturbed, identifying critical input parameters, and obtaining insight about the model behavior. Differential sensitivity analysis relies on calculating first-order partial derivatives of the model output with respect to its inputs; however, obtaining the derivatives of complex models can be challenging.
In this study, we show how differential sensitivity analysis of a complex framework such as PSHA can be carried out using algorithmic/automatic differentiation (AD). AD has already been successfully applied for sensitivity analyses in various domains such as oceanography and aerodynamics. First, we demonstrate the feasibility of the AD methodology by comparing AD-derived sensitivities with analytically derived sensitivities for a basic case of PSHA using a simple ground-motion prediction equation. Second, we derive sensitivities via AD for a more complex PSHA study using a stochastic simulation approach for the prediction of ground motions. The presented approach is general enough to accommodate more advanced PSHA studies of greater complexity.
Semi-automated detection of annual laminae (varves) in lake sediments using a fuzzy logic algorithm
(2015)
Annual laminae (varves) in lake sediments are typically visually identified, measured and counted, although numerous attempts have been made to automate this process. The reason for the failure of most of these automated algorithms for varve counting is the complexity of the seasonal laminations, typically rich in lateral fades variations and internal heterogeneities. In the manual counting of varves, the investigator acquired and interpreted flexible numbers of complex decision criteria to understand whether a particular simple lamination is a varve or not. Fuzzy systems simulate the flexible decision making process in a computer by introducing a smooth transition between true varve and false varve. In our investigation, we use an adaptive neuro fuzzy inference system (ANFIS) to detect varves on the basis of a digital image of the sediment. The results of the application of the ANFIS to laminated sediments from the Meerfelder Maar (Eifel, Germany) and from a landslide-dammed lake in the Quebrada de Cafayate of Argentina are compared with manual varve counts and possible reasons for the differences are discussed. (C) 2015 Elsevier B.V. All rights reserved.
Studies of seismic tomography have been highly successful at imaging the deep structure of subduction zones. In a study complementary to these tomographic studies, we use array seismology and reflected waves to image a stagnant slab in the mantle transition zone. Using P and S (SH) waves we find a steeply dipping reflector centred at ca. 400 km depth and ca. 550 km west of the present Mariana subduction zone (at 20N, 140E). The discovery of this anomaly in tomography and independently in array seismology (this paper) helps in understanding the evolution of the Mariana margin. The reflector/stagnant slab may be the remains of the hypothetical North New Guinea Plate, which was theorized to have subducted ca. 50 Ma.
The segmentation of major fault systems in subduction zones controls earthquake magnitude and location, but the causes for the existence of segment boundaries and the relationships between long-term deformation and the extent of earthquake rupture, are poorly understood. We compare permanent and seismic-cycle deformation patterns along the rupture zone of the 2010 Maule earthquake (M8.8), which ruptured 500 km of the Chile subduction margin. We analyzed the morphology of MIS-5 marine terraces using LiDAR topography and established their chronology and coeval origin with twelve luminescence ages, stratigraphy and geomorphic correlation, obtaining a virtually continuous distribution of uplift rates along the entire rupture zone. The mean uplift rate for these terraces is 0.5 m/ka. This value is exceeded in three areas, which have experienced rapid emergence of up to 1.6 m/ka; they are located at the northern, central, and southern sectors of the rupture zone, referred to as Topocalma, Carranza and Arauco, respectively. The three sectors correlate with boundaries of eight great earthquakes dating back to 1730. The Topocalma and Arauco sectors, located at the boundaries of the 2010 rupture, consist of broad zones of crustal warping with wavelengths of 60 and 90 km, respectively. These two regions coincide with the axes of oroclinal bending of the entire Andean margin and correlate with changes in curvature of the plate interface. Rapid uplift at Carranza, in turn, is of shorter wavelength and associated with footwall flexure of three crustal-scale normal faults. The uplift rate at Carranza is inversely correlated with plate coupling as well as with coseismic slip, suggesting permanent deformation may accumulate interseismically. We propose that the zones of upwarping at Arauco and Topocalma reflect changes in frictional properties of the megathrust resulting in barriers to the propagation of great earthquakes. Slip during the 1960 (M9.5) and 2010 events overlapped with the similar to 90-km-long zone of rapid uplift at Arauco; similarly, slip in 2010 and 1906 extended across the similar to 60-km-long section of the megathrust at Topocalma, but this area was completely breached by the 1730 (M similar to 9) event, which propagated southward until Carranza. Both Arauco and Topocalma show evidence of sustained rapid uplift since at least the middle Pleistocene. These two sectors might thus constitute discrete seismotectonic boundaries restraining most, but not all great earthquake ruptures. Based on our observations, such barriers might be breached during multi-segment super-cycle events. (C) 2015 Elsevier Ltd. All rights reserved.
Roads at risk
(2015)
Globalisation and interregional exchange of people, goods, and services has boosted the importance of and reliance on all kinds of transport networks. The linear structure of road networks is especially sensitive to natural hazards. In southern Norway, steep topography and extreme weather events promote frequent traffic disruption caused by debris flows. Topographic susceptibility and trigger frequency maps serve as input into a hazard appraisal at the scale of first-order catchments to quantify the impact of debris flows on the road network in terms of a failure likelihood of each link connecting two network vertices, e.g. road junctions. We compute total additional traffic loads as a function of traffic volume and excess distance, i.e. the extra length of an alternative path connecting two previously disrupted network vertices using a shortest-path algorithm. Our risk metric of link failure is the total additional annual traffic load, expressed as vehicle kilometres, because of debris-flow-related road closures. We present two scenarios demonstrating the impact of debris flows on the road network and quantify the associated path-failure likelihood between major cities in southern Norway. The scenarios indicate that major routes crossing the central and north-western part of the study area are associated with high link-failure risk. Yet options for detours on major routes are manifold and incur only little additional costs provided that drivers are sufficiently well informed about road closures. Our risk estimates may be of importance to road network managers and transport companies relying on speedy delivery of services and goods.
Interlayered blueschists and greenschists of the Coastal Cordillera (Chile) are part of a Late Palaeozoic accretionary complex. They represent metavolcanic rocks with oceanic affinities based on predominantly 01B-type REE patterns and immobile trace element ratios. Both rock types have similar mineralogies, albeit with different mineral modal abundances. Amphibole is the major mafic mineral and varies compositionally from glaucophane to actinolite. The presence of glaucophane relicts as cores in zoned amphiboles in both blueschists and greenschists is evidence for a pervasive high-pressure metamorphic stage, indicating that tectonic juxtaposition is an unlikely explanation for the cm-dm scale interlayering. During exhumation, a retrograde greenschist-facies overprint stabilized chlorite + albite + winchitic/actinolitic amphibole + phengitic white mica +/- epidote +/- K-feldspar at 0.4 +/- 0.1 GPa. Geochemical variability can be partly ascribed to primary magmatic and partly to secondary metasomatic processes that occurred under greenschist-facies conditions. Isocon diagrams of several adjacent blueschist-greenschist pairs with similar protolith geochemistry were used to evaluate metasomatic changes due to retrograde fluid-rock interaction. The most important geochemical changes are depletion of Si and Na and addition of water in the greenschists compared to the blueschists. Transition metals and LILE are mobilized to varying degrees. The unsystematic deviations from magmatic fractionation trends suggest open system conditions and influx of an external fluid. Pseudosection and water isopleth calculations show that the rocks were dehydrating during most of their exhumation history and remained at water-saturated conditions. The mineralogical changes, in particular breakdown of blue amphibole and replacement by chlorite, albite and calcic/sodic-calcic amphibole, are the prime cause for the distinct coloring. Pseudo-binary phase diagrams were used as a means to link bulk rock geochemical variability to modal and chemical changes in the mineralogy. The geochemical changes induced by fluid-rock interaction are important in two ways: First, the bulk rock chemistry is altered, leading to the stabilization of higher modal proportions of chlorite in the greenschists. Second, the retrograde overprint is a selective, layer-parallel fluid infiltration process, causing more intense greenschist-facies recrystallization in greenschist layers and therefore preferential preservation of blue amphibole in blueschist layers. Hence, the distinct colors were acquired by a combination of compositional variability, both primary magmatic and secondary metasomatic, and the different intensity of retrograde fluid infiltration. (C) 2014 Elsevier B.V. All rights reserved.
Vegetated dunes are recognized as important natural barriers that shelter inland ecosystems and coastlines suffering daily erosive impacts of the sea and extreme events, such as tsunamis. However, societal responses to erosion and shoreline retreat often result in man-made coastal defence structures that cover part of the intertidal and upper shore zones causing coastal squeeze and habitat loss, especially for upper shore biota, such as dune plants. Coseismic uplift of up to 2.0 m on the Peninsula de Arauco (South central Chile, ca. 37.5 degrees S) caused by the 2010 Maule earthquake drastically modified the coastal landscape, including major increases in the width of uplifted beaches and the immediate conversion of mid to low sandy intertidal habitat to supralittoral sandy habitat above the reach of average tides and waves. To investigate the early stage responses in species richness, cover and across-shore distribution of the hitherto absent dune plants, we surveyed two formerly intertidal armoured sites and a nearby intertidal unarmoured site on a sandy beach located on the uplifted coast of Llico (Peninsula de Arauco) over two years. Almost 2 years after the 2010 earthquake, dune plants began to recruit, then rapidly grew and produced dune hummocks in the new upper beach habitats created by uplift at the three sites. Initial vegetation responses were very similar among sites. However, over the course of the study, the emerging vegetated dunes of the armoured sites suffered a slowdown in the development of the spatial distribution process, and remained impoverished in species richness and cover compared to the unarmoured site. Our results suggest that when released from the effects of coastal squeeze, vegetated dunes can recover without restoration actions. However, subsequent human activities and management of newly created beach and dune habitats can significantly alter the trajectory of vegetated dune development. Management that integrates the effects of natural and human induced disturbances, and promotes the development of dune vegetation as natural barriers can provide societal and conservation benefits in coastal ecosystems.
Pollen productivity estimates (PPE) are used to quantitatively reconstruct variations in vegetation within a specific distance of the sampled pollen archive. Here, for the first time, PPEs from Siberia are presented. The study area (Khatanga region, Krasnoyarsk territory, Russia) is located in the Siberian Sub-arctic where Larix is the sole forest-line forming tree taxon. Pollen spectra from two different sedimentary environments, namely terrestrial mosses (n = 16) and lakes (n = 15, median radius similar to 100 m) and their surrounding vegetation were investigated to extract PPEs. Our results indicate some differences in pollen spectra between moss and lake pollen. Larix and Cyperaceae for example obtained higher representation in the lacustrine than in terrestrial moss samples. This highlights that in calibration studies, modem and fossil datasets should use archives of similar sedimentary origin. Results of an Extended R-Value model were applied to assess the relevant source area of pollen (RSAP) and to calculate the PPEs for both datasets. As expected, the RSAP of the moss samples was very small (about 10 m) compared to the lacustrine samples (about 25 km). Calculation of PPEs for the six most common taxa yielded generally similar results for both datasets. Relative to Poaceae (reference taxon, PPE = 1) Betula nana-type (PPEmoss: 1.8, PPElake: 1.8) and Alnus fruticosa-type (PPEmoss:6.4, PPElake:2.9) were overrepresented while Cyperaceae (PPEmoss:0.5, PPElake:0.1), Ericaceae (PPEmoss: 0.3, PPElake <0.01), Salix (PPEmoss:0.03, PPElake <00.1) and Larix (PPEmoss: <0.01, PPElake:0.2) were under-represented in the pollen spectra compared to the vegetation in the RSAP. The estimation for the dominant tree in the region, Larix gmelinii, is the first published result for this species, but needs to be considered very preliminary. The inferred sequence from over- to under-representation is mostly consistent with results from Europe; however, the absolute values show some differences. Gathering vegetation data was limited by the remoteness of our study area and a lack of high-resolute satellite imagery and vegetation maps. Our estimate may serve as a first reference to strengthen future vegetation reconstructions in this climate-sensitive region. (C) 2015 Elsevier B.V. All rights reserved.
It has been known for many decades that the lunar tidal influence in the equatorial electrojet (EEJ) is noticeably enhanced during Northern Hemisphere winters. Recent literature has discussed the role of stratospheric sudden warming (SSW) events behind the enhancement of lunar tides and the findings suggest a positive correlation between the lunar tidal amplitude and lower stratospheric parameters (zonal mean air temperature and zonal mean zonal wind) during SSW events. The positive correlation raises the question whether an inverse approach could also be developed which makes it possible to deduce the occurrence of SSW events before their direct observations (before 1952) from the amplitude of the lunar tides. This study presents an analysis technique based on the phase of the semi-monthly lunar tide to determine the lunar tidal modulation of the EEJ. A statistical approach using the superposed epoch analysis is also carried out to formulate a relation between the EEJ tidal amplitude and lower stratospheric parameters. Using these results, we have estimated a threshold value for the tidal wave power that could be used to identify years with SSW events from magnetic field observations.
Adjustment of empirically derived ground motion prediction equations (GMPEs), from a data- rich region/site where they have been derived to a data-poor region/site, is one of the major challenges associated with the current practice of seismic hazard analysis. Due to the fre- quent use in engineering design practices the GMPEs are often derived for response spectral ordinates (e.g., spectral acceleration) of a single degree of freedom (SDOF) oscillator. The functional forms of such GMPEs are based upon the concepts borrowed from the Fourier spectral representation of ground motion. This assumption regarding the validity of Fourier spectral concepts in the response spectral domain can lead to consequences which cannot be explained physically.
In this thesis, firstly results from an investigation that explores the relationship between Fourier and response spectra, and implications of this relationship on the adjustment issues of GMPEs, are presented. The relationship between the Fourier and response spectra is explored by using random vibration theory (RVT), a framework that has been extensively used in earthquake engineering, for instance within the stochastic simulation framework and in the site response analysis. For a 5% damped SDOF oscillator the RVT perspective of response spectra reveals that no one-to-one correspondence exists between Fourier and response spectral ordinates except in a limited range (i.e., below the peak of the response spectra) of oscillator frequencies. The high oscillator frequency response spectral ordinates are dominated by the contributions from the Fourier spectral ordinates that correspond to the frequencies well below a selected oscillator frequency. The peak ground acceleration (PGA) is found to be related with the integral over the entire Fourier spectrum of ground motion which is in contrast to the popularly held perception that PGA is a high-frequency phenomenon of ground motion.
This thesis presents a new perspective for developing a response spectral GMPE that takes the relationship between Fourier and response spectra into account. Essentially, this frame- work involves a two-step method for deriving a response spectral GMPE: in the first step two empirical models for the FAS and for a predetermined estimate of duration of ground motion are derived, in the next step, predictions from the two models are combined within the same RVT framework to obtain the response spectral ordinates. In addition to that, a stochastic model based scheme for extrapolating the individual acceleration spectra beyond the useable frequency limits is also presented. To that end, recorded acceleration traces were inverted to obtain the stochastic model parameters that allow making consistent extrapola- tion in individual (acceleration) Fourier spectra. Moreover an empirical model, for a dura- tion measure that is consistent within the RVT framework, is derived. As a next step, an oscillator-frequency-dependent empirical duration model is derived that allows obtaining the most reliable estimates of response spectral ordinates. The framework of deriving the response spectral GMPE presented herein becomes a self-adjusting model with the inclusion of stress parameter (∆σ) and kappa (κ0) as the predictor variables in the two empirical models. The entire analysis of developing the response spectral GMPE is performed on recently compiled RESORCE-2012 database that contains recordings made from Europe, the Mediterranean and the Middle East. The presented GMPE for response spectral ordinates should be considered valid in the magnitude range of 4 ≤ MW ≤ 7.6 at distances ≤ 200 km.
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.
Reconstructing climate from the Dead Sea sediment record using high-resolution micro-facies analyses
(2015)
The sedimentary record of the Dead Sea is a key archive for reconstructing climate in the eastern Mediterranean region, as it stores the environmental and tectonic history of the Levant for the entire Quaternary. Moreover, the lake is located at the boundary between Mediterranean sub-humid to semi-arid and Saharo-Arabian hyper-arid climates, so that even small shifts in atmospheric circulation are sensitively recorded in the sediments. This DFG-funded doctoral project was carried out within the ICDP Dead Sea Deep Drilling Project (DSDDP) that intended to gain the first long, continuous and high-resolution sediment core from the deep Dead Sea basin. The drilling campaign was performed in winter 2010-11 and more than 700 m of sediments were recovered. The main aim of this thesis was (1) to establish the lithostratigraphic framework for the ~455 m long sediment core from the deep Dead Sea basin and (2) to apply high-resolution micro-facies analyses for reconstructing and better understanding climate variability from the Dead Sea sediments.
Addressing the first aim, the sedimentary facies of the ~455 m long deep-basin core 5017-1 were described in great detail and characterised through continuous overview-XRF element scanning and magnetic susceptibility measurements. Three facies groups were classified: (1) the marl facies group, (2) the halite facies group and (3) a group involving different expressions of massive, graded and slumped deposits including coarse clastic detritus. Core 5017-1 encompasses a succession of four main lithological units. Based on first radiocarbon and U-Th ages and correlation of these units to on-shore stratigraphic sections, the record comprises the last ca 220 ka, i.e. the upper part of the Amora Formation (parts of or entire penultimate interglacial and glacial), the last interglacial Samra Fm. (~135-75 ka), the last glacial Lisan Fm. (~75-14 ka) and the Holocene Ze’elim Formation. A major advancement of this record is that, for the first time, also transitional intervals were recovered that are missing in the exposed formations and that can now be studied in great detail.
Micro-facies analyses involve a combination of high-resolution microscopic thin section analysis and µXRF element scanning supported by magnetic susceptibility measurements. This approach allows identifying and characterising micro-facies types, detecting event layers and reconstructing past climate variability with up to seasonal resolution, given that the analysed sediments are annually laminated. Within this thesis, micro-facies analyses, supported by further sedimentological and geochemical analyses (grain size, X-ray diffraction, total organic carbon and calcium carbonate contents) and palynology, were applied for two time intervals:
(1) The early last glacial period ~117-75 ka was investigated focusing on millennial-scale hydroclimatic variations and lake level changes recorded in the sediments. Thereby, distinguishing six different micro-facies types with distinct geochemical and sedimentological characteristics allowed estimating relative lake level and water balance changes of the lake. Comparison of the results to other records in the Mediterranean region suggests a close link of the hydroclimate in the Levant to North Atlantic and Mediterranean climates during the time of the build-up of Northern hemisphere ice sheets during the early last glacial period.
(2) A mostly annually laminated late Holocene section (~3700-1700 cal yr BP) was analysed in unprecedented detail through a multi-proxy, inter-site correlation approach of a shallow-water core (DSEn) and its deep-basin counterpart (5017-1). Within this study, a ca 1500 years comprising time series of erosion and dust deposition events was established and anchored to the absolute time-scale through 14C dating and age modelling. A particular focus of this study was the characterisation of two dry periods, from ~3500 to 3300 and from ~3000 to 2400 cal yr BP, respectively. Thereby, a major outcome was the coincidence of the latter dry period with a period of moist and cold climate in Europe related to a Grand Solar Minimum around 2800 cal yr BP and an increase in flood events despite overall dry conditions in the Dead Sea region during that time. These contrasting climate signatures in Europe and at the Dead Sea were likely linked through complex teleconnections of atmospheric circulation, causing a change in synoptic weather patterns in the eastern Mediterranean.
In summary, within this doctorate the lithostratigraphic framework of a unique long sediment core from the deep Dead Sea basin is established, which serves as a base for any further high-resolution investigations on this core. It is demonstrated in two case studies that micro-facies analyses are an invaluable tool to understand the depositional processes in the Dead Sea and to decipher past climate variability in the Levant on millennial to seasonal time-scales. Hence, this work adds important knowledge helping to establish the deep Dead Sea record as a key climate archive of supra-regional significance.
We present the results of our investigations on the radiocarbon dated core sediments from the Lake Tso Moriri, NW Himalaya aimed at reconstructing palaeohydrological changes in this climatically sensitive region. Based on the detailed geochemical, mineralogical and sedimentological analysis, we recognise several short-term fluctuations superimposed upon seven major palaeohydrological stages identified in this lake since similar to 26 cal ka. Stage I (>20.2 cal ka): shallow lake characterised by input of coarse-grained detrital sediments; Stage II (20.2-16.4 cal ka): lake deepening and intensification of this trend ca. 18 cal ka; Stage III (16.4-11.2 cal ka): rising lake levels with a short term wet phase (13.1-11.7 cal ka); Stage IV (11.2-8.5 cal ka): early Holocene hydrological maxima and highest lake levels inferred to have resulted from early Holocene Indian monsoon intensification, as records from central Asia indicate weaker westerlies during this interval; Stage V (8.5-5.5 cal ka): mid-Holocene climate deterioration; Stage VI (5.5-2.7 cal ka): progressive lowering of lake level; Stage VII (2.7-0 cal ka): onset of modern conditions. The reconstructed hydrological variability in Lake Tso Moriri is governed by temperature changes (meltwater inflow) and monsoon precipitation (increased runoff). A regional comparison shows considerable differences with other palaeorecords from peninsular India during late Holocene. (C) 2014 Elsevier Ltd and INQUA. All rights reserved.
In many procedures of seismic risk mitigation, ground motion simulations are needed to test systems or improve their effectiveness. For example they may be used to estimate the level of ground shaking caused by future earthquakes. Good physical models for ground motion simulation are also thought to be important for hazard assessment, as they could close gaps in the existing datasets. Since the observed ground motion in nature shows a certain variability, part of which cannot be explained by macroscopic parameters such as magnitude or position of an earthquake, it would be desirable that a good physical model is not only able to produce one single seismogram, but also to reveal this natural variability.
In this thesis, I develop a method to model realistic ground motions in a way that is computationally simple to handle, permitting multiple scenario simulations. I focus on two aspects of ground motion modelling. First, I use deterministic wave propagation for the whole frequency range – from static deformation to approximately 10 Hz – but account for source variability by implementing self-similar slip distributions and rough fault interfaces. Second, I scale the source spectrum so that the modelled waveforms represent the correct radiated seismic energy. With this scaling I verify whether the energy magnitude is suitable as an explanatory variable, which characterises the amount of energy radiated at high frequencies – the advantage of the energy magnitude being that it can be deduced from observations, even in real-time.
Applications of the developed method for the 2008 Wenchuan (China) earthquake, the 2003 Tokachi-Oki (Japan) earthquake and the 1994 Northridge (California, USA) earthquake show that the fine source discretisations combined with the small scale source variability ensure that high frequencies are satisfactorily introduced, justifying the deterministic wave propagation approach even at high frequencies. I demonstrate that the energy magnitude can be used to calibrate the high-frequency content in ground motion simulations.
Because deterministic wave propagation is applied to the whole frequency range, the simulation method permits the quantification of the variability in ground motion due to parametric uncertainties in the source description. A large number of scenario simulations for an M=6 earthquake show that the roughness of the source as well as the distribution of fault dislocations have a minor effect on the simulated variability by diminishing directivity effects, while hypocenter location and rupture velocity more strongly influence the variability. The uncertainty in energy magnitude, however, leads to the largest differences of ground motion amplitude between different events, resulting in a variability which is larger than the one observed.
For the presented approach, this dissertation shows (i) the verification of the computational correctness of the code, (ii) the ability to reproduce observed ground motions and (iii) the validation of the simulated ground motion variability. Those three steps are essential to evaluate the suitability of the method for means of seismic risk mitigation.
Based on high-spatiotemporal-resolution data, the authors perform a climatological study of strong rainfall events propagating from southeastern South America to the eastern slopes of the central Andes during the monsoon season. These events account for up to 70% of total seasonal rainfall in these areas. They are of societal relevance because of associated natural hazards in the form of floods and landslides, and they form an intriguing climatic phenomenon, because they propagate against the direction of the low-level moisture flow from the tropics. The responsible synoptic mechanism is analyzed using suitable composites of the relevant atmospheric variables with high temporal resolution. The results suggest that the low-level inflow from the tropics, while important for maintaining sufficient moisture in the area of rainfall, does not initiate the formation of rainfall clusters. Instead, alternating low and high pressure anomalies in midlatitudes, which are associated with an eastward-moving Rossby wave train, in combination with the northwestern Argentinean low, create favorable pressure and wind conditions for frontogenesis and subsequent precipitation events propagating from southeastern South America toward the Bolivian Andes.
Hydrothermal carbonization (HTC) produces carbon-rich nano-micro size particles. In this study, magnetic hydrochar (MHC) was prepared from model compound cellulose by simply adding ferrites during HTC. The effects of ferrites on HTC were evaluated by characterizing solid MHC and corresponding process liquid. Additionally, magnetic stability of MHC was tested by magnetic susceptibility method. Finally, MHC was used as support media for anaerobic films in anaerobic digestion (AD). Ash-free mass yield was around 50% less in MHC than hydrochar produced without ferrites at any certain HTC reaction condition, where organic part of MHC is mainly carbon. In fact, amorphous hydrochar was growing on the surface of inorganic ferrites. MHC maintained magnetic susceptibility regardless of reaction time at reaction temperature 250 degrees C. Pronounced inhibitory effects of magnetic hydrochar occurred during start-up of AD but diminished with prolong AD times. Visible biofilms were observed on the MHC by laser scanning microscope after AD. (C) 2015 Elsevier Ltd. All rights reserved.
Probing the core-mantle boundary beneath Europe and Western Eurasia: A detailed study using PcP
(2015)
We use PcP (the core reflected P phase) recordings of deep earthquakes and nuclear explosions from the Grafenberg (Germany) and NORSAR (Norway) arrays to investigate the core-mantle boundary region beneath Europe and western Eurasia. We find evidence for a previously unknown ultra-low velocity zone 600 km south-east of Moscow, located at the edge of a middle-size low shear- velocity region imaged in seismic tomography that is located beneath the Volga river region. The observed amplitude variations of PcP can be modelled by velocity reductions of P and S-waves of -5% and -15%, respectively, with a density increase of +15%. Travel time delays of pre-and postcursors are indicating a thickness of about 13 km for this ultra-low velocity region (ULVZ). However, our modelling also reveals highly ambiguous amplitude variations of PcP and a reflection off the top of the anomaly for various ULVZs and topography models. Accordingly, large velocity contrasts of up to -10% in V-P and -20% in Vs cannot be excluded. In general, the whole Volga river region shows a complex pattern of PcP amplitudes caused most likely by CMB undulations. Further PcP probes beneath Paris, Kiev and northern Italy indicate likely normal CMB conditions, whereas the samples below Finland and the Hungary-Slovakia border yield strongly amplified PcP signals suggesting strong CMB topography effects.
We evaluate the amplitude behaviour of PcP as a function of distance and several ULVZ models using the 1D reflectivity and the 2D Gauss beam method. The influence of the velocity and density perturbations is analysed as well as the anomaly thickness, the dominant period of the source wavelet and interface topographies. Strong variation of the PcP amplitude are obtained as a function of distance and of the impedance contrast. We also consider two types of topographies: undulations atop the CMB in the presence of flat ULVZs and vice versa. Where a broad range of CMB topography dimensions lead to large PcP amplitude variations, only large ULVZ undulations generate significant amplitude scattering. Consequently, this indicates that topography effects of anomalies may mask the true medium parameters as well as the ULVZ thickness. Moreover, there might be a possibility of misinterpreting the precursor as PcP, in particular for thin ULVZs. (C) 2015 Elsevier B.V. All rights reserved.
Remnants of hydrous melt formed at mantle depths have been identified and characterized within high-pressure leucogranulites of the Orlica-Snieznik Dome (Bohemian Massif, central Europe). They occur as nanogranites in garnet formed via partial melting of granitoids during the Variscan orogeny. Melt composition and H2O content have been investigated in situ after experimental re-homogenization of the nanogranites, and are consistent with melts produced experimentally from crustal lithologies at mantle depths. This is the first geochemical study of melt inclusions from natural crustal rocks equilibrated close to the stability field of coesite, shedding light on how continental crust melts during deep subduction. Whereas decompressional melting is commonly invoked for deeply subducted crustal lithologies, melting occurred near or at the metamorphic peak pressure in the Orlica-Snieznik granulites. Melting of deeply subducted crustal rocks significantly modifies the rheology and thus promotes fast exhumation: this process has a critical influence on the geodynamic evolution of subduction-collision orogens as well as crustal differentiation at depth.
Numerous studies investigated the influence of abiotic (meteorological conditions) and biotic factors (tree characteristics) on stemflow generation. Although these studies identified the variables that influence stemflow volumes in simply structured forests, the combination of tree characteristics that allows a robust prediction of stemflow volumes in species-rich forests is not well known. Many hydrological applications, however, require at least a rough estimate of stemflow volumes based on the characteristics of a forest stand. The need for robust predictions of stemflow motivated us to investigate the relationships between tree characteristics and stemflow volumes in a species-rich tropical forest located in central Panama. Based on a sampling setup consisting of ten rainfall collectors, 300 throughfall samplers and 60 stemflow collectors and cumulated data comprising 26 rain events, we derive three main findings. Firstly, stemflow represents a minor hydrological component in the studied 1-ha forest patch (1.0% of cumulated rainfall). Secondly, in the studied species-rich forest, single tree characteristics are only weakly related to stemflow volumes. The influence of multiple tree parameters (e.g. crown diameter, presence of large epiphytes and inclination of branches) and the dependencies among these parameters require a multivariate approach to understand the generation of stemflow. Thirdly, predicting stemflow in species-rich forests based on tree parameters is a difficult task. Although our best model can capture the variation in stemflow to some degree, a critical validation reveals that the model cannot provide robust predictions of stemflow. A reanalysis of data from previous studies in species-rich forests corroborates this finding. Based on these results and considering that for most hydrological applications, stemflow is only one parameter among others to estimate, we advocate using the base model, i.e. the mean of the stemflow data, to quantify stemflow volumes for a given study area. Studies in species-rich forests that wish to obtain predictions of stemflow based on tree parameters probably need to conduct a much more extensive sampling than currently implemented by most studies. Copyright (c) 2015 John Wiley & Sons, Ltd.
Body size has been widely recognised as a key factor determining community structure in ecosystems. We analysed size diversity patterns of phytoplankton, zooplankton and fish assemblages in 13 data sets from freshwater and marine sites with the aim to assess whether there is a general trend in the effect of predation and resource competition on body size distribution across a wide range of aquatic ecosystems. We used size diversity as a measure of the shape of size distribution. Size diversity was computed based on the Shannon-Wiener diversity expression, adapted to a continuous variable, i.e. as body size. Our results show that greater predation pressure was associated with reduced size diversity of prey at all trophic levels. In contrast, competition effects depended on the trophic level considered. At upper trophic levels (zooplankton and fish), size distributions were more diverse when potential resource availability was low, suggesting that competitive interactions for resources promote diversification of aquatic communities by size. This pattern was not found for phytoplankton size distributions where size diversity mostly increased with low zooplankton grazing and increasing nutrient availability. Relationships we found were weak, indicating that predation and competition are not the only determinants of size distribution. Our results suggest that predation pressure leads to accumulation of organisms in the less predated sizes, while resource competition tends to favour a wider size distribution.
The Central Pontides is an accretionary-type orogenic area within the Alpine-Himalayan orogenic belt characterized by pre-collisional tectonic continental growth. The region comprises Mesozoic subduction-accretionary complexes and an accreted intra-oceanic arc that are sandwiched between the Laurasian active continental margin and Gondwana-derived the Kırşehir Block. The subduction-accretion complexes mainly consist of an Albian-Turonian accretionary wedge representing the Laurasian active continental margin. To the north, the wedge consists of slate/phyllite and metasandstone intercalation with recrystallized limestone, Na-amphibole-bearing metabasite (PT= 7–12 kbar and 400 ± 70 ºC) and tectonic slices of serpentinite representing accreted distal part of a large Lower Cretaceous submarine turbidite fan deposited on the Laurasian active continental margin that was subsequently accreted and metamorphosed. Raman spectra of carbonaceous material (RSCM) of the metapelitic rocks revealed that the metaflysch sequence consists of metamorphic packets with distinct peak metamorphic temperatures. The majority of the metapelites are low-temperature (ca. 330 °C) slates characterized by lack of differentiation of the graphite (G) and D2 defect bands. They possibly represent offscraped distal turbidites along the toe of the Albian accretionary wedge. The rest are phyllites that are characterized by slightly pronounced G band with D2 defect band occurring on its shoulder. Peak metamorphic temperatures of these phyllites are constrained to 370-385 °C. The phyllites are associated with a strip of incipient blueschist facies metabasites which are found as slivers within the offscraped distal turbidites. They possibly represent underplated continental metasediments together with oceanic crustal basalt along the basal décollement. Tectonic emplacement of the underplated rocks into the offscraped distal turbidites was possibly achieved by out-of-sequence thrusting causing tectonic thickening and uplift of the wedge. 40Ar/39Ar phengite ages from the phyllites are ca. 100 Ma, indicating Albian subduction and regional HP metamorphism.
The accreted continental metasediments are underlain by HP/LT metamorphic rocks of oceanic origin along an extensional shear zone. The oceanic metamorphic sequence mainly comprises tectonically thickened deep-seated eclogite to blueschist facies metabasites and micaschists. In the studied area, metabasites are epidote-blueschists locally with garnet (PT= 17 ± 1 kbar and 500 ± 40 °C). Lawsonite-blueschists are exposed as blocks along the extensional shear zone (PT= 14 ± 2 kbar and 370–440 °C). They are possibly associated with low shear stress regime of the initial stage of convergence. Close to the shear zone, the footwall micaschists consist of quartz, phengite, paragonite, chlorite, rutile with syn-kinematic albite porphyroblast formed by pervasive shearing during exhumation. These types of micaschists are tourmaline-bearing and their retrograde nature suggests high-fluid flux along shear zones. Peak metamorphic mineral assemblages are partly preserved in the chloritoid-micaschist farther away from the shear zone representing the zero strain domains during exhumation. Three peak metamorphic assemblages are identified and their PT conditions are constrained by pseudosections produced by Theriak-Domino and by Raman spectra of carbonaceous material: 1) garnet-chloritoid-glaucophane with lawsonite pseudomorphs (P= 17.5 ± 1 kbar, T: 390-450 °C) 2) chloritoid with glaucophane pseudomorphs (P= 16-18 kbar, T: 475 ± 40 °C) and 3) relatively high-Mg chloritoid (17%) with jadeite pseudomorphs (P= 22-25 kbar; T: 440 ± 30 °C) in addition to phengite, paragonite, quartz, chlorite, rutile and apatite. The last mineral assemblage is interpreted as transformation of the chloritoid + glaucophane assemblage to chloritoid + jadeite paragenesis with increasing pressure. Absence of tourmaline suggests that the chloritoid-micaschist did not interact with B-rich fluids during zero strain exhumation. 40Ar/39Ar phengite age of a pervasively sheared footwall micaschist is constrained to 100.6 ± 1.3 Ma and that of a chloritoid-micaschist is constrained to 91.8 ± 1.8 Ma suggesting exhumation during on-going subduction with a southward younging of the basal accretion and the regional metamorphism. To the south, accretionary wedge consists of blueschist and greenschist facies metabasite, marble and volcanogenic metasediment intercalation. 40Ar/39Ar phengite dating reveals that this part of the wedge is of Middle Jurassic age partly overprinted during the Albian. Emplacement of the Middle Jurassic subduction-accretion complexes is possibly associated with obliquity of the Albian convergence.
Peak metamorphic assemblages and PT estimates of the deep-seated oceanic metamorphic sequence suggest tectonic stacking within wedge with different depths of burial. Coupling and exhumation of the distinct metamorphic slices are controlled by decompression of the wedge possibly along a retreating slab. Structurally, decompression of the wedge is evident by an extensional shear zone and the footwall micaschists with syn-kinematic albite porphyroblasts. Post-kinematic garnets with increasing grossular content and pseudomorphing minerals within the chloritoid-micaschists also support decompression model without an extra heating.
Thickening of subduction-accretionary complexes is attributed to i) significant amount of clastic sediment supply from the overriding continental domain and ii) deep level basal underplating by propagation of the décollement along a retreating slab. Underplating by basal décollement propagation and subsequent exhumation of the deep-seated subduction-accretion complexes are connected and controlled by slab rollback creating a necessary space for progressive basal accretion along the plate interface and extension of the wedge above for exhumation of the tectonically thickened metamorphic sequences. This might be the most common mechanism of the tectonic thickening and subsequent exhumation of deep-seated HP/LT subduction-accretion complexes.
To the south, the Albian-Turonian accretionary wedge structurally overlies a low-grade volcanic arc sequence consisting of low-grade metavolcanic rocks and overlying metasedimentary succession is exposed north of the İzmir-Ankara-Erzincan suture (İAES), separating Laurasia from Gondwana-derived terranes. The metavolcanic rocks mainly consist of basaltic andesite/andesite and mafic cognate xenolith-bearing rhyolite with their pyroclastic equivalents, which are interbedded with recrystallized pelagic limestone and chert. The metavolcanic rocks are stratigraphically overlain by recrystallized micritic limestone with rare volcanogenic metaclastic rocks. Two groups can be identified based on trace and rare earth element characteristics. The first group consists of basaltic andesite/andesite (BA1) and rhyolite with abundant cognate gabbroic xenoliths. It is characterized by relative enrichment of LREE with respect to HREE. The rocks are enriched in fluid mobile LILE, and strongly depleted in Ti and P reflecting fractionation of Fe-Ti oxides and apatite, which are found in the mafic cognate xenoliths. Abundant cognate gabbroic xenoliths and identical trace and rare earth elements compositions suggest that rhyolites and basaltic andesites/andesites (BA1) are cogenetic and felsic rocks were derived from a common mafic parental magma by fractional crystallization and accumulation processes. The second group consists only of basaltic andesites (BA2) with flat REE pattern resembling island arc tholeiites. Although enriched in LILE, this group is not depleted in Ti or P.
Geochemistry of the metavolcanic rocks indicates supra-subduction volcanism evidenced by depletion of HFSE and enrichment of LILE. The arc sequence is sandwiched between an Albian-Turonian subduction-accretionary complex representing the Laurasian active margin and an ophiolitic mélange. Absence of continent derived detritus in the arc sequence and its tectonic setting in a wide Cretaceous accretionary complex suggest that the Kösdağ Arc was intra-oceanic. This is in accordance with basaltic andesites (BA2) with island arc tholeiite REE pattern.
Zircons from two metarhyolite samples give Late Cretaceous (93.8 ± 1.9 and 94.4 ± 1.9 Ma) U/Pb ages. Low-grade regional metamorphism of the intra-oceanic arc sequence is constrained 69.9 ± 0.4 Ma by 40Ar/39Ar dating on metamorphic muscovite from a metarhyolite indicating that the arc sequence became part of a wide Tethyan Cretaceous accretionary complex by the latest Cretaceous. The youngest 40Ar/39Ar phengite age from the overlying subduction-accretion complexes is 92 Ma confirming southward younging of an accretionary-type orogenic belt. Hence, the arc sequence represents an intra-oceanic paleo-arc that formed above the sinking Tethyan slab and finally accreted to Laurasian active continental margin. Abrupt non-collisional termination of arc volcanism was possibly associated with southward migration of the arc volcanism similar to the Izu-Bonin-Mariana arc system.
The intra-oceanic Kösdağ Arc is coeval with the obducted supra-subduction ophiolites in NW Turkey suggesting that it represents part of the presumed but missing incipient intra-oceanic arc associated with the generation of the regional supra-subduction ophiolites. Remnants of a Late Cretaceous intra-oceanic paleo-arc and supra-subduction ophiolites can be traced eastward within the Alp-Himalayan orogenic belt. This reveals that Late Cretaceous intra-oceanic subduction occurred as connected event above the sinking Tethyan slab. It resulted as arc accretion to Laurasian active margin and supra-subduction ophiolite obduction on Gondwana-derived terranes.
A detailed architectural analysis has been performed on the porcelaneous species Taberina bingistani Henson, 1948 from diverse localities, including its type locality in Iran. Characterized by a subcylindrical to slightly conical shape, with a cribbed apertural face, the Cenomanian Middle-East species bingistani is removed from the American early-Paleogene genus Taberina Keijzer, 1945 (type species: T. cubana) and ascribed to the new genus Praetaberina. In the early stages of growth, the chambers are planispiral-involute, becoming uncoiled in the later stages. The interior of the chambers shows numerous septula and pillars. The septula are radial and distributed in the marginal area of the chamber lumen, while the pillars occupy the central area and are regularly distributed between the intercameral foramina. Pillars are in alignment from one chamber to the next. A new species, from southeast Italy, Praetaberina apula, is also erected. It differs from P. bingistani in its smaller size and significantly reduced structural elements. Both species characterize upper Cenomanian shallow-water deposits.
The Upper Devonian Munster Basin of southern Ireland has traditionally been viewed as a post-orogenic molasse deposit that was sourced from the Caledonides of central Ireland and subsequently deformed by the end Carboniferous Variscan orogenic event. The basin fill is composed of super-mature quartz arenite sandstone that clearly represents a second cycle of deposition. The source of this detritus is now recognized as Lower Devonian Dingle Basin red bed sequences to the north. This genetic link is based on the degree of similarity in the detrital mica chemistry in both of these units; micas plot in identical fields and define the same trends. In addition, the two sequences show increased textural and chemical maturity up-sequence and define indistinguishable Ar-40/Ar-39 age ranges for the detrital mica grains. Partial resetting of the Ar ages can be attributed to elevated heat flow in the region caused by Munster Basin extension and subsequent Variscan deformation. The combined evidence from southwest Ireland therefore points to a Caledonian or possibly Taconian primary source area that initially shed detritus into the Lower Devonian Dingle Basin which was subsequently recycled into the Upper Devonian Munster Basin following mid-Devonian Acadian basin inversion. (C) 2014 International Association for Gondwana Research. Published by Elsevier B.V. All rights reserved.
Scientific theories of how subduction and plate tectonics began on Earth-and what the tectonic structure of Earth was before this-remain enigmatic and contentious(1). Understanding viable scenarios for the onset of subduction and plate tectonics(2,3) is hampered by the fact that subduction initiation processes must have been markedly different before the onset of global plate tectonics because most present-day subduction initiation mechanisms require acting plate forces and existing zones of lithospheric weakness, which are both consequences of plate tectonics(4). However, plume-induced subduction initiation(5-9) could have started the first subduction zone without the help of plate tectonics. Here, we test this mechanism using high-resolution three-dimensional numerical thermomechanical modelling. We demonstrate that three key physical factors combine to trigger self-sustained subduction: (1) a strong, negatively buoyant oceanic lithosphere; (2) focused magmatic weakening and thinning of lithosphere above the plume; and (3) lubrication of the slab interface by hydrated crust. We also show that plume-induced subduction could only have been feasible in the hotter early Earth for old oceanic plates. In contrast, younger plates favoured episodic lithospheric drips rather than self-sustained subduction and global plate tectonics.
Plan oblique relief shows terrain with a side view on a two-dimensional map, resulting in visualizations where the third dimension of the terrain is more explicit than on traditional two-dimensional maps. Existing plan oblique maps are static: the angle of terrain inclination is not adjustable and the orientation of plan oblique inclination does not change with the orientation of the map. This article introduces two complementary methods that address these issues by using the 3D graphics pipeline to render plan oblique relief for tile-based web maps. The goal is to allow users to adjust the terrain inclination and map rotation angles to better visualize the third dimension of the terrain. The first method pre-renders plan oblique tiles with a server-side application. The tiles are visualized with a standard web mapping framework. The second method renders plan oblique relief on-the-fly in a web browser using WebGL and a customized version of OpenLayers 3, which enables users to select arbitrary terrain inclination and map rotation angles. The second method uses a tiled digital terrain model that is loaded by the web browser. The browser applies the plan oblique transformation, computes a shaded relief, and texturizes the terrain with tiled map layers.
Endosulfan - an agricultural insecticide and banned by Stockholm Convention - is produced as a 2:1 to 7:3 mixture of isomers endosulfan I (ESI) and endosulfan II (ESII). Endosulfan is transformed under aerobic conditions into endosulfan sulfate (ESS). The study shows for 76 sampling locations in German forests that endosulfan is abundant in all samples with an opposite ratio between the ESI and ESII than the technical product, where the main metabolite ESS is found with even higher abundance. The ratio between ESI/ESII and ESS show clear dependence on the type of stands (coniferous vs. deciduous) and humus type and increases from deciduous via mixed to coniferous forest stands. The study argues for a systematic monitoring of ESI, ESII, and ESS and underlines the need for further research, specifically on the fate of endosulfan including biomagnifications and bioaccumulation in soil. (C) 2015 Elsevier Ltd. All rights reserved.
Several solutions have been proposed to explain the long-standing kinematic observation that postcollisional upper crustal shortening within the Himalaya and Asia is much less than the magnitude of India-Asia convergence. Here we implement these hypotheses in global plate reconstructions and test paleolatitudes predicted by the global apparent polar wander path against independent, and the most robust paleomagnetic data. Our tests demonstrate that (1) reconstructed 600-750km postcollisional intra-Asian shortening is a minimum value; (2) a 52Ma collision age is only consistent with paleomagnetic data if intra-Asian shortening was 900km; a 56-58Ma collision age requires greater intra-Asian shortening; (3) collision ages of 34 or 65Ma incorrectly predict Late Cretaceous and Paleogene paleolatitudes of the Tibetan Himalaya (TH); and (4) Cretaceous counterclockwise rotation of India cannot explain the paleolatitudinal divergence between the TH and India. All hypotheses, regardless of collision age, require major Cretaceous extension within Greater India.
The Tibetan Himalaya represents the northernmost continental unit of the Indian plate that collided with Asia in the Cenozoic. Paleomagnetic studies on the Tibetan Himalaya can help constrain the dimension and paleogeography of "Greater India,' the Indian plate lithosphere that subducted and underthrusted below Asia after initial collision. Here we present a paleomagnetic investigation of a Jurassic (limestones) and Lower Cretaceous (volcaniclastic sandstones) section of the Tibetan Himalaya. The limestones yielded positive fold test, showing a prefolding origin of the isolated remanent magnetizations. Detailed paleomagnetic analyses, rock magnetic tests, end-member modeling of acquisition curves of isothermal remanent magnetization, and petrographic investigation reveal that the magnetic carrier of the Jurassic limestones is authigenic magnetite, whereas the dominant magnetic carrier of the Lower Cretaceous volcaniclastic sandstones is detrital magnetite. Our observations lead us to conclude that the Jurassic limestones record a prefolding remagnetization, whereas the Lower Cretaceous volcaniclastic sandstones retain a primary remanence. The volcaniclastic sandstones yield an Early Cretaceous paleolatitude of 55.5 degrees S [52.5 degrees S, 58.6 degrees S] for the Tibetan Himalaya, suggesting it was part of the Indian continent at that time. The size of "Greater India' during Jurassic time cannot be estimated from these limestones. Instead, a paleolatitude of the Tibetan Himalaya of 23.8 degrees S [21.8 degrees S, 26.1 degrees S] during the remagnetization process is suggested. It is likely that the remagnetization, caused by the oxidation of early diagenetic pyrite to magnetite, was induced during 103-83 or 77-67 Ma. The inferred paleolatitudes at these two time intervals imply very different tectonic consequences for the Tibetan Himalaya.
Sedimentological, palaeontological and mineralogical analyses of sediments from the endorheic Al Jafr Basin were conducted to better understand the depositional and hydrological conditions on the southern Jordan Plateau in the late Quaternary. Surficially exposed carbonate-rich sediments in the western part of the basin contain ostracod (micro-crustacean) shells of Ilyocypris cf. bradyi, Candona neglecta, Heterocypris salina, Fabaeformiscandona fabaeformis, Pseudocandona sp. and Herpetocypris brevicaudata. The shells of these and other more rare species, and charophyte and mollusc remains indicate that the sediments were formed in a wetland setting of shallow freshwater to slightly oligohaline ponds, streams and swamps. The present more northern distribution of some of the recorded taxa implies that climate conditions were probably cooler during the wetland formation. Radiocarbon age data for biogenic carbonate from two locations suggest that the wetland setting existed during the second half of Marine Isotope Stage (MIS) 3 or possibly earlier. A significantly higher water table must have existed in the basin during wetland formation; and wetter climate conditions are inferred for the catchment or at least for its highest and most humid westernmost part. Deflation and local sediment accumulation by wind and occasional sheet-wash events apparently prevailed in the region since MIS 2. Our newly presented data and inferences do not support the reconstruction of a previously reported large and relatively deep Pleistocene lake in the Al Jafr Basin. However, more extensive studies are certainly required for a detailed assessment of the Quaternary hydrological conditions in southern Jordan. (C) 2014 Elsevier Ltd and INQUA. All rights reserved.
The organic-carbon (OC) pool accumulated in Arctic permafrost (perennially frozen ground) equals the carbon stored in the modern atmosphere. To give an idea of how Yedoma region permafrost could respond under future climatic warming, we conducted a study to quantify the organic-matter quality (here defined as the intrinsic potential to be further transformed, decomposed, and mineralized) of late Pleistocene (Yedoma) and Holocene (thermokarst) deposits on the Buor-Khaya Peninsula, northeast Siberia. The objective of this study was to develop a stratigraphic classified organic-matter quality characterization. For this purpose the degree of organic-matter decomposition was estimated by using a multiproxy approach. We applied sedimentological (grain-size analyses, bulk density, ice content) and geochemical parameters (total OC, stable carbon isotopes (delta C-13),total organic carbon : nitrogen (C / N) ratios) as well as lipid biomarkers (n-alkanes, n-fatty acids, hopanes, triterpenoids, and biomarker indices, i.e., average chain length, carbon preference index (CPI), and higher-plant fatty-acid index (HPFA)). Our results show that the Yedoma and thermokarst organic-matter qualities for further decomposition exhibit no obvious degradation-depth trend. Relatively, the C / N and delta C-13 values and the HPFA index show a significantly better preservation of the organic matter stored in thermokarst deposits compared to Yedoma deposits. The CPI data suggest less degradation of the organic matter from both deposits, with a higher value for Yedoma organic matter. As the interquartile ranges of the proxies mostly over-lap, we interpret this as indicating comparable quality for further decomposition for both kinds of deposits with likely better thermokarst organic-matter quality. Supported by principal component analyses, the sediment parameters and quality proxies of Yedoma and thermokarst deposits could not be unambiguously separated from each other. This revealed that the organic-matter vulnerability is heterogeneous and depends on different decomposition trajectories and the previous decomposition and preservation history. Elucidating this was one of the major new contributions of our multiproxy study. With the addition of biomarker data, it was possible to show that permafrost organic-matter degradation likely occurs via a combination of (uncompleted) degradation cycles or a cascade of degradation steps rather than as a linear function of age or sediment facies. We conclude that the amount of organic matter in the studied sediments is high for mineral soils and of good quality and therefore susceptible to future decomposition. The lack of depth trends shows that permafrost acts like a giant freezer, preserving the constant quality of ancient organic matter. When undecomposed Yedoma organic matter is mobilized via thermokarst processes, the fate of this carbon depends largely on the environmental conditions; the carbon could be preserved in an undecomposed state till refreezing occurs. If modern input has occurred, thermokarst organic matter could be of a better quality for future microbial decomposition than that found in Yedoma deposits.
One of the most critical problems affecting geophysical data acquisition procedures is related to the misorientation of multicomponent sensors with respect to a common reference system (e.g., geographic north). In many applications, misoriented sensors affect data analysis procedures, leading to errors in results and interpretations. These problems generally occur in applications where the orientation of the sensor cannot be actively controlled and is not known a priori, e.g., geophysical sensors deployed in borehole installations or on the seafloor. We have developed a quaternion-based method for the optimal reorientation of multicomponent geophysical sensors. In contrast to other approaches, we took into account the full time-series record from all sensor components. Therefore, our method could be applied to all time-series data and was not restricted to a certain type of geophysical sensor. Our method allows the robust calculation of relative reorientations between two-component or three-component sensors. By using a reference sensor in an iterative process, this result can be extended to the estimation of absolute sensor orientations. In addition to finding an optimal solution for a full 3D sensor rotation, we have established a rigorous scheme for the estimation of uncertainties of the resulting orientation parameters. We tested the feasibility and applicability of our method using synthetic data examples for a vertical seismic profile and an ocean bottom seismometer array. We noted that the quaternion-based reorientation method is superior to the standard approach of a single-parameter estimation of rotation angles.
The use of topographic metrics for estimating the susceptibility to, and reconstructing the characteristics of, debris flows has a long research tradition, although largely devoted to humid mountainous terrain. The exceptional 2010 monsoonal rainstorms in the high-altitude mountain desert of Ladakh and Zanskar, NW India, were a painful reminder of how susceptible arid regions are to rainfall-triggered flash floods, landslides, and debris flows. The rainstorms of August 4-6 triggered numerous debris flows, killing 182 people, devastating 607 houses, and more than 10 bridges around Ladakh's capital of Leh. The lessons from this disaster motivated us to revisit methods of predicting (a) flow parameters such as peak discharge and maximum velocity from field and remote sensing data, and (b) the susceptibility to debris flows from catchment morphometry. We focus on quantifying uncertainties tied to these approaches. Comparison of high-resolution satellite images pre- and post-dating the 2010 rainstorm reveals the extent of damage and catastrophic channel widening. Computations based on these geomorphic markers indicate maximum flow velocities of 1.6-6.7 m s(-1) with runout of up to similar to 10 km on several alluvial fans that sustain most of the region's settlements. We estimate median peak discharges of 310-610 m(3) s(-1), which are largely consistent with previous estimates. Monte Carlo-based error propagation for a single given flow-reconstruction method returns a variance in discharge similar to one derived from juxtaposing several different flow reconstruction methods. We further compare discriminant analysis, classification tree modelling, and Bayesian logistic regression to predict debris-flow susceptibility from morphometric variables of 171 catchments in the Ladakh Range. These methods distinguish between fluvial and debris flow-prone catchments at similar success rates, but Bayesian logistic regression allows quantifying uncertainties and relationships between potential predictors. We conclude that, in order to be robust and reliable, morphometric reconstruction of debris-flow properties and susceptibility requires careful assessment and reporting of errors and uncertainties. (C) 2015 Elsevier B.V. All rights reserved.
For near surface geophysical surveys, small-fixed offset loop-loop electromagnetic induction (EMI) sensors are usually placed parallel to the ground surface (i.e., both loops are at the same height above ground). In this study, we evaluate the potential of making measurements with a system that is not parallel to the ground; i.e., by positioning the system at different inclinations with respect to ground surface. First, we present the Maxwell theory for inclined magnetic dipoles over a homogeneous half space. By analyzing the sensitivities of such configurations, we,show that varying the angle of the system would result in improved imaging capabilities. For example, we show that acquiring data with a vertical system allows detection of a conductive body with a better lateral resolution compared to data acquired using standard horizontal configurations. The synthetic responses are presented for a heterogeneous medium and compared to field data acquired in the historical Park Sanssouci in Potsdam, Germany. After presenting a detailed sensitivity analysis and synthetic examples of such ground conductivity measurements, we suggest a new strategy of acquisition that allows to better estimate the true distribution of electrical conductivity using instruments with a fixed, small offset between the loops. This strategy is evaluated using field data collected at a well-constrained test-site in Horstwalde (Germany). Here, the target buried utility pipes are best imaged using vertical system configurations demonstrating the potential of our approach for typical applications. (C) 2015 Elsevier B.V. Pill rights reserved.
Non-invasive imaging techniques to study O-2 micro-patterns around pesticide treated lupine roots
(2015)
The soil root interface is a highly heterogeneous system, e.g. in terms of O-2 and pH distribution. The destructive character of conventional methods disturbs the natural conditions of those biogeochemical gradients. Therefore, experiments aiming to control these influences and study pesticide kinetics under given O-2 and pH conditions suffer from a large uncertainty of the "real" O-2/pH at a certain position. Our approach with two different imaging techniques will examine the soil-root interface as well as the dissipation of the applied pesticide at a high spatial resolution.
The obtained outcomes show directly that the pH has an influence on enantioselective dissipation of the acetanilide fungicide metalaxyl. In areas with high pH from an applied racemic mixture, the R-enantiomer dissipates faster than the S-enantiomer. Moreover, we found significantly reduced oxygen values in the bulk soil and vicinity of metalaxyl treated roots compared to control plant roots. The combination of matrix-assisted laser desorption/ionization mass spectrometry (MALDI) and fluorescence imaging indicated the oxygen-dependent behavior of metalaxyl at the root surface.
The results presented here underline the great potential of combining different imaging methods to examine the soil-root interfaces as well as the dissipation of organic pollutants in small soil compartments. (C) 2014 Elsevier B.V. All rights reserved.
Previous studies suggested that electric and/or magnetic field fluctuations observed in the nighttime topside ionosphere at midlatitudes generally originate from quiet time nocturnal medium-scale traveling ionospheric disturbances (MSTIDs). However, decisive evidences for the connection between the two have been missing. In this study we make use of the multispacecraft observations of midlatitude magnetic fluctuations (MMFs) in the nighttime topside ionosphere by the Swarm constellation. The analysis results show that the area hosting MMFs is elongated in the NW-SE (NE-SW) direction in the Northern (Southern) Hemisphere. The elongation direction and the magnetic field polarization support that the area hosting MMFs is nearly field aligned. All these properties of MMFs suggest that they have close relationship with MSTIDs. Expectation values of root-mean-square field-aligned currents associated with MMFs are up to about 4nA/m(2). MMF coherency significantly drops for longitudinal distances of 1 degrees.
Modern imaging spectrometers produce an ever-growing amount of data, which increases the need for automated analysis techniques. The algorithms employed, such as the United States Geological Survey (USGS) Tetracorder and the Mineral Identification and Characterization Algorithm (MICA), use a standardized spectral library and expert knowledge for the detection of surface cover types. Correct absorption feature definition and isolation are key to successful material identification using these algorithms. Here, a new continuum removal and feature isolation technique is presented, named the 'Geometric Hull Technique'. It is compared to the well-established, knowledge-based Tetracorder feature database together with the adapted state of the art techniques scale-space filtering, alpha shapes and convex hull.
The results show that the geometric hull technique yields the smallest deviations from the feature definitions of the MICA Group 2 library with a median difference of only 8nm for the position of the features and a median difference of only 15% for the feature shapes. The modified scale-space filtering hull technique performs second best with a median feature position difference of 16nm and a median difference of 25% for the feature shapes. The scale-space alpha hull technique shows a 23nm median position difference and a median deviation of 77% for the feature shapes. The geometric hull technique proposed here performs best amongst the four feature isolation techniques and may be an important building block for next generation automatic mapping algorithms.
Earth's life-support systems are in rapid decline, yet we have few metrics or indicators with which to track these changes. The world's governments are calling for biodiversity and ecosystem-service monitoring to guide and evaluate international conservation policy as well as to incorporate natural capital into their national accounts. The Group on Earth Observations Biodiversity Observation Network (GEO BON) has been tasked with setting up this monitoring system. Here we explore the immediate feasibility of creating a global ecosystem-service monitoring platform under the GEO BON framework through combining data from national statistics, global vegetation models, and production function models. We found that nine ecosystem services could be annually reported at a national scale in the short term: carbon sequestration, water supply for hydropower, and non-fisheries marine products, crop, livestock, game meat, fisheries, mariculture, and timber production. Reported changes in service delivery over time reflected ecological shocks (e.g., droughts and disease outbreaks), highlighting the immediate utility of this monitoring system. Our work also identified three opportunities for creating a more comprehensive monitoring system. First, investing in input data for ecological process models (e.g., global land-use maps) would allow many more regulating services to be monitored. Currently, only 1 of 9 services that can be reported is a regulating service. Second, household surveys and censuses could help evaluate how nature affects people and provides non-monetary benefits. Finally, to forecast the sustainability of service delivery, research efforts could focus on calculating the total remaining biophysical stocks of provisioning services. Regardless, we demonstrated that a preliminary ecosystem-service monitoring platform is immediately feasible. With sufficient international investment, the platform could evolve further into a much-needed system to track changes in our planet's life-support systems. (C) 2015 Elsevier Ltd. All rights reserved.
Multi-scale analysis of electrical conductivity of peatlands for the assessment of peat properties
(2015)
Peatlands store large amounts of carbon. This storage function has been reduced through intensive drainage, which leads to the decomposition of peat, resulting in a loss of carbon. Measurements of the real (sigma) and imaginary part (sigma) of electrical conductivity can deliver information on peat properties, such as the pore fluid conductivity (sigma(w)), cation exchange capacity (CEC), bulk density ((b)), water content (WC) and soil organic matter (SOM) content. These properties change with the peat's degree of decomposition (DD). To explore the relationships between the peat properties, sigma, sigma and DD, we focused on three different types of survey and scales. First, point measurements were made with a conductivity probe at various locations over a large area of northeast Germany to determine the degree of correlation between sigma and DD. Second, nine of these locations were selected for sampling to determine which of the properties sigma(w), CEC, (b), WC and SOM predominantly influence sigma and sigma. This multisite dataset includes the entire range of DD and was analysed in the laboratory. Third, one site was selected for a survey of sigma including sampling, to identify which properties mainly control sigma in a single-site approach. Statistical analysis revealed that for the multisite laboratory dataset, sigma(w) has the strongest effect on sigma, followed by CEC, whereas sigma is mainly determined by CEC. In a single-site approach, WC followed by CEC had a dominant effect on sigma. No clear correlation could be observed between (i) DD and peat properties and (ii) DD and sigma or sigma. This is because of the complex changes in properties with increasing DD.