Refine
Has Fulltext
- no (140)
Year of publication
- 2013 (140) (remove)
Document Type
- Article (140) (remove)
Language
- English (140)
Is part of the Bibliography
- yes (140)
Keywords
- Climate change (4)
- Land use change (4)
- Iran (3)
- exhumation (3)
- Alborz Mountains (2)
- Asian monsoon (2)
- Earthquake source observations (2)
- Geostatistics (2)
- Holocene (2)
- Mountain basins (2)
Institute
- Institut für Geowissenschaften (140) (remove)
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.