Refine
Year of publication
- 2013 (187) (remove)
Document Type
- Article (140)
- Doctoral Thesis (35)
- Review (5)
- Preprint (4)
- Other (2)
- Monograph/Edited Volume (1)
Is part of the Bibliography
- yes (187)
Keywords
- Climate change (4)
- Land use change (4)
- Iran (3)
- Northeast German Basin (3)
- Soil hydrology (3)
- climate change (3)
- erosion (3)
- exhumation (3)
- remote sensing (3)
- thermochronology (3)
Institute
- Institut für Geowissenschaften (187) (remove)
A new sedimentary sequence from Lago di Venere on Pantelleria Island, located in the Strait of Sicily between Tunisia and Sicily was recovered. The lake is located in the coastal infra-Mediterranean vegetation belt at 2 m a.s.l. Pollen, charcoal and sedimentological analyses are used to explore linkages among vegetation, fire and climate at a decadal scale over the past 1200 years. A dry period from ad 800 to 1000 that corresponds to the Medieval Warm Period' (WMP) is inferred from sedimentological analysis. The high content of carbonate recorded in this period suggests a dry phase, when the ratio of evaporation/precipitation was high. During this period the island was dominated by thermophilous and drought-tolerant taxa, such as Quercus ilex, Olea, Pistacia and Juniperus. A marked shift in the sediment properties is recorded at ad 1000, when carbonate content became very low suggesting wetter conditions until ad 1850-1900. Broadly, this period coincides with the Little Ice Age' (LIA), which was characterized by wetter and colder conditions in Europe. During this time rather mesic conifers (i.e. Pinus pinaster), shrubs and herbs (e.g. Erica arborea and Selaginella denticulata) expanded, whereas more drought-adapted species (e.g. Q. ilex) declined. Charcoal data suggest enhanced fire activity during the LIA probably as a consequence of anthropogenic burning and/or more flammable fuel (e.g. resinous Pinus biomass). The last century was characterized by a shift to high carbonate content, indicating a change towards drier conditions, and re-expansion of Q. ilex and Olea. The post-LIA warming is in agreement with historical documents and meteorological time series. Vegetation dynamics were co-determined by agricultural activities on the island. Anthropogenic indicators (e.g. Cerealia-type, Sporormiella) reveal the importance of crops and grazing on the island. Our pollen data suggest that extensive logging caused the local extinction of deciduous Quercus pubescens around ad1750.
The study of outcrop modeling is located at the interface between two fields of expertise, Sedimentology and Computing Geoscience, which respectively investigates and simulates geological heterogeneity observed in the sedimentary record. During the last past years, modeling tools and techniques were constantly improved. In parallel, the study of Phanerozoic carbonate deposits emphasized the common occurrence of a random facies distribution along single depositional domain. Although both fields of expertise are intrinsically linked during outcrop simulation, their respective advances have not been combined in literature to enhance carbonate modeling studies. The present study re-examines the modeling strategy adapted to the simulation of shallow-water carbonate systems, based on a close relationship between field sedimentology and modeling capabilities. In the present study, the evaluation of three commonly used algorithms Truncated Gaussian Simulation (TGSim), Sequential Indicator Simulation (SISim), and Indicator Kriging (IK), were performed for the first time using visual and quantitative comparisons on an ideally suited carbonate outcrop. The results show that the heterogeneity of carbonate rocks cannot be fully simulated using one single algorithm. The operating mode of each algorithm involves capabilities as well as drawbacks that are not capable to match all field observations carried out across the modeling area. Two end members in the spectrum of carbonate depositional settings, a low-angle Jurassic ramp (High Atlas, Morocco) and a Triassic isolated platform (Dolomites, Italy), were investigated to obtain a complete overview of the geological heterogeneity in shallow-water carbonate systems. Field sedimentology and statistical analysis performed on the type, morphology, distribution, and association of carbonate bodies and combined with palaeodepositional reconstructions, emphasize similar results. At the basin scale (x 1 km), facies association, composed of facies recording similar depositional conditions, displays linear and ordered transitions between depositional domains. Contrarily, at the bedding scale (x 0.1 km), individual lithofacies type shows a mosaic-like distribution consisting of an arrangement of spatially independent lithofacies bodies along the depositional profile. The increase of spatial disorder from the basin to bedding scale results from the influence of autocyclic factors on the transport and deposition of carbonate sediments. Scale-dependent types of carbonate heterogeneity are linked with the evaluation of algorithms in order to establish a modeling strategy that considers both the sedimentary characteristics of the outcrop and the modeling capabilities. A surface-based modeling approach was used to model depositional sequences. Facies associations were populated using TGSim to preserve ordered trends between depositional domains. At the lithofacies scale, a fully stochastic approach with SISim was applied to simulate a mosaic-like lithofacies distribution. This new workflow is designed to improve the simulation of carbonate rocks, based on the modeling of each scale of heterogeneity individually. Contrarily to simulation methods applied in literature, the present study considers that the use of one single simulation technique is unlikely to correctly model the natural patterns and variability of carbonate rocks. The implementation of different techniques customized for each level of the stratigraphic hierarchy provides the essential computing flexibility to model carbonate systems. Closer feedback between advances carried out in the field of Sedimentology and Computing Geoscience should be promoted during future outcrop simulations for the enhancement of 3-D geological models.
Based on a numerical model of the Northeast German Basin (NEGB), we investigate the sensitivity of the calculated thermal field as resulting from heat conduction, forced and free convection in response to consecutive horizontal and vertical mesh refinements. Our results suggest that computational findings are more sensitive to consecutive horizontal mesh refinements than to changes in the vertical resolution. In addition, the degree of mesh sensitivity depends strongly on the type of the process being investigated, whether heat conduction, forced convection or free thermal convection represents the active heat driver. In this regard, heat conduction exhibits to be relative robust to imposed changes in the spatial discretization. A systematic mesh sensitivity is observed in areas where forced convection promotes an effective role in shorten the background conductive thermal field. In contrast, free thermal convection is to be regarded as the most sensitive heat transport process as demonstrated by non-systematic changes in the temperature field with respect to imposed changes in the model resolution.
Reliable hydrological monitoring is the basis for sound water management in drained wetlands. Since statistical methods cannot be employed for unobserved or sparsely monitored areas, the primary design (first set-up) may be arbitrary in most instances. The objective of this paper is therefore to provide a guideline for designing the initial hydrological monitoring network. A scheme is developed that handles different parts of monitoring and hydrometry in wetlands, focusing on the positioning of surface water and groundwater gauges. For placement of the former, control units are used which correspond to areas whose water levels can be regulated separately. The latter are arranged depending on hydrological response units, defined by combinations of soil type and land use, and the chosen surface water monitoring sites. A practical application of the approach is shown for an investigation area in the Spreewald region in north-east Germany. The presented scheme leaves a certain degree of freedom to its user, allowing the inclusion of expert knowledge or special concerns. Based on easily obtainable data, the developed hydrological network serves as a first step in the iterative procedure of monitoring network optimisation. Copyright (c) 2013 John Wiley & Sons, Ltd.
A total of 271 pollen records were selected from a large collection of both raw and digitized pollen spectra from eastern continental Asia (70 degrees-135 degrees E and 18 degrees-55 degrees N). Following pollen percentage recalculations, taxonomic homogenization, and age-depth model revision, the pollen spectra were interpolated at a 500-year resolution and a taxonomically harmonized and temporally standardized fossil pollen dataset established with 226 pollen taxa, covering the last 22 cal lea. Of the 271 pollen records, 85% were published since 1990, with reliable chronologies and high temporal resolutions; of these, 50% have raw data with complete pollen assemblages, ensuring the quality of this dataset The pollen records available for each 500-year time slice are well distributed over all main vegetation types and climatic zones of the study area, making their pollen spectra suitable for paleovegetation and paleoclimate research. Such a dataset can be used as an example for the development of similar datasets for other regions of the world.
The Pleistocene archeological record in East Africa has revealed unusual accumulations of Acheulean handaxes at prehistoric sites. In particular, there has been intensive debate concerning whether the artifact accumulation at the Middle Pleistocene Olorgesailie (Southern Kenya Rift) and Kariandusi (Central Kenya Rift) sites were a result of fluvial reworking or of in situ deposition by hominids. We used a two-step approach to test the hypothesis of fluvial reworking. Firstly, the behavior of handaxes in water currents was investigated in a current flume and the flow threshold required to reorientate the handaxes was determined. The results of these experiments suggested that, in relatively high energy and non-steady flow conditions, handaxes will reorientate themselves perpendicular to the current direction. Secondly, an automated image analysis routine was developed and applied to archeological plans from three Acheulean sites, two at Olorgesailie and one at Kariandusi, in order to determine the orientations of the handaxes. A Rayleigh test was then applied to the orientation data to test for a preferred orientation. The results revealed that the handaxes at the Upper Kariandusi Site and the Olorgesailie Main Site Mid Trench had a preferential orientation, suggesting reworking by a paleocurrent. The handaxes from the Olorgesailie Main Site H/6A, however, appeared to be randomly oriented and in situ deposition by the producers therefore remains a possibility.
In soils and sediments there is a strong coupling between local biogeochemical processes and the distribution of water, electron acceptors, acids and nutrients. Both sides are closely related and affect each other from small scale to larger scales. Soil structures such as aggregates, roots, layers or macropores enhance the patchiness of these distributions. At the same time it is difficult to access the spatial distribution and temporal dynamics of these parameter. Noninvasive imaging techniques with high spatial and temporal resolution overcome these limitations. And new non-invasive techniques are needed to study the dynamic interaction of plant roots with the surrounding soil, but also the complex physical and chemical processes in structured soils. In this study we developed an efficient non-destructive in-situ method to determine biogeochemical parameters relevant to plant roots growing in soil. This is a quantitative fluorescence imaging method suitable for visualizing the spatial and temporal pH changes around roots. We adapted the fluorescence imaging set-up and coupled it with neutron radiography to study simultaneously root growth, oxygen depletion by respiration activity and root water uptake. The combined set up was subsequently applied to a structured soil system to map the patchy structure of oxic and anoxic zones induced by a chemical oxygen consumption reaction for spatially varying water contents. Moreover, results from a similar fluorescence imaging technique for nitrate detection were complemented by a numerical modeling study where we used imaging data, aiming to simulate biodegradation under anaerobic, nitrate reducing conditions.
This study aims to establish, evaluate, and apply a modern pollen-climate transfer function from the transition zone between arctic tundra and light-needled taiga in Arctic Siberia. Lacustrine samples (n = 96) from the northern Siberian lowlands of Yakutia were collected along four north-to-south transects crossing the arctic forest line. Samples span a broad temperature and precipitation gradient (mean July temperature, T-July: 7.5-18.7 degrees C; mean annual precipitation, P-ann: 114-315 mm/yr). Redundancy analyses are used to examine the relationship between the modern pollen signal and corresponding vegetation types and climate. Performance of transfer functions for T-July and P-ann were cross-validated and tested for spatial autocorrelation effects. The root mean square errors of prediction are 1.67 degrees C for T-July and 40 mm/yr for P-ann. A climate reconstruction based on fossil pollen spectra from a Siberian Arctic lake sediment core spanning the Holocene yielded cold conditions for the Late Glacial (1-2 degrees C below present T-July). Warm and moist conditions were reconstructed for the early to mid Holocene (2 degrees C higher T-July than present), and climate conditions similar to modern ones were reconstructed for the last 4000 years. In conclusion, our modern pollen data set fills the gap of existing regional calibration sets with regard to the underrepresented Siberian tundra-taiga transition zone. The Holocene climate reconstruction indicates that the temperature deviation from modern values was only moderate despite the assumed Arctic sensitivity to present climate change.
Pseudotachylyte veins frequently associated with mylonites and ultramylonites occur within migmatitic paragneisses, metamonzodiorites, as well as felsic and mafic granulites at the base of the section of the Hercynian lower crust exposed in Calabria (Southern Italy). The crustal section is tectonically superposed on lower grade units. Ultramylonites and pseudotachylytes are particularly well developed in migmatitic paragneisses, whereas sparse fault-related pseudotachylytes and thin mylonite/ultramylonite bands occur in granulite-facies rocks. The presence of sillimanite and clinopyroxene in ultramylonites and mylonites indicates that relatively high-temperature conditions preceded the formation of pseudotachylytes. We have analysed pseudotachylytes from different rock types to ascertain their deep crustal origin and to better understand the relationships between brittle and ductile processes during deformation of the deeper crust. Different protoliths were selected to test how lithology controls pseudotachylyte composition and textures. In migmatites and felsic granulites, euhedral or cauliflower-shaped garnets directly crystallized from pseudotachylyte melts of near andesitic composition. This indicates that pseudotachylytes originated at deep crustal conditions (> 0.75 GPa). In mafic protoliths, quenched needle-to-feather-shaped high-alumina orthopyroxene occurs in contact with newly crystallized plagioclase. The pyroxene crystallizes in garnet-free and garnet-bearing veins. The simultaneous growth of orthopyroxene and plagioclase as well as almandine, suggests lower crustal origin, with pressures in excess of 0.85 GPa. The existence of melts of different composition in the same vein indicates the stepwise, non-equilibrium conditions of frictional melting. Melt formed and intruded into pre-existing anisotropies. In mafic granulites, brittle faulting is localized in a previously formed thin high-temperature mylonite bands. migmatitic gneisses are deformed into ultramylonite domains characterized by s-c fabric. Small grain size and fluids lowered the effective stress on the c planes favouring a seismic event and the consequent melt generation. Microstructures and ductile deformation of pseudotachylytes suggest continuous ductile flow punctuated by episodes of high-strain rate, leading to seismic events and melting.
Environmental parameters such as rainfall, temperature and relative humidity can affect the composition of higher plant leaf wax. The abundance and distribution of leaf wax biomarkers, such as long chain n-alkanes, in sedimentary archives have therefore been proposed as proxies reflecting climate change. However, a robust palaeoclimatic interpretation requires a thorough understanding of how environmental changes affect leaf wax n-alkane distributions in living plants. We have analysed the concentration and chain length distribution of leaf wax n-alkanes in Acacia and Eucalyptus species along a 1500 km climatic gradient in northern Australia that ranges from subtropical to arid. We show that aridity affected the concentration and distribution of n-alkanes for plants in both genera. For both Acacia and Eucalyptus n-alkane concentration increased by a factor of ten to the dry centre of Australia, reflecting the purpose of the wax in preventing water loss from the leaf. Furthermore, Acacian-alkanes decreased in average chain length (ACL) towards the arid centre of Australia, whereas Eucalyptus ACL increased under arid conditions. Our observations demonstrate that n-alkane concentration and distribution in leaf wax are sensitive to hydroclimatic conditions. These parameters could therefore potentially be employed in palaeorecords to estimate past environmental change. However, our finding of a distinct response of n-alkane ACL values to hydrological changes in different taxa also implies that the often assumed increase in ACL under drier conditions is not a robust feature for all plant species and genera and as such additional information about the prevalent vegetation are required when ACL values are used as a palaeoclimate proxy.
Convexities in the longitudinal profiles of actively incising rivers are typically considered to represent the morphologic signal of a transient response to external perturbations in tectonic or climatic forcing. Distinguishing such knickzones from those that may be anchored to the channel network by spatial variations in rock uplift, however, can be challenging. Here, we combine stream profile analysis, Be-10 watershed-averaged erosion rates, and numerical modeling of stream profile evolution to evaluate whether knickzones in the Abukuma massif of northeast Japan represent a temporal or spatial change in rock uplift rate in relation to forearc shortening. Knickzones in channels that drain the eastern flank of the Abukuma massif are characterized by breaks in slope-area scaling and separate low-gradient, alluvial upper-channel segments from high-gradient, deeply-incised lower channel segments. Average erosion rates inferred from Be-10 concentrations in modern sediment below knickzones exceed erosion rates above knickzones by 20-50%. Although profile convexities could be interpreted as a transient response to an increase in rock uplift rate associated with slip on the range-bounding fault, geologic constraints on the initiation of fault slip and the magnitude of displacement cannot be reconciled with a recent, spatially uniform increase in slip rate. Rather, we find that knickzone position, stream profile gradients, and basin averaged erosion rates are best explained by a relatively abrupt spatial increase in uplift rate localized above a flat-ramp transition in the fault system. These analyses highlight the importance of considering spatially non-uniform uplift in the interpretation of stream profile evolution and demonstrate that the adjustment of river profiles to fault displacement can provide constraints on fault geometry in actively eroding landscapes. (C) 2013 Elsevier B.V. All rights reserved.
We review the geodynamic evolution of the Aegean-Anatolia region and discuss strain localisation there over geological times. From Late Eocene to Present, crustal deformation in the Aegean backarc has localised progressively during slab retreat. Extension started with the formation of the Rhodope Metamorphic Core Complex (Eocene) and migrated to the Cyclades and the northern Menderes Massif (Oligocene and Miocene), accommodated by crustal-scale detachments and a first series of core complexes (MCCs). Extension then localised in Western Turkey, the Corinth Rift and the external Hellenic arc after Messinian times, while the North Anatolian Fault penetrated the Aegean Sea. Through time the direction and style of extension have not changed significantly except in terms of localisation. The contributions of progressive slab retreat and tearing, basal drag, extrusion tectonics and tectonic inheritance are discussed and we favour a model (I) where slab retreat is the main driving engine, (2) successive slab tearing episodes are the main causes of this stepwise strain localisation and (3) the inherited heterogeneity of the crust is a major factor for localising detachments. The continental crust has an inherited strong heterogeneity and crustal-scale contacts such as major thrust planes act as weak zones or as zones of contrast of resistance and viscosity that can localise later deformation. The dynamics of slabs at depth and the asthenospheric flow due to slab retreat also have influence strain localisation in the upper plate. Successive slab ruptures from the Middle Miocene to the late Miocene have isolated a narrow strip of lithosphere, still attached to the African lithosphere below Crete. The formation of the North Anatolian Fault is partly a consequence of this evolution. The extrusion of Anatolia and the Aegean extension are partly driven from below (asthenospheric flow) and from above (extrusion of a lid of rigid crust).
The occurrence of neritic microbial carbonates is often related to ecological refuges, where grazers and other competitors are reduced by environmental conditions, or to post-extinction events (e.g. in the Late Devonian, Early Triassic). Here, we present evidence for Middle Jurassic (Bajocian) microbial mounds formed in the normal marine, shallow neritic setting of an inner, ramp system from the High Atlas of Morocco. The microbial mounds are embedded in cross-bedded oolitic facies. Individual mounds show low relief domal geometries (up to 3 m high and 4.5 m across), but occasionally a second generation of mounds exhibits tabular geometries (<1 m high). The domes are circular in plan view and have intact tops, lacking evidence of current influence on mound preferred growth direction or distribution patterns, or truncation. The mound fades consists almost entirely of non-laminated, micritic thrombolites with branching morphologies and fine-grained, clotted and peloidal fabrics. Normal marine biota are present but infrequent. Several lines of evidence document that microbial mound growth alternates with time intervals of active ooid shoal deposition. This notion is of general significance when compared with modern Bahamian microbialites that co-exist with active sub-aquatic dunes. Furthermore, the lack of detailed studies of Middle Jurassic, normal marine shallow neritic microbial mounds adds a strong motivation for the present study. Specifically, Bajocian mounds formed on a firmground substratum during transgressive phases under condensed sedimentation. Furthermore, a transient increase in nutrient supply in the prevailing mesotrophic setting, as suggested by the heterotrophic-dominated biota, may have controlled microbial mound stages.
Agricultural soils have a high potential for sequestration of atmospheric carbon due to their volume and several promising management options. However, there is a remarkable lack of information about the status quo of organic carbon in agricultural soils. In this study a comprehensive data set of 384 cropland soils and 333 grassland soils within the state of Bavaria in southeast Germany was analyzed in order to provide representative information on total amount, regional distribution and driving parameters of soil organic carbon (SOC) and nitrogen (N) in agricultural soils of central Europe. The results showed that grassland soils stored higher amounts of SOC (11.8 kg m(-2)) and N (0.92 kg m(-2)) than cropland soils (9.0 and 0.66 kg m(-2), respectively) due to moisture-induced accumulation of soil organic matter (SOM) in B horizons. Surprisingly, no distinct differences were found for the A horizons since tillage led to a relocation of SOM with depth in cropland soils. Statistical analyses of driving factors for SOM storage revealed soil moisture, represented by the topographic wetness index (TWI), as the most important parameter for both cropland and grassland soils. Climate effects (mean annual temperature and precipitation) were of minor importance in agricultural soils because management options counteracted them to a certain extent, particularly in cropland soils. The distribution of SOC and N stocks within Bavaria based on agricultural regions confirmed the importance of soil moisture since the highest cropland SOC and N stocks were found for tertiary hills and loess regions, which exhibited large areas with potentially high soil moisture content in extant floodplains. Grassland soils showed the highest accumulation of SOC and N in the Alps and Pre-Alps as a result of low temperatures, high amounts of precipitation and high soil moisture content in areas of glacial denudation. Soil class was identified as a further driving parameter for SOC and N storage in cropland soils. In total, cropland and grassland soils in Bavaria store 242 and 134 Mt SOC as well as 19 and 12 Mt N down to a soil depth of 1 m or the parent material, respectively.
In this paper a change-point detection method is proposed by extending the singular spectrum transformation (SST) developed as one of the capabilities of singular spectrum analysis (SSA). The method uncovers change points related with trends and periodicities. The potential of the proposed method is demonstrated by analysing simple model time series including linear functions and sine functions as well as real world data (precipitation data in Kenya). A statistical test of the results is proposed based on a Monte Carlo simulation with surrogate methods. As a result, the successful estimation of change points as inherent properties in the representative time series of both trend and harmonics is shown. With regards to the application, we find change points in the precipitation data of Kenyan towns (Nakuru, Naivasha, Narok, and Kisumu) which coincide with the variability of the Indian Ocean Dipole (IOD) suggesting its impact of extreme climate in East Africa.
In probabilistic seismic-hazard analysis, epistemic uncertainties are commonly treated within a logic-tree framework in which the branch weights express the degree of belief of an expert in a set of models. For the calculation of the distribution of hazard curves, these branch weights represent subjective probabilities. A major challenge for experts is to provide logically consistent weight estimates (in the sense of Kolmogorovs axioms), to be aware of the multitude of heuristics, and to minimize the biases which affect human judgment under uncertainty. We introduce a platform-independent, interactive program enabling us to quantify, elicit, and transfer expert knowledge into a set of subjective probabilities by applying experimental design theory, following the approach of Curtis and Wood (2004). Instead of determining the set of probabilities for all models in a single step, the computer-driven elicitation process is performed as a sequence of evaluations of relative weights for small subsets of models. From these, the probabilities for the whole model set are determined as a solution of an optimization problem. The result of this process is a set of logically consistent probabilities together with a measure of confidence determined from the amount of conflicting information which is provided by the expert during the relative weighting process. We experiment with different scenarios simulating likely expert behaviors in the context of knowledge elicitation and show the impact this has on the results. The overall aim is to provide a smart elicitation technique, and our findings serve as a guide for practical applications.
Water management and environmental protection is vulnerable to extreme low flows during streamflow droughts. During the last decades, in most rivers of Central Europe summer runoff and low flows have decreased. Discharge projections agree that future decrease in runoff is likely for catchments in Brandenburg, Germany. Depending on the first-order controls on low flows, different adaption measures are expected to be appropriate. Small catchments were analyzed because they are expected to be more vulnerable to a changing climate than larger rivers. They are mainly headwater catchments with smaller ground water storage. Local characteristics are more important at this scale and can increase vulnerability. This thesis mutually evaluates potential adaption measures to sustain minimum runoff in small catchments of Brandenburg, Germany, and similarities of these catchments regarding low flows. The following guiding questions are addressed: (i) Which first-order controls on low flows and related time scales exist? (ii) Which are the differences between small catchments regarding low flow vulnerability? (iii) Which adaption measures to sustain minimum runoff in small catchments of Brandenburg are appropriate considering regional low flow patterns? Potential adaption measures to sustain minimum runoff during periods of low flows can be classified into three categories: (i) increase of groundwater recharge and subsequent baseflow by land use change, land management and artificial ground water recharge, (ii) increase of water storage with regulated outflow by reservoirs, lakes and wetland water management and (iii) regional low flow patterns have to be considered during planning of measures with multiple purposes (urban water management, waste water recycling and inter-basin water transfer). The question remained whether water management of areas with shallow groundwater tables can efficiently sustain minimum runoff. Exemplary, water management scenarios of a ditch irrigated area were evaluated using the model Hydrus-2D. Increasing antecedent water levels and stopping ditch irrigation during periods of low flows increased fluxes from the pasture to the stream, but storage was depleted faster during the summer months due to higher evapotranspiration. Fluxes from this approx. 1 km long pasture with an area of approx. 13 ha ranged from 0.3 to 0.7 l\s depending on scenario. This demonstrates that numerous of such small decentralized measures are necessary to sustain minimum runoff in meso-scale catchments. Differences in the low flow risk of catchments and meteorological low flow predictors were analyzed. A principal component analysis was applied on daily discharge of 37 catchments between 1991 and 2006. Flows decreased more in Southeast Brandenburg according to meteorological forcing. Low flow risk was highest in a region east of Berlin because of intersection of a more continental climate and the specific geohydrology. In these catchments, flows decreased faster during summer and the low flow period was prolonged. A non-linear support vector machine regression was applied to iteratively select meteorological predictors for annual 30-day minimum runoff in 16 catchments between 1965 and 2006. The potential evapotranspiration sum of the previous 48 months was the most important predictor (r²=0.28). The potential evapotranspiration of the previous 3 months and the precipitation of the previous 3 months and last year increased model performance (r²=0.49, including all four predictors). Model performance was higher for catchments with low yield and more damped runoff. In catchments with high low flow risk, explanatory power of long term potential evapotranspiration was high. Catchments with a high low flow risk as well as catchments with a considerable decrease in flows in southeast Brandenburg have the highest demand for adaption. Measures increasing groundwater recharge are to be preferred. Catchments with high low flow risk showed relatively deep and decreasing groundwater heads allowing increased groundwater recharge at recharge areas with higher altitude away from the streams. Low flows are expected to stay low or decrease even further because long term potential evapotranspiration was the most important low flow predictor and is projected to increase during climate change. Differences in low flow risk and runoff dynamics between catchments have to be considered for management and planning of measures which do not only have the task to sustain minimum runoff.
Scarcity of hydrological data, especially streamflow discharge and groundwater level series, restricts the understanding of channel transmission losses (TL) in drylands. Furthermore, the lack of information on spatial river dynamics encompasses high uncertainty on TL analysis in large rivers. The objective of this study was to combine the information from streamflow and groundwater level series with multi-temporal satellite data to derive a hydrological concept of TL for a reach of the Middle Jaguaribe River (MJR) in semi-arid north-eastern Brazil. Based on this analysis, we proposed strategies for its modelling and simulation. TL take place in an alluvium, where river and groundwater can be considered to be hydraulically connected. Most losses certainly infiltrated only through streambed and levees and not through the flood plains, as could be shown by satellite image analysis. TL events whose input river flows were smaller than a threshold did not reach the outlet of the MJR. TL events whose input flows were higher than this threshold reached the outlet losing on average 30% of their input. During the dry seasons (DS) and at the beginning of rainy seasons (DS/BRS), no river flow is expected for pre-events, and events have vertical infiltration into the alluvium. At the middle and the end of the rainy seasons (MRS/ERS), river flow sustained by base flow occurs before/after events, and lateral infiltration into the alluvium plays a major role. Thus, the MJR shifts from being a losing river at DS/BRS to become a losing/gaining (mostly losing) river at MRS/ERS. A model of this system has to include the coupling of river and groundwater flow processes linked by a leakage approach.
Even though quite different in occurrence and consequences, from a modeling perspective many natural hazards share similar properties and challenges. Their complex nature as well as lacking knowledge about their driving forces and potential effects make their analysis demanding: uncertainty about the modeling framework, inaccurate or incomplete event observations and the intrinsic randomness of the natural phenomenon add up to different interacting layers of uncertainty, which require a careful handling. Nevertheless deterministic approaches are still widely used in natural hazard assessments, holding the risk of underestimating the hazard with disastrous effects. The all-round probabilistic framework of Bayesian networks constitutes an attractive alternative. In contrast to deterministic proceedings, it treats response variables as well as explanatory variables as random variables making no difference between input and output variables. Using a graphical representation Bayesian networks encode the dependency relations between the variables in a directed acyclic graph: variables are represented as nodes and (in-)dependencies between variables as (missing) edges between the nodes. The joint distribution of all variables can thus be described by decomposing it, according to the depicted independences, into a product of local conditional probability distributions, which are defined by the parameters of the Bayesian network. In the framework of this thesis the Bayesian network approach is applied to different natural hazard domains (i.e. seismic hazard, flood damage and landslide assessments). Learning the network structure and parameters from data, Bayesian networks reveal relevant dependency relations between the included variables and help to gain knowledge about the underlying processes. The problem of Bayesian network learning is cast in a Bayesian framework, considering the network structure and parameters as random variables itself and searching for the most likely combination of both, which corresponds to the maximum a posteriori (MAP score) of their joint distribution given the observed data. Although well studied in theory the learning of Bayesian networks based on real-world data is usually not straight forward and requires an adoption of existing algorithms. Typically arising problems are the handling of continuous variables, incomplete observations and the interaction of both. Working with continuous distributions requires assumptions about the allowed families of distributions. To "let the data speak" and avoid wrong assumptions, continuous variables are instead discretized here, thus allowing for a completely data-driven and distribution-free learning. An extension of the MAP score, considering the discretization as random variable as well, is developed for an automatic multivariate discretization, that takes interactions between the variables into account. The discretization process is nested into the network learning and requires several iterations. Having to face incomplete observations on top, this may pose a computational burden. Iterative proceedings for missing value estimation become quickly infeasible. A more efficient albeit approximate method is used instead, estimating the missing values based only on the observations of variables directly interacting with the missing variable. Moreover natural hazard assessments often have a primary interest in a certain target variable. The discretization learned for this variable does not always have the required resolution for a good prediction performance. Finer resolutions for (conditional) continuous distributions are achieved with continuous approximations subsequent to the Bayesian network learning, using kernel density estimations or mixtures of truncated exponential functions. All our proceedings are completely data-driven. We thus avoid assumptions that require expert knowledge and instead provide domain independent solutions, that are applicable not only in other natural hazard assessments, but in a variety of domains struggling with uncertainties.
Seismicity induced by coal mining in the Ruhr region, Germany, has been monitored continuously over the last 25 yr. In 2006, a dense temporary network (HAMNET) was deployed to locally monitor seismicity induced by longwall mining close to the town of Hamm. Between 2006 July and 2007 July, more than 7000 events with magnitudes M-L from -1.7 to 2.0 were detected. The spatiotemporal distribution of seismicity shows high correlation with the mining activity. In order to monitor rupture processes, we set up an automated source inversion routine and successfully perform double couple and full moment tensor (MT) inversions for more than 1000 events with magnitudes above M-L -0.5. The source inversion is based on a full waveform approach, both in the frequency and in the time domain, providing information about the centroid location, focal mechanism, scalar moment and full MT. Inversion results indicate a strong dominance of normal faulting focal mechanisms, with a steeper plane and a subhorizontal one. Fault planes are oriented parallel to the mining stopes. We classify the focal mechanisms based on their orientation and observe different frequency-magnitude distributions for families of events with different focal mechanisms; the overall frequency-magnitude distribution is not fitting the Gutenberg-Richter relation. Full MTs indicate that non-negligible opening tensile components accompanied normal faulting source mechanisms. Finally, extended source models are investigated for largest events. Results suggest that the rupture processes mostly occurred along the subvertical planes.
Automated seismic event location by travel-time stacking an application to mining induced seismicity
(2013)
Antarctic glacier forfields are extreme environments and pioneer sites for ecological succession. The Antarctic continent shows microbial community development as a natural laboratory because of its special environment, geographic isolation and little anthropogenic influence. Increasing temperatures due to global warming lead to enhanced deglaciation processes in cold-affected habitats and new terrain is becoming exposed to soil formation and accessible for microbial colonisation. This study aims to understand the structure and development of glacier forefield bacterial communities, especially how soil parameters impact the microorganisms and how those are adapted to the extreme conditions of the habitat. To this effect, a combination of cultivation experiments, molecular, geophysical and geochemical analysis was applied to examine two glacier forfields of the Larsemann Hills, East Antarctica. Culture-independent molecular tools such as terminal restriction length polymorphism (T-RFLP), clone libraries and quantitative real-time PCR (qPCR) were used to determine bacterial diversity and distribution. Cultivation of yet unknown species was carried out to get insights in the physiology and adaptation of the microorganisms. Adaptation strategies of the microorganisms were studied by determining changes of the cell membrane phospholipid fatty acid (PLFA) inventory of an isolated bacterium in response to temperature and pH fluctuations and by measuring enzyme activity at low temperature in environmental soil samples. The two studied glacier forefields are extreme habitats characterised by low temperatures, low water availability and small oligotrophic nutrient pools and represent sites of different bacterial succession in relation to soil parameters. The investigated sites showed microbial succession at an early step of soil formation near the ice tongue in comparison to closely located but rather older and more developed soil from the forefield. At the early step the succession is influenced by a deglaciation-dependent areal shift of soil parameters followed by a variable and prevalently depth-related distribution of the soil parameters that is driven by the extreme Antarctic conditions. The dominant taxa in the glacier forefields are Actinobacteria, Acidobacteria, Proteobacteria, Bacteroidetes, Cyanobacteria and Chloroflexi. The connection of soil characteristics with bacterial community structure showed that soil parameter and soil formation along the glacier forefield influence the distribution of certain phyla. In the early step of succession the relative undifferentiated bacterial diversity reflects the undifferentiated soil development and has a high potential to shift according to past and present environmental conditions. With progressing development environmental constraints such as water or carbon limitation have a greater influence. Adapting the culturing conditions to the cold and oligotrophic environment, the number of culturable heterotrophic bacteria reached up to 108 colony forming units per gram soil and 148 isolates were obtained. Two new psychrotolerant bacteria, Herbaspirillum psychrotolerans PB1T and Chryseobacterium frigidisoli PB4T, were characterised in detail and described as novel species in the family of Oxalobacteraceae and Flavobacteriaceae, respectively. The isolates are able to grow at low temperatures tolerating temperature fluctuations and they are not specialised to a certain substrate, therefore they are well-adapted to the cold and oligotrophic environment. The adaptation strategies of the microorganisms were analysed in environmental samples and cultures focussing on extracellular enzyme activity at low temperature and PLFA analyses. Extracellular phosphatases (pH 11 and pH 6.5), β-glucosidase, invertase and urease activity were detected in the glacier forefield soils at low temperature (14°C) catalysing the conversion of various compounds providing necessary substrates and may further play a role in the soil formation and total carbon turnover of the habitat. The PLFA analysis of the newly isolated species C. frigidisoli showed that the cold-adapted strain develops different strategies to maintain the cell membrane function under changing environmental conditions by altering the PLFA inventory at different temperatures and pH values. A newly discovered fatty acid, which was not found in any other microorganism so far, significantly increased at decreasing temperature and low pH and thus plays an important role in the adaption of C. frigidisoli. This work gives insights into the diversity, distribution and adaptation mechanisms of microbial communities in oligotrophic cold-affected soils and shows that Antarctic glacier forefields are suitable model systems to study bacterial colonisation in connection to soil formation.
Die Ziele der Energiewende sind ehrgeizig. Der Vortrag zeigt, welche Nutzungsoptionen der Untergrund bietet - z.B. geothermische Energiegewinnung oder geologische Speicherung. Für eine gesellschaftsweite, sachliche Diskussion werden konkrete Zahlen nicht nur zu den Chancen, sondern auch zu den Risiken benötigt.
In undisturbed tropical montane rainforests massive organic layers accommodate the majority of roots and only a small fraction of roots penetrate the mineral soil. We investigated the contribution of vegetation to slope stability in such environments by modifying a standard model for slope stability to include an organic layer with distinct mechanical properties. The importance of individual model parameters was evaluated using detailed measurements of soil and vegetation properties to reproduce the observed depth of 11 shallow landslides in the Andes of southern Ecuador. By distinguishing mineral soil, organic layer and above-ground biomass, it is shown that in this environment vegetation provides a destabilizing effect mainly due to its contribution to the mass of the organic layer (up to 973 t ha-1 under wet conditions). Sensitivity analysis shows that the destabilizing effect of the mass of soil and vegetation can only be effective on slopes steeper than 37.9 degrees. This situation applies to 36% of the study area. Thus, on the steep slopes of this megadiverse ecosystem, the mass of the growing forest promotes landsliding, which in turn promotes a new cycle of succession. This feedback mechanism is worth consideration in further investigations of the impact of landslides on plant diversity in similar environments.
The Upper Cretaceous shallow-water carbonates of the Pyrenean Basin (NE Spain) host rich and diverse larger foraminiferal associations which witness the recovery of this group of protozoans after the dramatic extinction of the Cenomanian-Turonian boundary interval. In this paper a new, large discoidal porcelaneous foraminifer, Broeckina gassoensis sp. nov., is described from the middle Coniacian shallow-water deposits of the Collada Gasso Formation, in the Boixols Thrust Sheet. This is the first complex porcelaneous larger foraminifer of the Late Cretaceous global community maturation cycle recorded in the Pyrenean bioprovince. It differs from the late Santonian-early Campanian B. dufrenoyi for its smaller size in A and B generations and the less developed endoskeleton, which shows short septula. Broeckina gassoensis sp. nov. has been widely employed as a stratigraphic marker in the regional geological literature, under the name of "Broeckina", but its age was so far controversial. Its middle Coniacian age (lowermost part of the Peroniceras tridorsatum ammonite zone), established in this paper by strontium isotope stratigraphy, indicates that it took about 5 My after the Cenomanian-Turonian boundary crisis to re-evolve the complex test architecture of larger foraminifera, which is functional to their relation with photosymbiotic algae and K-strategy.
Surface lake sediments, 28 from Hoh Xil, 24 from northeastern China, 99 from Lake Bosten, 31 from Ulungur and 26 from Heihai were collected to determine C-13 and O-18 values. Considering the impact factors, conductivity, alkalinity, pH, TOC, C/N and carbonate-content in the sediments, Cl, P, S, and metal element ratios of Mg/Ca, Sr/Ca, Fe/Mn of bulk sediments as environmental variables enable evaluation of their influences on C-13 and O-18 using principal component analysis (PCA) method. The closure and residence time of lakes can influence the correlation between C-13 and O-18. Lake water will change from fresh to brackish with increasing reduction and eutrophication effects. Mg/Ca in the bulk sediment indicates the characteristic of residence time, Sr/Ca and Fe/Mn infer the salinity of lakes. Carbonate formation processes and types can influence the C-13-O-18 correlation. O-18 will be heavier from Mg-calcite and aragonite formed in a high-salinity water body than calcite formed in freshwater conditions. When carbonate content is less than 30%, there is no relationship with either C-13 or O-18, and also none between C-13 and O-18. More than 30%, carbonate content, however, co-varies highly to C-13 and O-18, and there is also a high correlation between C-13 and O-18. Vegetation conditions and primary productivity of lakes can influence the characteristics of C-13 and O-18, and their co-variance. Total organic matter content (TOC) in the sediments is higher with more terrestrial and submerged plants infilling. In northeastern and northwestern China, when organic matter in the lake sediments comes from endogenous floating organisms and algae, the C-13 value is high. C-13 is in the range of -4%o to 0 parts per thousand when organic matter comes mainly from floating organisms (C/N<6); in the range of -4 parts per thousand to 8 parts per thousand when organic matter comes from diatoms (C/N=6 to 8); and -8 parts per thousand to -4 parts per thousand when organic matter comes from aquatic and terrestrial plants (C/N>8).
Natural and human-induced erosion supplies high amounts of soil organic carbon (OC) to terrestrial drainage networks. Yet OC fluxes in rivers were considered in global budgets only recently. Modern estimates of annual carbon burial in inland river sediments of 0.6 Gt C, or 22% of C transferred from terrestrial ecosystems to river channels, consider only lakes and reservoirs and disregard any long-term carbon burial in hillslope or floodplain sediments. Here we present the first assessment of sediment-bound OC storage in Central Europe from a synthesis of similar to 1500 Holocene hillslope and floodplain sedimentary archives. We show that sediment storage increases with drainage-basin size due to more extensive floodplains in larger river basins. However, hillslopes retain hitherto unrecognized high amounts of eroded soils at the scale of large river basins such that average agricultural erosion rates during the Holocene would have been at least twice as high as reported previously. This anthropogenic hillslope sediment storage exceeds floodplain storage in drainage basins <10(5) km(2), challenging the notion that floodplains are the dominant sedimentary sinks. In terms of carbon burial, OC concentrations in floodplains exceed those on hillslopes, and net OC accumulation rates in floodplains (0.70.2 g C m(-2)a(-1)) surpass those on hillslopes (0.40.1 g C m(-2)a(-1)) over the last 7500 years. We conclude that carbon burial in floodplains and on hillslopes in Central Europe exceeds terrestrial carbon storage in lakes and reservoirs by at least 2 orders of magnitude and should thus be considered in continental carbon budgets.
Within a research project about future sustainable water management options in the Elbe River basin, quasi-natural discharge scenarios had to be provided. The semi-distributed eco-hydrological model SWIM was utilised for this task. According to scenario simulations driven by the stochastical climate model STAR, the region would get distinctly drier. However, this thesis focuses on the challenge of meeting the requirement of high model fidelity even for smaller sub-basins. Usually, the quality of the simulations is lower at inner points than at the outlet. Four research paper chapters and the discussion chapter deal with the reasons for local model deviations and the problem of optimal spatial calibration. Besides other assessments, the Markov Chain Monte Carlo method is applied to show whether evapotranspiration or precipitation should be corrected to minimise runoff deviations, principal component analysis is used in an unusual way to evaluate local precipitation alterations by land cover changes, and remotely sensed surface temperatures allow for an independent view on the evapotranspiration landscape. The overall insight is that spatially explicit hydrological modelling of such a large river basin requires a lot of local knowledge. It probably needs more time to obtain such knowledge as is usually provided for hydrological modelling studies.
Changes in rainfall interception along a secondary forest succession gradient in lowland Panama
(2013)
Secondary forests are rapidly expanding in tropical regions. Yet, despite the importance of understanding the hydrological consequences of land-cover dynamics, the relationship between forest succession and canopy interception is poorly understood. This lack of knowledge is unfortunate because rainfall interception plays an important role in regional water cycles and needs to be quantified for many modeling purposes. To help close this knowledge gap, we designed a throughfall monitoring study along a secondary succession gradient in a tropical forest region of Panama. The investigated gradient comprised 20 forest patches 3 to 130 yr old. We sampled each patch with a minimum of 20 funnel-type throughfall collectors over a continuous 2month period that had nearly 900 mm of rain. During the same period, we acquired forest inventory data and derived several forest structural attributes. We then applied simple and multiple regression models (Bayesian model averaging, BMA) and identified those vegetation parameters that had the strongest influence on the variation of canopy interception. Our analyses yielded three main findings. First, canopy interception changed rapidly during forest succession. After only a decade, throughfall volumes approached levels that are typical for mature forests. Second, a parsimonious (simple linear regression) model based on the ratio of the basal area of small stems to the total basal area outperformed more complex multivariate models (BMA approach). Third, based on complementary forest inventory data, we show that the influence of young secondary forests on interception in realworld fragmented landscapes might be detectable only in regions with a substantial fraction of young forests. Our re-sults suggest that where entire catchments undergo forest regrowth, initial stages of succession may be associated with a substantial decrease of streamflow generation. Our results further highlight the need to study hydrological processes in all forest succession stages, including early ones.
Nowadays, an increasing amount of seismic data is collected by daily observatory routines. The basic step for successfully analyzing those data is the correct detection of various event types. However, the visually scanning process is a time-consuming task. Applying standard techniques for detection like the STA/LTAtrigger still requires the manual control for classification. Here, we present a useful alternative. The incoming data stream is scanned automatically for events of interest. A stochastic classifier, called hidden Markov model, is learned for each class of interest enabling the recognition of highly variable waveforms. In contrast to other automatic techniques as neural networks or support vector machines the algorithm allows to start the classification from scratch as soon as interesting events are identified. Neither the tedious process of collecting training samples nor a time-consuming configuration of the classifier is required. An approach originally introduced for the volcanic task force action allows to learn classifier properties from a single waveform example and some hours of background recording. Besides a reduction of required workload this also enables to detect very rare events. Especially the latter feature provides a milestone point for the use of seismic devices in alpine warning systems. Furthermore, the system offers the opportunity to flag new signal classes that have not been defined before. We demonstrate the application of the classification system using a data set from the Swiss Seismological Survey achieving very high recognition rates. In detail we document all refinements of the classifier providing a step-by-step guide for the fast set up of a well-working classification system.
Die Anpassung von Sektoren an veränderte klimatische Bedingungen erfordert ein Verständnis von regionalen Vulnerabilitäten. Vulnerabilität ist als Funktion von Sensitivität und Exposition, welche potentielle Auswirkungen des Klimawandels darstellen, und der Anpassungsfähigkeit von Systemen definiert. Vulnerabilitätsstudien, die diese Komponenten quantifizieren, sind zu einem wichtigen Werkzeug in der Klimawissenschaft geworden. Allerdings besteht von der wissenschaftlichen Perspektive aus gesehen Uneinigkeit darüber, wie diese Definition in Studien umgesetzt werden soll. Ausdiesem Konflikt ergeben sich viele Herausforderungen, vor allem bezüglich der Quantifizierung und Aggregierung der einzelnen Komponenten und deren angemessenen Komplexitätsniveaus. Die vorliegende Dissertation hat daher zum Ziel die Anwendbarkeit des Vulnerabilitätskonzepts voranzubringen, indem es in eine systematische Struktur übersetzt wird. Dies beinhaltet alle Komponenten und schlägt für jede Klimaauswirkung (z.B. Sturzfluten) eine Beschreibung des vulnerablen Systems vor (z.B. Siedlungen), welches direkt mit einer bestimmten Richtung eines relevanten klimatischen Stimulus in Verbindung gebracht wird (z.B. stärkere Auswirkungen bei Zunahme der Starkregentage). Bezüglich der herausfordernden Prozedur der Aggregierung werden zwei alternative Methoden, die einen sektorübergreifenden Überblick ermöglichen, vorgestellt und deren Vor- und Nachteile diskutiert. Anschließend wird die entwickelte Struktur einer Vulnerabilitätsstudie mittels eines indikatorbasierten und deduktiven Ansatzes beispielhaft für Gemeinden in Nordrhein-Westfalen in Deutschland angewandt. Eine Übertragbarkeit auf andere Regionen ist dennoch möglich. Die Quantifizierung für die Gemeinden stützt sich dabei auf Informationen aus der Literatur. Da für viele Sektoren keine geeigneten Indikatoren vorhanden waren, werden in dieser Arbeit neue Indikatoren entwickelt und angewandt, beispielsweise für den Forst- oder Gesundheitssektor. Allerdings stellen fehlende empirische Daten bezüglich relevanter Schwellenwerte eine Lücke dar, beispielsweise welche Stärke von Klimaänderungen eine signifikante Auswirkung hervorruft. Dies führt dazu, dass die Studie nur relative Aussagen zum Grad der Vulnerabilität jeder Gemeinde im Vergleich zum Rest des Bundeslandes machen kann. Um diese Lücke zu füllen, wird für den Forstsektor beispielhaft die heutige und zukünftige Sturmwurfgefahr von Wäldern berechnet. Zu diesem Zweck werden die Eigenschaften der Wälder mit empirischen Schadensdaten eines vergangenen Sturmereignisses in Verbindung gebracht. Der sich daraus ergebende Sensitivitätswert wird anschließend mit den Windverhältnissen verknüpft. Sektorübergreifende Vulnerabilitätsstudien erfordern beträchtliche Ressourcen, was oft deren Anwendbarkeit erschwert. In einem nächsten Schritt wird daher das Potential einer Vereinfachung der Komplexität anhand zweier sektoraler Beispiele untersucht. Um das Auftreten von Waldbränden vorherzusagen, stehen zahlreiche meteorologische Indices zur Verfügung, welche eine Spannbreite unterschiedlicher Komplexitäten aufweisen. Bezüglich der Anzahl monatlicher Waldbrände weist die relative Luftfeuchtigkeit für die meisten deutschen Bundesländer eine bessere Vorhersagekraft als komplexere Indices auf. Dies ist er Fall, obgleich sie selbst als Eingangsvariable für die komplexeren Indices verwendet wird. Mit Hilfe dieses einzelnen meteorologischen Faktors kann also die Waldbrandgefahr in deutschen Region ausreichend genau ausgedrückt werden, was die Ressourceneffizienz von Studien erhöht. Die Methodenkomplexität wird auf ähnliche Weise hinsichtlich der Anwendung des ökohydrologischen Modells SWIM für die Region Brandenburg untersucht. Die interannuellen Bodenwasserwerte, welche durch dieses Modell simuliert werden, können nur unzureichend durch ein einfacheres statistisches Modell, welches auf denselben Eingangsdaten aufbaut, abgebildet werden. Innerhalb eines Zeithorizonts von Jahrzehnten, kann der statistische Ansatz jedoch das Bodenwasser zufriedenstellend abbilden und zeigt eine Dominanz der Bodeneigenschaft Feldkapazität. Dies deutet darauf hin, dass die Komplexität im Hinblick auf die Anzahl der Eingangsvariablen für langfristige Berechnungen reduziert werden kann. Allerdings sind die Aussagen durch fehlende beobachtete Bodenwasserwerte zur Validierung beschränkt. Die vorliegenden Studien zur Vulnerabilität und ihren Komponenten haben gezeigt, dass eine Anwendung noch immer wissenschaftlich herausfordernd ist. Folgt man der hier verwendeten Vulnerabilitätsdefinition, treten zahlreiche Probleme bei der Implementierung in regionalen Studien auf. Mit dieser Dissertation wurden Fortschritte bezüglich der aufgezeigten Lücken bisheriger Studien erzielt, indem eine systematische Struktur für die Beschreibung und Aggregierung von Vulnerabilitätskomponenten erarbeitet wurde. Hierfür wurden mehrere Ansätze diskutiert, die jedoch Vor- und Nachteile besitzen. Diese sollten vor der Anwendung von zukünftigen Studien daher ebenfalls sorgfältig abgewogen werden. Darüber hinaus hat sich gezeigt, dass ein Potential besteht einige Ansätze zu vereinfachen, jedoch sind hierfür weitere Untersuchungen nötig. Insgesamt konnte die Dissertation die Anwendung von Vulnerabilitätsstudien als Werkzeug zur Unterstützung von Anpassungsmaßnahmen stärken.
The terrestrial biosphere is a key component of the global carbon cycle and its carbon balance is strongly influenced by climate. Continuing environmental changes are thought to increase global terrestrial carbon uptake. But evidence is mounting that climate extremes such as droughts or storms can lead to a decrease in regional ecosystem carbon stocks and therefore have the potential to negate an expected increase in terrestrial carbon uptake. Here we explore the mechanisms and impacts of climate extremes on the terrestrial carbon cycle, and propose a pathway to improve our understanding of present and future impacts of climate extremes on the terrestrial carbon budget.
The geomorphic evolution of the Makran Range of SE-Iran and SW-Pakistan has been controlled by the prevailing SW-Asian monsoon and Mediterranean winter rainfall climate and the surface uplift processes resulting from the Arabia-Eurasia collision. The impact of climate on Quaternary fluvial and alluvial sequence formation and their regional correlation has been little investigated due to limited age control of these sequences. Using Be-10 cosmogenic nuclide exposure ages we established a Middle to Late Pleistocene terrace chronology. Our record tentatively indicates that terrace levels were abandoned towards the transition to or during warmer/pluvial periods (interglacials and/or interstadials) back to Marine Isotope Stage (MIS) 7, but abandoned ages show a large spread. It is hypothesized that pluvial phases correspond with times of enhanced SW-monsoons and a northward shift of the Intertropical Convergence Zone (ITCZ). Furthermore, orbital periodidties can be deduced on frequencies related to obliquity and precession cycles. Overall, caution has to be placed in sampling and interpreting alluvial deposits, which may have complex inheritance patterns and spatially and temporarily variable catchment erosion histories and terrace-channel dynamics.
Beside the dominant climate control on terrace formation, elevated channel steepness indices around major thrusts and numerous knickpoints indicate an additionally tectonic influence on terrace formation. Local incision rates (mean similar to 0.6-0.8 min.a(-1)) are variable in space and time but are similar to uplift rates obtained from coastal terraces and thus suggest a regional surface uplift. (C) 2013 Elsevier B.V. All rights reserved.
Black shales are sedimentary rocks with a high content of organic carbon, which leads to a dark grayish to black color. Due to their potential to contain oil or gas, black shales are of great interest for the support of the worldwide energy supply. An integrated seismic investigation of the Lower Palaeozoic black shales was carried out at the Danish island Bornholm to locate the shallow-lying Alum Shale layer and its surrounding formations and to characterize its potential as a source rock. Therefore, two seismic experiments at a total of three crossing profiles were carried out in October 2010 and in June 2012 in the southern part of the island. Two different active measurements were conducted with either a weight drop source or a minivibrator. Additionally, the ambient noise field was recorded at the study location over a time interval of about one day, and also a laboratory analysis of borehole samples was carried out. The seismic profiles were positioned as close as possible to two scientific boreholes which were used for comparative purposes. The seismic field data was analyzed with traveltime tomography, surface wave inversion and seismic interferometry to obtain the P-wave and S-wave velocity models of the subsurface. The P-wave velocity models which were determined for all three profiles clearly locate the Alum Shale layer between the Komstad Limestone layer on top and the Læså Sandstone Formation at the base of the models. The black shale layer has P-wave velocities around 3 km/s which are lower compared to the adjacent formations. Due to a very good agreement of the sonic log and the vertical velocity profiles of the two seismic lines, which are directly crossing the borehole where the sonic log was conducted, the reliability of the traveltime tomography is proven. A correlation of the seismic velocities with the content of organic carbon is an important task for the characterization of the reservoir properties of a black shale formation. It is not possible without calibration but in combination with a full 2D tomographic image of the subsurface it gives the subsurface distribution of the organic material. The S-wave model obtained with surface wave inversion of the vibroseis data of one of the profiles images the Alum Shale layer also very well with S-wave velocities around 2 km/s. Although individual 1D velocity models for each of the source positions were determined, the subsurface S-wave velocity distribution is very uniform with a good match between the single models. A really new approach described here is the application of seismic interferometry to a really small study area and a quite short time interval. Also new is the selective procedure of only using time windows with the best crosscorrelation signals to achieve the final interferograms. Due to the small scale of the interferometry even P-wave signals can be observed in the final crosscorrelations. In the laboratory measurements the seismic body waves were recorded for different pressure and temperature stages. Therefore, samples of different depths of the Alum Shale were available from one of the scientific boreholes at the study location. The measured velocities have a high variance with changing pressure or temperature. Recordings with wave propagation both parallel and perpendicular to the bedding of the samples reveal a great amount of anisotropy for the P-wave velocity, whereas the S-wave velocity is almost independent of the wave direction. The calculated velocity ratio is also highly anisotropic with very low values for the perpendicular samples and very high values for the parallel ones. Interestingly, the laboratory velocities of the perpendicular samples are comparable to the velocities of the field experiments indicating that the field measurements are sensitive to wave propagation in vertical direction. The velocity ratio is also calculated with the P-wave and S-wave velocity models of the field experiments. Again, the Alum Shale can be clearly separated from the adjacent formations because it shows overall very low vP/vS ratios around 1.4. The very low velocity ratio indicates the content of gas in the black shale formation. With the combination of all the different methods described here, a comprehensive interpretation of the seismic response of the black shale layer can be made and the hydrocarbon source rock potential can be estimated.
Retention and transformation of nutrients within a river catchment are important mechanisms influencing water quality measured at the watershed outlet. Nutrient storage and reduction can occur in soils as well as in the river and should be considered in water quality modelling. Consideration is possible using various methods at several points during modelling cascade. The study compares the effects of five different equation sets implemented into the Soil and Water Integrated Model (SWIM), one describing terrestrial and four in-stream retention with a rising complexity (including algal growth and death at the highest complexity level). The influences of the different methods alone and in combinations on water quality model outputs (NO3-N, NH4-N, PO4-P) were analyzed for the outlet of the large-scale Saale basin in Germany. Experiments revealed that nutrient forms coming primarily from diffuse sources are mostly influenced by retention processes in the soils of the catchment, and river processes are less important. Nutrients introduced to the river mainly by point sources are more subject to retention by in-stream processes, but both nutrient retention and transformation processes in soils and rivers have to be included. Although the best overall results could be achieved at the highest complexity level, the calibration efforts for this case are extremely high, and only minor improvements of overall model performance with the highest complexity were detected. Therefore, it could be reasoned that for some research questions also less complex model approaches would be sufficient, which could help to reduce unnecessary complexity and diminish high uncertainty in water quality modelling at the catchment scale. (C) 2013 Elsevier B.V. All rights reserved.
Finding evidence for ultrahigh-pressure (UHP) metamorphism in crustal rocks is far from straightforward. The index minerals coesite and diamond are incredibly inconspicuous and are therefore difficult to use as UHP prospecting tools. Consequently, petrographers rely on recognizing subtle breakdown microstructures that result from pressure release during the return to the surface of the once deeply buried rock. Similarly, many other UHP minerals are first suspected on the basis of typical reaction or exsolution microstructures. Thus, the painstaking use of microscopic techniques has been fundamental to the tremendous advances in characterizing, quantifying, and understanding macroscopic-scale, deep continental subduction, rapid exhumation, and mountain-building processes.
The deep thermal field in sedimentary basins can be affected by convection, conduction or both resulting from the structural inventory, physical properties of geological layers and physical processes taking place therein. For geothermal energy extraction, the controlling factors of the deep thermal field need to be understood to delineate favorable drill sites and exploitation compartments. We use geologically based 3-D finite element simulations to figure out the geologic controls on the thermal field of the geothermal research site Gro Schonebeck located in the E part of the North German Basin. Its target reservoir consists of Permian Rotliegend clastics that compose the lower part of a succession of Late Carboniferous to Cenozoic sediments, subdivided into several aquifers and aquicludes. The sedimentary succession includes a layer of mobilized Upper Permian Zechstein salt which plays a special role for the thermal field due to its high thermal conductivity. Furthermore, the salt is impermeable and due to its rheology decouples the fault systems in the suprasalt units from subsalt layers. Conductive and coupled fluid and heat transport simulations are carried out to assess the relative impact of different heat transfer mechanisms on the temperature distribution. The measured temperatures in 7 wells are used for model validation and show a better fit with models considering fluid and heat transport than with a purely conductive model. Our results suggest that advective and convective heat transport are important heat transfer processes in the suprasalt sediments. In contrast, thermal conduction mainly controls the subsalt layers. With a third simulation, we investigate the influence of a major permeable and of three impermeable faults dissecting the subsalt target reservoir and compare the results to the coupled model where no faults are integrated. The permeable fault may have a local, strong impact on the thermal, pressure and velocity fields whereas the impermeable faults only cause deviations of the pressure field.
This cumulative dissertation explored the use of the detection of natural background of fast neutrons, the so-called cosmic-ray neutron sensing (CRS) approach to measure field-scale soil moisture in cropped fields. Primary cosmic rays penetrate the top atmosphere and interact with atmospheric particles. Such interaction results on a cascade of high-energy neutrons, which continue traveling through the atmospheric column. Finally, neutrons penetrate the soil surface and a second cascade is produced with the so-called secondary cosmic-ray neutrons (fast neutrons). Partly, fast neutrons are absorbed by hydrogen (soil moisture). Remaining neutrons scatter back to the atmosphere, where its flux is inversely correlated to the soil moisture content, therefore allowing a non-invasive indirect measurement of soil moisture. The CRS methodology is mainly evaluated based on a field study carried out on a farmland in Potsdam (Brandenburg, Germany) along three crop seasons with corn, sunflower and winter rye; a bare soil period; and two winter periods. Also, field monitoring was carried out in the Schaefertal catchment (Harz, Germany) for long-term testing of CRS against ancillary data. In the first experimental site, the CRS method was calibrated and validated using different approaches of soil moisture measurements. In a period with corn, soil moisture measurement at the local scale was performed at near-surface only, and in subsequent periods (sunflower and winter rye) sensors were placed in three depths (5 cm, 20 cm and 40 cm). The direct transfer of CRS calibration parameters between two vegetation periods led to a large overestimation of soil moisture by the CRS. Part of this soil moisture overestimation was attributed to an underestimation of the CRS observation depth during the corn period ( 5-10 cm), which was later recalculated to values between 20-40 cm in other crop periods (sunflower and winter rye). According to results from these monitoring periods with different crops, vegetation played an important role on the CRS measurements. Water contained also in crop biomass, above and below ground, produces important neutron moderation. This effect was accounted for by a simple model for neutron corrections due to vegetation. It followed crop development and reduced overall CRS soil moisture error for periods of sunflower and winter rye. In Potsdam farmland also inversely-estimated soil hydraulic parameters were determined at the field scale, using CRS soil moisture from the sunflower period. A modelling framework coupling HYDRUS-1D and PEST was applied. Subsequently, field-scale soil hydraulic properties were compared against local scale soil properties (modelling and measurements). Successful results were obtained here, despite large difference in support volume. Simple modelling framework emphasizes future research directions with CRS soil moisture to parameterize field scale models. In Schaefertal catchment, CRS measurements were verified using precipitation and evapotranspiration data. At the monthly resolution, CRS soil water storage was well correlated to these two weather variables. Also clearly, water balance could not be closed due to missing information from other compartments such as groundwater, catchment discharge, etc. In the catchment, the snow influence to natural neutrons was also evaluated. As also observed in Potsdam farmland, CRS signal was strongly influenced by snow fall and snow accumulation. A simple strategy to measure snow was presented for Schaefertal case. Concluding remarks of this dissertation showed that (a) the cosmic-ray neutron sensing (CRS) has a strong potential to provide feasible measurement of mean soil moisture at the field scale in cropped fields; (b) CRS soil moisture is strongly influenced by other environmental water pools such as vegetation and snow, therefore these should be considered in analysis; (c) CRS water storage can be used for soil hydrology modelling for determination of soil hydraulic parameters; and (d) CRS approach has strong potential for long term monitoring of soil moisture and for addressing studies of water balance.
The last known eruption at Cerro Machin Volcano (CMV) in the Central Cordillera of Colombia occurred 900 years BP and ended with the formation of a dacitic lava dome. The dome rocks contain both normally and reversely zoned plagioclase (An(24-54)), unzoned and reversely zoned amphiboles of dominantly tschermakite and pargasite/magnesio-hastingsite composition and olivine xenocrysts (Fo = 85-88) with amphibole/clinopyroxene overgrowth, all suggesting interaction with mafic magma at depth. Plagioclase additionally exhibits complex oscillatory zoning patterns reflecting repeated replenishment, fractionation and changes in intrinsic conditions in the magma reservoir. Unzoned amphiboles and cores of the reversely zoned amphiboles give identical crystallization conditions of 910 +/- 30 degrees C and 360 +/- 70 MPa, corresponding to a depth of about 13 +/- 2 km, at moderately oxidized conditions (f(O2) = +0.5 +/- 0.2 Delta NNO). The water content in the melt, calculated based on amphibole chemistry, is 7.1 +/- 0.4 wt.%. Rims of the reversely zoned amphiboles are relatively enriched in MgO and yield higher crystallization temperatures (T = 970 +/- 25 degrees C), slightly lower melt H2O contents (6.1 +/- 0.7 wt.%) and overlapping pressures (410 +/- 100 MPa). We suggest that these rims crystallized following an influx of mafic melt into a resident magma reservoir at mid-crustal depths, further supported by the occurrence of xenocrystic olivine. Crystallization of biotite, albite-rich plagioclase and quartz occurred at comparatively low temperatures (probably <800 degrees C) during early stages of ascent or storage at shallower levels. Based on amphibole mineral chemistry, the felsic resident melt had a rhyolitic composition (71 +/- 2 wt.% SiO2), whereas the hybrid magma, from which the amphibole rims crystallized, was dacitic (64 +/- 3 wt.% SiO2). The bulk rock chemistry of the CMV lava dome dacites is homogenous. They have elevated (La/Nb)(N) ratios of 3.8-4.5, typical for convergent margin magmas, and display several geochemical characteristics, of adakites. Both Sr and Nd isotope compositions (Sr-87/Sr-86 similar to 0.70497, Nd-143/Nd-144 similar to 0.51267) are among the most radiogenic observed for the Northern Volcanic Zone of the Andes. They are distinct from oceanic crust that has been subducted in the region, pointing to a continental crustal control on the isotope composition and hence the adakitic signature, possibly in a crustal "hot zone".
Sophisticated methods have been developed and become standard in analysing floods as well as for assessing flood risk. However, increasingly critique of the current standards and scientific practice can be found both in the flood hydrology community as well as in the risk community who argue that the considerable amount of information already available on natural disasters has not been adequately deployed and brought to effective use. We describe this phenomenon as a failure to synthesize knowledge that results from barriers and ignorance in awareness, use and management of the entire spectrum of relevant content, that is, data, information and knowledge. In this paper we argue that the scientific community in flood risk research ignores event-specific analysis and documentations as another source of data. We present results from a systematic search that includes an intensive study on sources and ways of information dissemination of flood-relevant publications. We obtain 186 documents that contain information on the sources, pathways, receptors and/or consequences for any of the 40 strongest trans-basin floods in Germany in the period 1952-2002. This study therefore provides the most comprehensive metadata collection of flood documentations for the considered geographical space and period. A total of 87.5% of all events have been documented, and especially the most severe floods have received extensive coverage. Only 30% of the material has been produced in the scientific/academic environment, and the majority of all documents (about 80%) can be considered grey literature (i.e. literature not controlled by commercial publishers). Therefore, ignoring grey sources in flood research also means ignoring the largest part of knowledge available on single flood events (in Germany). Further, the results of this study underpin the rapid changes in information dissemination of flood event literature over the last decade. We discuss the options and obstacles of incorporating this data into the knowledge-building process in light of the current technological developments and international, interdisciplinary debates for data curation.
The northeastern margin of the Tibetan Plateau is a tectonically active region consisting of a series of faults with bounded intermountain basins and is located in the transition zone between the Tibetan Plateau and the Loess Plateau. Active deformation that may affect the topography in this region can be quantified using geomorphic indices. Therefore, we applied geomorphic indices such as the hypsometric integral and the stream length gradient index to infer neo-tectonics in the northeastern margin of the Tibetan Plateau. Different time-scaled geodetic leveling data and river incision rates were also integrated into the investigation. The results show that the hypsometric integrals are not significantly affected by lithology but spatially correspond to the hanging walls of thrust faults. The hypsometric integrals are also positively correlated with the leveling data. Although the stream length gradient index is influenced by lithology, its most pronounced anomalies of the stream length gradient are associated with the thrust faults. Consequently, the uplift in the northeast margin of the Tibetan Plateau appeared to be concentrated along the hanging walls of the thrust faults.
Facies bodies in geostatistical models of deep-water depositional environments generally represent channel-levee-overbank-lobe morphologies. Such models adequately capture one set of the erosional and depositional processes resulting from turbidity currents traveling downslope to the ocean basin floor. However, depositional morphologies diverge from the straight forward channel-levee-overbank-lobe paradigm when the topography of the slope or the shape of the basin impacts the timing and magnitude of turbidity current deposition. Subaqueous mass-transport-deposits (MTDs) present the need for an exception to the channel-levee-overbank-lobe archetype. Irregular surface topography of subaqueous MTDs can play a primary role in controlling sand deposition from turbidity currents. MTD topography creates mini-basins in which sand accumulates in irregularly-shaped deposits. These accumulations are difficult to laterally correlate using well-log data due to their variable and unpredictable shape and size. Prediction is further complicated because sandstone bodies typical of this setting are difficult to resolve in seismic-reflection data. An event-based model is presented, called DFTopoSim, which simulates debris flows and turbidity currents. The accommodation space on top of and between debris flow lobes is filled in by sand from turbidity currents. When applied to a subsurface case in the Molasse Basin of Upper Austria, DFTopoSim predicts sand packages consistent with observations from core, well, and seismic data and the interpretation of the sedimentologic processes. DFTopoSim expands the set of available geostatistical deep-water depositional models beyond the standard channel-levee-overbank-lobe model.
Emerged marine terraces and paleoshorelines along plate margins are prominent geomorphic markers that can be used to quantify the rates and patterns of crustal deformation. The northern margin of the Central Anatolian Plateau has been interpreted as an actively deforming orogenic wedge between the North Anatolian Fault and the Black Sea. Here we use uplifted marine terraces across principal faults on the Sinop Peninsula at the central northern side of the Pontide orogenic wedge to unravel patterns of Quaternary faulting and orogenic wedge behavior. We leveled the present-day elevations of paleoshorelines and dated marine terrace deposits using optically stimulated luminescence (OSL) to determine coastal uplift. The elevations of the paleoshorelines vary between 4 +/- 0.2 and 67 +/- 1.4 m above sea level and OSL ages suggest terrace formation episodes during interglacial periods at ca 125, 190, 400 and 570 ka, corresponding to marine isotopic stages (MIS) 5e, 7a, 11 and 15. Mean apparent vertical displacement rates (without eustatic correction) deduced from these terraces range between 0.02 and 0.18 mm/a, with intermittent faster rates of up to 0.26 mm/a. We obtained higher rates at the eastern and southern parts of the peninsula, toward the hinterland, indicating non-uniform uplift across the different morphotectonic segments of the peninsula. Our data are consistent with active on- and offshore faulting across the Sinop Peninsula. When integrated with regional tectonic observations, the faulting pattern reflects shortening distributed over a broad region of the northern margin of the Central Anatolian Plateau during the Quaternary.
The tectonically driven growth of mountains reflects the characteristics of the underlying fault systems and the applied tectonic forces. Over time, fault networks might be relatively static, but stress conditions could change and result in variations in fault slip orientation. Such a tectonic landscape would transition from a simple to a composite state: the topography of simple landscapes is correlated with a single set of tectonic boundary conditions, while composite landscapes contain inherited topography due to earlier deformation under different boundary conditions. We use fault interaction modeling to compare vertical displacement fields with topographic metrics to differentiate the two types of landscapes. By successively rotating the axis of maximum horizontal stress, we produce a suite of vertical displacement fields for comparison with real landscapes. We apply this model to a transpressional duplex in the south central Alborz Mountains of Iran, where NW oriented compression was superseded by neotectonic NE compression. The consistency between the modeled displacement field and real landforms indicates that the duplex topography is mostly compatible with the modern boundary conditions, but might include a small remnant from the earlier deformation phase. Our approach is applicable for various tectonic settings and represents an approach to identify the changing boundary conditions that produce composite landscapes. It may be particularly useful for identifying changes that occurred in regions where river profiles may no longer record a signal of the change or where the spatial pattern of uplift is complex.
Increases in animal products consumption and the associated environmental consequences have been a matter of scientific debate for decades. Consequences of such increases include rises in greenhouse gas emissions, growth of consumptive water use, and perturbation of global nutrients cycles. These consequences vary spatially depending on livestock types, their densities and their production system. In this letter, we investigate the spatial distribution of embodied crop calories in animal products. On a global scale, about 40% of the global crop calories are used as livestock feed (we refer to this ratio as crop balance for livestock) and about 4 kcal of crop products are used to generate 1 kcal of animal products (embodied crop calories of around 4). However, these values vary greatly around the world. In some regions, more than 100% of the crops produced is required to feed livestock requiring national or international trade to meet the deficit in livestock feed. Embodied crop calories vary between less than 1 for 20% of the livestock raising areas worldwide and greater than 10 for another 20% of the regions. Low values of embodied crop calories are related to production systems for ruminants based on fodder and forage, while large values are usually associated with production systems for non-ruminants fed on crop products. Additionally, we project the future feed demand considering three scenarios: (a) population growth, (b) population growth and changes in human dietary patterns and (c) changes in population, dietary patterns and feed conversion efficiency. When considering dietary changes, we project the global feed demand to be almost doubled (1.8-2.3 times) by 2050 compared to 2000, which would force us to produce almost equal or even more crops to raise our livestock than to directly nourish ourselves in the future. Feed demand is expected to increase over proportionally in Africa, South-Eastern Asia and Southern Asia, putting additional stress on these regions.
Changing food consumption patterns and associated greenhouse gas (GHG) emissions have been a matter of scientific debate for decades. The agricultural sector is one of the major GHG emitters and thus holds a large potential for climate change mitigation through optimal management and dietary changes. We assess this potential, project emissions, and investigate dietary patterns and their changes globally on a per country basis between 1961 and 2007. Sixteen representative and spatially differentiated patterns with a per capita calorie intake ranging from 1,870 to >3,400 kcal/day were derived. Detailed analyses show that low calorie diets are decreasing worldwide, while in parallel diet composition is changing as well: a discernable shift towards more balanced diets in developing countries can be observed and steps towards more meat rich diets as a typical characteristics in developed countries. Low calorie diets which are mainly observable in developing countries show a similar emission burden than moderate and high calorie diets. This can be explained by a less efficient calorie production per unit of GHG emissions in developing countries. Very high calorie diets are common in the developed world and exhibit high total per capita emissions of 3.7-6.1 kg CO2eq./day due to high carbon intensity and high intake of animal products. In case of an unbridled demographic growth and changing dietary patterns the projected emissions from agriculture will approach 20 Gt CO2eq./yr by 2050.
A high resolution multi proxy (pollen, grain size, total organic carbon) record from a small mountain lake (Lake Khuisiin; 46.6 degrees N, 101.8 degrees E; 2270 m a.s.l.) in the south eastern Khangai Mountains of central Mongolia has been used to explore changes in vegetation and climate over the last 1200 years. The pollen data indicates that the vegetation changed from dry steppe dominated by Poaceae and Artemisia (ca AD 760-950), to Larix forest steppe (ca AD 950-1170), Larix Betula forest steppe (ca AD 1170-1380), meadow dominated by Cyperaceae and Poaceae (ca AD 1380-1830), and Larix Betula forest steppe (after similar to AD 1830). The cold-wet period between AD 1380 and 1830 may relate to the Little Ice Age. Environmental changes were generally subtle and climate change seems to have been the major driver of variations in vegetation until at least the early part of the 20th century, suggesting that either the level of human activity was generally low, or the relationship between human activity and vegetation did not alter substantially between AD 760 and 1830. A review of centennial scale moisture records from China and Mongolia revealed that most areas experienced major changes at ca AD 1500 and AD 1900. However, the moisture availability since AD 1500 varied between sites, with no clear regional pattern or relationship to present day conditions. Both the reconstructions and the moisture levels simulation on a millennium scale performed in the MPI Earth System Model indicate that the monsoon-westerlies transition area shows a greater climate variability than those areas influenced by the westerlies, or by the summer monsoon only.
Purpose: Dry land vegetation is expected to respond sensitively to climate change and the projected variability of rainfall events. Rainfall as a water source is an obvious factor for the water supply of vegetation. However, the interaction of water and surface on rocky desert slopes with a patchy soil cover is also vital for vegetation in drylands. In particular, runoff on rocky surfaces and infiltration capacity of soil patches determine plant available water. Process-based studies into rock-soil interaction benefit from rainfall simulation, but require an approach accounting for the micro-scale heterogeneity of the slope surfaces. This study therefore aims at developing a suitable procedure for examining rock-soil interaction and the relevance of soil volume for storing plant available water in the northern Negev, Israel.
Materials and methods: To determine the amount of rainfall required to fill the available soil water storage capacity rainfall simulation experiments were conducted. The design of the rainfall-simulator and the selection of the plots aimed specifically at observing infiltration into small soil patches on a micro-scale relevant for the prevalent vegetation cover.
Results and discussion: The preliminary results of the study in the Negev Desert indicate that the ratio between soil volume and frequency of rainfall events determine the effect of climate change on plant available water and thus ultimately vegetation cover.
Conclusions: Based on the experiments examining runoff and soil moisture the qualitative understanding of hillslope ecohydrology in a rocky desert environment can be expanded into a quantitative assessment of the potential impact of varying rainfall conditions. The study also illustrates the contribution of rainfall simulation experiments for studies on the impact of climate change.
Explaning change in flood hazard in the Mekong river : the hypothesis of nonstationary variance
(2013)
Digital flow networks derived from digital elevation models (DEMs) sensitively react to errors due to measurement, data processing and data representation. Since high-resolution DEMs are increasingly used in geomorphological and hydrological research, automated and semi-automated procedures to reduce the impact of such errors on flow networks are required. One such technique is stream-carving, a hydrological conditioning technique to ensure drainage connectivity in DEMs towards the DEM edges. Here we test and modify a state-of-the-art carving algorithm for flow network derivation in a low-relief, agricultural landscape characterized by a large number of spurious, topographic depressions. Our results show that the investigated algorithm reconstructs a benchmark network insufficiently in terms of carving energy, distance and a topological network measure. The modification to the algorithm that performed best, combines the least-cost auxiliary topography (LCAT) carving with a constrained breaching algorithm that explicitly takes automatically identified channel locations into account. We applied our methods to a low relief landscape, but the results can be transferred to flow network derivation of DEMs in moderate to mountainous relief in situations where the valley bottom is broad and flat and precise derivations of the flow networks are needed.
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.
Through their relevance for sediment budgets and the sensitivity of geomorphic systems, geomorphic coupling and (sediment) connectivity represent important topics in geomorphology. Since the introduction of the systems perspective to physical geography by Chorley and Kennedy (1971), a catchment has been perceived as consisting of landscape elements (e.g. landforms, subcatchments) that are coupled by geomorphic processes through sediment transport. In this study, we present a novel application of mathematical graph theory to explore the network structure of coarse sediment pathways in a central alpine catchment. Numerical simulation models for rockfall, debris flows, and (hillslope and channel) fluvial processes are used to establish a spatially explicit graph model of sediment sources, pathways and sinks. The raster cells of a digital elevation model form the nodes of this graph, and simulated sediment trajectories represent the corresponding edges. Model results are validated by visual comparison with the field situation and aerial photos. The interaction of sediment pathways, i.e. where the deposits of a geomorphic process form the sources of another process, forms sediment cascades, represented by paths (a succession of edges) in the graph model. We show how this graph can be used to explore upslope (contributing area) and downslope (source to sink) functional connectivity by analysing its nodes, edges and paths. The analysis of the spatial distribution, composition and frequency of sediment cascades yields information on the relative importance of geomorphic processes and their interaction (however regardless of their transport capacity). In the study area, the analysis stresses the importance of mass movements and their interaction, e.g. the linkage of large rockfall source areas to debris flows that potentially enter the channel network. Moreover, it is shown that only a small percentage of the study area is coupled to the channel network which itself is longitudinally disconnected by natural and anthropogenic barriers. Besides the case study, we discuss the methodological framework and alternatives for node and edge representations of graph models in geomorphology. We conclude that graph theory provides an excellent methodological framework for the analysis of geomorphic systems, especially for the exploration of quantitative approaches towards sediment connectivity.
Soilscapes of the post-glacial morainic regions of the youngest glaciation are characterized by small hydrological kettle hole catchments forming hummocky soil landscapes. The spatial heterogeneity of subsurface structures as well as erosion-controlled pedogenesis under arable land use may complicate hydrological modeling. Our aim was to generate a soil landscape model for a small representative kettle hole catchment based on geoelectrical exploration and soil profile information. For a 1-ha catchment located in the northeastern German lowlands near the town of Prenzlau, electrical resistivity transects were determined by a multi electrode system (IMPETUS 12 Fs) and electrical conductivity (ECa) was mapped by using the electromagnetic induction (EMI) device EM38DD in both the vertical and horizontal modes. The 1-m digital elevation model (DEM) was obtained by kriging from high resolution manual elevation data determined with a leveling device (ZEISS Ni 40). Soil profile data from 26 boreholes distributed radially around the central pond were used to identify boundaries between soil horizons. The soil is characterized by varying topography and morphology of diagnostic horizons such as M- (colluvium), Bt- (clay illuviation), and C- (parent glacial till). By EMI mapping we identified (i) the boundary between erosive and colluvial areas around the kettle hole, and modeled (ii) the subsurface morphology of loamy horizons. Electrical resistivity tomography results coincide with these findings and allow for distinguishing between sandy and loamy dominated areas both in vertical and horizontal direction, respectively. This soil model of soil textural properties could be used for hydrological modeling.
The knowledge of the local soil structure is important for the assessment of seismic hazards. A widespread, but time-consuming technique to retrieve the parameters of the local underground is the drilling of boreholes. Another way to obtain the shear wave velocity profile at a given location is the inversion of surface wave dispersion curves. To ensure a good resolution for both superficial and deeper layers, the used dispersion curves need to cover a wide frequency range. This wide frequency range can be obtained using several arrays of seismic sensors or a single array comprising a large number of sensors. Consequently, these measurements are time-consuming. A simpler alternative is provided by the use of the ellipticity of Rayleigh waves. The frequency dependence of the ellipticity is tightly linked to the shear wave velocity profile. Furthermore, it can be measured using a single seismic sensor. As soil structures obtained by scaling of a given model exhibit the same ellipticity curve, any inversion of the ellipticity curve alone will be ambiguous. Therefore, additional measurements which fix the absolute value of the shear wave velocity profile at some points have to be included in the inversion process. Small-scale spatial autocorrelation measurements or MASW measurements can provide the needed data. Using a theoretical soil structure, we show which parts of the ellipticity curve have to be included in the inversion process to get a reliable result and which parts can be omitted. Furthermore, the use of autocorrelation or high-frequency dispersion curves will be highlighted. The resulting guidelines for inversions including ellipticity data are then applied to real data measurements collected at 14 different sites during the European NERIES project. It is found that the results are in good agreement with dispersion curve measurements. Furthermore, the method can help in identifying the mode of Rayleigh waves in dispersion curve measurements.
The Arctic tundra, covering approx. 5.5 % of the Earth’s land surface, is one of the last ecosystems remaining closest to its untouched condition. Remote sensing is able to provide information at regular time intervals and large spatial scales on the structure and function of Arctic ecosystems. But almost all natural surfaces reveal individual anisotropic reflectance behaviors, which can be described by the bidirectional reflectance distribution function (BRDF). This effect can cause significant changes in the measured surface reflectance depending on solar illumination and sensor viewing geometries. The aim of this thesis is the hyperspectral and spectro-directional reflectance characterization of important Arctic tundra vegetation communities at representative Siberian and Alaskan tundra sites as basis for the extraction of vegetation parameters, and the normalization of BRDF effects in off-nadir and multi-temporal remote sensing data. Moreover, in preparation for the upcoming German EnMAP (Environmental Mapping and Analysis Program) satellite mission, the understanding of BRDF effects in Arctic tundra is essential for the retrieval of high quality, consistent and therefore comparable datasets. The research in this doctoral thesis is based on field spectroscopic and field spectro-goniometric investigations of representative Siberian and Alaskan measurement grids. The first objective of this thesis was the development of a lightweight, transportable, and easily managed field spectro-goniometer system which nevertheless provides reliable spectro-directional data. I developed the Manual Transportable Instrument platform for ground-based Spectro-directional observations (ManTIS). The outcome of the field spectro-radiometrical measurements at the Low Arctic study sites along important environmental gradients (regional climate, soil pH, toposequence, and soil moisture) show that the different plant communities can be distinguished by their nadir-view reflectance spectra. The results especially reveal separation possibilities between the different tundra vegetation communities in the visible (VIS) blue and red wavelength regions. Additionally, the near-infrared (NIR) shoulder and NIR reflectance plateau, despite their relatively low values due to the low structure of tundra vegetation, are still valuable information sources and can separate communities according to their biomass and vegetation structure. In general, all different tundra plant communities show: (i) low maximum NIR reflectance; (ii) a weakly or nonexistent visible green reflectance peak in the VIS spectrum; (iii) a narrow “red-edge” region between the red and NIR wavelength regions; and (iv) no distinct NIR reflectance plateau. These common nadir-view reflectance characteristics are essential for the understanding of the variability of BRDF effects in Arctic tundra. None of the analyzed tundra communities showed an even closely isotropic reflectance behavior. In general, tundra vegetation communities: (i) usually show the highest BRDF effects in the solar principal plane; (ii) usually show the reflectance maximum in the backward viewing directions, and the reflectance minimum in the nadir to forward viewing directions; (iii) usually have a higher degree of reflectance anisotropy in the VIS wavelength region than in the NIR wavelength region; and (iv) show a more bowl-shaped reflectance distribution in longer wavelength bands (>700 nm). The results of the analysis of the influence of high sun zenith angles on the reflectance anisotropy show that with increasing sun zenith angles, the reflectance anisotropy changes to azimuthally symmetrical, bowl-shaped reflectance distributions with the lowest reflectance values in the nadir view position. The spectro-directional analyses also show that remote sensing products such as the NDVI or relative absorption depth products are strongly influenced by BRDF effects, and that the anisotropic characteristics of the remote sensing products can significantly differ from the observed BRDF effects in the original reflectance data. But the results further show that the NDVI can minimize view angle effects relative to the contrary spectro-directional effects in the red and NIR bands. For the researched tundra plant communities, the overall difference of the off-nadir NDVI values compared to the nadir value increases with increasing sensor viewing angles, but on average never exceeds 10 %. In conclusion, this study shows that changes in the illumination-target-viewing geometry directly lead to an altering of the reflectance spectra of Arctic tundra communities according to their object-specific BRDFs. Since the different tundra communities show only small, but nonetheless significant differences in the surface reflectance, it is important to include spectro-directional reflectance characteristics in the algorithm development for remote sensing products.
In this study, we evaluated the potential of the hydrogen isotopic composition of algal lipid biomarkers as a proxy for past hydroclimatic variability in hypersaline Isabel Lake, Mexico (Eastern Pacific). We compared rainfall variability recorded in the region over the last 65 years with changes in delta D values of the most abundant compounds preserved in the uppermost 16 cm of lake sediment. Changes in delta D values of the 1,15-C-32 diol (delta D-diol), a specific biomarker of algal populations, were related to rainfall variability; specifically, n-alkyl diols were more deuterium-enriched (depleted) during wetter (drier) periods. Strikingly, neither the magnitude of lipid biomarker isotopic changes over interannual timescales (of up to 70-80 parts per thousand) nor the direction of that variability can be explained by changes in delta D values of the water source or salinity fluctuations (approximately 30 on the practical salinity scale) controlled by seasonal rainfall. However, changes in sedimentary biomarker composition, higher total organic carbon content and less negative delta C-13 values of the 1,15-C-32 diol indicate enhanced algal growth during wetter periods. We find that these conditions result in less negative delta D values of n-alkyl diols. We hypothesize that due to higher lipid demand during enhanced algal growth, an increasing proportion of hydrogen for lipid synthesis is derived from the cytosol via oxidation of polysaccharides, which may cause a deuterium enrichment of the acetogenic compounds. This study has significant implications for paleohydrological reconstructions using algal lipid delta D values, particularly in highly seasonal environments such as Isabel Lake. In such environments, delta D values of specific algal lipid biomarkers may not record the full seasonal cycle in source water delta D but appear to be mainly controlled by the physiological state of algal populations. Our data provide the first evidence that changes in D/H fractionation due to algal growth conditions can be recorded in sediments. For paleoclimate reconstructions in such environments, algal growth conditions should be constrained with additional proxy data (delta C-13 values of the same biomarkers), as the net D/H fractionation between water and lipids may not have been constant over time.
North of Naran in the Kaghan Valley (NE Pakistan), the metabasites of the melange units lying below the Kohistan Arc, contain glaucophane. Typically they reflect blueschist-metamorphic conditions (0.7 GPa, 400 degrees C). Associated graphite-rich metapelites with quartz veins underwent upper greenschist to amphibolite conditions. Near Naran we observed in quartz grains of type 3 veins first minute relics of Fe-Mg carpholite indicating earlier blueschist metamorphic conditions. P-T estimates indicate 1.2-1.6 GPa at 380-410 degrees C, pressure and temperature values typical for blueschist metamorphic conditions. Changes in mineral assemblages and abundant sudoite component in associated chlorite point to a pressure drop after peak I conditions. We assign the observed changes to peak I conditions occurring during a Cretaceous subduction event. Temperatures estimated with Raman graphite-thermometry clearly indicate a significant subsequent rise of post-peak I temperatures up to 500 degrees C. This is compatible with the amphibolite peak II assigned to the Tertiary continental collision that leads to subduction of the Indian Plate and ultra-high-pressure metamorphism. During subduction the blueschist metamorphic metapelites underwent dehydration, which caused alteration in the overlying lithospheric mantle. In a hydrated lithospheric mantle density is significantly reduced which enhanced subduction of continental crust in the Higher Himalaya. This P-T evolution is typical for a collision orogen with a high plateau but remarkably contrasting findings from Eastern Anatolia, where plateau building is in "statu nascendi" (e.g., Oberhänsli et al., 2010).
Extremely rare veinlets and reaction textures composed of symplectites of olivine (similar to Fo(43-55)) + plagioclase +/- spinel +/- ilmenite, associated with more common pyroxene + plagioclase and amphibole + plagioclase varieties, are preserved within eclogites and garnet pyroxenites in the Moldanubian Zone of the Bohemian Massif. Thermodynamic modelling integrated with conventional geothermometry conducted on an eclogite reveals that the symplectite-forming stage occurred at high T (similar to 850 degrees C) and low P (< 6 and > 2 center dot 5 kbar). The development of the different symplectite types reflects reactions that took place in micro-scale domains. The breakdown of high-P garnet controlled the formation of olivine-bearing and amphibole + plagioclase symplectites, whereas breakdown of high-P omphacite led to formation of pyroxene + plagioclase symplectites. In addition, post-eclogite facies but pre-symplectite stage porphyroblastic amphibole and phlogopite were also replaced by olivine-bearing symplectites. Material transfer calculations and thermodynamic modelling indicate that the formation of different symplectite types was linked despite their different bulk compositions. For example, the olivine-bearing symplectites gained Fe +/- Mg, whereas adjacent amphibole + plagioclase and pyroxene + plagioclase symplectites show losses in Fe and Mg; Al, Si and Ca were also variably exchanged. The olivine-bearing symplectites were particularly sensitive to Na despite the small concentration of this element. In eclogites where Na was readily available, the plagioclase composition in the olivine-bearing symplectites shifted from pure anorthite to bytownite, with the less calcic feldspar partitioning Si and inhibiting the formation of orthopyroxene. This regional high-T, low-P granulite-facies symplectite overprint may have been caused by advective heat loss from rapidly exhumed high-T, high-P granulitic bodies (Gfohl Unit) that were emplaced into and over the middle crust (Monotonous and Varied Series) during Carboniferous continent-continent collision.
In the recent past, the Alpine Lech valley (Austria) experienced three damaging flood events within 6 years despite the various structural flood protection measures in place. For an improved flood risk management, the analysis of flood damage potentials is a crucial component. Since the expansion of built-up areas and their associated values is seen as one of the main drivers of rising flood losses, the goal of this study is to analyze the spatial development of the assets at risk, particularly of residential areas, due to land use changes over a historic period (since 1971) and up to possible shifts in future (until 2030). The analysis revealed that the alpine study area was faced to remarkable land use changes like urbanization and the decline of agriculturally used grassland areas. Although the major agglomeration of residential areas inside the flood plains took place before 1971, a steady growth of values at risk can still be observed until now. Even for the future, the trend is ongoing, but depends very much on the assumed land use scenario and the underlying land use policy. Between 1971 and 2006, the annual growth rate of the damage potential of residential areas amounted to 1.1 % ('constant values,' i.e., asset values at constant prices of reference year 2006) or 3.0 % ('adjusted values,' i.e., asset values adjusted by GDP increase at constant prices of reference year 2006) for three flood scenarios. For the projected time span between 2006 and 2030, a further annual increase by 1.0 % ('constant values') or even 4.2 % ('adjusted values') may be possible when the most extreme urbanization scenario 'Overall Growth' is considered. Although socio-economic development is regarded as the main driver for increasing flood losses, our analysis shows that settlement development does not preferably take place within flood prone areas.
Sediments of a thermokarst system on the north-eastern Tibetan Plateau were studied to infer changes in the lacustrine depositional environment related to climatic changes since the early Holocene. The thermokarst pond with a length of 360 m is situated in a 14.5 x 6 km tectonically unaffected intermontane basin, which is underlain by discontinuous permafrost.
A lake sediment core and bankside lacustrine onshore deposits were analysed. Additionally, fossil lake sediments were investigated, which document a former lake-level high stand. The sediments are mainly composed of marls with variable amounts of silt carbonate micrite, and organic matter.
On the basis of sedimentological (grain size data), geochemical (XRF), mineralogical (XRD) and micropaleontological data (ostracods and chironomide assemblages) a reconstruction of a paleolake environment was achieved.
Lacustrine sediments with endogenic carbonate precipitation suggest a lacustrine environment since at least 19.0 cal ka BP. However, because of relocation and reworking processes in the lake, the sediments did not provide distinct information about the ultimate formation of the lake. The high amount of endogenic carbonate suggests prolonged still-water conditions at about 9.3 cal ka BP. Ostracod shells and chironomid head capsules in fossil lake sediments indicate at least one former lake-level high stand, which were developed between the early and middle Holocene. From the late Holocene the area was possibly characterized by a lake-level decline, documented by a hiatus between lacustrine sediments and a reworked loess or loess-like horizon. After the lake-level decline and the following warming period, the area was affected by thermally-induced subsidence and a re-flooding of the basin because of thawing permafrost.
Improving our understanding of biodiversity and ecosystem functioning and our capacity to inform ecosystem management requires an integrated framework for functional biodiversity research (FBR). However, adequate integration among empirical approaches (monitoring and experimental) and modelling has rarely been achieved in FBR. We offer an appraisal of the issues involved and chart a course towards enhanced integration. A major element of this path is the joint orientation towards the continuous refinement of a theoretical framework for FBR that links theory testing and generalization with applied research oriented towards the conservation of biodiversity and ecosystem functioning. We further emphasize existing decision-making frameworks as suitable instruments to practically merge these different aims of FBR and bring them into application. This integrated framework requires joint research planning, and should improve communication and stimulate collaboration between modellers and empiricists, thereby overcoming existing reservations and prejudices. The implementation of this integrative research agenda for FBR requires an adaptation in most national and international funding schemes in order to accommodate such joint teams and their more complex structures and data needs.
While sea level rise is one of the most likely consequences of climate change, the provoked costs remain highly uncertain. Based on a block-maxima approach, we provide a stochastic framework to estimate the increase of expected damages with sea level rise as well as with meteorological changes and demonstrate the application to two case studies. In addition, the uncertainty of the damage estimations due to the stochastic nature of extreme events is studied. Starting with the probability distribution of extreme flood levels, we calculate the distribution of implied damages in a specific region employing stage-damage functions. Universal relations of the expected damages and their standard deviation, which demonstrate the importance of the shape of the damage function, are provided. We also calculate how flood protection reduces the damages leading to a more complex picture, where the extreme value behavior plays a fundamental role. Citation: Boettle, M., D. Rybski, and J. P. Kropp (2013), How changing sea level extremes and protection measures alter coastal flood damages, Water Resour. Res., 49, 1199-1210, doi: 10.1002/wrcr.20108.
Logging and large earthquakes are disturbances that may significantly affect hydrological and erosional processes and process rates, although in decisively different ways. Despite numerous studies that have documented the impacts of both deforestation and earthquakes on water and sediment fluxes, a number of details regarding the timing and type of de- and reforestation; seismic impacts on subsurface water fluxes; or the overall geomorphic work involved have remained unresolved. The main objective of this thesis is to address these shortcomings and to better understand and compare the hydrological and erosional process responses to such natural and man-made disturbances. To this end, south-central Chile provides an excellent natural laboratory owing to its high seismicity and the ongoing conversion of land into highly productive plantation forests. In this dissertation I combine paired catchment experiments, data analysis techniques, and physics-based modelling to investigate: 1) the effect of plantation forests on water resources, 2) the source and sink behavior of timber harvest areas in terms of overland flow generation and sediment fluxes, 3) geomorphic work and its efficiency as a function of seasonal logging, 4) possible hydrologic responses of the saturated zone to the 2010 Maule earthquake and 5) responses of the vadose zone to this earthquake. Re 1) In order to quantify the hydrologic impact of plantation forests, it is fundamental to first establish their water balances. I show that tree species is not significant in this regard, i.e. Pinus radiata and Eucalyptus globulus do not trigger any decisive different hydrologic response. Instead, water consumption is more sensitive to soil-water supply for the local hydro-climatic conditions. Re 2) Contradictory opinions exist about whether timber harvest areas (THA) generate or capture overland flow and sediment. Although THAs contribute significantly to hydrology and sediment transport because of their spatial extent, little is known about the hydrological and erosional processes occurring on them. I show that THAs may act as both sources and sinks for overland flow, which in turn intensifies surface erosion. Above a rainfall intensity of ~20 mm/h, which corresponds to <10% of all rainfall, THAs may generate runoff whereas below that threshold they remain sinks. The overall contribution of Hortonian runoff is thus secondary considering the local rainfall regime. The bulk of both runoff and sediment is generated by Dunne, saturation excess, overland flow. I also show that logging may increase infiltrability on THAs which may cause an initial decrease in streamflow followed by an increase after the groundwater storage has been refilled. Re 3) I present changes in frequency-magnitude distributions following seasonal logging by applying Quantile Regression Forests at hitherto unprecedented detail. It is clearly the season that controls the hydro-geomorphic work efficiency of clear cutting. Logging, particularly dry seasonal logging, caused a shift of work efficiency towards less flashy and mere but more frequent moderate rainfall-runoff events. The sediment transport is dominated by Dunne overland flow which is consistent with physics-based modelling using WASA-SED. Re 4) It is well accepted that earthquakes may affect hydrological processes in the saturated zone. Assuming such flow conditions, consolidation of saturated saprolitic material is one possible response. Consolidation raises the hydraulic gradients which may explain the observed increase in discharge following earthquakes. By doing so, squeezed water saturates the soil which in turn increases the water accessible for plant transpiration. Post-seismic enhanced transpiration is reflected in the intensification of diurnal cycling. Re 5) Assuming unsaturated conditions, I present the first evidence that the vadose zone may also respond to seismic waves by releasing pore water which in turn feeds groundwater reservoirs. By doing so, water tables along the valley bottoms are elevated thus providing additional water resources to the riparian vegetation. By inverse modelling, the transient increase in transpiration is found to be 30-60%. Based on the data available, both hypotheses, are not testable. Finally, when comparing the hydrological and erosional effects of the Maule earthquake with the impact of planting exotic plantation forests, the overall observed earthquake effects are comparably small, and limited to short time scales.
During the African Humid Period (AHP, 15-5 ka BP) an almost 300 m deep paleo-lake covering 2200 km(2) developed in the Suguta Valley, in the Northern Kenya Rift Data from lacustrine sediments and paleo-shorelines indicate that a large paleo-lake already existed by 13.9 ka BP, and record rapid water level fluctuations of up to 100 m within periods of 100 years or less, and a final lowstand at the end of the AHP (5 ka BP). We used a hydro-balance model to assess the abruptness of these water level fluctuations and identify their causes. We observed that fluctuations within the AHP were caused by abrupt changes in precipitation of 26-40%. Despite the absence of continuous lacustrine data documenting the onset of the AHP in the Suguta Valley, we conclude from the hydro-balance model that only an abrupt onset to the AHP, prior to 14.8 ka BP, could have led to high water levels recorded. The modeling results suggest that the sudden increase in rainfall was the direct consequence of an eastward migration of the Congo Air Boundary (CAB), caused by an enhanced atmospheric pressure gradient between East Africa and southern Asia during a northern hemisphere (NH) summer insolation maximum. In contrast the end of the AHP must have been gradual despite an abrupt change in the source of precipitation when a decreasing pressure gradient between Asia and Africa prevented the CAB from reaching the study area. This abruptness was probably buffered by a contemporaneous change in precession producing an insolation maximum at the equator during September-October. This change would have meant that the only rain source was the Intertropical Convergence Zone (IT CZ), which would have carried a greater amount of moisture during the short rainy season thus slowing the fall in water level over a period of about 1000 years in association with the reduction in insolation. The results of this study provide an indication of the amount of time available for humans in north-eastern Africa to adapt in response to a changing climate, from hunting and gathering to farming and herding.
Several mechanisms are proposed to be part of the earthquake triggering process, including static stress interactions and dynamic stress transfer. Significant differences of these mechanisms are particularly expected in the spatial distribution of aftershocks. However, testing the different hypotheses is challenging because it requires the consideration of the large uncertainties involved in stress calculations as well as the appropriate consideration of secondary aftershock triggering which is related to stress changes induced by smaller pre- and aftershocks. In order to evaluate the forecast capability of different mechanisms, I take the effect of smaller--magnitude earthquakes into account by using the epidemic type aftershock sequence (ETAS) model where the spatial probability distribution of direct aftershocks, if available, is correlated to alternative source information and mechanisms. Surface shaking, rupture geometry, and slip distributions are tested. As an approximation of the shaking level, ShakeMaps are used which are available in near real-time after a mainshock and thus could be used for first-order forecasts of the spatial aftershock distribution. Alternatively, the use of empirical decay laws related to minimum fault distance is tested and Coulomb stress change calculations based on published and random slip models. For comparison, the likelihood values of the different model combinations are analyzed in the case of several well-known aftershock sequences (1992 Landers, 1999 Hector Mine, 2004 Parkfield). The tests show that the fault geometry is the most valuable information for improving aftershock forecasts. Furthermore, they reveal that static stress maps can additionally improve the forecasts of off--fault aftershock locations, while the integration of ground shaking data could not upgrade the results significantly. In the second part of this work, I focused on a procedure to test the information content of inverted slip models. This allows to quantify the information gain if this kind of data is included in aftershock forecasts. For this purpose, the ETAS model based on static stress changes, which is introduced in part one, is applied. The forecast ability of the models is systematically tested for several earthquake sequences and compared to models using random slip distributions. The influence of subfault resolution and segment strike and dip is tested. Some of the tested slip models perform very good, in that cases almost no random slip models are found to perform better. Contrastingly, for some of the published slip models, almost all random slip models perform better than the published slip model. Choosing a different subfault resolution hardly influences the result, as long the general slip pattern is still reproducible. Whereas different strike and dip values strongly influence the results depending on the standard deviation chosen, which is applied in the process of randomly selecting the strike and dip values.
Reliable estimations of magnitude of completeness (M-c) are essential for a correct interpretation of seismic catalogues. The spatial distribution of M-c may be strongly variable and difficult to assess in mining environments, owing to the presence of galleries, cavities, fractured regions, porous media and different mineralogical bodies, as well as in consequence of inhomogeneous spatial distribution of the seismicity. We apply a 3-D modification of the probabilistic magnitude of completeness (PMC) method, which relies on the analysis of network detection capabilities. In our approach, the probability to detect an event depends on its magnitude, source receiver Euclidian distance and source receiver direction. The suggested method is proposed for study of the spatial distribution of the magnitude of completeness in a mining environment and here is applied to a 2-months acoustic emission (AE) data set recorded at the Morsleben salt mine, Germany. The dense seismic network and the large data set, which includes more than one million events, enable a detailed testing of the method. This method is proposed specifically for strongly heterogeneous media. Besides, it can also be used for specific network installations, with sensors with a sensitivity, dependent on the direction of the incoming wave (e.g. some piezoelectric sensors). In absence of strong heterogeneities, the standards PMC approach should be used. We show that the PMC estimations in mines strongly depend on the source receiver direction, and cannot be correctly accounted using a standard PMC approach. However, results can be improved, when adopting the proposed 3-D modification of the PMC method. Our analysis of one central horizontal and vertical section yields a magnitude of completeness of about M-c approximate to 1 (AE magnitude) at the centre of the network, which increases up to M-c approximate to 4 at further distances outside the network; the best detection performance is estimated for a NNE-SSE elongated region, which corresponds to the strike direction of the low-attenuating salt body. Our approach provides us with small-scale details about the capability of sensors to detect an earthquake, which can be linked to the presence of heterogeneities in specific directions. Reduced detection performance in presence of strong structural heterogeneities (cavities) is confirmed by synthetic waveform modelling in heterogeneous media.
The coastal stretch of north-eastern Mediterranean Morocco holds vitally important ecological, social, and economic functions. The implementation of large-scale luxury tourism resorts shall push socio-economic development and facilitate the shift from a mainly agrarian to a service economy. Sufficient water availability and intact beaches are among the key requirements for the successful realization of regional development plans. The water situation is already critical, additional water-intense sectors could overstrain the capacity of water resources. Further, coastal erosion caused by sea-level rise is projected. Regional climate change is observable, and must be included in regional water management. Long-term climate trends are assessed for the larger region (Moulouya basin) and for the near-coastal zone at Saidia. The amount of additional water demand is assessed for the large-dimensioned Saidia resort; including the monthly, seasonal and annual tourist per capita water need under inclusion of irrigated golf courses and garden areas. A shift of climate patterns is observed, a lengthening of the dry summer season, and as well a significant decline of annual precipitation. Thus, current water scarcity is mainly human-induced; however, climate change will aggravate the situation. As a consequence, severe environmental damage due to water scarcity is likely and could impinge on the quality of local tourism. The re-adjustment of current management routines is therefore essential. Possible adjustments are discussed and the analysis concludes with management recommendations for innovative regional water management of tourism facilities.
Da geologische Störungen können als Grundwasserleiter, -Barrieren oder als gemischte leitende /stauende Fluidsysteme wirken. Aufgrund dessen können Störungen maßgeblich den Grundwasserfluss im Untergrund beeinflussen, welcher deutliche Veränderungen des tiefen thermischen Feldes bewirken kann. Grundwasserdynamik und Temperaturveränderungen sind wiederum entscheidende Faktoren für die Exploration geothermischer Energie. Diese Studie untersuchte den Einfluss von Störungen auf das Fluidsystem und das thermische Feld im Untergrund. Sie erforschte die physikalischen Prozesse, welche das Fluidverhalten und die Temperaturverteilung in Störungen und in den umgebenden Gesteinen. Dazu wurden 3D Finite Elemente Simulationen des gekoppelten Fluid und Wärmetransports für synthetische sowie reale Modelszenarien auf unterschiedlichen Skalen durchgeführt. Um den Einfluss einer schräg einfallenden Störung systematisch durch die schrittweise Veränderung der hydraulischen Öffnungsweite und der Permeabilität, zu untersuchen, wurde ein klein-skaliges synthetisches Modell entwickelt. Ein inverser linearer Zusammenhang wurde festgestellt, welcher zeigt, dass sich die Fluidgeschwindigkeit in der Störung jeweils um ~1e-01 m/s verringert, wenn die Öffnungsweite der Störung um jeweils eine Magnitude vergrößert wird. Ein hoher Permeabilitätskontrast zwischen Störung und umgebender Matrix begünstigt die Fluidadvektion hin zur Störung und führt zu ausgeprägten Druck- und Temperaturveränderungen innerhalb und um die Störung herum. Bei geringem Permeabilitätskontrast zwischen Störung und umgebendem Gestein findet hingegen kein Fluidfluss in der Störung statt, wobei das hydrostatische Druck- sowie das Temperaturfeld unverändert bleiben. Auf Grundlage der synthetischen Modellierungsergebnisse wurde der Einfluss von Störungen auf einer größeren Skala anhand eines komplexeren (realen) geologischen Systems analysiert. Dabei handelt es sich um ein 3D Modell des Geothermiestandortes Groß Schönebeck, der ca. 40 km nördlich von Berlin liegt. Die Integration von einer permeablen und drei impermeablen Hauptstörungen, zeigte unterschiedlich starke Einflüsse auf Fluidzirkulation, Temperatur – und Druckfeld. Die modellierte konvektive Zirkulation in der permeablen Störung verändert das thermische Feld stark (bis zu 15 K). In den gering durchlässigen Störungen wird die Wärme ausschließlich durch Diffusion geleitet. Der konduktive Wärmetransport beeinflusst das thermische Feld nicht, bewirkt jedoch lokale Veränderungen des hydrostatischen Druckfeldes. Um den Einfluss großer Störungszonen mit kilometerweitem vertikalen Versatz auf das geothermische Feld der Beckenskala zu untersuchen, wurden gekoppelte Fluid- und Wärmetransportsimulationen für ein 3D Strukturmodell des Gebietes Brandenburg durchgeführt (Noack et al. 2010; 2013). Bezüglich der Störungspermeabilität wurden verschiedene geologische Szenarien modelliert, von denen zwei Endgliedermodelle ausgewertet wurden. Die Ergebnisse zeigten, dass die undurchlässigen Störungen den Fluidfluss nur lokal beeinflussen. Da sie als hydraulische Barrieren wirken, wird der Fluidfluss mir sehr geringen Geschwindigkeiten entlang der Störungen innerhalb eines Bereichs von ~ 1 km auf jeder Seite umgelenkt. Die modellierten lokalen Veränderungen des Grundwasserzirkulationssystems haben keinen beobachtbaren Effekt auf das Temperaturfeld. Hingegen erzeugen permeable Störungszonen eine ausgeprägte thermische Signatur innerhalb eines Einflussbereichs von ~ 2.4-8.8 km in -1000 m Tiefe und ~6-12 km in -3000 m Tiefe. Diese thermische Signatur, in der sich kältere und wärmere Temperaturbereiche abwechseln, wird durch auf- und abwärts gerichteten Fluidfluss innerhalb der Störung verursacht, der grundsätzlich durch existierende Gradienten in der hydraulischen Druckhöhe angetrieben wird. Alle Studien haben gezeigt, dass Störungen einen beachtlichen Einfluss auf den Fluid-, und Wärmefluss haben. Es stellte sich heraus, dass die Permeabilität in der Störung und in den umgebenden geologischen Schichten so wie der spezifische geologische Rahmen entscheidende Faktoren in der Ausbildung verschiedener Wärmetransportmechanismen sind, die sich in Störungen entwickeln können. Die von permeablen Störungen verursachten Temperaturveränderungen können lokal, jedoch groß sein, genauso wie die durch hydraulisch leitende und nichtleitende Störungen hervorgerufenen Veränderungen des Fluidystems. Letztlich haben die Simulationen für die unterschiedlich skalierten Modelle gezeigt, dass die Ergebnisse sich nicht aufeinander übertragen lassen und dass es notwendig ist, jeden geologischen Rahmen hinsichtlich Konfiguration und Größenskala gesondert zu betrachten. Abschließend hat diese Studie demonstriert, dass die Betrachtung von Störungen in 3D Finiten Elementen Modellen für die Simulation von gekoppeltem Fluid- und Wärmetransport auf unterschiedlichen Skalen möglich ist. Da diese Art von numerischen Simulationen sowohl die geologische Struktur des Untergrunds sowie die im Erdinnern ablaufenden physikalischen Prozesse integriert, können sie einen wertvollen Beitrag leisten, indem sie Feld- und Laborgestützte Untersuchungen vervollständigen.
We integrate the information of multiple tomographic models acquired from the earth's surface by modifying a statistical approach recently developed for the integration of cross-borehole tomographic models. In doing so, we introduce spectral cluster analysis as the new core of the model integration procedure to capture the spatial heterogeneity present in all considered tomographic models and describe this heterogeneity in a fuzzy sense. Because spectral cluster algorithms analyze model structure locally, they are considered relatively robust with regard to systematically and spatially varying imaging capabilities typical for geophysical tomographic surveys conducted on the earth's surface. Using a synthetic aquifer example, a fuzzy spectral cluster algorithm can be used to integrate the information provided by 2D tomographic refraction seismic and DC resistivity surveys. The integrated information in the fuzzy membership domain is then used to derive an integrated zonal geophysical model outlining the major structural units present in both input models. We also explain how the fuzzy membership information can be used to identify optimal locations for sparse logging of additional target parameters, i.e., porosity information in our synthetic example. We demonstrate how this sparse porosity information can be extrapolated based on all tomographic input models. The resultant 2D porosity model matches the original porosity distribution reasonably well within the spatial resolution limits of the underlying tomographic models. Consecutively, we apply this approach to a field data base acquired over a former river channel. Sparse information about natural gamma radiation is available and extrapolated on the basis of the fuzzy membership information obtained by spectral cluster analysis of 2D P-wave velocity and electrical resistivity models. This field data shows that the presented parameter extrapolation procedure is robust, even if the locations of target parameter acquisition have not been optimized with regard to the fuzzy membership information.
We have undertaken a high resolution palaeoclimate reconstruction on radiocarbon dated palaeolake sediments from the Spiti valley, NW Himalaya. This site lies in the climatically sensitive winter westerlies and Indian Summer Monsoon (ISM) transitional regime and provides an opportunity to reconstruct the precipitation seasonality, and extreme precipitation events that are characterised by intensified erosion. The lake sediments reveal distinct lithofacies that provide evidence of changes in depositional environment and climate during early to mid Holocene (8.7-6.1 cal ka BP). We have identified three stages during the period of lake's existence: the Stage I (8.7-7.6 cal ka BP) is marked by lake establishment; Stage II (similar to 7.6-6.8 cal ka BP) by sustained cooler periods and weakened summer monsoon, and Stage III (similar to 6.8-6.1 cal ka BP) by a shift from colder to warmer climate with stronger ISM. We have identified several short term cooler periods at ca. 8.7, 8.5, 8.3 and 7.2-6.9 cal ka BP. Based on an overview of regional climate records we show that there is an abrupt switch in precipitation seasonality ca. 6.8 cal ka BP that is followed by the onset of the intensified monsoon in the NW Himalaya. (C) 2013 Elsevier Ltd and INQUA. All rights reserved.
The spatial distribution of soil types is controlled by a set of environmental factors such as climate, organisms, parent material and topography as well as time and space. A change of these factors will lead to a change in the spatial distribution of soil types. In this study, we use a digital soil mapping approach to improve our knowledge about major soil type distributing factors in the steppe regions of Inner Mongolia (China) which currently undergo tremendous environmental change, e.g. climate and land use change. We use Random Forests in an effort to map Reference Soil Groups according to the World Reference Base for Soil Resources (WRB) in the Xilin River catchment. We benefit from the superior prediction capabilities of RF and additional interpretive results in order to identify the major environmental factors that control spatial patterns of soil types. The nine WRB soil groups that were identified and spatially predicted for the study area are Arenosol, Calcisol, Cambisol, Chernozem, Cryosol, Gleysol, Kastanozem, Phaeozem and Regosol.
Model and prediction performances of the RF model are high with an Out-of-Bag error of 51.6% for the model and a misclassification error for the predicted map of 28.9%. The main controlling factors of soil type distribution are land use, a set of topographic variables, geology and climate. However, land use and climate are of major importance and topography and geology are of minor importance. The visualizations of the predictions, the variable importance measures as result of RF and the comparisons of these with the spatial distribution of the environmental factors delivered additional, quantitative information of these controlling factors and revealed that intensively grazed areas are subjected to soil degradation. However, most of the area is still governed by natural soil forming processes which are driven by climate, topography and geology. Most importantly though, our study revealed that a shift towards warmer temperatures and lower precipitation regimes will lead to a change of the spatial distribution of RSGs towards steppe soils that store less carbon, i.e. a decrease of spatial extent of Phaeozems and an increase of spatial extent of Chernozems and Kastanozems.
Landslides are one of the biggest natural hazards in Georgia, a mountainous country in the Caucasus. So far, no systematic monitoring and analysis of the dynamics of landslides in Georgia has been made. Especially as landslides are triggered by extrinsic processes, the analysis of landslides together with precipitation and earthquakes is challenging. In this thesis I describe the advantages and limits of remote sensing to detect and better understand the nature of landslide in Georgia. The thesis is written in a cumulative form, composing a general introduction, three manuscripts and a summary and outlook chapter. In the present work, I measure the surface displacement due to active landslides with different interferometric synthetic aperture radar (InSAR) methods. The slow landslides (several cm per year) are well detectable with two-pass interferometry. In same time, the extremely slow landslides (several mm per year) could be detected only with time series InSAR techniques. I exemplify the success of InSAR techniques by showing hitherto unknown landslides, located in the central part of Georgia. Both, the landslide extent and displacement rate is quantified. Further, to determine a possible depth and position of potential sliding planes, inverse models were developed. Inverse modeling searches for parameters of source which can create observed displacement distribution. I also empirically estimate the volume of the investigated landslide using displacement distributions as derived from InSAR combined with morphology from an aerial photography. I adapted a volume formula for our case, and also combined available seismicity and precipitation data to analyze potential triggering factors. A governing question was: What causes landslide acceleration as observed in the InSAR data? The investigated area (central Georgia) is seismically highly active. As an additional product of the InSAR data analysis, a deformation area associated with the 7th September Mw=6.0 earthquake was found. Evidences of surface ruptures directly associated with the earthquake could not be found in the field, however, during and after the earthquake new landslides were observed. The thesis highlights that deformation from InSAR may help to map area prone landslides triggering by earthquake, potentially providing a technique that is of relevance for country wide landslide monitoring, especially as new satellite sensors will emerge in the coming years.
In the last decade growing interest has emerged in quantifying the spatial and temporal variations in mountain building. Until recently, insufficient data have been available to attempt such a task at the scale of large orogens such as the Himalaya. The Himalaya accommodates ongoing convergence between India and Eurasia and is a focal point for studying orogen evolution and hypothesized interactions between tectonics and climate. Here we integrate 1126 published bedrock mineral cooling ages with a transient 1D Monte-Carlo thermal-kinematic erosion model to quantify the denudation histories along similar to 2700 km of the Himalaya. The model free parameter is a temporally variable denudation rate from 50 Ma to present. Thermophysical material properties and boundary conditions were tuned to individual study areas. Monte-Carlo simulations were conducted to identify the range of denudation histories that can reproduce the observed cooling ages. Results indicate large temporal and spatial variations in denudation and these are resolvable across different tectonic units of the Himalaya. More specifically, across > 1000 km of the southern Greater Himalaya denudation rates were highest (similar to 1.5-3 mm/yr) between similar to 10 and 2 Ma and lower (0.5-2.6 mm/yr) over the last 2 My. These differences are best determined in the NW-Himalaya. In contrast to this, across the similar to 2500 km length of the northern Greater Himalaya denudation rates vary over length scales of similar to 300-1700 km. Slower denudation (<1 mm/yr) occurred between 10 and 4 Ma followed by a large increase (1.2-2.6 mm/yr) in the last similar to 4 Ma. We find that only the southern Greater Himalayan Sequence clearly supports a continuous co-evolution of tectonics, climate and denudation. Results from the higher elevation northern Greater Himalaya suggest either tectonic driven variations in denudation due to a ramp-flat geometry in the main decollement and/or recent glacially enhanced denudation.
The northward motion of the Pamir indenter with respect to Eurasia has resulted in coeval thrusting, strike-slip faulting, and normal faulting. The eastern Pamir is currently deformed by east-west oriented extension, accompanied by uplift and exhumation of the Kongur Shan (7719m) and Muztagh Ata (7546m) gneiss domes. Both domes are an integral part of the footwall of the Kongur Shan extensional fault system (KES), a 250 km long, north-south oriented graben. Why active normal faulting within the Pamir is primarily localized along the KES and not distributed more widely throughout the orogen has remained unclear. In addition, relatively little is known about how deformation has evolved throughout the Cenozoic, despite refined estimates on present-day crustal deformation rates and microseismicity, which indicate where crustal deformation is presently being accommodated. To better constrain the spatiotemporal evolution of faulting along the KES, we present 39 new apatite fission track, zircon U-Th-Sm/He, and Ar-40/Ar-39 cooling ages from a series of footwall transects along the KES graben shoulder. Combining these data with present-day topographic relief, 1-D thermokinematic and exhumational modeling documents successive stages, rather than synchronous deformation and gneiss dome exhumation. While the exhumation of the Kongur Shan commenced during the late Miocene, extensional processes in the Muztagh Ata massif began earlier and have slowed down since the late Miocene. We present a new model of synorogenic extension suggesting that thermal and density effects associated with a lithospheric tear fault along the eastern margin of the subducting Alai slab localize extensional upper plate deformation along the KES and decouple crustal motion between the central/western Pamir and eastern Pamir/Tarim basin.
Architectural analysis of the Late Cretaceous alveolinaceans of the Caribbean palaeobioprovince has made it possible to separate four genera: Praechubbina, Chubbinella gen. nov., Chubbina and Caribalveolina. The first three genera belong to the family Rhapydioninidae, while the fourth is placed in the family Alveolinidae. Two species, Praechubbina breviclaustra and P. oxchucensis sp. nov., represent the primitive genus Praechubbina, while the species cardenasensis and obesa, previously ascribed to this genus, must be reassigned respectively to Chubbinella gen. nov. and Caribalveolina. The species Chubbina jamaicensis, C. macgillavryi and C. fourcadei sp. nov. complete the inventory of Chubbina. The alveolinid genus Caribalveolina comprises two species, C. obesa and C. michaudi. Caribbean alveolinaceans include two successive assemblages. The lower assemblage is characterized by Praechubbina oxchucensis, P. brevisclaustra, Chubbinella cardenasensis and Caribalveolina obesa. The upper assemblage is represented by the genus Chubbina, with C. fourcadei, C. jamaicensis and C. macgillavryi, and Caribalveolina michaudi. The age of the lower assemblage is uncertain (probably Late CampanianEarly Maastrichtian), while the upper assemblage has been dated by strontium isotope stratigraphy as Late Maastrichtian.
A new occurrence of eclogites was found in the Kesandere valley in the eastern most part of the Bitlis complex, SE Anatolia. These high-pressure (HP) relics were preserved in calc-arenitic metasediments within the high-grade metamorphic basement of the Bitlis complex. The eclogitic parageneses were strongly overprinted during decompression and heating. These new eclogites locality complements the evidence of blueschist-facies metamorphism documented recently in the meta-sedimentary cover sequence of this part of the Bitlis complex. Thermodynamic calculations suggest peak conditions of ca. 480-540 degrees C/1.9-2.4GPa. New U/Pb dates of 84.4 +/-.9 and 82.4 +/-.9Ma were obtained on zircons from two Kesandere eclogite samples. On the basis of geochemical criteria, these dates are interpreted to represent zircon crystallization during the eclogitic peak stage. Kesandere eclogites differ from those previously described in the western Bitlis complex (Mt. Gablor locality) in terms of lithologic association, protolithic origin, and peak P-T conditions (600-650 degrees C/1.0-2.0GPa, respectively). On the other hand, eclogitic metamorphism of Kesandere metasediments occurred shortly before blueschist-facies metamorphism of the sedimentary cover (79-74Ma Ar-40/Ar-39 white mica). Therefore, the exhumation of Kesandere eclogites started between ca. 82 and 79Ma, while the meta-sedimentary cover was being buried. During this short time span, Kesandere eclogite were likely uplifted from similar to 65 to 35km depth, indicating a syn-subduction exhumation rate of similar to 4.3mm/a. Subsequently, eclogite- and blueschist-facies rocks were likely retrogressed contemporarily during collision-type metamorphism (around 72-69Ma). The Bitlis HP rocks thus sample a subduction zone that separated the Bitlis-Puturge (Bistun?) block from the South-Armenian block, further north. To the south, Eocene metasediments of the Urse formation are imbricated below the Bitlis complex. They contain (post Eocene) blueschists, testifying separation from the Arabian plate and southward migration of the subduction zone. The HT overprint of Kesandere eclogites can be related to the asthenospheric flow provoked by subducting slab retreat or break off.
This study presents a reconstruction of the Late Holocene climate in Kamchatka based on chironomid remains from a 332 cm long composite sediment core recovered from Dvuyurtochnoe Lake (Two-Yurts Lake, TYL) in central Kamchatka. The oldest recovered sediments date to about 4500 cal years BP. Chironomid head capsules from TYL reflect a rich and diverse fauna. An unknown morphotype of Tanytarsini, Tanytarsus type klein, was found in the lake sediments. Our analysis reveals four chironomid assemblage zones reflecting four different climatic periods in the Late Holocene. Between 4500 and 4000 cal years BP, the chironomid composition indicates a high lake level, well-oxygenated lake water conditions and close to modern temperatures (similar to 13 degrees C). From 4000 to 1000 cal years BP, two consecutive warm intervals were recorded, with the highest reconstructed temperature reaching 16.8 degrees C between 3700 and 2800 cal years BP. Cooling trend, started around 1100 cal years BP led to low temperatures during the last stage of the Holocene. Comparison with other regional studies has shown that termination of cooling at the beginning of late Holocene is relatively synchronous in central Kamchatka, South Kurile, Bering and Japanese Islands and take place around 3700 cal years BP. From ca 3700 cal years BP to the last millennium, a newly strengthened climate continentality accompanied by general warming trend with minor cool excursions led to apparent spatial heterogeneity of climatic patterns in the region. Some timing differences in climatic changes reconstructed from chironomid record of TYL sediments and late Holocene events reconstructed from other sites and other proxies might be linked to differences in local forcing mechanisms or caused by the different degree of dating precision, the different temporal resolution, and the different sensitive responses of climate proxies to the climate variations. Further high-resolution stratigraphic studies in this region are needed to understand the spatially complex pattern of climate change in Holocene in Kamchatka and the surrounding region.
Leaf wax n-alkanes of terrestrial plants are long-chain hydrocarbons that can persist in sedimentary records over geologic timescales. Since meteoric water is the primary source of hydrogen used in leaf wax synthesis, the hydrogen isotope composition (delta D value) of these biomarkers contains information on hydrological processes. Consequently, leaf wax n-alkane delta D values have been advocated as powerful tools for paleohydrological research. The exact kind of hydrological information that is recorded in leaf wax n-alkanes remains, however, unclear because critical processes that determine their delta D values have not yet been resolved. In particular the effects of evaporative deuterium (D)-enrichment of leaf water on the delta D values of leaf wax n-alkanes have not yet been directly assessed and quantified. Here we present the results of a study where we experimentally tested if and by what magnitude evaporative D-enrichment of leaf water affects the delta D of leaf wax n-alkanes in angiosperm C3 and C4 plants. Our study revealed that n-alkane delta D values of all plants that we investigated were affected by evaporative D-enrichment of leaf water. For dicotyledonous plants we found that the full extent of leaf water evaporative D-enrichment is recorded in leaf wax n-alkane delta D values. For monocotyledonous plants we found that between 18% and 68% of the D-enrichment in leaf water was recorded in the delta D values of their n-alkanes. We hypothesize that the different magnitudes by which evaporative D-enrichment of leaf water affects the delta D values of leaf wax n-alkanes in monocotyledonous and dicotyledonous plants is the result of differences in leaf growth and development between these plant groups. Our finding that the evaporative D-enrichment of leaf water affects the delta D values of leaf wax n-alkanes in monocotyledonous and dicotyledonous plants albeit at different magnitudes - has important implications for the interpretation of leaf wax n-alkane delta D values from paleohydrological records. In addition, our finding opens the door to employ delta D values of leaf wax n-alkanes as new ecohydrological proxies for evapotranspiration that can be applied in contemporary plant and ecosystem research.
Low molecular weight organic acids (LMWOAs) are important nutrients for microbes. However, most LMWOAs do not exist freely in the environment but are bound to macromolecular organic matter, e.g. kerogen, lignite and coal. During burial and geological maturation of sedimentary macromolecular organic matter biological and abiological processes promote the liberation of LMWOAs into the surrounding sediment. Through this process, microbes in sedimentary subsurface environments are supplied with essential nutrients. To estimate the feedstock potential of buried macromolecular organic matter to many environments it is important to determine the amount of LMWOAs that are bound to such a matrix. However, high-pressure and high temperature are a key feature of deep subsurface environments, and these physical parameters have a profound influence on chemical reaction kinetics. Therefore it is essential for the estimation of the feedstock potential to generate high-pressure and high temperature for the liberation of LMWOAs to recreate true in-situ conditions. This work presents a newly developed, inexpensive incubation system for biological and geological samples. It allows the application of high-pressure and high temperature as well as a subsampling of the liquid phase without loss of pressure, thereby not disturbing the on-going processes. When simulating the liberation of LMWOAs from sedimentary organic matter, the newly developed incubation system produces more realistic results than other extraction systems like Soxhlet. The extraction products remain in the extraction medium throughout the extraction, influencing the chemical conditions of the extraction medium. Sub-bituminous coal samples from New Zealand as well as lignite samples from Germany were extracted at elevated temperature (90˚C) and pressure (5 MPa). The main LMWOAs released from these low rank coals were formate, acetate and oxalate. Extraction efficiency was increased by two to four times for formate, acetate and oxalate in comparison to existing extraction methods without pressurisation and with demineralised water. This shows the importance of pressure for the simulation of true in-situ conditions and suggests that the amount of bioavailable LMWOAs is higher than previously thought. With the increase in carbon capture and storage (CCS) and the enhanced recovery of oil and gas (EOR/EGR), more and more CO2 becomes injected into the underground. However, the effects of elevated concentrations of carbon dioxide on sedimentary organic matter are rarely investigated. As the incuabtion system allows the manipulation of the composition and partial pressure of dissolved gasses, the effect of highly gas-enriched (CO2, CO2/SO2, CO2/NO2; to simulate flue gas conditions) waters on the extraction yield of LMWOAs from macromolecular organic matter was evaluated. For sub-bituminous coal the concentrations of all LMWAOs decreased upon the addition of gas, irrespective of its composition, whereas for lignite formate always and acetate mostly increased, while oxalate decreased. This suggests an positive effect on the nutrient supply for the subsurface microbiota of lignite layers, as formate and acetate are the most common LMWOAs used for microbial metabolism. In terrestrial mud volcanoes (TMVs), sedimentary material is rapidly ascending from great depth to the surface. Therefore LMWOAs that were produced from buried macromolecular organic matter at depth are also brought up to the surface, and fuel heterotrophic microbial ecosystems at the surface. TMVs represent geochemically and microbiologically diverse habitats, which are supplied with organic substrates and electron acceptors from deep-seated hydrocarbon-generating systems and intersected shallow aquifers, respectively. The main electron donor in TMVs in Azerbaijan is sulphate, and microbial sulphate reduction leads to the production of a wide range of reduced sulphur species that are key players in several biological processes. In our study we estimated the effect of LMWOAs on the sulphur metabolising activity of microorganims in TMVs from Azerbaijan. The addition of a mixture of volatile fatty acids containing acetate and other LMWOAs showed significant positive response to the sulphate reduction rate (SRR) of samples of several mud volcanoes. Further investigations on the temperature dependency of the SRR and the characterisation of thermophilic sulphate-reducing bacteria (SRB) showed a connection between the deep hot subsurface and the surface.
During the last glacial period, the North Atlantic region experienced pronounced, millennial-scale alternations between cold, stadial conditions and milder interstadial conditions-commonly referred to as Dansgaard-Oeschger oscillations-as well as periods of massive iceberg discharge known as Heinrich events(1). Changes in Northern Hemisphere temperature, as recorded in Greenland(2-4), are thought to have affected the location of the Atlantic intertropical convergence zone(5,6) and the strength of the Indian summer monsoon(7,8). Here we use high-resolution records of sediment colour-a measure of terrigenous versus biogenic content-from the Cariaco Basin off the coast of Venezuela and the Arabian Sea to assess teleconnections with the North Atlantic climate system during the last glacial period. The Cariaco record indicates that the intertropical convergence zone migrated seasonally over the site during mild stadial conditions, but was permanently displaced south of the basin during peak stadials and Heinrich events. In the Arabian Sea, we find evidence of a weak Indian summer monsoon during the stadial events. The tropical records show a more variable response to North Atlantic cooling than the Greenland temperature records. We therefore suggest that Greenland climate is especially sensitive to variations in the North Atlantic system-in particular sea-ice extent-whereas the intertropical convergence zone and Indian monsoon system respond primarily to variations in mean Northern Hemisphere temperature.
This study examines patterns of climate variability by mapping shoreline changes between 1973 and 2008 for two shallow lakes in Malawi, East Africa. Multi-temporal LANDSAT data covering eight different years within the period of investigation were utilised for lake-area mapping, using image classification techniques. The approach was verified using simple comparisons with local rainfall data and satellite altimetry data. Results indicated that the lake areas varied between the different years investigated and that, although the lakes were affected by the same climate-forcing mechanisms, the individual basin characteristics had a critical effect on their responses. The most likely drivers for the climate-related shoreline variations are the El Nino/Southern Oscillation (ENSO) and the Indian Ocean Dipole (IOD), causing synchronous changes to the water levels in both lakes, albeit showing differences in climate signal amplification.
The Shanderman eclogites and related metamorphosed oceanic rocks mark the site of closure of the Palaeotethys ocean in northern Iran. The protolith of the eclogites was an oceanic tholeiitic basalt with MORB composition. Eclogite occurs within a serpentinite matrix, accompanied by mafic rocks resembling a dismembered ophiolite. The eclogitic mafic rocks record different stages of metamorphism during subduction and exhumation. Minerals formed during the prograde stages are preserved as inclusions in peak metamorphic garnet and omphacite. The rocks experienced blueschist facies metamorphism on their prograde path and were metamorphosed in eclogite facies at the peak of metamorphism. The peak metamorphic mineral paragenesis of the rocks is omphacite, garnet (pyrope-rich), glaucophane, paragonite, zoisite and rutile. Based on textural relations, post-peak stages can be divided into amphibolite and greenschist facies. Pressure and temperature estimates for eclogite facies minerals (peak of metamorphism) indicate 15-20kbar at similar to 600 degrees C. The pre-peak blueschist facies assemblage yields <11kbar and 400-460 degrees C. The average pressure and temperature of the post-peak amphibolite stage was 5-6kbar, similar to 470 degrees C. The Shanderman eclogites were formed by subduction of Palaeotethys oceanic crust to a depth of no more than 75km. Subduction was followed by collision between the Central Iran and Turan blocks, and then exhumation of the high pressure rocks in northern Iran.
Climatic variations and human activity now and increasingly in the future cause land cover changes and introduce perturbations in the terrestrial carbon reservoirs in vegetation, soil and detritus. Optical remote sensing and in particular Imaging Spectroscopy has shown the potential to quantify land surface parameters over large areas, which is accomplished by taking advantage of the characteristic interactions of incident radiation and the physico-chemical properties of a material. The objective of this thesis is to quantify key soil parameters, including soil organic carbon, using field and Imaging Spectroscopy. Organic carbon, iron oxides and clay content are selected to be analyzed to provide indicators for ecosystem function in relation to land degradation, and additionally to facilitate a quantification of carbon inventories in semiarid soils. The semiarid Albany Thicket Biome in the Eastern Cape Province of South Africa is chosen as study site. It provides a regional example for a semiarid ecosystem that currently undergoes land changes due to unadapted management practices and furthermore has to face climate change induced land changes in the future. The thesis is divided in three methodical steps. Based on reflectance spectra measured in the field and chemically determined constituents of the upper topsoil, physically based models are developed to quantify soil organic carbon, iron oxides and clay content. Taking account of the benefits limitations of existing methods, the approach is based on the direct application of known diagnostic spectral features and their combination with multivariate statistical approaches. It benefits from the collinearity of several diagnostic features and a number of their properties to reduce signal disturbances by influences of other spectral features. In a following step, the acquired hyperspectral image data are prepared for an analysis of soil constituents. The data show a large spatial heterogeneity that is caused by the patchiness of the natural vegetation in the study area that is inherent to most semiarid landscapes. Spectral mixture analysis is performed and used to deconvolve non-homogenous pixels into their constituent components. For soil dominated pixels, the subpixel information is used to remove the spectral influence of vegetation and to approximate the pure spectral signature coming from the soil. This step is an integral part when working in natural non-agricultural areas where pure bare soil pixels are rare. It is identified as the largest benefit within the multi-stage methodology, providing the basis for a successful and unbiased prediction of soil constituents from hyperspectral imagery. With the proposed approach it is possible (1) to significantly increase the spatial extent of derived information of soil constituents to areas with about 40 % vegetation coverage and (2) to reduce the influence of materials such as vegetation on the quantification of soil constituents to a minimum. Subsequently, soil parameter quantities are predicted by the application of the feature-based soil prediction models to the maps of locally approximated soil signatures. Thematic maps showing the spatial distribution of the three considered soil parameters in October 2009 are produced for the Albany Thicket Biome of South Africa. The maps are evaluated for their potential to detect erosion affected areas as effects of land changes and to identify degradation hot spots in regard to support local restoration efforts. A regional validation, carried out using available ground truth sites, suggests remaining factors disturbing the correlation of spectral characteristics and chemical soil constituents. The approach is developed for semiarid areas in general and not adapted to specific conditions in the study area. All processing steps of the developed methodology are implemented in software modules, where crucial steps of the workflow are fully automated. The transferability of the methodology is shown for simulated data of the future EnMAP hyperspectral satellite. Soil parameters are successfully predicted from these data despite intense spectral mixing within the lower spatial resolution EnMAP pixels. This study shows an innovative approach to use Imaging Spectroscopy for mapping of key soil constituents, including soil organic carbon, for large areas in a non-agricultural ecosystem and under consideration of a partially vegetation coverage. It can contribute to a better assessment of soil constituents that describe ecosystem processes relevant to detect and monitor land changes. The maps further provide an assessment of the current carbon inventory in soils, valuable for carbon balances and carbon mitigation products.
The Arctic is considered as a focal region in the ongoing climate change debate. The currently observed and predicted climate warming is particularly pronounced in the high northern latitudes. Rising temperatures in the Arctic cause progressive deepening and duration of permafrost thawing during the arctic summer, creating an ‘active layer’ with high bioavailability of nutrients and labile carbon for microbial consumption. The microbial mineralization of permafrost carbon creates large amounts of greenhouse gases, including carbon dioxide and methane, which can be released to the atmosphere, creating a positive feedback to global warming. However, to date, the microbial communities that drive the overall carbon cycle and specifically methane production in the Arctic are poorly constrained. To assess how these microbial communities will respond to the predicted climate changes, such as an increase in atmospheric and soil temperatures causing increased bioavailability of organic carbon, it is necessary to investigate the current status of this environment, but also how these microbial communities reacted to climate changes in the past. This PhD thesis investigated three records from two different study sites in the Russian Arctic, including permafrost, lake shore and lake deposits from Siberia and Chukotka. A combined stratigraphic approach of microbial and molecular organic geochemical techniques were used to identify and quantify characteristic microbial gene and lipid biomarkers. Based on this data it was possible to characterize and identify the climate response of microbial communities involved in past carbon cycling during the Middle Pleistocene and the Late Pleistocene to Holocene. It is shown that previous warmer periods were associated with an expansion of bacterial and archaeal communities throughout the Russian Arctic, similar to present day conditions. Different from this situation, past glacial and stadial periods experienced a substantial decrease in the abundance of Bacteria and Archaea. This trend can also be confirmed for the community of methanogenic archaea that were highly abundant and diverse during warm and particularly wet conditions. For the terrestrial permafrost, a direct effect of the temperature on the microbial communities is likely. In contrast, it is suggested that the temperature rise in scope of the glacial-interglacial climate variations led to an increase of the primary production in the Arctic lake setting, as can be seen in the corresponding biogenic silica distribution. The availability of this algae-derived carbon is suggested to be a driver for the observed pattern in the microbial abundance. This work demonstrates the effect of climate changes on the community composition of methanogenic archae. Methanosarcina-related species were abundant throughout the Russian Arctic and were able to adapt to changing environmental conditions. In contrast, members of Methanocellales and Methanomicrobiales were not able to adapt to past climate changes. This PhD thesis provides first evidence that past climatic warming led to an increased abundance of microbial communities in the Arctic, closely linked to the cycling of carbon and methane production. With the predicted climate warming, it may, therefore, be anticipated that extensive amounts of microbial communities will develop. Increasing temperatures in the Arctic will affect the temperature sensitive parts of the current microbiological communities, possibly leading to a suppression of cold adapted species and the prevalence of methanogenic archaea that tolerate or adapt to increasing temperatures. These changes in the composition of methanogenic archaea will likely increase the methane production potential of high latitude terrestrial regions, changing the Arctic from a carbon sink to a source.
Annually laminated (varved) lake sediments with intercalated detrital layers resulting from sedimentary input by runoff events are ideal archives to establish precisely dated records of past extreme runoff events. In this study, the mid- to late Holocene varved sediments of Lake Mondsee (Upper Austria) were analysed by combining sedimentological, geophysical and geochemical methods. This approach allows to distinguish two types of detrital layers related to different types of extreme runoff events (floods and debris flows) and to detect changes in flood activity during the last 7100 years. In total, 271 flood and 47 debris flow layers, deposited during spring and summer, were identified, which cluster in 18 main flood episodes (FE 1-18) with durations of 30-50 years each. These main flood periods occurred during the Neolithic (7100-7050 vyr BP and 6470-4450 vyr BP), the late Bronze Age and the early Iron Age (3300-3250 and 2800-2750 vyr BP), the late Iron Age (2050-2000 vyr BP), throughout the Dark Ages Cold Period (1500-1200 vyr BP), and at the end of the Medieval Warm Period and the Little Ice Age (810-430 vyr BP).
Summer flood episodes in Lake Mondsee are generally more abundant during the last 1500 years, often coinciding with major advances of Alpine glaciers. Prior to 1500 vyr BP, spring/summer floods and debris flows are generally less frequent, indicating a lower number of intense rainfall events that triggered erosion. In comparison with the increase of late Holocene flood activity in western and northwestern (NW) Europe, commencing already as early as 2800 yr BP, the hydro-meteorological shift in the Lake Mondsee region occurred much later. These time lags in the onset of increased hydrological activity might be either due to regional differences in atmospheric circulation pattern or to the sensitivity of the individual flood archives. The Lake Mondsee sediments represent the first precisely dated and several millennia long summer flood record for the northeastern (NE) Alps, a key region at the climatic boundary of Atlantic, Mediterranean and East European air masses, aiding a better understanding of regional and seasonal peculiarities of flood occurrence under changing climate conditions. (C) 2013 Elsevier Ltd. All rights reserved.
Any understanding of sediment routing from mountain belts to their forelands and offshore sinks remains incomplete without estimates of intermediate storage that decisively buffers sediment yields from erosion rates, attenuates water and sediment fluxes, and protects underlying bedrock from incision. We quantify for the first time the sediment stored in > 38000 mainly postglacial Himalayan valley fills, based on an empirical volume-area scaling of valley-fill outlines automatically extracted from digital topographic data. The estimated total volume of 690(+452/-242) km(3) is mostly contained in few large valley fills > 1 km(3), while catastrophic mass wasting adds another 177(31) km(3). Sediment storage volumes are highly disparate along the strike of the orogen. Much of the Himalaya's stock of sediment is sequestered in glacially scoured valleys that provide accommodation space for similar to 44% of the total volume upstream of the rapidly exhuming and incising syntaxes. Conversely, the step-like long-wave topography of the central Himalayas limits glacier extent, and thus any significant glacier-derived storage of sediment away from tectonic basins. We show that exclusive removal of Himalayan valley fills could nourish contemporary sediment flux from the Indus and Brahmaputra basins for > 1 kyr, though individual fills may attain residence times of > 100 kyr. These millennial lag times in the Himalayan sediment routing system may sufficiently buffer signals of short-term seismic as well as climatic disturbances, thus complicating simple correlation and interpretation of sedimentary archives from the Himalayan orogen, its foreland, and its submarine fan systems. (C) 2013 Elsevier B.V. All rights reserved.
Species distribution models are useful for identifying driving environmental factors that determine earthworm distributions as well as for predicting earthworm distribution patterns and abundances at different scales. However, due to large efforts in data acquisition, studies on larger scales are rare and often focus on single species or earthworms in general. In this study, we use boosted regression tree models (BRTs) for predicting the distribution of the three functional earthworm types, i.e. anecics, endogeics and epigeics, in an agricultural area in Baden-Wurttemberg (Southwest Germany).
First, we predicted presence and absence and later earthworm abundances, considering predictors depicting land management, topography, and soil conditions as well as biotic interaction by using the abundance of the other functional earthworm types. The final presence-absence models performed reasonably well, with explained deviances between 24 and 51% after crossvalidation. Models for abundances of anecics and endogeics were less successful, since the high small-scale variability and patchiness in earthworm abundance influenced the representativeness of the field measurements. This resulted in a significant model uncertainty, which is practically very difficult to overcome with earthworm sampling campaigns at the catchment scale.
Results showed that management practices (i.e. disturbances), topography, soil conditions, and biotic interactions with other earthworm groups are the most relevant predictors for spatial distribution (incidence) patterns of all three functional groups. The response curves and contributions of predictors differ for the three functional earthworm types. Epigeics are also controlled by topographic features, endogeics by soil parameters.
Fluid flow in low-permeable carbonate rocks depends on the density of fractures, their interconnectivity and on the formation of fault damage zones. The present-day stress field influences the aperture hence the transmissivity of fractures whereas paleostress fields are responsible for the formation of faults and fractures. In low-permeable reservoir rocks, fault zones belong to the major targets. Before drilling, an estimate for reservoir productivity of wells drilled into the damage zone of faults is therefore required. Due to limitations in available data, a characterization of such reservoirs usually relies on the use of numerical techniques. The requirements of these mathematical models encompass a full integration of the actual fault geometry, comprising the dimension of the fault damage zone and of the fault core, and the individual population with properties of fault zones in the hanging and foot wall and the host rock. The paper presents both the technical approach to develop such a model and the property definition of heterogeneous fault zones and host rock with respect to the current stress field. The case study describes a deep geothermal reservoir in the western central Molasse Basin in southern Bavaria, Germany. Results from numerical simulations indicate that the well productivity can be enhanced along compressional fault zones if the interconnectivity of fractures is lateral caused by crossing synthetic and antithetic fractures. The model allows a deeper understanding of production tests and reservoir properties of faulted rocks.
Multi-proxy dating of Holocene maar lakes and Pleistocene dry maar sediments in the Eifel, Germany
(2013)
During the last twelve years the ELSA Project (Eifel Laminated Sediment Archive) at Mainz University has drilled a total of about 52 cores from 27 maar lakes and filled-in maar basins in the Eifel/Germany. Dating has been completed for the Holocene cores using 6 different methods (Pb-210 and Cs-137 activities, palynostratigraphy, event markers, varve counting, C-14) In general, the different methods consistently complement one another within error margins. Event correlation was used for relating typical lithological changes with historically known events such as the two major Holocene flood events at 1342 AD and ca 800 BC. Dating of MIS2-MIS3 core sections is based on greyscale tuning, radiocarbon and OSL dating, magnetostratigraphy and tephrochronology. The lithological changes in the sediment cores demonstrate a sequence of events similar to the North Atlantic rapid climate variability of the Last Glacial Cycle. The warmest of the MIS3 interstadials was GI14, when a forest with abundant spruce covered the Eifel area from 55 to 48 ka BP, i.e. during a time when also other climate archives in Europe suggested very warm conditions. The forest of this "Early Stage 3 warm phase" developed subsequently into a steppe with scattered birch and pine, and finally into a glacial desert at around 25 ka BP. Evidence for Mono Lake and Laschamp geomagnetic excursions is found in two long cores. Several large eruptions during Middle and Late Pleistocene (Ulmener Maar - 11,000 varve years BP, Laacher See - 12,900 varve years BP, Mosenberg volcanoes/Meerfelder Maar 41-45 cal ka BP, Dumpel Maar 116 ka BP, Glees Maar - 151 ka BP) produced distinct ash-layers crucial for inter-core and inter-site correlations. The oldest investigated maar of the Eifel is Ar-40/Ar-39 dated to the time older than 520 ka BP.
The persistence of topography within ancient orogens remains one of the outstanding questions in landscape evolution. In the eastern North American Appalachians, this question is manifest in the outstanding problem of whether topographic relief is in a quasi-equilibrium state, decaying slowly over many millennia, or whether relief has increased during the late Cenozoic. Here we present quantitative geomorphic data from the nonglaciated portion of the Susquehanna River drainage basin that provide insight into these end-member models. Analysis of channel profiles draining upland catchments in the northern Valley and Ridge, Appalachian Plateau, Blue Ridge, and Piedmont provinces reveals that a large number of streams have well defined knickpoints clustered at 300-600 m elevation but not systematically associated with transitions from weak to resistant substrate. Cosmogenic Be-10 inventories of modern stream sediment indicate that erosion rates are spatially variable, ranging from similar to 5-30 m/Myr above knickpoints to similar to 50-100 m/Myr below knickpoints. Overall, channel gradients, normalized for drainage area, scale linearly with catchment-averaged erosion rates. Collectively, regionally consistent spatial relationships among erosion rate, channel steepness, and knickpoints reveal an ongoing wave of transient channel adjustment to a change in relative base level. Reconstructions of relict channel profiles above knickpoints suggest that higher rates of incision are associated with similar to 100-150 m of relative base level fall that accompanied epierogenic rock uplift rather than a change to a more erosive climate or drainage reorganization. Channel response timescales imply that the onset of relative base level change predates similar to 3.5 Ma and may have begun as early as similar to 15 Ma. We suggest that adjustment of the channel network was likely driven by changes in mantle dynamics along the eastern seaboard of North America during the Neogene.
The multiple high-pressure (HP), low-temperature (LT) metamorphic units of Western and Central Anatolia offer a great opportunity to investigate the subduction-and continental accretion-related evolution of the eastern limb of the long-lived Aegean subduction system. Recent reports of the HP-LT index mineral Fe-Mg-carpholite in three metasedimentary units of the Gondwana-derived Anatolide-Tauride continental block (namely the Afyon Zone, the Oren Unit and the southern Menderes Massif) suggest a more complicated scenario than the single-continental accretion model generally put forward in previous studies. This study presents the first isotopic dates (white mica Ar-40-Ar-39 geochronology), and where possible are combined with P-T estimates (chlorite thermometry, phengite barometry, multi-equilibrium thermobarometry), on carpholite-bearing rocks from these three HP-LT metasedimentary units. It is shown that, in the Afyon Zone, carpholite-bearing assemblages were retrogressed through greenschist-facies conditions at c. 67-62 Ma. Early retrograde stages in the Oren Unit are dated to 63-59 Ma. In the Kurudere-Nebiler Unit (HP Mesozoic cover of the southern Menderes Massif), HP retrograde stages are dated to c. 45 Ma, and post-collisional cooling to c. 26 Ma. These new results support that the Oren Unit represents the westernmost continuation of the Afyon Zone, whereas the Kurudere-Nebiler Unit correlates with the Cycladic Blueschist Unit of the Aegean Domain. In Western Anatolia, three successive HP-LT metamorphic belts thus formed: the northernmost Tavsanli Zone (c. 88-82 Ma), the Oren-Afyon Zone (between 70 and 65 Ma), and the Kurudere-Nebiler Unit (c. 52-45 Ma). The southward younging trend of the HP-LT metamorphism from the upper and internal to the deeper and more external structural units, as in the Aegean Domain, points to the persistence of subduction in Western Anatolia between 93-90 and c. 35 Ma. After the accretion of the Menderes-Tauride terrane, in Eocene times, subduction stopped, leading to continental collision and associated Barrovian-type metamorphism. Because, by contrast, the Aegean subduction did remain active due to slab roll-back and trench migration, the eastern limb (below Southwestern Anatolia) of the Hellenic slab was dramatically curved and consequently teared. It therefore is suggested that the possibility for subduction to continue after the accretion of buoyant (e.g. continental) terranes probably depends much on palaeogeography.
From the Middle Eocene to Early Oligocene, the Earth experienced the most significant climatic cooling of the Cenozoic era. The Eocene-Oligocene transition (EOT) represents the culmination of this climatic cooling, leading to the onset of the Antarctic glaciation and, consequently, to the beginning of the present-day icehouse world. Whereas the response of deep-sea systems to this climate transition has been widely studied, its impact on the shallow-water carbonate realm is poorly constrained. Here, the sedimentary expression of the EOT in two shallow-marine carbonate successions (Nago and San Valentino, northern Italy) belonging to the Calcare di Nago Formation is presented. The chronostratigraphic framework was constructed by integrating litho-, bio-, and isotope-stratigraphic data (C and Sr isotopes), allowing to correlate these shallow-marine successions with pelagic sections in central Italy (Massignano), Tanzania (TOP Sites 12 and 17), and the Indian Ocean (ODP Site 744). Within several sections in northern Italy, including Nago and San Valentino, a Priabonian (Late Eocene) transgression is recorded. Oxygen isotopes of ODP Site 744 show a coeval negative shift of 0.4 parts per thousand., suggesting a glacio-eustatic origin for this transgression. In the Nago and San Valentino sections, no prominent sequence boundary has been detected that would indicate a rapid sea-level drop occurring together with the positive shift in delta O-18 defining the EOT-1 cooling event. Instead, a gradual shallowing of the depositional environment is observed. At TDP Sites 12 and 17, the EOT-1 is followed by a negative shift in delta O-18 of around 0.4 parts per thousand, which correlates with a relative deepening of the environment in the studied sections and suggests a melting pulse between EOT-1 and the Oligocene isotope event 1 (Oi-1). The positive delta O-18 shift related to the Oi-1 translates in San Valentino into a change in carbonate factory from a photozoan association dominated by larger benthic foraminifera, corals, and red algae to a heterozoan association dominated by bryozoans. The same bryozoan fades occurs in several Italian localities near the Eocene-Oligocene boundary. This fades is interpreted to represent an analogue of modern cool-water carbonates and results from a cooling pulse of at least regional scale, associated to the Oi-1 event.
Vertical flow systems filled with porous medium have been shown to efficiently remove volatile organic contaminants (VOCs) from contaminated groundwater. To apply this semi-natural remediation strategy it is however necessary to distinguish between removal due to biodegradation and due to volatile losses to the atmosphere. Especially for (potentially) toxic VOCs, the latter needs to be minimized to limit atmospheric emissions. In this study, numerical simulation was used to investigate quantitatively the removal of volatile organic compounds in two pilot-scale water treatment systems: an unplanted vertical flow filter and a planted one, which could also be called a vertical flow constructed wetland, both used for the treatment of contaminated groundwater. These systems were intermittently loaded with contaminated water containing benzene and MTBE as main VOCs. The highly dynamic but permanently unsaturated conditions in the porous medium facilitated aerobic biodegradation but could lead to volatile emissions of the contaminants. Experimental data from porous material analyses, flow rate measurements, solute tracer and gas tracer test, as well as contaminant concentration measurements at the boundaries of the systems were used to constrain a numerical reactive transport modeling approach. Numerical simulations considered unsaturated water flow, transport of species in the aqueous and the gas phase as well as aerobic degradation processes, which made it possible to quantify the rates of biodegradation and volatile emissions and calculating their contribution to total contaminant removal. A range of degradation rates was determined using experimental results of both systems under two operation modes and validated by field data obtained at different operation modes applied to the filters. For both filters, simulations and experimental data point to high biodegradation rates, if the flow filters have had time to build up their removal capacity. For this case volatile emissions are negligible and total removal can be attributed to biodegradation, only. The simulation study thus supports the use of both of these vertical flow systems for the treatment of groundwater contaminated with VOCs and the use of reactive transport modeling for the assessment of VOCs removal and operation modes in these high performance treatment systems.