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) (remove)
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)
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.
Metabolically active microbial communities are present in a wide range of subsurface environments. Techniques like enumeration of microbial cells, activity measurements with radiotracer assays and the analysis of porewater constituents are currently being used to explore the subsurface biosphere, alongside with molecular biological analyses. However, many of these techniques reach their detection limits due to low microbial activity and abundance. Direct measurements of microbial turnover not just face issues of insufficient sensitivity, they only provide information about a single specific process but in sediments many different process can occur simultaneously. Therefore, the development of a new technique to measure total microbial activity would be a major improvement. A new tritium-based hydrogenase-enzyme assay appeared to be a promising tool to quantify total living biomass, even in low activity subsurface environments. In this PhD project total microbial biomass and microbial activity was quantified in different subsurface sediments using established techniques (cell enumeration and pore water geochemistry) as well as a new tritium-based hydrogenase enzyme assay. By using a large database of our own cell enumeration data from equatorial Pacific and north Pacific sediments and published data it was shown that the global geographic distribution of subseafloor sedimentary microbes varies between sites by 5 to 6 orders of magnitude and correlates with the sedimentation rate and distance from land. Based on these correlations, global subseafloor biomass was estimated to be 4.1 petagram-C and ~0.6 % of Earth's total living biomass, which is significantly lower than previous estimates. Despite the massive reduction in biomass the subseafloor biosphere is still an important player in global biogeochemical cycles. To understand the relationship between microbial activity, abundance and organic matter flux into the sediment an expedition to the equatorial Pacific upwelling area and the north Pacific Gyre was carried out. Oxygen respiration rates in subseafloor sediments from the north Pacific Gyre, which are deposited at sedimentation rates of 1 mm per 1000 years, showed that microbial communities could survive for millions of years without fresh supply of organic carbon. Contrary to the north Pacific Gyre oxygen was completely depleted within the upper few millimeters to centimeters in sediments of the equatorial upwelling region due to a higher supply of organic matter and higher metabolic activity. So occurrence and variability of electron acceptors over depth and sites make the subsurface a complex environment for the quantification of total microbial activity. Recent studies showed that electron acceptor processes, which were previously thought to thermodynamically exclude each other can occur simultaneously. So in many cases a simple measure of the total microbial activity would be a better and more robust solution than assays for several specific processes, for example sulfate reduction rates or methanogenesis. Enzyme or molecular assays provide a more general approach as they target key metabolic compounds. Since hydrogenase enzymes are ubiquitous in microbes, the recently developed tritium-based hydrogenase radiotracer assay is applied to quantify hydrogenase enzyme activity as a parameter of total living cell activity. Hydrogenase enzyme activity was measured in sediments from different locations (Lake Van, Barents Sea, Equatorial Pacific and Gulf of Mexico). In sediment samples that contained nitrate, we found the lowest cell specific enzyme activity around 10^(-5) nmol H_(2) cell^(-1) d^(-1). With decreasing energy yield of the electron acceptor used, cell-specific hydrogenase activity increased and maximum values of up to 1 nmol H_(2) cell^(-1) d^(-1) were found in samples with methane concentrations of >10 ppm. Although hydrogenase activity cannot be converted directly into a turnover rate of a specific process, cell-specific activity factors can be used to identify specific metabolism and to quantify the metabolically active microbial population. In another study on sediments from the Nankai Trough microbial abundance and hydrogenase activity data show that both the habitat and the activity of subseafloor sedimentary microbial communities have been impacted by seismic activities. An increase in hydrogenase activity near the fault zone revealed that the microbial community was supplied with hydrogen as an energy source and that the microbes were specialized to hydrogen metabolism.
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.
Considerable effort has been devoted to the development of simulation algorithms for facies modeling, whereas a discussion of how to combine those techniques has not existed. The integration of multiple geologic data into a three-dimensional model, which requires the combination of simulation techniques, is yet a current challenge for reservoir modeling. This article presents a thought process that guides the acquisition and modeling of geologic data at various scales. Our work is based on outcrop data collected from a Jurassic carbonate ramp located in the High Atlas mountain range of Morocco. The study window is 1 km (0.6 mi) wide and 100 m (328.1 ft) thick. We describe and model the spatial and hierarchical arrangement of carbonate bodies spanning from largest to smallest: (1) stacking pattern of high-frequency depositional sequences, (2) facies association, and (3) lithofacies. Five sequence boundaries were modeled using differential global position system mapping and light detection and ranging data. The surface-based model shows a low-angle profile with modest paleotopographic relief at the inner-to-middle ramp transition. Facies associations were populated using truncated Gaussian simulation to preserve ordered trends between the inner, middle, and outer ramps. At the lithofacies scale, field observations and statistical analysis show a mosaiclike distribution that was simulated using a fully stochastic approach with sequential indicator simulation.
This study observes that the use of one single simulation technique is unlikely to correctly model the natural patterns and variability of carbonate rocks. The selection and implementation of different techniques customized for each level of the stratigraphic hierarchy will provide the essential computing flexibility to model carbonate settings. This study demonstrates that a scale-dependent modeling approach should be a common procedure when building subsurface and outcrop models.
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.
Palaeoenvironmental implications of evaporative gaylussite crystals from Lonar Lake, central India
(2013)
We have undertaken petrographic, mineralogical, geochemical and isotopic investigations on carbonate minerals found within a 10-m-long core from Lonar Lake, central India, with the aim of evaluating their potential as palaeoenvironmental proxies. The core encompasses the entire Holocene and is the first well-dated high-resolution record from central India. While calcite and/or aragonite were found throughout the core, the mineral gaylussite was found only in two specific intervals (46303890 and 2040560 cal a BP). Hydrochemical and isotope data from inflowing streams and lake waters indicate that evaporitic processes play a dominant role in the precipitation of carbonates within this lake. Isotopic (18O and 13C) studies on the evaporative gaylussite crystals and residual bulk carbonates (calcite) from the long core show that evaporation is the major control on 18O enrichment in both the minerals. However, in case of 13C additional mechanisms, for example methanogenesis (gaylussite) and phytoplankton productivity (calcium carbonate), play an additional important role in some intervals. We also discuss the relevance of our investigation for palaeoclimate reconstruction and late Holocene monsoon variability.
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.
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.
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.
With this paper, we assess the present-day conductive thermal field of the Glueckstadt Graben in NW Germany that is characterized by large salt walls and diapirs structuring the graben fill. We use a finite element method to calculate the 3D steady-state conductive thermal field based on a lithosphere-scale 3D structural model that resolves the first-order structural characteristics of the graben and its underlying lithosphere. Model predictions are validated against measured temperatures in six deep wells. Our investigations show that the interaction of thickness distributions and thermal rock properties of the different geological layers is of major importance for the distribution of temperatures in the deep subsurface of the Glueckstadt Graben. However, the local temperatures may result from the superposed effects of different controlling factors. Especially, the upper sedimentary part of the model exhibits huge lateral temperature variations, which correlate spatially with the shape of the thermally highly conductive Permian salt layer. Variations in thickness and geometry of the salt cause two major effects, which provoke considerable lateral temperature variations for a given depth. (1) The "chimney effect" causes more efficient heat transport within salt diapirs. As a consequence positive thermal anomalies develop in the upper part and above salt structures, where the latter are covered by much less conductive sediments. In contrast, negative thermal anomalies are noticeable underneath salt structures. (2) The "thermal blanketing effect" is caused by thermally low conductive sediments that provoke the local storage of heat where these insulating sediments are present. The latter effect leads to both local and regional thermal anomalies. Locally, this translates to higher temperatures where salt margin synclines are filled with thick insulating clastic sediments. For the regional anomalies the cumulative insulating effects of the entire sediment fill results in a long-wavelength variation of temperatures in response to heat refraction effects caused by the contrast between insulating sediments and highly conductive crystalline crust. Finally, the longest wavelength of temperature variations is caused by the depth position of the isothermal lithosphere-asthenosphere boundary defining the regional variations of the overall geothermal gradient. We find that a conductive thermal model predicts observed temperatures reasonably well for five of the six available wells, whereas the steady-state conductive approach appears not to be valid for the sixth well.
Soil moisture dynamics are affected by complex interactions among several factors. Understanding the relative importance of these factors is still an important challenge in the study of water fluxes and solute transport in unsaturated media. In this study, the spatio-temporal variability of surface soil moisture was investigated in a 10 ha flat cropped field located in northern Italy. Soil moisture was measured on a regular 50 x 50 m grid on seven dates during the growing season. For each measurement campaign, the spatial variability of the soil moisture was compared with the spatial variability of the soil texture and crop properties. In particular, to better understand the role of the vegetation, the spatio-temporal variability of two different parameters - leaf area index and crop height - was monitored on eight dates at different crop development stages. Statistical and geostatistical analysis was then applied to explore the interactions between these variables. In agreement with other studies, the results show that the soil moisture variability changes according to the average value within the field, with the standard deviation reaching a maximum value under intermediate mean soil moisture conditions and the coefficient of variation decreasing exponentially with increasing mean soil moisture. The controls of soil moisture variability change according to the average soil moisture within the field. Under wet conditions, the spatial distribution of the soil moisture reflects the variability of the soil texture. Under dry conditions, the spatial distribution of the soil moisture is affected mostly by the spatial variability of the vegetation. The interaction between these two factors is more important under intermediate soil moisture conditions. These results confirm the importance of considering the average soil moisture conditions within a field when investigating the controls affecting the spatial variability of soil moisture. This study highlights the importance of considering the spatio-temporal variability of the vegetation in investigating soil moisture dynamics, especially under intermediate and dry soil moisture conditions. The results of this study have important implications in different hydrological applications, such as for sampling design, ranking stability application, indirect measurements of soil properties and model parameterisation.
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.
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.
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.
Beta diversity is a conceptual link between diversity at local and regional scales. Various additional methodologies of quantifying this and related phenomena have been applied. Among them, measures of pairwise (dis)similarity of sites are particularly popular. Undersampling, i.e. not recording all taxa present at a site, is a common situation in ecological data. Bias in many metrics related to beta diversity must be expected, but only few studies have explicitly investigated the properties of various measures under undersampling conditions. On the basis of an empirical data set, representing near-complete local inventories of the Lepidoptera from an isolated Pacific island, as well as simulated communities with varying properties, we mimicked different levels of undersampling. We used 14 different approaches to quantify beta diversity, among them dataset-wide multiplicative partitioning (i.e. true beta diversity') and pairwise site x site dissimilarities. We compared their values from incomplete samples to true results from the full data. We used these comparisons to quantify undersampling bias and we calculated correlations of the dissimilarity measures of undersampled data with complete data of sites. Almost all tested metrics showed bias and low correlations under moderate to severe undersampling conditions (as well as deteriorating precision, i.e. large chance effects on results). Measures that used only species incidence were very sensitive to undersampling, while abundance-based metrics with high dependency on the distribution of the most common taxa were particularly robust. Simulated data showed sensitivity of results to the abundance distribution, confirming that data sets of high evenness and/or the application of metrics that are strongly affected by rare species are particularly sensitive to undersampling. The class of beta measure to be used should depend on the research question being asked as different metrics can lead to quite different conclusions even without undersampling effects. For each class of metric, there is a trade-off between robustness to undersampling and sensitivity to rare species. In consequence, using incidence-based metrics carries a particular risk of false conclusions when undersampled data are involved. Developing bias corrections for such metrics would be desirable.
The Acheulean technological tradition, characterized by a large (>10 cm) flake-based component, represents a significant technological advance over the Oldowan. Although stone tool assemblages attributed to the Acheulean have been reported from as early as circa 1.6-1.75 Ma, the characteristics of these earliest occurrences and comparisons with later assemblages have not been reported in detail. Here, we provide a newly established chronometric calibration for the Acheulean assemblages of the Konso Formation, southern Ethiopia, which span the time period similar to 1.75 to <1.0 Ma. The earliest Konso Acheulean is chronologically indistinguishable from the assemblage recently published as the world's earliest with an age of similar to 1.75 Ma at Kokiselei, west of Lake Turkana, Kenya. This Konso assemblage is characterized by a combination of large picks and crude bifaces/unifaces made predominantly on large flake blanks. An increase in the number of flake scars was observed within the Konso Formation handaxe assemblages through time, but this was less so with picks. The Konso evidence suggests that both picks and handaxes were essential components of the Acheulean from its initial stages and that the two probably differed in function. The temporal refinement seen, especially in the handaxe forms at Konso, implies enhanced function through time, perhaps in processing carcasses with long and stable cutting edges. The documentation of the earliest Acheulean at similar to 1.75 Ma in both northern Kenya and southern Ethiopia suggests that behavioral novelties were being established in a regional scale at that time, paralleling the emergence of Homo erectus-like hominid morphology.
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.
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.
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.
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.
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.
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.
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.
Flood risk is expected to increase in many regions of the world in the next decades with rising flood losses as a consequence. First and foremost, it can be attributed to the expansion of settlement and industrial areas into flood plains and the resulting accumulation of assets. For a future-oriented and a more robust flood risk management, it is therefore of importance not only to estimate potential impacts of climate change on the flood hazard, but also to analyze the spatio-temporal dynamics of flood exposure due to land use changes. In this study, carried out in the Alpine Lech Valley in Tyrol (Austria), various land use scenarios until 2030 were developed by means of a spatially explicit land use model, national spatial planning scenarios and current spatial policies. The combination of the simulated land use patterns with different inundation scenarios enabled us to derive statements about possible future changes in flood-exposed built-up areas. The results indicate that the potential assets at risk depend very much on the selected socioeconomic scenario. The important conditions affecting the potential assets at risk that differ between the scenarios are the demand for new built-up areas as well as on the types of conversions allowed to provide the necessary areas at certain locations. The range of potential changes in flood-exposed residential areas varies from no further change in the most moderate scenario 'Overall Risk' to 119 % increase in the most extreme scenario 'Overall Growth' (under current spatial policy) and 159 % increase when disregarding current building restrictions.
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 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.
On 2012 May 20 and 29, two damaging earthquakes with magnitudes M-w 6.1 and 5.9, respectively, struck the Emilia-Romagna region in the sedimentary Po Plain, Northern Italy, causing 26 fatalities, significant damage to historical buildings and substantial impact to the economy of the region. The earthquake sequence included four more aftershocks with M-w, >= 5.0, all at shallow depths (about 7-9 km), with similar WNW-ESE striking reverse mechanism. The timeline of the sequence suggests significant static stress interaction between the largest events. We perform here a detailed source inversion, first adopting a point source approximation and considering pure double couple and full moment tensor source models. We compare different extended source inversion approaches for the two largest events, and find that the rupture occurred in both cases along a subhorizontal plane, dipping towards SSW Directivity is well detected for the May 20 main shock, indicating that the rupture propagated unilaterally towards SE. Based on the focal mechanism solution, we further estimate the co-seismic static stress change induced by the May 20 event. By using the rate-and-state model and a Poissonian earthquake occurrence, we infer that the second largest event of May 29 was induced with a probability in the range 0.2-0.4. This suggests that the segment of fault was already prone to rupture. Finally, we estimate peak ground accelerations for the two main events as occurred separately or simultaneously. For the scenario involving hypothetical rupture areas of both main events, we estimate M-w = 6.3 and an increase of ground acceleration by 50 per cent. The approach we propose may help to quantify rapidly which regions are invested by a significant increase of the hazard, bearing the potential for large aftershocks or even a second main shock.
Preface to the special issue "Triggered and induced seismicity: probabilities and discrimination"
(2013)
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.
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.
An important tool for understanding deformation occurring within a subduction zone is the measurement of seismic anisotropy through observations of shear wave splitting (SWS). In Sumatra, two temporary seismic networks were deployed between December 2007 and February 2009, covering the fore arc between the fore-arc islands to the back arc. We use SKS and local SWS measurements to determine the type, amount, and location of anisotropy. Local SWS measurements from the fore-arc islands exhibit trench-parallel fast directions which can be attributed to shape preferred orientation of cracks/fractures in the overriding sediments. In the Sumatran Fault region, the predominant fast direction is fault/trench parallel, while in the back-arc region it is trench perpendicular. The trench-perpendicular measurements exhibit a positive correlation between delay time and raypath length in the mantle wedge, while the fault-parallel measurements are similar to the fault-parallel fast directions observed for two crustal events at the Sumatran Fault. This suggests that there are two layers of anisotropy: one due to entrained flow within the mantle wedge and a second layer within the overriding crust due to the shear strain caused by the Sumatran Fault. SKS splitting results show a NNW-SSE fast direction with delay times of 0.8-3.0s. The fast directions are approximately parallel to the absolute plate motion of the subducting Indo-Australian Plate. The small delay times exhibited by the local SWS (0.05-0.45s), in combination with the large SKS delay times, suggest that the anisotropy generating the teleseismic SWS is dominated by entrained flow in the asthenosphere below the slab.
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.
A problem encountered by many distributed hydrological modelling studies is high simulation errors at interior gauges when the model is only globally calibrated at the outlet. We simulated river runoff in the Elbe River basin in central Europe (148 268 km(2)) with the semi-distributed eco-hydrological model SWIM (Soil and Water Integrated Model). While global parameter optimisation led to Nash-Sutcliffe efficiencies of 0.9 at the main outlet gauge, comparisons with measured runoff series at interior points revealed large deviations. Therefore, we compared three different strategies for deriving sub-basin evapotranspiration: (1) modelled by SWIM without any spatial calibration, (2) derived from remotely sensed surface temperatures, and (3) calculated from long-term precipitation and discharge data. The results show certain consistencies between the modelled and the remote sensing based evapotranspiration rates, but there seems to be no correlation between remote sensing and water balance based estimations. Subsequent analyses for single sub-basins identify amongst others input weather data and systematic error amplification in inter-gauge discharge calculations as sources of uncertainty. The results encourage careful utilisation of different data sources for enhancements in distributed hydrological modelling.
Astronomically tuned cyclic sedimentary successions provide unprecedented insight into the temporal evolution of depositional systems and major geologic events. However, placing astronomically calibrated records into an absolute time frame with confidence requires independent and precise geochronologic constraints. Astronomical tuning of the precessionally modulated sedimentary cycles of the Mediterranean Basin deposited during the Messinian Salinity Crisis (5.96-5.33 Ma) has indicated an similar to 90 k.y. "Messinian gap", corresponding to the evaporative drawdown of the Mediterranean following the closure of the Mediterranean-Atlantic gateway. In the Messinian deposits, a volcanic ash dated by Ar-40/Ar-39 geochronology was used to anchor the sedimentary cycles to the insolation curve. However, the uncertainty of the Ar-40/Ar-39 date introduces a potential two-cycle (similar to 40 k.y.) uncertainty in the tuning. Using high-precision chemical abrasion-thermal ionization mass spectrometry (CA-TIMS) U-Pb geochronology on single zircon grains from two Messinian ash layers in Italy, we obtained dates of 5.5320 +/- 0.0046 Ma and 5.5320 +/- 0.0074 Ma with sub-precessional resolution. Combined with our astronomical tuning of the Messinian Lower Evaporites, the results refine the duration of the "Messinian gap" to at most 28 or 58 +/- 9.6 k.y., which correlates with either the TG12 glacial interval alone, or both TG12 and TG14 glacial intervals, supporting the hypothesis of a glacio-eustatic contribution in fully isolating the Mediterranean from the Atlantic Ocean. Our new U-Pb dates also allow us to infer a precessionally modulated cyclicity for the post-evaporitic deposits, and hence enable us to tune those successions to the insolation curve.
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.
Various techniques are utilized by the seismological community, extractive industries, energy and geoengineering companies to identify earthquake nucleation processes in close proximity to engineering operation points. These operations may comprise fluid extraction or injections, artificial water reservoir impoundments, open pit and deep mining, deep geothermal power generations or carbon sequestration. In this letter to the editor, we outline several lines of investigation that we suggest to follow to address the discrimination problem between natural seismicity and seismic events induced or triggered by geoengineering activities. These suggestions have been developed by a group of experts during several meetings and workshops, and we feel that their publication as a summary report is helpful for the geoscientific community. Specific investigation procedures and discrimination approaches, on which our recommendations are based, are also published in this Special Issue (SI) of Journal of Seismology.
This study provides a detailed analysis of the mid-Holocene to present-day precipitation change in the Asian monsoon region. We compare for the first time results of high resolution climate model simulations with a standardised set of mid-Holocene moisture reconstructions. Changes in the simulated summer monsoon characteristics (onset, withdrawal, length and associated rainfall) and the mechanisms causing the Holocene precipitation changes are investigated. According to the model, most parts of the Indian subcontinent received more precipitation (up to 5 mm/day) at mid-Holocene than at present-day. This is related to a stronger Indian summer monsoon accompanied by an intensified vertically integrated moisture flux convergence. The East Asian monsoon region exhibits local inhomogeneities in the simulated annual precipitation signal. The sign of this signal depends on the balance of decreased pre-monsoon and increased monsoon precipitation at mid-Holocene compared to present-day. Hence, rainfall changes in the East Asian monsoon domain are not solely associated with modifications in the summer monsoon circulation but also depend on changes in the mid-latitudinal westerly wind system that dominates the circulation during the pre-monsoon season. The proxy-based climate reconstructions confirm the regional dissimilarities in the annual precipitation signal and agree well with the model results. Our results highlight the importance of including the pre-monsoon season in climate studies of the Asian monsoon system and point out the complex response of this system to the Holocene insolation forcing. The comparison with a coarse climate model simulation reveals that this complex response can only be resolved in high resolution simulations.
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.
Explaning change in flood hazard in the Mekong river : the hypothesis of nonstationary variance
(2013)
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.
The M (w) 6.2 Baladeh earthquake occurred on 28 May 2004 in the Alborz Mountains, northern Iran. This earthquake was the first strong shock in this intracontinental orogen for which digital regional broadband data are available. The Baladeh event provides a rare opportunity to study fault geometry and ongoing deformation processes using modern seismological methods. A joint inversion for hypocentres and a velocity model plus a surface-wave group dispersion curve analysis were used to obtain an adapted velocity model, customised for mid- and long-period waveform modelling. Based on the new velocity model, regional waveform data of the mainshock and larger aftershocks (M (w) a parts per thousand yen3.3) were inverted for moment tensors. For the Baladeh mainshock, this included inversion for kinematic parameters. All analysed earthquakes show dominant thrust mechanisms at depths between 14 and 26 km, with NW-SE striking fault planes. The mainshock ruptured a 28A degrees south-dipping area of 24 x 21 km along a north-easterly direction. The rupture plane of the mainshock does not coincide with the aftershock distribution, neither in map view nor with respect to depth. The considered aftershocks form two main clusters. The eastern cluster is associated with the mainshock. The western cluster does not appear to be connected with the rupture plane of the mainshock but, instead, indicates a second activated fault plane dipping at 85A degrees towards the north.
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.
This study presents the theory, applicability, and merits of the new THERIAK_D add-on for the open source Theriak/Domino software package. The add-on works as an interface between Theriak and user-generated scripts, providing the opportunity to process phase equilibrium computation parameters in a programming environment (e. g., C or MATLABV (R)). THERIAK_D supports a wide range of features such as calculating the solid rock density or testing the stability of mineral phases along any pressure-temperature (P-T) path and P-T grid. To demonstrate applicability, an example is given in which the solid rock density of a 2-D-temperature-pressure field is calculated, portraying a simplified subduction zone. Consequently, the add-on effectively combines thermodynamics and geodynamic modeling. The carefully documented examples could be easily adapted for a broad range of applications. THERIAK_D is free, and the program, user manual, and source codes may be downloaded from http://www.min.unikiel.de/similar to ed/theriakd/.
Background: Short lived, iteroparous animals in seasonal environments experience variable social and environmental conditions over their lifetime. Animals can be divided into those with a "young-of-the-year" life history (YY, reproducing and dying in the summer of birth) and an "overwinter" life history (OW, overwintering in a subadult state before reproducing next spring).
We investigated how behavioural patterns across the population were affected by season and sex, and whether variation in behaviour reflects the variation in life history patterns of each season. Applications of pace-of-life (POL) theory would suggest that long-lived OW animals are shyer in order to increase survival, and YY are bolder in order to increase reproduction. Therefore, we expected that in winter and spring samples, when only OW can be sampled, the animals should be shyer than in summer and autumn, when both OW and YY animals can be sampled. We studied common vole (Microtus arvalis) populations, which express typical, intra-annual density fluctuation. We captured a total of 492 voles at different months over 3 years and examined boldness and activity level with two standardised behavioural experiments.
Results: Behavioural variables of the two tests were correlated with each other. Boldness, measured as short latencies in both tests, was extremely high in spring compared to other seasons. Activity level was highest in spring and summer, and higher in males than in females.
Conclusion: Being bold in laboratory tests may translate into higher risk-taking in nature by being more mobile while seeking out partners or valuable territories. Possible explanations include asset-protection, with OW animals being rather old with low residual reproductive value in spring. Therefore, OW may take higher risks during this season. Offspring born in spring encounter a lower population density and may have higher reproductive value than offspring of later cohorts. A constant connection between life history and animal personality, as suggested by the POL theory, however, was not found. Nevertheless, correlations of traits suggest the existence of animal personalities. In conclusion, complex patterns of population dynamics, seasonal variation in life histories, and variability of behaviour due to asset-protection may cause complex seasonal behavioural dynamics in a population.
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.
In sedimentary basins, rock thermal conductivity can vary both laterally and vertically, thus altering the basin’s thermal structure locally and regionally. Knowledge of the thermal conductivity of geological formations and its spatial variations is essential, not only for quantifying basin evolution and hydrocarbon maturation processes, but also for understanding geothermal conditions in a geological setting. In conjunction with the temperature gradient, thermal conductivity represents the basic input parameter for the determination of the heat-flow density; which, in turn, is applied as a major input parameter in thermal modeling at different scales. Drill-core samples, which are necessary to determine thermal properties by laboratory measurements, are rarely available and often limited to previously explored reservoir formations. Thus, thermal conductivities of Mesozoic rocks in the North German Basin (NGB) are largely unknown. In contrast, geophysical borehole measurements are often available for the entire drilled sequence. Therefore, prediction equations to determine thermal conductivity based on well-log data are desirable. In this study rock thermal conductivity was investigated on different scales by (1) providing thermal-conductivity measurements on Mesozoic rocks, (2) evaluating and improving commonly applied mixing models which were used to estimate matrix and pore-filled rock thermal conductivities, and (3) developing new well-log based equations to predict thermal conductivity in boreholes without core control. Laboratory measurements are performed on sedimentary rock of major geothermal reservoirs in the Northeast German Basin (NEGB) (Aalenian, Rhaethian-Liassic, Stuttgart Fm., and Middle Buntsandstein). Samples are obtained from eight deep geothermal wells that approach depths of up to 2,500 m. Bulk thermal conductivities of Mesozoic sandstones range between 2.1 and 3.9 W/(m∙K), while matrix thermal conductivity ranges between 3.4 and 7.4 W/(m∙K). Local heat flow for the Stralsund location averages 76 mW/m², which is in good agreement to values reported previously for the NEGB. For the first time, in-situ bulk thermal conductivity is indirectly calculated for entire borehole profiles in the NEGB using the determined surface heat flow and measured temperature data. Average bulk thermal conductivity, derived for geological formations within the Mesozoic section, ranges between 1.5 and 3.1 W/(m∙K). The measurement of both dry- and water-saturated thermal conductivities allow further evaluation of different two-component mixing models which are often applied in geothermal calculations (e.g., arithmetic mean, geometric mean, harmonic mean, Hashin-Shtrikman mean, and effective-medium theory mean). It is found that the geometric-mean model shows the best correlation between calculated and measured bulk thermal conductivity. However, by applying new model-dependent correction, equations the quality of fit could be significantly improved and the error diffusion of each model reduced. The ‘corrected’ geometric mean provides the most satisfying results and constitutes a universally applicable model for sedimentary rocks. Furthermore, lithotype-specific and model-independent conversion equations are developed permitting a calculation of water-saturated thermal conductivity from dry-measured thermal conductivity and porosity within an error range of 5 to 10%. The limited availability of core samples and the expensive core-based laboratory measurements make it worthwhile to use petrophysical well logs to determine thermal conductivity for sedimentary rocks. The approach followed in this study is based on the detailed analyses of the relationships between thermal conductivity of rock-forming minerals, which are most abundant in sedimentary rocks, and the properties measured by standard logging tools. By using multivariate statistics separately for clastic, carbonate and evaporite rocks, the findings from these analyses allow the development of prediction equations from large artificial data sets that predict matrix thermal conductivity within an error of 4 to 11%. These equations are validated successfully on a comprehensive subsurface data set from the NGB. In comparison to the application of earlier published approaches formation-dependent developed for certain areas, the new developed equations show a significant error reduction of up to 50%. These results are used to infer rock thermal conductivity for entire borehole profiles. By inversion of corrected in-situ thermal-conductivity profiles, temperature profiles are calculated and compared to measured high-precision temperature logs. The resulting uncertainty in temperature prediction averages < 5%, which reveals the excellent temperature prediction capabilities using the presented approach. In conclusion, data and methods are provided to achieve a much more detailed parameterization of thermal models.
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.
As part of the International Continental Drilling Program deep lake drilling project Paleo Van, we investigated sulfate reduction (SR) in deep sediment cores of the saline, alkaline (salinity 21.4%0, alkalinity 155 m mEq-1, pH 9.81) Lake Van, Turkey. The cores were retrieved in the Northern Basin (NB) and at Ahlat Ridge (AR) and reached a maximum depth of 220 m. Additionally, 65-75 cm long gravity cores were taken at both sites. SR rates (SRR) were low (<22 nmol cm-3 day-1) compared to lakes with higher salinity and alkalinity, indicating that salinity and alkalinity are not limiting SR in Lake Van. Both sites differ significantly in rates and depth distribution of SR. In NB, SRR are up to 10 times higher than at AR. SR could be detected down to 19 mblf (meters below lake floor) at NB and down to 13 mblf at AR. Although SRR were lower at AR than at NB, organic matter (OM) concentrations were higher. In contrast, dissolved OM in the pore water at AR contained more macromolecular OM and less low molecular weight OM.VVe thus suggest, that OM content alone cannot be used to infer microbial activity at Lake Van but that quality of OM has an important impact as well. These differences suggest that biogeochemical processes in lacustrine sediments are reacting very sensitively to small variations in geological, physical, or chemical parameters over relatively short distances.
Silicon (Si), although not considered essential, has beneficial effects on plant growth which are mostly associated with the ability to accumulate amorphous (phytogenic) Si, e.g., as phytoliths. Phytogenic Si is the most active Si pool in the soil-plant system because of its great surface-to-volume ratio, amorphous structure, and high water solubility. Despite the high abundance of Si in terrestrial biogeosystems and its importance, e.g., for the global C cycle, little is known about Si fluxes between soil and plants and Si pools used by plants. This study aims at elucidating the contribution of various soil Si pools to Si uptake by wheat. As pH affects dissolution of Si pools and Si uptake by plants, the effect of pH (4.5 and 7) was evaluated. Wheat was grown on Si-free pellets mixed with one of the following Si pools: quartz sand (crystalline), anorthite powder (crystalline), or silica gel (amorphous). Silicon content was measured in aboveground biomass, roots, and soil solution 4 times in intervals of 7 d. At pH 4.5, plants grew best on anorthite, but pH did not significantly affect Si-uptake rates. Total Si contents in plant biomass were significantly higher in the silica-gel treatment compared to all other treatments, with up to 26 mg g(-1) in aboveground biomass and up to 17 mg g(-1) in roots. Thus, Si uptake depends on the conversion of Si into plant-available silicic acid. This conversion occurs too slowly for crystalline Si phases, therefore Si uptake from treatments with quartz sand and anorthite did not differ from the control. For plants grown on silica gel, real Si-uptake rates were higher than the theoretical value calculated based on water transpiration. This implies that Si uptake by wheat is driven not only by passive water flux but also by active transporters, depending on Si concentration in the aqueous phase, thus on type of Si pool. These results show that Si uptake by plants as well as plant growth are significantly affected by the type of Si pool and factors controlling its solubility.
Automated seismic event location by travel-time stacking an application to mining induced seismicity
(2013)
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.
In many hydrological applications, ground-wave velocity measurements are increasingly used to map and monitor shallow soil water content. In this study, we propose an automated spectral velocity analysis method to determine the direct ground-wave (DGW) velocity from common midpoint (CMP) or multi-offset ground-penetrating radar (GPR) data. The method introduced in this paper is a variation of the well-known spectral velocity analysis for seismic and GPR reflection events where velocity spectra are computed using different coherency measures along hyperbolas following the normal moveout model. Here, the unnormalized cross-correlation is computed between waveforms across data gathers that are corrected with a linear moveout equation using a predefined range of velocities. Peaks in the resulting velocity spectra identify linear events in the GPR data gathers like DGW events and allow for estimating the corresponding velocities. In addition to obtaining a DGW velocity measurement, we propose a robust method to estimate the associated velocity uncertainties based on the width of the peak in the calculated velocity spectrum. Our proposed method is tested on synthetic data examples to evaluate the influence of subsurface velocity, surveying geometry and signal frequency on the accuracy of estimated ground-wave velocities. In addition, we investigate the influence of such velocity uncertainties on subsequent soil water content estimates using an established petrophysical relationship. Furthermore, we apply our approach to analyse field data, which were collected across a test site in Canada to monitor a wide range of seasonal soil moisture variations. A comparison between our spectral velocity estimates and results derived from manually picked ground-wave arrivals shows good agreement, which illustrates that our spectral velocity analysis is a feasible tool to analyse DGW arrivals in multi-offset GPR data gathers in an objective and more automated manner.
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.
The Late Permian Zechstein Group in northeastern Germany is characterized by shelf and slope carbonates that rimmed a basin extending from eastern England through the Netherlands and Germany to Poland. Conventional reservoirs are found in grainstones rimming islands created by pre-existing paleohighs and platform-rimming shoals that compose steep margins in the north and ramp deposits in the southern part. The slope and basin deposits are characterized by debris flows and organic-rich mudstones. Lagoonal and basinal evaporites formed the seal for these carbonate and underlying sandstone reservoirs. The objective of this investigation is to evaluate potential unconventional reservoirs in organic-rich, fine-grained and/or tight mudrocks in slope and basin as well as platform carbonates occurring in this stratigraphic interval. Therefore, a comprehensive study was conducted that included sedimentology, sequence stratigraphy, petrography, and geochemistry. Sequence stratigraphic correlations from shelf to basin are crucial in establishing a framework that allows correlation of potential productive facies in fine-grained, organic-rich basinal siliceous and calcareous mudstones or interfingering tight carbonates and siltstones, ranging from the lagoon, to slope to basin, which might be candidates for forming an unconventional reservoir. Most organic-rich shales worldwide are associated with eustatic transgressions. The basal Zechstein cycles, Z1 and Z2, contain organic-rich siliceous and calcareous mudstones and carbonates that form major transgressive deposits in the basin. Maturities range from over-mature (gas) in the basin to oil-generation on the slope with variable TOC contents. This sequence stratigraphic and sedimentologic evaluation of the transgressive facies in the Z1 and Z2 assesses the potential for shale-gas/oil and hybrid unconventional plays. Potential unconventional reservoirs might be explored in laminated organic-rich mudstones within the oil window along the northern and southern slopes of the basin. Although the Zechstein Z1 and Z2 cycles might have limited shale-gas potential because of low thickness and deep burial depth to be economic at this point, unconventional reservoir opportunities that include hybrid and shale-oil potential are possible in the study area.
Quantifying lithogenic inputs to the North Pacific Ocean using the long-lived thorium isotopes
(2013)
Dissolved Th-232 is added to the ocean though the partial dissolution of lithogenic materials such as aerosol dust in the same way as other lithogenically sourced and more biologically important trace metals such as Fe. Oceanic Th-230, on the other hand, is sourced primarily from the highly predictable decay of dissolved U-234. The rate at which dissolved Th-232 is released by mineral dissolution can be constrained by a Th removal rate derived from Th-230:U-234 disequilibria, assuming steady-state. Calculated fluxes of dissolved Th-232 can in turn be used to estimate fluxes of other lithogenically sourced dissolved metals as well as the original lithogenic supplies, such as aerosol dust deposition, given the concentration and fractional solubility of Th (or other metals) in the lithogenic material. This method is applied to 7 water column profiles from the Innovative North Pacific Experiment (INOPEX) cruise of 2009 and 2 sites from the subtropical North Pacific. The structure of shallow depth profiles suggests rapid scavenging at the surface and at least partial regeneration of dissolved Th-232 at 100-200 m depth. This rapid cycling could involve colloidal Th generated during mineral dissolution, which may not be subject to the same removal rates as the more truly dissolved Th-230. An additional deep source of Th-232 was revealed in deep waters, most likely dissolution of seafloor sediments, and offers a constraint on dissolved trace element supply due to boundary exchange. (C) 2013 Elsevier B.V. All rights reserved.
Large areas in the humid tropics are currently undergoing land-use change. The decrease of tropical rainforest, which is felled for land clearing and timber production, is countered by increasing areas of tree plantations and secondary forests. These changes are known to affect the regional water cycle as a result of plant-specific water demand and by influencing key soil properties which determine hydrological flow paths. One of these key properties sensitive to land-use change is the saturated hydraulic conductivity (Ks) as it governs vertical percolation of water within the soil profile. Low values of Ks in a certain soil depth can form an impeding layer and lead to perched water tables and the development of predominantly lateral flow paths such as overland flow. These processes can induce nutrient redistribution, erosion and soil degradation and thus affect ecosystem services and human livelihoods. Due to its sensitivity to land-use change, Ks is commonly used to assess the associated changes in hydrological flow paths. The objective of this dissertation was to assess the effect of land-use change on hydrological flow paths by analysing Ks as indicator variable. Sources of Ks variability, their implications for Ks monitoring and the relationship between Ks and near-surface hydrological flow paths in the context of land-use change were studied. The research area was located in central Panama, a country widely experiencing the abovementioned changes in land use. Ks is dependent on both static, soil-inherent properties such as particle size and clay mineralogy and dynamic, land use-dependent properties such as organic carbon content. By conducting a pair of studies with one of these influences held constant in each, the importance of static and dynamic properties for Ks was assessed. Applying a space-for-time approach to sample Ks under secondary forests of different age classes on comparable soils, a recovery of Ks from the former pasture use was shown to require more than eight years. The process was limited to the 0−6 cm sampling depth and showed large variability among replicates. A wavelet analysis of a Ks transect crossing different soil map units under comparable land cover, old-growth tropical rainforest, showed large small-scale variability, which was attributed to biotic influences, as well as a possible but non-conclusive influence of soil types. The two results highlight the importance of dynamic, land use-dependent influences on Ks. Monitoring studies can help to quantify land use-induced change of Ks, but there is a variety of sampling designs which differ in efficiency of estimating mean Ks. A comparative study of four designs and their suitability for Ks monitoring is used to give recommendations about designing a Ks monitoring scheme. Quantifying changes in spatial means of Ks for small catchments with a rotational stratified sampling design did not prove to be more efficient than Simple Random Sampling. The lack of large-scale spatial structure prevented benefits of stratification, and large small-scale variability resulting from local biotic processes and artificial effects of destructive sampling caused a lack of temporal consistency in the re-sampling of locations, which is part of the rotational design. The relationship between Ks and near-surface hydrological flow paths is of critical importance when assessing the consequences of land-use change in the humid tropics. The last part of this dissertation aimed at disclosing spatial relationships between Ks and overland flow as influenced by different land cover types. The effects of Ks on overland-flow generation were spatially variable, different between planar plots and incised flowlines and strongly influenced by land-cover characteristics. A simple comparison of Ks values and rainfall intensities was insufficient to describe the observed pattern of overland flow. Likewise, event flow in the stream was apparently not directly related to overland flow response patterns within the catchments. The study emphasises the importance of combining pedological, hydrological, meteorological and botanical measurements to comprehensively understand the land use-driven change in hydrological flow paths. In summary, Ks proved to be a suitable parameter for assessing the influence of land-use change on soils and hydrological processes. The results illustrated the importance of land cover and spatial variability of Ks for decisions on sampling designs and for interpreting overland-flow generation. As relationships between Ks and overland flow were shown to be complex and dependent on land cover, an interdisciplinary approach is required to comprehensively understand the effects of land-use change on soils and near-surface hydrological flow paths in the humid tropics.
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.
The potential of weather radar observations for hydrological and meteorological research and applications is undisputed, particularly with increasing world-wide radar coverage. However, several barriers impede the use of weather radar data. These barriers are of both scientific and technical nature. The former refers to inherent measurement errors and artefacts, the latter to aspects such as reading specific data formats, geo-referencing, visualisation. The radar processing library wradlib is intended to lower these barriers by providing a free and open source tool for the most important steps in processing weather radar data for hydro-meteorological and hydrological applications. Moreover, the community-based development approach of wradlib allows scientists to share their knowledge about efficient processing algorithms and to make this knowledge available to the weather radar community in a transparent, structured and well-documented way.
Siberian larch forests and the ion content of thaw lakes form a geochemically functional entity
(2013)
Siberian larch forests growing on shallow permafrost soils have not, until now, been considered to be controlling the abiotic and biotic characteristics of the vast number of thaw-lake ecosystems. Here we show, using four independent data sets (a modern data set from 201 lakes from the tundra to taiga, and three lake-core records), that lake-water geochemistry in Yakutia is highly correlated with vegetation. Alkalinity increases with catchment forest density. We postulate that in this arid area, higher evapotranspiration in larch forests compared with that in the tundra vegetation leads to local salt accumulation in soils. Solutes are transported to nearby thaw lakes during rain events and snow melt, but are not fully transported into rivers, because there is no continuous groundwater flow within permafrost soils. This implies that potentially large shifts in the chemical characteristics of aquatic ecosystems to known warming are absent because of the slow response of catchment forests to climate change.
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.
Silicate weathering - initiated by major mineralogical transformations at the base of ten meters of clay-rich saprolite generates the exceptionally low weathering flux found in streams draining the crystalline rocks of the mountainous and humid tropical Highlands of Sri Lanka. This conclusion is reached from a thorough investigation of the mineralogical, chemical, and Sr isotope compositions of samples within a regolith profile extending >10 m from surface soil through the weathering front in charnockite bedrock (a high-grade metamorphic rock), corestones formed at the weathering front, as well as from the chemical composition of the dissolved loads in nearby streams. Weatherable minerals and soluble elements are fully depleted at the top of the profile, showing that the system is supply-limited, such that weathering fluxes are controlled directly by the supply of fresh minerals. We determine the weathering rates using two independent means: (1) in situ-produced cosmogenic nuclides in surface soil and creek sediments in the close vicinity of the regolith combined with immobile element mass balance across the regolith and (2) river dissolved loads. Silicate weathering rates determined from both approaches range from 16 to 36 t km(-2) y(-1), corresponding to a weathering front advance rate of 6-14 mm ky(-1). These rates agree across the 10(1) to 10(4) - y time scales over which our rate metrics integrate, suggesting that the weathering system operates at steady state. Within error these rates are furthermore compatible with those obtained by modeling the advance rate of the weathering front from chemical gradients and mineral dissolution rates. The silicate weathering flux out of the weathering profile, measured on small creeks, amounts to 84% of the profile's export flux; the remaining 16% is contributed by non-silicate, atmospheric-derived input. The silicate weathering flux, as measured by dissolved loads in large catchments, amounts to ca. 50% of the total dissolved flux; the remainder being contributed by dust, rain, and weathering of local marble bands. Spheroidal weathering is the key processes of converting the fresh bedrock into saprolite at the weathering front. The mineralogical composition of weathering rinds shows that the sequence of mineral decomposition is: pyroxene; plagioclase; biotite; K-feldspar. Observable biotite alteration does not appear to initiate spheroidal weathering within corestones; therefore, we infer that other processes than biotite oxidation, like pyroxene oxidation, clay formation from pyroxene and plagioclase decomposition, the development of secondary porosity by plagioclase dissolution, or even microbiologic processes at depth enable the coupling between slow advance of the weathering front and slow erosion at the surface. The comparison to tectonically more active tropical landscapes lets us conclude that the combination of hard rock with tightly interlocked mineral grains and slow erosion in the absence of tectonically-induced landscape rejuvenation lead to these exceptionally low weathering rates. (C) 2013 Elsevier Ltd. All rights reserved.
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.
The evolution of most orogens typically records cogenetic shortening and extension. Pervasive normal faulting in an orogen, however, has been related to late syn- and post-collisional stages of mountain building with shortening focused along the peripheral sectors of the orogen. While extensional processes constitute an integral part of orogenic evolution, the spatiotemporal characteristics and the kinematic linkage of structures related to shortening and extension in the core regions of the orogen are often not well known. Related to the India-Eurasia collision, the Himalaya forms the southern margin of the Tibetan Plateau and constitutes the most prominent Cenozoic type example of a collisional orogen. While thrusting is presently observed along the foothills of the orogen, several generations of extensional structures have been detected in the internal, high-elevation regions, both oriented either parallel or perpendicular to the strike of the orogen. In the NW Indian Himalaya, earthquake focal mechanisms, seismites and ubiquitous normal faulting in Quaternary deposits, and regional GPS measurements reveal ongoing E-W extension. In contrast to other extensional structures observed in the Himalaya, this extension direction is neither parallel nor perpendicular to the NE-SW regional shortening direction. In this study, I took advantage of this obliquity between the trend of the orogen and structures related to E-W oriented extension in order to address the question of the driving forces of different extension directions. Thus, extension might be triggered triggered by processes within the Tibetan Plateau or originates from the curvature of the Himalayan orogen. In order to elaborate on this topic, I present new fault-kinematic data based on systematic measurements of approximately 2000 outcrop-scale brittle fault planes with displacements of up to several centimeters that cover a large area of the NW Indian Himalaya. This new data set together with field observations relevant for relative chronology allows me to distinguish six different deformation styles. One of the main results are that the overall strain pattern derived from this data reflects the regionally important contractional deformation pattern very well, but also reveals significant extensional deformation. In total, I was able to identify six deformation styles, most of which are temporally and spatially linked and represent protracted shortening, but also significant extensional directions. For example, this is the first data set where a succession of both, arc-normal and E-W extension have been documented in the Himalaya. My observations also furnish the basis for a detailed overview of the younger extensional deformation history in the NW Indian Himalaya. Field and remote-sensing based geomorphic analyses, and geochronologic 40Ar/39Ar data on synkinematic muscovites along normal faults help elucidate widespread E-W extension in the NW Indian Himalaya which must have started at approximately 14-16 Ma, if not earlier. In addition, I documented and mapped fault scarps in Quaternary sedimentary deposits using satellite imagery and field inspection. Furthermore, I made field observations of regional normal faults, compiled structures from geological maps and put them in a regional context. Finally, I documented seismites in lake sediments close to the currently most active normal fault in the study area in order to extend the (paleo) seismic record of this particular fault. Taken together, this data sets document that E-W extension is the dominant active deformation style in the internal parts of the orogen. In addition, the combined field, geomorphic and remote-sensing data sets prove that E-W extension occurs in a much more larger region toward the south and west than the seismicity data have suggested. In conclusion, the data presented here reveal the importance of extension in a region, which is still dominated by ongoing collision and shortening. The regional fault distribution and cross-cutting relationships suggest that extension parallel and perpendicular to the strike of the orogen are an integral part of the southward propagation of the active thrust front and the associated lateral growth of the Himalayan arc. In the light of a wide range of models proposed for extension in the Himalaya and the Tibetan plateau, I propose that E-W extension in the NW Indian Himalaya is transferred from the Tibetan Plateau due the inability of the Karakorum fault (KF) to adequately accommodate ongoing E-W extension on the Tibetan Plateau. Furthermore, in line with other observations from Tibet, the onset of E-W normal faulting in the NW Himalaya may also reflect the attainment of high topography in this region, which generated crustal stresses conducive to spatially extensive extension.
Slow fourier transform
(2013)
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.
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.
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.
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.
A general increase in precipitation has been observed in Germany in the last century, and potential changes in flood generation and intensity are now at the focus of interest. The aim of the paper is twofold: a) to project the future flood conditions in Germany accounting for various river regimes (from pluvial to nival-pluvial regimes) and under different climate scenarios (the high, A2, low, B1, and medium, A1B, emission scenarios) and b) to investigate sources of uncertainty generated by climate input data and regional climate models. Data of two dynamical Regional Climate Models (RCMs), REMO (REgional Model) and CCLM (Cosmo-Climate Local Model), and one statistical-empirical RCM, Wettreg (Wetterlagenbasierte Regionalisierungsmethode: weather-type based regionalization method), were applied to drive the eco-hydrological model SWIM (Soil and Water Integrated Model), which was previously validated for 15 gauges in Germany. At most of the gauges, the 95 and 99 percentiles of the simulated discharge using SWIM with observed climate data had a good agreement with the observed discharge for 1961-2000 (deviation within +/- 10 %). However, the simulated discharge had a bias when using RCM climate as input for the same period. Generalized Extreme Value (GEV) distributions were fitted to the annual maximum series of river runoff for each realization for the control and scenario periods, and the changes in flood generation over the whole simulation time were analyzed. The 50-year flood values estimated for two scenario periods (2021-2060, 2061-2100) were compared to the ones derived from the control period using the same climate models. The results driven by the statistical-empirical model show a declining trend in the flood level for most rivers, and under all climate scenarios. The simulations driven by dynamical models give various change directions depending on region, scenario and time period. The uncertainty in estimating high flows and, in particular, extreme floods remains high, due to differences in regional climate models, emission scenarios and multi-realizations generated by RCMs.
Sustainable land use in Mountain Regions under global change synthesis across scales and disciplines
(2013)
Mountain regions provide essential ecosystem goods and services (EGS) for both mountain dwellers and people living outside these areas. Global change endangers the capacity of mountain ecosystems to provide key services. The Mountland project focused on three case study regions in the Swiss Alps and aimed to propose land-use practices and alternative policy solutions to ensure the provision of key EGS under climate and land-use changes. We summarized and synthesized the results of the project and provide insights into the ecological, socioeconomic, and political processes relevant for analyzing global change impacts on a European mountain region. In Mountland, an integrative approach was applied, combining methods from economics and the political and natural sciences to analyze ecosystem functioning from a holistic human-environment system perspective. In general, surveys, experiments, and model results revealed that climate and socioeconomic changes are likely to increase the vulnerability of the EGS analyzed. We regard the following key characteristics of coupled human-environment systems as central to our case study areas in mountain regions: thresholds, heterogeneity, trade-offs, and feedback. Our results suggest that the institutional framework should be strengthened in a way that better addresses these characteristics, allowing for (1) more integrative approaches, (2) a more network-oriented management and steering of political processes that integrate local stakeholders, and (3) enhanced capacity building to decrease the identified vulnerability as central elements in the policy process. Further, to maintain and support the future provision of EGS in mountain regions, policy making should also focus on project-oriented, cross-sectoral policies and spatial planning as a coordination instrument for land use in general.
High applications of P fertilizers and manure are general practice in intensive agriculture and may cause eutrophication in adjacent streams. Bioavailability of P can be estimated by sequential extractions commonly used for soil or sediment. A single combined method may facilitate more effective comparisons of topsoils and adjoining stream sediments, and enhance management decisions. In this study, the suitability of an established soil P sequential extraction was tested on stream bed sediments. The study was conducted in the Sumas River watershed in the agricultural Lower Fraser Valley, Canada. Sediment samples with differing land use (forest, low and high intensity agriculture) from 1993, 1994, 2008, and 2009 from 14 sites along the Sumas River and tributaries were used. Total sequential extraction concentrations were in agreement with aqua regia digestion (Rs=0.96) and showed consistency over the study time sequence. P fractions released by 0.5 M NaHCO3 (median 14 %), 0.1 M NaOH (33 %), and 1.0 M HCl (38 %) were significantly (alpha=0.05) higher than P released by other extractants. These three extraction steps provide a practical and time-effective assessment of P lability in stream sediments and may be used as a combined scheme for sediment and soil. Analytical results further revealed that land use has a major and characteristic impact on P lability. With a land use change from forest to intensive agriculture, results showed an increase in total P concentrations (30 to 4,000 ppm) and in P lability, in particular for the moderately labile NaOH-P fraction (20 to 50 %).
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.
Soil moisture estimation under low vegetation cover using a multi-angular polarimetric decomposition
(2013)
The estimation of volumetric soil moisture under low agricultural vegetation from fully polarimetric synthetic aperture radar (SAR) data at L-band using a multi-angular polarimetric decomposition is investigated. Radar polarimetry provides the framework to decompose the backscattered signal into different canonical scattering mechanisms referring to scattering contributions from the underlying soil and the vegetation cover. Multiangular observation diversity further increases the information space for soil moisture inversion enabling higher inversion rates and a stable inversion performance. The developed approach was applied on the multi-angular L-band data set acquired by German Aerospace Center's ESAR sensor as part of the OPAQUE campaign in 2008. The obtained results are compared against ground measurements collected by the OPAQUE team over a variety of vegetated agricultural fields. The validation of the estimated against ground measured soil moisture results in an root mean square error level of 6-8 vol.% including all test fields with a variety of crop types.
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.
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.
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).
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.
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.
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.
The internal geological structure of the Northeast German Basin (NEGB) is affected by intense salt diapirism and by the presence of several stratified aquifer complexes of regional relevance. The shallow Quaternary to late Tertiary freshwater aquifer is separated from the underlying Mesozoic saline aquifers by an embedded Tertiary clay enriched aquitard (Rupelian Aquitard). An important feature of this aquitard is that hydraulic connections between the upper and lower aquifers do exist in areas where the Rupelian Aquitard is missing (hydrogeological windows). Three-dimensional thermohaline numerical simulations are carried out to investigate the effects of such hydrogeological windows in the Rupelian Aquitard on the resulting groundwater, temperature, and salinity distributions. Numerical results suggest that hydrogeological windows act as preferential domains of hydraulic interconnectivity between the different aquifers at depth and enable vigorous heat and mass transport which causes a mixing of warm and saline groundwater with cold and less saline groundwater within both aquifers. In areas where the Rupelian Aquitard confines the Mesozoic aquifer, dissolved solutes from major salt structures are transported laterally giving rise to plumes of variable salinity content ranging from few hundreds of meters to several tens of kilometers. Furthermore, destabilizing thermal buoyancy forces may overwhelm counteracting stabilizing salinity induced forces offside of salt domes. This may result in buoyant upward groundwater flow transporting heat and mass to shallower levels within the same Mesozoic Aquifer.
The hypersaline crater lake and its catchment on seabird island Isabel (Pacific, off Mexico) was studied to explore the influence of strong seasonal variations in rainfall/evaporation and guano contribution on its limnology. The hypersaline lake water (HSW, 78 %) is up to 2.2-times enriched in inert ions relative to mean seawater. Rainfall during summer dilutes the HSW to form a less saline rainwater body (RWB) above a chemolimnion between 2 and 4 m water depth. The RWB is inhabited first by diatoms and ostracods followed later on by cyanobacteria and ciliates. Evaporation of > 1.5 m depth of lake water over the dry season increases the salinity of the RWB until the water column becomes isohaline at HSW concentrations in the late dry season. Differences in the stable isotope composition of water and primary producers in RWB and HSW reflect this development. Introduction of seabird guano and the decrease of salinity fuel a high primary production in the RWB with higher delta(CDIC)-C-13 and delta(13)Corg of particulate organic matter than in the HSW. The high N supply leads to high delta N-15 NH4 values (+ 39 % in the HSW) as the consequence of ammonia volatilization that is strongest during guano maturation and with evaporative salinity increase from the HSW. Precipitation of carbonate (calcite and aragonite) from the RWB and the HSW is hindered by the high concentration of guano-derived P. This inhibition may be overcome with evaporative supersaturation during particularly dry conditions. Carbonate may also precipitate during particularly wet conditions from the dilute RWB, where the P-concentration is reduced during an active phytoplankton production that raises the pH. Differences in the stable isotope signatures of carbon and oxygen in HSW and RWB (+ 5 % delta(CDIC)-C-13 and -3 % d18OH2O) suggest the processes of carbonate precipitation can be distinguished based on the isotope signature of the carbonates deposited. Changes in the lake system are indicated when lower temperatures and higher rainfall in the 2006 wet season introduced more and less mature guano to the lake. The lower pH was accompanied by lower ammonia volatilization and carbonate precipitation as indicated by an increased concentration of NH4, Ca, Sr and DIC, while delta H-2, delta(NNH4)-N-15, and salinity were lower. According to our results, the observed sediment laminations should reflect the introduction of catchment material (including guano) with runoff, the RWB plankton production, and the carbonate precipitation in relation to its origin and seasonality.
The Armutlu peninsula, located in the eastern Marmara Sea, coincides with the western end of the rupture of the 17 August 1999, Izmit M-W 7.6 earthquake which is the penultimate event of an apparently westward migrating series of strong and disastrous earthquakes along the NAFZ during the past century. We present new seismotectonic data of this key region in order to evaluate previous seismotectonic models and their implications for seismic hazard assessment in the eastern Marmara Sea. Long term kinematics were investigated by performing paleo strain reconstruction from geological field investigations by morphotectonic and kinematic analysis of exposed brittle faults. Short term kinematics were investigated by inverting for the moment tensor of 13 small to moderate recent earthquakes using surface wave amplitude spectra. Our results confirm previous models interpreting the eastern Marmara Sea Region as an active transtensional pull-apart environment associated with significant NNE-SSW extension and vertical displacement. At the northern peninsula, long term deformation pattern did not change significantly since Pliocene times contradicting regional tectonic models which postulate a newly formed single dextral strike slip fault in the Marmara Sea Region. This area is interpreted as a horsetail splay fault structure associated with a major normal fault segment that we call the Waterfall Fault. Apart from the Waterfall Fault, the stress strain relation appears complex associated with a complicated internal fault geometry, strain partitioning, and reactivation of pre-existing plane structures. At the southern peninsula, recent deformation indicates active pull-apart tectonics constituted by NE-SW trending dextral strike slip faults. Earthquakes generated by stress release along large rupture zones seem to be less probable at the northern, but more probable at the southern peninsula. Additionally, regional seismicity appears predominantly driven by plate boundary stresses as transtensional faulting is consistent with the southwest directed far field deformation of the Anatolian plate. (C) 2013 Elsevier B.V. All rights reserved.
The electricity system is particularly susceptible to climate change due to the close interconnectedness between electricity production, consumption and climate. This study provides a country based relative analysis of 21 European countries' electricity system susceptibility to climate change. Taking into account 14 quantitative influencing factors, the susceptibility of each country is examined both for the current and projected system with the result being a relative ranked index. Luxembourg and Greece are the most susceptible relatively due in part to their inability to meet their own electricity consumption demand with inland production, and the fact that the majority of their production is from more susceptible sources, primarily combustible fuels. Greece experiences relatively warm mean temperatures, which are expected to increase in the future leading to greater summer electricity consumption, increasing susceptibility. Norway was found to be the least susceptible, relatively, due to its consistent production surplus, which is primarily from hydro (a less susceptible source) and a likely decrease of winter electricity consumption as temperatures rise due to climate change. The findings of this study enable countries to identify the main factors that increase their electricity system susceptibility and proceed with adaptation measures that are the most effective in decreasing susceptibility.
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.
Cratons with their thick lithospheric roots can influence the thermal structure, and thus the convective flow, in the surrounding mantle. As mantle temperatures are hard to measure directly, depth variations in the mantle transition zone (MTZ) discontinuities are often employed as a proxy. Here, we use a large new data set of P-receiver functions to map the 410 km and 660 km discontinuities beneath the western edge of the East European Craton and adjacent Phanerozoic Europe across the most fundamental lithospheric boundary in Europe, the Trans-European Suture Zone (TESZ). We observe significantly shorter travel times for conversions from both MTZ discontinuities within the craton, caused by the high velocities of the cratonic root. By contrast, the differential travel time across the MTZ is normal to only slightly raised. This implies that any insulating effect of the cratonic keel does not reach the MTZ. In contrast to earlier observations in Siberia, we do not find any trace of a discontinuity at 520 km depth, which indicates a rather dry MTZ beneath the western edge of the craton. Within most of covered Phanerozoic Europe, the MTZ differential travel time is remarkably uniform and in agreement with standard Earth models. No widespread thermal effects of the various episodes of Caledonian and Variscan subduction that took place during the amalgamation of the continent remain. Only more recent tectonic events, related to Alpine subduction and Quarternary volcanism in the Eifel area, can be traced. While the East European craton shows no distinct imprint into the MTZ, we discover the signature of the TESZ in the MTZ in the form of a linear region of about 350 km width with a 1.5 s increase in differential travel time, which could either be caused by high water content or decreased temperature. Taking into account results of recent S-wave tomographies, raised water content in the MTZ cannot be the main cause for this observation. Accordingly, we explain the increase, equivalent to a 15 km thicker MTZ, by a temperature decrease of about 80 K. We discuss two alternative models for this temperature reduction, either a remnant of subduction or an indication of downwelling due to small-scale, edge-driven convection caused by the contrast in lithospheric thickness across the TESZ. Any subducted lithosphere found in the MTZ at this location is unlikely to be related to Variscan subduction along the TESZ, though, as Eurasia has moved significantly northward since the Variscan orogeny.
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.
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.
AimAdvancement in ecological methods predicting species distributions is a crucial precondition for deriving sound management actions. Maximum entropy (MaxEnt) models are a popular tool to predict species distributions, as they are considered able to cope well with sparse, irregularly sampled data and minor location errors. Although a fundamental assumption of MaxEnt is that the entire area of interest has been systematically sampled, in practice, MaxEnt models are usually built from occurrence records that are spatially biased towards better-surveyed areas. Two common, yet not compared, strategies to cope with uneven sampling effort are spatial filtering of occurrence data and background manipulation using environmental data with the same spatial bias as occurrence data. We tested these strategies using simulated data and a recently collated dataset on Malay civet Viverra tangalunga in Borneo.
LocationBorneo, Southeast Asia.
MethodsWe collated 504 occurrence records of Malay civets from Borneo of which 291 records were from 2001 to 2011 and used them in the MaxEnt analysis (baseline scenario) together with 25 environmental input variables. We simulated datasets for two virtual species (similar to a range-restricted highland and a lowland species) using the same number of records for model building. As occurrence records were biased towards north-eastern Borneo, we investigated the efficacy of spatial filtering versus background manipulation to reduce overprediction or underprediction in specific areas.
ResultsSpatial filtering minimized omission errors (false negatives) and commission errors (false positives). We recommend that when sample size is insufficient to allow spatial filtering, manipulation of the background dataset is preferable to not correcting for sampling bias, although predictions were comparatively weak and commission errors increased.
Main ConclusionsWe conclude that a substantial improvement in the quality of model predictions can be achieved if uneven sampling effort is taken into account, thereby improving the efficacy of species conservation planning.
Biogeochemical turnover in hyporheic zones is known to have the potential to affect the chemical signature of surface water cycling through shallow streambed sediments. This study investigates the impact of streambed physical properties on the fate of nitrate and dissolved oxygen in groundwater upwelling through the streambed of a lowland river. For analyzing depth-dependent patterns and zonation of nitrogen concentrations, diffuse gel probes in shallow (top 15 cm) streambed sediments have been deployed in a nested setup together with multilevel minipiezometers for streambed sediments of 15-150 cm. Spatial heterogeneity of groundwater upwelling was controlled by patterns of low-conductivity peat and clay strata that caused locally confined conditions, suggesting increased streambed residence times. Nitrate concentrations in the upwelling groundwater changed by up to 68.06 mg L-1 within the top 15 cm of streambed sediments and by up to 107.47 mg L-1 at depths of 15-150 cm, indicating that significant nitrogen turnover was not restricted to shallow streambed sediments. Intensive reduction of nitrate concentrations was found, in particular, in vicinity of low-conductivity streambed strata. The coincidence of confined groundwater upwelling and reduced oxygen concentrations at these locations suggests that increased residence times and associated depletion of dissolved oxygen create conditions favorable for nitrate reduction. Our results highlight that increased nitrogen turnover at aquifer-river interfaces is not necessarily limited to shallow streambed zones, where surface water is mixing with groundwater, but can affect upwelling groundwater in reactive hot spots that extend to greater streambed depths and beyond hyporheic mixing zones. Citation: Krause, S., C. Tecklenburg, M. Munz, and E. Naden (2013), Streambed nitrogen cycling beyond the hyporheic zone: Flow controls on horizontal patterns and depth distribution of nitrate and dissolved oxygen in the upwelling groundwater of a lowland river,
The impact of land use changes on sustainable development is of increasing interest in many regions of the world. This study aimed to test the transferability of the Framework for Participatory Impact Assessment (FoPIA), which was originally developed in the European context, to developing countries, in which lack of data often prevents the use of data-driven impact assessment methods. The core aspect of FoPIA is the stakeholder-based assessment of alternative land use scenarios. Scenario impacts on regional sustainability are assessed by using a set of nine regional land use functions (LUFs), which equally cover the economic, social and environmental dimensions of sustainability. The cases analysed in this study include (1) the alternative spatial planning policies around the Merapi volcano and surrounding areas of Yogyakarta City, Indonesia; (2) the large-scale afforestation of agricultural areas to reduce soil erosion in Guyuan, China; (3) the expansion of soil and water conservation measures in the Oum Zessar watershed, Tunisia; (4) the agricultural intensification and the potential for organic agriculture in Bijapur, India; and (5) the land degradation and land conflicts resulting from land division and privatisation in Narok, Kenya. All five regions are characterised by population growth, partially combined with considerable economic development, environmental degradation problems and social conflicts. Implications of the regional scenario impacts as well as methodological aspects are discussed. Overall, FoPIA proved to be a useful tool for diagnosing regional human-environment interactions and for supporting the communication and social learning process among different stakeholder groups.
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.