Refine
Year of publication
- 2015 (256) (remove)
Document Type
- Article (203)
- Doctoral Thesis (33)
- Postprint (4)
- Review (4)
- Habilitation Thesis (3)
- Other (3)
- Monograph/Edited Volume (2)
- Preprint (2)
- Bachelor Thesis (1)
- Conference Proceeding (1)
Is part of the Bibliography
- yes (256)
Keywords
- Holocene (6)
- climate change (4)
- erosion (4)
- India-Asia collision (3)
- Ostracoda (3)
- Pamir (3)
- Seismicity and tectonics (3)
- ecohydrology (3)
- landscape evolution (3)
- Cerrado (2)
Institute
- Institut für Geowissenschaften (256) (remove)
Previous studies suggested that electric and/or magnetic field fluctuations observed in the nighttime topside ionosphere at midlatitudes generally originate from quiet time nocturnal medium-scale traveling ionospheric disturbances (MSTIDs). However, decisive evidences for the connection between the two have been missing. In this study we make use of the multispacecraft observations of midlatitude magnetic fluctuations (MMFs) in the nighttime topside ionosphere by the Swarm constellation. The analysis results show that the area hosting MMFs is elongated in the NW-SE (NE-SW) direction in the Northern (Southern) Hemisphere. The elongation direction and the magnetic field polarization support that the area hosting MMFs is nearly field aligned. All these properties of MMFs suggest that they have close relationship with MSTIDs. Expectation values of root-mean-square field-aligned currents associated with MMFs are up to about 4nA/m(2). MMF coherency significantly drops for longitudinal distances of 1 degrees.
Permafrost inundated since the last glacial maximum is degrading, potentially releasing trapped or stabilized greenhouse gases, but few observations of the depth of ice-bonded permafrost (IBP) below the seafloor exist for most of the arctic continental shelf. We use spectral ratios of the ambient vibration seismic wavefield, together with estimated shear wave velocity from the dispersion curves of surface waves, for estimating the thickness of the sediment overlying the IBP. Peaks in spectral ratios modeled for three-layered 1-D systems correspond with varying thickness of the unfrozen sediment. Seismic receivers were deployed on the seabed around Muostakh Island in the central Laptev Sea, Siberia. We derive depths of the IBP between 3.7 and 20.7m15%, increasing with distance from the shoreline. Correspondence between expected permafrost distribution, modeled response, and observational data suggests that the method is promising for the determination of the thickness of unfrozen sediment.
Earthquake clustering has proven the most useful tool to forecast changes in seismicity rates in the short and medium term (hours to months), and efforts are currently being made to extend the scope of such models to operational earthquake forecasting. The overarching goal of the research presented in this thesis is to improve physics-based earthquake forecasts, with a focus on aftershock sequences. Physical models of triggered seismicity are based on the redistribution of stresses in the crust, coupled with the rate-and-state constitutive law proposed by Dieterich to calculate changes in seismicity rate. This type of models are known as Coulomb- rate and-state (CRS) models. In spite of the success of the Coulomb hypothesis, CRS models typically performed poorly in comparison to statistical ones, and they have been underepresented in the operational forecasting context. In this thesis, I address some of these issues, and in particular these questions: (1) How can we realistically model the uncertainties and heterogeneity of the mainshock stress field? (2) What is the effect of time dependent stresses in the postseismic phase on seismicity? I focus on two case studies from different tectonic settings: the Mw 9.0 Tohoku megathrust and the Mw 6.0 Parkfield strike slip earthquake. I study aleatoric uncertainties using a Monte Carlo method. I find that the existence of multiple receiver faults is the most important source of intrinsic stress heterogeneity, and CRS models perform better when this variability is taken into account. Epistemic uncertainties inherited from the slip models also have a significant impact on the forecast, and I find that an ensemble model based on several slip distributions outperforms most individual models. I address the role of postseismic stresses due to aseismic slip on the mainshock fault (afterslip) and to the redistribution of stresses by previous aftershocks (secondary triggering). I find that modeling secondary triggering improves model performance. The effect of afterslip is less clear, and difficult to assess for near-fault aftershocks due to the large uncertainties of the afterslip models. Off-fault events, on the other hand, are less sensitive to the details of the slip distribution: I find that following the Tohoku earthquake, afterslip promotes seismicity in the Fukushima region. To evaluate the performance of the improved CRS models in a pseudo-operational context, I submitted them for independent testing to a collaborative experiment carried out by CSEP for the 2010-2012 Canterbury sequence. Preliminary results indicate that physical models generally perform well compared to statistical ones, suggesting that CRS models may have a role to play in the future of operational forecasting. To facilitate efforts in this direction, and to enable future studies of earthquake triggering by time dependent processes, I have made the code open source. In the final part of this thesis I summarize the capabilities of the program and outline technical aspects regarding performance and parallelization strategies.
Information on extreme precipitation for future climate is needed to assess the changes in the frequency and intensity of flooding. The primary source of information in climate change impact studies is climate model projections. However, due to the coarse resolution and biases of these models, they cannot be directly used in hydrological models. Hence, statistical downscaling is necessary to address climate change impacts at the catchment scale.
This study compares eight statistical downscaling methods (SDMs) often used in climate change impact studies. Four methods are based on change factors (CFs), three are bias correction (BC) methods, and one is a perfect prognosis method. The eight methods are used to downscale precipitation output from 15 regional climate models (RCMs) from the ENSEMBLES project for 11 catchments in Europe. The overall results point to an increase in extreme precipitation in most catchments in both winter and summer. For individual catchments, the downscaled time series tend to agree on the direction of the change but differ in the magnitude. Differences between the SDMs vary between the catchments and depend on the season analysed. Similarly, general conclusions cannot be drawn regarding the differences between CFs and BC methods. The performance of the BC methods during the control period also depends on the catchment, but in most cases they represent an improvement compared to RCM outputs. Analysis of the variance in the ensemble of RCMs and SDMs indicates that at least 30% and up to approximately half of the total variance is derived from the SDMs. This study illustrates the large variability in the expected changes in extreme precipitation and highlights the need for considering an ensemble of both SDMs and climate models. Recommendations are provided for the selection of the most suitable SDMs to include in the analysis.
In a recent BAMS article, it is argued that community-based Open Source Software (OSS) could foster scientific progress in weather radar research, and make weather radar software more affordable, flexible, transparent, sustainable, and interoperable.
Nevertheless, it can be challenging for potential developers and users to realize these benefits: tools are often cumbersome to install; different operating systems may have particular issues, or may not be supported at all; and many tools have steep learning curves.
To overcome some of these barriers, we present an open, community-based virtual machine (VM). This VM can be run on any operating system, and guarantees reproducibility of results across platforms. It contains a suite of independent OSS weather radar tools (BALTRAD, Py-ART, wradlib, RSL, and Radx), and a scientific Python stack. Furthermore, it features a suite of recipes that work out of the box and provide guidance on how to use the different OSS tools alone and together. The code to build the VM from source is hosted on GitHub, which allows the VM to grow with its community.
We argue that the VM presents another step toward Open (Weather Radar) Science. It can be used as a quick way to get started, for teaching, or for benchmarking and combining different tools. It can foster the idea of reproducible research in scientific publishing. Being scalable and extendable, it might even allow for real-time data processing.
We expect the VM to catalyze progress toward interoperability, and to lower the barrier for new users and developers, thus extending the weather radar community and user base.
The Central Pontides is an accretionary-type orogenic area within the Alpine-Himalayan orogenic belt characterized by pre-collisional tectonic continental growth. The region comprises Mesozoic subduction-accretionary complexes and an accreted intra-oceanic arc that are sandwiched between the Laurasian active continental margin and Gondwana-derived the Kırşehir Block. The subduction-accretion complexes mainly consist of an Albian-Turonian accretionary wedge representing the Laurasian active continental margin. To the north, the wedge consists of slate/phyllite and metasandstone intercalation with recrystallized limestone, Na-amphibole-bearing metabasite (PT= 7–12 kbar and 400 ± 70 ºC) and tectonic slices of serpentinite representing accreted distal part of a large Lower Cretaceous submarine turbidite fan deposited on the Laurasian active continental margin that was subsequently accreted and metamorphosed. Raman spectra of carbonaceous material (RSCM) of the metapelitic rocks revealed that the metaflysch sequence consists of metamorphic packets with distinct peak metamorphic temperatures. The majority of the metapelites are low-temperature (ca. 330 °C) slates characterized by lack of differentiation of the graphite (G) and D2 defect bands. They possibly represent offscraped distal turbidites along the toe of the Albian accretionary wedge. The rest are phyllites that are characterized by slightly pronounced G band with D2 defect band occurring on its shoulder. Peak metamorphic temperatures of these phyllites are constrained to 370-385 °C. The phyllites are associated with a strip of incipient blueschist facies metabasites which are found as slivers within the offscraped distal turbidites. They possibly represent underplated continental metasediments together with oceanic crustal basalt along the basal décollement. Tectonic emplacement of the underplated rocks into the offscraped distal turbidites was possibly achieved by out-of-sequence thrusting causing tectonic thickening and uplift of the wedge. 40Ar/39Ar phengite ages from the phyllites are ca. 100 Ma, indicating Albian subduction and regional HP metamorphism.
The accreted continental metasediments are underlain by HP/LT metamorphic rocks of oceanic origin along an extensional shear zone. The oceanic metamorphic sequence mainly comprises tectonically thickened deep-seated eclogite to blueschist facies metabasites and micaschists. In the studied area, metabasites are epidote-blueschists locally with garnet (PT= 17 ± 1 kbar and 500 ± 40 °C). Lawsonite-blueschists are exposed as blocks along the extensional shear zone (PT= 14 ± 2 kbar and 370–440 °C). They are possibly associated with low shear stress regime of the initial stage of convergence. Close to the shear zone, the footwall micaschists consist of quartz, phengite, paragonite, chlorite, rutile with syn-kinematic albite porphyroblast formed by pervasive shearing during exhumation. These types of micaschists are tourmaline-bearing and their retrograde nature suggests high-fluid flux along shear zones. Peak metamorphic mineral assemblages are partly preserved in the chloritoid-micaschist farther away from the shear zone representing the zero strain domains during exhumation. Three peak metamorphic assemblages are identified and their PT conditions are constrained by pseudosections produced by Theriak-Domino and by Raman spectra of carbonaceous material: 1) garnet-chloritoid-glaucophane with lawsonite pseudomorphs (P= 17.5 ± 1 kbar, T: 390-450 °C) 2) chloritoid with glaucophane pseudomorphs (P= 16-18 kbar, T: 475 ± 40 °C) and 3) relatively high-Mg chloritoid (17%) with jadeite pseudomorphs (P= 22-25 kbar; T: 440 ± 30 °C) in addition to phengite, paragonite, quartz, chlorite, rutile and apatite. The last mineral assemblage is interpreted as transformation of the chloritoid + glaucophane assemblage to chloritoid + jadeite paragenesis with increasing pressure. Absence of tourmaline suggests that the chloritoid-micaschist did not interact with B-rich fluids during zero strain exhumation. 40Ar/39Ar phengite age of a pervasively sheared footwall micaschist is constrained to 100.6 ± 1.3 Ma and that of a chloritoid-micaschist is constrained to 91.8 ± 1.8 Ma suggesting exhumation during on-going subduction with a southward younging of the basal accretion and the regional metamorphism. To the south, accretionary wedge consists of blueschist and greenschist facies metabasite, marble and volcanogenic metasediment intercalation. 40Ar/39Ar phengite dating reveals that this part of the wedge is of Middle Jurassic age partly overprinted during the Albian. Emplacement of the Middle Jurassic subduction-accretion complexes is possibly associated with obliquity of the Albian convergence.
Peak metamorphic assemblages and PT estimates of the deep-seated oceanic metamorphic sequence suggest tectonic stacking within wedge with different depths of burial. Coupling and exhumation of the distinct metamorphic slices are controlled by decompression of the wedge possibly along a retreating slab. Structurally, decompression of the wedge is evident by an extensional shear zone and the footwall micaschists with syn-kinematic albite porphyroblasts. Post-kinematic garnets with increasing grossular content and pseudomorphing minerals within the chloritoid-micaschists also support decompression model without an extra heating.
Thickening of subduction-accretionary complexes is attributed to i) significant amount of clastic sediment supply from the overriding continental domain and ii) deep level basal underplating by propagation of the décollement along a retreating slab. Underplating by basal décollement propagation and subsequent exhumation of the deep-seated subduction-accretion complexes are connected and controlled by slab rollback creating a necessary space for progressive basal accretion along the plate interface and extension of the wedge above for exhumation of the tectonically thickened metamorphic sequences. This might be the most common mechanism of the tectonic thickening and subsequent exhumation of deep-seated HP/LT subduction-accretion complexes.
To the south, the Albian-Turonian accretionary wedge structurally overlies a low-grade volcanic arc sequence consisting of low-grade metavolcanic rocks and overlying metasedimentary succession is exposed north of the İzmir-Ankara-Erzincan suture (İAES), separating Laurasia from Gondwana-derived terranes. The metavolcanic rocks mainly consist of basaltic andesite/andesite and mafic cognate xenolith-bearing rhyolite with their pyroclastic equivalents, which are interbedded with recrystallized pelagic limestone and chert. The metavolcanic rocks are stratigraphically overlain by recrystallized micritic limestone with rare volcanogenic metaclastic rocks. Two groups can be identified based on trace and rare earth element characteristics. The first group consists of basaltic andesite/andesite (BA1) and rhyolite with abundant cognate gabbroic xenoliths. It is characterized by relative enrichment of LREE with respect to HREE. The rocks are enriched in fluid mobile LILE, and strongly depleted in Ti and P reflecting fractionation of Fe-Ti oxides and apatite, which are found in the mafic cognate xenoliths. Abundant cognate gabbroic xenoliths and identical trace and rare earth elements compositions suggest that rhyolites and basaltic andesites/andesites (BA1) are cogenetic and felsic rocks were derived from a common mafic parental magma by fractional crystallization and accumulation processes. The second group consists only of basaltic andesites (BA2) with flat REE pattern resembling island arc tholeiites. Although enriched in LILE, this group is not depleted in Ti or P.
Geochemistry of the metavolcanic rocks indicates supra-subduction volcanism evidenced by depletion of HFSE and enrichment of LILE. The arc sequence is sandwiched between an Albian-Turonian subduction-accretionary complex representing the Laurasian active margin and an ophiolitic mélange. Absence of continent derived detritus in the arc sequence and its tectonic setting in a wide Cretaceous accretionary complex suggest that the Kösdağ Arc was intra-oceanic. This is in accordance with basaltic andesites (BA2) with island arc tholeiite REE pattern.
Zircons from two metarhyolite samples give Late Cretaceous (93.8 ± 1.9 and 94.4 ± 1.9 Ma) U/Pb ages. Low-grade regional metamorphism of the intra-oceanic arc sequence is constrained 69.9 ± 0.4 Ma by 40Ar/39Ar dating on metamorphic muscovite from a metarhyolite indicating that the arc sequence became part of a wide Tethyan Cretaceous accretionary complex by the latest Cretaceous. The youngest 40Ar/39Ar phengite age from the overlying subduction-accretion complexes is 92 Ma confirming southward younging of an accretionary-type orogenic belt. Hence, the arc sequence represents an intra-oceanic paleo-arc that formed above the sinking Tethyan slab and finally accreted to Laurasian active continental margin. Abrupt non-collisional termination of arc volcanism was possibly associated with southward migration of the arc volcanism similar to the Izu-Bonin-Mariana arc system.
The intra-oceanic Kösdağ Arc is coeval with the obducted supra-subduction ophiolites in NW Turkey suggesting that it represents part of the presumed but missing incipient intra-oceanic arc associated with the generation of the regional supra-subduction ophiolites. Remnants of a Late Cretaceous intra-oceanic paleo-arc and supra-subduction ophiolites can be traced eastward within the Alp-Himalayan orogenic belt. This reveals that Late Cretaceous intra-oceanic subduction occurred as connected event above the sinking Tethyan slab. It resulted as arc accretion to Laurasian active margin and supra-subduction ophiolite obduction on Gondwana-derived terranes.
Development of geophysical methods to characterize methane hydrate reservoirs on a laboratory scale
(2015)
Gas hydrates are crystalline solids composed of water and gas molecules. They are stable at elevated pressure and low temperatures. Therefore, natural gas hydrate deposits occur at continental margins, permafrost areas, deep lakes, and deep inland seas. During hydrate formation, the water molecules rearrange to form cavities which host gas molecules. Due to the high pressure during hydrate formation, significant amounts of gas can be stored in hydrate structures. The water-gas ratio hereby can reach up to 1:172 at 0°C and atmospheric pressure. Natural gas hydrates predominantly contain methane. Because methane constitutes both a fuel and a greenhouse gas, gas hydrates are a potential energy resource as well as a potential source for greenhouse gas.
This study investigates the physical properties of methane hydrate bearing sediments on a laboratory scale. To do so, an electrical resistivity tomography (ERT) array was developed and mounted in a large reservoir simulator (LARS). For the first time, the ERT array was applied to hydrate saturated sediment samples under controlled temperature, pressure, and hydrate saturation conditions on a laboratory scale. Typically, the pore space of (marine) sediments is filled with electrically well conductive brine. Because hydrates constitute an electrical isolator, significant contrasts regarding the electrical properties of the pore space emerge during hydrate formation and dissociation. Frequent measurements during hydrate formation experiments permit the recordings of the spatial resistivity distribution inside LARS. Those data sets are used as input for a new data processing routine which transfers the spatial resistivity distribution into the spatial distribution of hydrate saturation. Thus, the changes of local hydrate saturation can be monitored with respect to space and time.
This study shows that the developed tomography yielded good data quality and resolved even small amounts of hydrate saturation inside the sediment sample. The conversion algorithm transforming the spatial resistivity distribution into local hydrate saturation values yielded the best results using the Archie-var-phi relation. This approach considers the increasing hydrate phase as part of the sediment frame, metaphorically reducing the sample’s porosity. In addition, the tomographical measurements showed that fast lab based hydrate formation processes cause small crystallites to form which tend to recrystallize.
Furthermore, hydrate dissociation experiments via depressurization were conducted in order to mimic the 2007/2008 Mallik field trial. It was observed that some patterns in gas and water flow could be reproduced, even though some setup related limitations arose.
In two additional long-term experiments the feasibility and performance of CO2-CH4 hydrate exchange reactions were studied in LARS. The tomographical system was used to monitor the spatial hydrate distribution during the hydrate formation stage. During the subsequent CO2 injection, the tomographical array allowed to follow the CO2 migration front inside the sediment sample and helped to identify the CO2 breakthrough.
Methanhydrate sind besonders in Verbindung mit den steigenden Weltmarktpreisen für Öl und Gas in den vergangenen Jahren mehr und mehr in den Fokus der Energiewirtschaft geraten, was zu einer starken Zunahme der angewandten Forschungsprojekte auf diesem Gebiet führte. Da Methanhydrat nur unter hohem Druck und niedrigen Temperaturen stabil ist, ist die Gewinnung natürlicher Proben für Laboruntersuchungen technisch sehr aufwendig und vor allem teuer. Zur Charakterisierung der Eigenschaften hydratführender Reservoire ist man häufig auf die Herstellung synthetischer Proben angewiesen. Die Eigenschaften der synthetisierten Proben sind dabei abhängig von der Herstellungsmethode und man ist noch immer auf der Suche nach Verfahren, mit denen sich möglichst „naturnahe“ Proben mit vertretbarem Aufwand erzeugen lassen.
In der vorliegenden Arbeit wurde eine neue, relativ schnell durchführbare Methode getestet, die im Porenraum von Sedimenten schwimmende bzw. gefügestützende Hydrate bildet, wie sie in der Natur vorkommen. Gleichzeitig erzeugt sie eine gleichmäßige Verteilung des Hydrats über die Probe und bietet gute Kontrolle über den Hydratgehalt. Sie funktioniert wie folgt: Eine mit einer KCl-Lösung gesättigte Sedimentprobe wird zu einem bestimmten Teil ausgefroren und das übrige Wasser mit Methan verdrängt. Durch Anlegen eines Methandrucks im Stabilitätsbereich wird das Eis zu Methanhydrat umgesetzt. Im Anschluss wird die Probe erneut mit einer KCl-Lösung gesättigt. Anhand seismischer Messungen konnte bestätigt werden, dass Hydrat mit dem gewünschten Hydrathabitus erzeugt wurde. Des Weiteren wurde gezeigt, dass die eishaltigen Proben aufgrund ähnlicher physikalischer Eigenschaften bereits vor der Umsetzung des Eises zu Methanhydrat als Näherung für Proben mit Porenraumhydrat verwendet werden können.
This thesis contains three experimental studies addressing the interplay between deformation and the mineral reaction between natural calcite and magnesite. The solid-solid mineral reaction between the two carbonates causes the formation of a magnesio-calcite precursor layer and a dolomite reaction rim in every experiment at isostatic annealing and deformation conditions.
CHAPTER 1 briefly introduces general aspects concerning mineral reactions in nature and diffusion pathways for mass transport. Moreover, results of previous laboratory studies on the influence of deformation on mineral reactions are summarized. In addition, the main goals of this study are pointed out.
In CHAPTER 2, the reaction between calcite and magnesite single crystals is examined at isostatic annealing conditions. Time series performed at a fixed temperature revealed a diffusion-controlled dolomite rim growth. Two microstructural domains could be identified characterized by palisade-shaped dolomite grains growing into the magnesite and granular dolomite growing towards calcite. A model was provided for the dolomite rim growth based on the counter-diffusion of CaO and MgO. All reaction products exhibited a characteristic crystallographic relationship with respect to the calcite reactant. Moreover, kinetic parameters of the mineral reaction were determined out of a temperature series at a fixed time. The main goal of the isostatic test series was to gain information about the microstructure evolution, kinetic parameters, chemical composition and texture development of the reaction products. The results were used as a reference to quantify the influence of deformation on the mineral reaction.
CHAPTER 3 deals with the influence of non-isostatic deformation on dolomite and magnesio-calcite layer production between calcite and magnesite single crystals. Deformation was achieved by triaxial compression and by torsion. Triaxial compression up to 38 MPa axial stress at a fixed time showed no significant influence of stress and strain on dolomite formation. Time series conducted at a fixed stress yield no change in growth rates for dolomite and magnesio-calcite at low strains. Slightly larger magnesio-calcite growth rates were observed at strains above >0.1. High strains at similar stresses were caused by the activation of additional glide systems in the calcite single crystal and more mobile dislocations in the magnesio-calcite grains, providing fast diffusion pathways. In torsion experiments a gradual decrease in dolomite and magnesio-calcite layer thickness was observed at a critical shear strain. During deformation, crystallographic orientations of reaction products rearranged with respect to the external framework. A direct effect of the mineral reaction on deformation could not be recognized due to the relatively small reaction product widths.
In CHAPTER 4, the influence of starting material microfabrics and the presence of water on the reaction kinetics was evaluated. In these experimental series polycrystalline material was in contact with single crystals or two polycrystalline materials were used as reactants. Isostatic annealing resulted in different dolomite and magnesio-calcite layer thicknesses, depending on starting material microfabrics. The reaction progress at the magnesite interface was faster with smaller magnesite grain size, because grain boundaries provided fast pathways for diffusion and multiple nucleation sites for dolomite formation. Deformation by triaxial compression and torsion yield lower dolomite rim thicknesses compared to annealed samples for the same time. This was caused by grain coarsening of polycrystalline magnesite during deformation. In contrast, magnesio-calcite layers tended to be larger during deformation, which triggered enhanced diffusion along grain boundaries. The presence of excess water had no significant influence on the reaction kinetics, at least if the reactants were single crystals.
In CHAPTER 5 general conclusions about the interplay between deformation and the mineral reaction in the carbonate system are presented.
Finally, CHAPTER 6 highlights possible future work in the carbonate system based on the results of this study.
Scientific theories of how subduction and plate tectonics began on Earth-and what the tectonic structure of Earth was before this-remain enigmatic and contentious(1). Understanding viable scenarios for the onset of subduction and plate tectonics(2,3) is hampered by the fact that subduction initiation processes must have been markedly different before the onset of global plate tectonics because most present-day subduction initiation mechanisms require acting plate forces and existing zones of lithospheric weakness, which are both consequences of plate tectonics(4). However, plume-induced subduction initiation(5-9) could have started the first subduction zone without the help of plate tectonics. Here, we test this mechanism using high-resolution three-dimensional numerical thermomechanical modelling. We demonstrate that three key physical factors combine to trigger self-sustained subduction: (1) a strong, negatively buoyant oceanic lithosphere; (2) focused magmatic weakening and thinning of lithosphere above the plume; and (3) lubrication of the slab interface by hydrated crust. We also show that plume-induced subduction could only have been feasible in the hotter early Earth for old oceanic plates. In contrast, younger plates favoured episodic lithospheric drips rather than self-sustained subduction and global plate tectonics.
We present an experimental approach to study the three-dimensional microstructure of gas diffusion layer (GDL) materials under realistic compression conditions. A dedicated compression device was designed that allows for synchrotron-tomographic investigation of circular samples under well-defined compression conditions. The tomographic data provide the experimental basis for stochastic modeling of nonwoven GDL materials. A plain compression tool is used to study the fiber courses in the material at different compression stages. Transport relevant geometrical parameters, such as porosity, pore size, and tortuosity distributions, are exemplarily evaluated for a GDL sample in the uncompressed state and for a compression of 30 vol.%. To mimic the geometry of the flow-field, we employed a compression punch with an integrated channel-rib-profile. It turned out that the GDL material is homogeneously compressed under the ribs, however, much less compressed underneath the channel. GDL fibers extend far into the channel volume where they might interfere with the convective gas transport and the removal of liquid water from the cell. (C) 2015 AIP Publishing LLC.
Effect of mass wasting on soil organic carbon storage and coastal erosion in permafrost environments
(2015)
Accelerated permafrost thaw under the warming Arctic climate can have a significant impact on Arctic landscapes. Areas underlain by permafrost store high amounts of soil organic carbon (SOC). Permafrost disturbances may contribute to increased release of carbon dioxide and methane to the atmosphere. Coastal erosion, amplified through a decrease in Arctic sea-ice extent, may also mobilise SOC from permafrost. Large expanses of permafrost affected land are characterised by intense mass-wasting processes such as solifluction, active-layer detachments and retrogressive thaw slumping. Our aim is to assess the influence of mass wasting on SOC storage and coastal erosion.
We studied SOC storage on Herschel Island by analysing active-layer and permafrost samples, and compared non-disturbed sites to those characterised by mass wasting. Mass-wasting sites showed decreased SOC storage and material compaction, whereas sites characterised by material accumulation showed increased storage. The SOC storage on Herschel Island is also significantly correlated to catenary position and other slope characteristics. We estimated SOC storage on Herschel Island to be 34.8 kg C m-2. This is comparable to similar environments in northwest Canada and Alaska.
Coastal erosion was analysed using high resolution digital elevation models (DEMs). Two LIDAR scanning of the Yukon Coast were done in 2012 and 2013. Two DEMs with 1 m horizontal resolution were generated and used to analyse elevation changes along the coast. The results indicate considerable spatial variability in short-term coastline erosion and progradation. The high variability was related to the presence of mass-wasting processes. Erosion and deposition extremes were recorded where the retrogressive thaw slump (RTS) activity was most pronounced. Released sediment can be transported by longshore drift and affects not only the coastal processes in situ but also along adjacent coasts.
We also calculated volumetric coastal erosion for Herschel Island by comparing a stereo-photogrammetrically derived DEM from 2004 with LIDAR DEMs. We compared this volumetric erosion to planimetric erosion, which was based on coastlines digitised from satellite imagery. We found a complex relationship between planimetric and volumetric coastal erosion, which we attribute to frequent occurrence of mass-wasting processes along the coasts. Our results suggest that volumetric erosion corresponds better with environmental forcing and is more suitable for the estimation of organic carbon fluxes than planimetric erosion.
Mass wasting can decrease SOC storage by several mechanisms. Increased aeration following disturbance may increase microbial activity, which accelerates organic matter decomposition. New hydrological conditions that follow the mass wasting event can cause leaching of freshly exposed material. Organic rich material can also be directly removed into the sea or into a lake. On the other hand the accumulation of mobilised material can result in increased SOC storage. Mass-wasting related accumulations of mobilised material can significantly impact coastal erosion in situ or along the adjacent coast by longshore drift. Therefore, the coastline movement observations cannot completely resolve the actual sediment loss due to these temporary accumulations. The predicted increase of mass-wasting activity in the course of Arctic warming may increase SOC mobilisation and coastal erosion induced carbon fluxes.
The long-term temporal and spatial changes in statistical, source, and stress characteristics of one cluster of induced seismicity recorded at The Geysers geothermal field (U.S.) are analyzed in relation to the field operations, fluid migration, and constraints on the maximum likely magnitude. Two injection wells, Prati-9 and Prati-29, located in the northwestern part of the field and their associated seismicity composed of 1776 events recorded throughout a 7year period were analyzed. The seismicity catalog was relocated, and the source characteristics including focal mechanisms and static source parameters were refined using first-motion polarity, spectral fitting, and mesh spectral ratio analysis techniques. The source characteristics together with statistical parameters (b value) and cluster dynamics were used to investigate and understand the details of fluid migration scheme in the vicinity of injection wells. The observed temporal, spatial, and source characteristics were clearly attributed to fluid injection and fluid migration toward greater depths, involving increasing pore pressure in the reservoir. The seasonal changes of injection rates were found to directly impact the shape and spatial extent of the seismic cloud. A tendency of larger seismic events to occur closer to injection wells and a correlation between the spatial extent of the seismic cloud and source sizes of the largest events was observed suggesting geometrical constraints on the maximum likely magnitude and its correlation to the average injection rate and volume of fluids present in the reservoir.
Barite concretions and bands are widely distributed in black shale-chert horizons in the Yurtus Formation of Lower Cambrian in Aksu area, northwestern Tarim Basin, NW China. They mainly consist of coarse-grained anhedral to euhedral barite crystals with minor dolomites and pyrites. Petrological features indicate these concretions grew from the porewater in unconsolidated sediments at shallow burial below sediment-water interface. The slight deviation of Sr-87/Sr-86 ratios (0.7083 to 0.7090) and significant elevated delta S-34 values (56.8-76.4 parts per thousand CDT) of barite samples with respect to those of the Early Cambrian seawater further support that barite deposits precipitated from the enclosed porewater in sediment column, which evolved from the penecontemporaneous seawater with weak interaction with the host fine-grained siliciclastic sediments and highly-depleted sulfate in response to prolonged strong bacterial sulfate reduction without necessary renewal. The abundant organic matters in the basal Yurtus Formation should have facilitated developing sulfate-depleted methanogenesis zone and sulfate-methane transition zone (SMTZ) slightly after deposition. Therefore, barite deposits in the Yurtus Formation most likely resulted from diagenetic barium cycling and persistently grew from the porewater in the static SMTZ with a low sedimentation rate in the Early Cambrian. In comparison with the distribution of sedimentary barites in geological records, we tentatively proposed that a transition in diagenetic barium cycling and associated mineralization may have occurred from the Precambrian to Cambrian periods; this scenario may be causally linked to the changes in marine ecology (the advent of mesozooplankton and associated faecal pellet) and geochemistry (the increase of seawater sulfate concentration). Thus, the occurrence of diagenetic barite deposits in the Yurtus Formation implies that diagenetic barium cycling and more effective scavenging of barium from CH4- and Ba-rich porewaters within sediments might have become an nonnegligible process in continental margin areas, at least, since the earliest Cambrian, which could have significantly impacted the marine barium cycling. (C) 2015 Elsevier B.V. All rights reserved.
Limits to lichenometry
(2015)
Lichenometry is a straightforward and inexpensive method for dating Holocene rock surfaces. The rationale is that the diameter of the largest lichen scales with the age of the originally fresh rock surface that it colonised. The success of the method depends on finding the largest lichen diameters, a suitable lichen-growth model, and a robust calibration curve. Recent critique of the method motivates us to revisit the accuracy and uncertainties of lichenometry. Specifically, we test how well lichenometry is capable of resolving the ages of different lobes of large active rock glaciers in the Kyrgyz Tien Shan. We use a bootstrapped quantile regression to calibrate local growth curves of Xanthoria elegans, Aspicilia tianshanica, and Rhizocarpon geographicum, and report a nonlinear decrease in dating accuracy with increasing lichen diameter. A Bayesian type of an analysis of variance demonstrates that our calibration allows discriminating credibly between rock-glacier lobes of different ages despite the uncertainties tied to sample size and correctly identifying the largest lichen thalli. Our results also show that calibration error grows with lichen size, so that the separability of rock-glacier lobes of different ages decreases, while the tendency to assign coeval ages increases. The abundant young (<200 yr) specimen of fast-growing X elegans are in contrast with the fewer, slow-growing, but older (200-1500 yr) R. geographicum and A. tianshanica, and record either a regional reactivation of lobes in the past 200 years, or simply a censoring effect of lichen mortality during early phases of colonisation. The high variance of lichen sizes captures the activity of rock-glacier lobes, which is difficult to explain by regional climatic cooling or earthquake triggers alone. Therefore, we caution against inferring palaeoclimatic conditions from the topographic position of rock-glacier lobes. We conclude that lichenometry works better as a tool for establishing a relative, rather than an absolute, chronology of rock-glacier lobes in the northern Tien Shan. (C) 2015 Elsevier Ltd. All rights reserved.
The 2008 eruption of Chaiten volcano in southern Chile severely impacted several densely forested river catchments by supplying excess pyroclastic sediment to the channel networks. Our aim is to substantiate whether and how channel geometry and forest stands changed in the Rayas River following the sudden input of pyroclastic sediment. We measured the resulting changes to channel geometry and riparian forest stands along 17.6 km of the impacted gravel-bed Rayas River (294 km(2)) from multiple high-resolution satellite images, aerial photographs, and fieldwork to quantify yield volume characteristics of the forest stands. Limited channel changes during the last 60 years before the eruption reflect a dynamic equilibrium condition of the river corridor, despite the high annual precipitation and the sediment supply from Chaiten and Michinmahuida volcanoes in the headwaters. Images taken in 1945, 2004, and 2005 show that total size of the vegetated channel islands nearly doubled between 1945 and 2004 and remained unchanged between 2004 and 2005. Pyroclastic sediment entering the Rayas River after the 2008 eruption caused only minor average channel widening (7%), but killed all island vegetation in the study reach. Substantial shifts in the size distribution of in-channel vegetation patches reflect losses in total island area of 46% from 2005 to 2009 and an additional 34% from 2009 to 2012. The estimated pulsed release of organic carbon into the channel, mainly in the form of large wood from obliterated island and floodplain forests, was 78-400 tC/km/y and surpasses most documented yields from small mountainous catchments with similar rainfall, forest cover, and disturbance history, while making up between 20% and 60% of the annual carbon burial rate of fluvial sediments in the northern Patagonian fjords. We conclude that the carbon footprint of the 2008 Chaiten eruption on the Rayas River was more significant than the measured geomorphic impacts on channel geometry for the first five years following disturbance. The modest post-eruptive geomorphic response in this river is a poor indicator of its biogeochemical response. (C) 2015 Elsevier B.V. All rights reserved.
Mass wasting is an important process for denuding hillslopes and lowering ridge crests in active mountain belts such as the Himalaya-Karakoram ranges (HKR). Such a high-relief landscape is likely to be at its mechanical threshold, maintained by competing rapid rock uplift, river incision, and pervasive slope failure. We introduce excess topography, Z(E), for quantifying potentially unstable rock-mass volumes inclined at angles greater than a specified threshold angle. We find that Z(E) peaks along major fluvial and glacial inner gorges, which is also where the majority of 492 large (>0.1 km(2)) rock-slope failures occur in the Himalaya's largest cluster of documented Pleistocene to Holocene bedrock landslides. Our data reveal that bedrock landslides in the HKR chiefly detached from near or below the median elevation, whereas glaciers and rock glaciers occupy higher-elevation bands almost exclusively. Less than 10% of the area of the HKR is upslope of glaciers, such that possible censoring of evidence of large bedrock landslides above the permanent snow line barely affects this finding. Bedrock landslides appear to preferentially undermine topographic relief in response to fluvial and glacial incision along inner gorges, unless more frequent and smaller undetected failures, or rigorous (peri-)glacial erosion, compensate for this role at higher elevation. Either way, the distinct patterns of excess topography and large bedrock landsliding in the HKR juxtapose two stacked domains of landslide and (peri-)glacial erosion that may respond to different time scales of perturbation. Our findings call for more detailed analysis of vertical erosional domains and their geomorphic coupling in active mountain belts.
In this study, we investigate how immersive 3D geovisualization can be used in higher education. Based on MacEachren and Kraak's geovisualization cube, we examine the usage of immersive 3D geovisualization and its usefulness in a research-based learning module on flood risk, called GEOSimulator. Results of a survey among participating students reveal benefits, such as better orientation in the study area, higher interactivity with the data, improved discourse among students and enhanced motivation through immersive 3D geovisualization. This suggests that immersive 3D visualization can effectively be used in higher education and that 3D CAVE settings enhance interactive learning between students.
Roads at risk
(2015)
Globalisation and interregional exchange of people, goods, and services has boosted the importance of and reliance on all kinds of transport networks. The linear structure of road networks is especially sensitive to natural hazards. In southern Norway, steep topography and extreme weather events promote frequent traffic disruption caused by debris flows. Topographic susceptibility and trigger frequency maps serve as input into a hazard appraisal at the scale of first-order catchments to quantify the impact of debris flows on the road network in terms of a failure likelihood of each link connecting two network vertices, e.g. road junctions. We compute total additional traffic loads as a function of traffic volume and excess distance, i.e. the extra length of an alternative path connecting two previously disrupted network vertices using a shortest-path algorithm. Our risk metric of link failure is the total additional annual traffic load, expressed as vehicle kilometres, because of debris-flow-related road closures. We present two scenarios demonstrating the impact of debris flows on the road network and quantify the associated path-failure likelihood between major cities in southern Norway. The scenarios indicate that major routes crossing the central and north-western part of the study area are associated with high link-failure risk. Yet options for detours on major routes are manifold and incur only little additional costs provided that drivers are sufficiently well informed about road closures. Our risk estimates may be of importance to road network managers and transport companies relying on speedy delivery of services and goods.
The use of topographic metrics for estimating the susceptibility to, and reconstructing the characteristics of, debris flows has a long research tradition, although largely devoted to humid mountainous terrain. The exceptional 2010 monsoonal rainstorms in the high-altitude mountain desert of Ladakh and Zanskar, NW India, were a painful reminder of how susceptible arid regions are to rainfall-triggered flash floods, landslides, and debris flows. The rainstorms of August 4-6 triggered numerous debris flows, killing 182 people, devastating 607 houses, and more than 10 bridges around Ladakh's capital of Leh. The lessons from this disaster motivated us to revisit methods of predicting (a) flow parameters such as peak discharge and maximum velocity from field and remote sensing data, and (b) the susceptibility to debris flows from catchment morphometry. We focus on quantifying uncertainties tied to these approaches. Comparison of high-resolution satellite images pre- and post-dating the 2010 rainstorm reveals the extent of damage and catastrophic channel widening. Computations based on these geomorphic markers indicate maximum flow velocities of 1.6-6.7 m s(-1) with runout of up to similar to 10 km on several alluvial fans that sustain most of the region's settlements. We estimate median peak discharges of 310-610 m(3) s(-1), which are largely consistent with previous estimates. Monte Carlo-based error propagation for a single given flow-reconstruction method returns a variance in discharge similar to one derived from juxtaposing several different flow reconstruction methods. We further compare discriminant analysis, classification tree modelling, and Bayesian logistic regression to predict debris-flow susceptibility from morphometric variables of 171 catchments in the Ladakh Range. These methods distinguish between fluvial and debris flow-prone catchments at similar success rates, but Bayesian logistic regression allows quantifying uncertainties and relationships between potential predictors. We conclude that, in order to be robust and reliable, morphometric reconstruction of debris-flow properties and susceptibility requires careful assessment and reporting of errors and uncertainties. (C) 2015 Elsevier B.V. All rights reserved.
Geochemistry and Ar-40/Ar-39 age of Early Carboniferous dolerite sills in the southern Baltic Sea
(2015)
The Early Carboniferous magmatic event in the southern Baltic Sea is manifested by dolerite intrusions. The presumable area in which the dolerite intrusions occur ranges from 30 to 60 km in east-west direction, and is about 100 km in north-south direction. The dolerites were sampled in well D1-1 and investigated by applying chemical analysis and Ar-40/Ar-39 step-heating dating. Dolerites are classified as alkali and sodic, characterized by high TiO2 (3.92, 3.99 wt%) and P2O5 (1.67, 1.77 wt%) and low MgO (4.89, 4.91 wt%) concentrations, enriched in light rare earth elements, originated from an enriched mantle magma source and emplaced in a continental rift tectonic setting. The 351 +/- 11 Ma Ar-40/Ar-39 plateau age for groundmass plagioclase indicates a considerable age gap with the 310-250 Ma magmatism in southern Scandinavia and northern Germany. The magmatic rocks in the Baltic Sedimentary Basin are coeval with alkaline intrusions of NE Poland. Both magmatic provinces lie in the northwestward prolongation of the Pripyat-Dnieper-Donetsk Rift (370-359 Ma) and may constitute a later phase of magmatic activity of this propagating rift system.
Damage due to floods has increased during the last few decades, and further increases are expected in several regions due to climate change and growing vulnerability. To address the projected increase in flood risk, a combination of structural and non-structural flood risk mitigation measures is considered as a promising adaptation strategy. Such a combination takes into account that flood defence systems may fail, and prepares for unexpected crisis situations via land-use planning and private damage reduction, e.g. via building precautionary measures, and disaster response. However, knowledge about damage-reducing measures is scarce and often fragmented since based on case studies. For instance, it is believed that private precautionary measures, like shielding with water shutters or building fortification, are especially effective in areas with frequent flood events and low flood water levels. However, some of these measures showed a significant damage-reducing effect also during the extreme flood event in 2002 in Germany. This review analyses potentials of land-use planning and private flood precautionary measures as components of adaptation strategies for global change. Focus is on their implementation, their damage-reducing effects and their potential contribution to address projected changes in flood risk, particularly in developed countries.
One of the most critical problems affecting geophysical data acquisition procedures is related to the misorientation of multicomponent sensors with respect to a common reference system (e.g., geographic north). In many applications, misoriented sensors affect data analysis procedures, leading to errors in results and interpretations. These problems generally occur in applications where the orientation of the sensor cannot be actively controlled and is not known a priori, e.g., geophysical sensors deployed in borehole installations or on the seafloor. We have developed a quaternion-based method for the optimal reorientation of multicomponent geophysical sensors. In contrast to other approaches, we took into account the full time-series record from all sensor components. Therefore, our method could be applied to all time-series data and was not restricted to a certain type of geophysical sensor. Our method allows the robust calculation of relative reorientations between two-component or three-component sensors. By using a reference sensor in an iterative process, this result can be extended to the estimation of absolute sensor orientations. In addition to finding an optimal solution for a full 3D sensor rotation, we have established a rigorous scheme for the estimation of uncertainties of the resulting orientation parameters. We tested the feasibility and applicability of our method using synthetic data examples for a vertical seismic profile and an ocean bottom seismometer array. We noted that the quaternion-based reorientation method is superior to the standard approach of a single-parameter estimation of rotation angles.
Sea-surface temperatures (SSTs) based on radiolarian assemblage changes are estimated for the last 160 kyr, from a sediment core (Y9) recovered from Pukaki Saddle, northeast of Campbell Plateau. Site Y9 lies beneath Subantarctic Surface Water (SAW) immediately to the north of the Subantarctic Front (SAF), which in this region is bathymetrically constrained by the edges of Campbell Plateau and defines the northern boundary of the Antarctic Circumpolar Current (ACC). Radiolarian assemblages are characterised by an exceptionally high abundance of the Antarctic to subantarctic species Antarctissa spp. (up to 68%), especially during glacial intervals. SST estimates are derived using Factor Analysis and the Modern Analog Technique. Both methods capture the glacial-interglacial (G-I) pattern. The SST reconstructions show the changing relative influence of distinct water masses during the past G-I cycle, with major temperature variations of the order of 7-9 degrees C at glacial Terminations. Glacials (marine isotope stages (MIS) 6 and 2) are associated with particularly cool SSTs that are indicative of a more vigorous SAF/ACC and an enhancement of the inflow through Pukaki Saddle and/or frequent development of cold-core eddies at the SAF. By contrast, the influence of warmer waters and relaxation of the ACC during interglacials can be inferred from temperatures slightly warmer (e.g., mid-Holocene) and/or comparable to present day (e.g., MIS 5e). During these intervals, relatively warmer temperatures most likely indicate a higher warmcore eddy activity due to a strengthened Subtropical Front and/or a weakened inflow of cool water through Pukaki Saddle and/or an increased stratification in the Campbell Plateau region. Furthermore, the SST record is characterised by an abrupt warming at ca. 10 kyr (i.e., Termination l), the occurrence of a reversal at Termination I, and a warming event at the end of MIS 4, coinciding with the A4 event in the Byrd ice core. These characteristics, together with the pronounced G-I cycle shown by the SST estimates, suggest that Site Y9 is influenced by major oceanographic changes in the SW Pacific and responds to thermal changes at high southern latitudes. (C) 2014 Elsevier B.V. All rights reserved.
High resolution isotopic (delta O-18 and delta C-13) investigations on endogenic carbonates (calcite/aragonite) from Tso Moriri Lake, NW Himalaya show dramatic fluctuations during the late glacial and the early Holocene, and a persistent enrichment trend during the late Holocene. Changes in this lake are largely governed by the [input (meltwater + monsoon precipitation)/evaporationj (WE) ratio, also reflected in changes in the carbonate mineralogy with aragonite being formed during periods of lowest I/E. Using new isotopic data on endogenic carbonates in combination with the available data on geochemistry, mineralogy, and reconstructed mean annual precipitation, we demonstrate that the late glacial and early Holocene carbonate delta O-18 variability resulted from fluctuating Indian summer monsoon (ISM) precipitation in NW Himalaya. This region experienced increasing ISM precipitation between ca. 13.1 and 11.7 cal ka and highest ISM precipitation during the early Holocene (11.2-8.5 cal ka). However, during the late Holocene, evaporation was the dominant control on the carbonate delta O-18. Regional comparison of reconstructed hydrological changes from Tso Moriri Lake with other archives from the Asian summer monsoon and westerlies domain shows that the intensified westerly influence that resulted in higher lake levels (after 8 cal ka) in central Asia was not strongly felt in NW Himalaya. (C) 2015 Elsevier B.V. All rights reserved.
Combinations of runoff characteristics are commonly used to represent distinct conceptual models of stormflow generation. In this study, three runoff characteristics: hydrograph response, time source of runoff water, and flow path are used to classify catchments. Published data from the scientific literature are used to provide evidence from small, forested catchments. Each catchment was assigned to one of the eight conceptual models, depending on the combination of quick/slow response, old/new water, and overland/subsurface flow. A standard procedure was developed to objectively diagnose the predominant conceptual model of stormflow generation for each catchment and assess its temporal and spatial support. The literature survey yielded 42 catchments, of which 30 catchments provide a complete set of qualitative runoff characteristics resulting in one of the eight conceptual models. The majority of these catchments classify as subsurface flow path dominated. No catchments were found for conceptual models representing combinations of quick response-new water-subsurface flow (SSF), slow-new-SSF, slow-old-overland flow (OF) nor new-slow-OF. Of the 30 qualitatively classified catchments, 24 provide a complete set of quantitative measures. In summary, the field support is strong for 19 subsurface-dominated catchments and is weak for 5 surface flow path dominated catchments (six catchments had insufficient quantitative data). Two alternative explanations exist for the imbalance of field support between the two flow path classes: (1) the selection of research catchments in past field studies was mainly to explain quick hydrograph response in subsurface dominated catchments; (2) catchments with prevailing subsurface flow paths are more common in nature. We conclude that the selection of research catchments needs to cover a wider variety of environmental conditions which should lead to a broader, and more widely applicable, spectrum of resulting conceptual models and process mechanisms. This is a prerequisite in studies where catchment organization and similarity approaches are used to develop catchment classification systems in order to regionalize stormflow.
Remnants of hydrous melt formed at mantle depths have been identified and characterized within high-pressure leucogranulites of the Orlica-Snieznik Dome (Bohemian Massif, central Europe). They occur as nanogranites in garnet formed via partial melting of granitoids during the Variscan orogeny. Melt composition and H2O content have been investigated in situ after experimental re-homogenization of the nanogranites, and are consistent with melts produced experimentally from crustal lithologies at mantle depths. This is the first geochemical study of melt inclusions from natural crustal rocks equilibrated close to the stability field of coesite, shedding light on how continental crust melts during deep subduction. Whereas decompressional melting is commonly invoked for deeply subducted crustal lithologies, melting occurred near or at the metamorphic peak pressure in the Orlica-Snieznik granulites. Melting of deeply subducted crustal rocks significantly modifies the rheology and thus promotes fast exhumation: this process has a critical influence on the geodynamic evolution of subduction-collision orogens as well as crustal differentiation at depth.
Hyperspectral REE (Rare Earth Element) Mapping of Outcrops-Applications for Neodymium Detection
(2015)
In this study, an in situ application for identifying neodymium (Nd) enriched surface materials that uses multitemporal hyperspectral images is presented (HySpex sensor). Because of the narrow shape and shallow absorption depth of the neodymium absorption feature, a method was developed for enhancing and extracting the necessary information for neodymium from image spectra, even under illumination conditions that are not optimal. For this purpose, the two following approaches were developed: (1) reducing noise and analyzing changing illumination conditions by averaging multitemporal image scenes and (2) enhancing the depth of the desired absorption band by deconvolving every image spectrum with a Gaussian curve while the rest of the spectrum remains unchanged (Richardson-Lucy deconvolution). To evaluate these findings, nine field samples from the Fen complex in Norway were analyzed using handheld X-ray fluorescence devices and by conducting detailed laboratory-based geochemical rare earth element determinations. The result is a qualitative outcrop map that highlights zones that are enriched in neodymium. To reduce the influences of non-optimal illumination, particularly at the studied site, a minimum of seven single acquisitions is required. Sharpening the neodymium absorption band allows for robust mapping, even at the outer zones of enrichment. From the geochemical investigations, we found that iron oxides decrease the applicability of the method. However, iron-related absorption bands can be used as secondary indicators for sulfidic ore zones that are mainly enriched with rare earth elements. In summary, we found that hyperspectral spectroscopy is a noninvasive, fast and cost-saving method for determining neodymium at outcrop surfaces.
We study segregation of the subducted oceanic crust (OC) at the core-mantle boundary and its ability to accumulate and form large thermochemical piles (such as the seismically observed Large Low Shear Velocity Provinces (LLSVPs)). Our high-resolution numerical simulations of thermochemical mantle convection suggest that the longevity of LLSVPs for up to three billion years, and possibly longer, can be ensured by a balance in the rate of segregation of high-density OC material to the core-mantle boundary (CMB) and the rate of its entrainment away from the CMB by mantle upwellings. For a range of parameters tested in this study, a large-scale compositional anomaly forms at the CMB, similar in shape and size to the LLSVPs. Neutrally buoyant thermochemical piles formed by mechanical stirringwhere thermally induced negative density anomaly is balanced by the presence of a fraction of dense anomalous materialbest resemble the geometry of LLSVPs. Such neutrally buoyant piles tend to emerge and survive for at least 3Gyr in simulations with quite different parameters. We conclude that for a plausible range of values of density anomaly of OC material in the lower mantleit is likely that it segregates to the CMB, gets mechanically mixed with the ambient material, and forms neutrally buoyant large-scale compositional anomalies similar in shape to the LLSVPs.
A SSHAC Level 3 Probabilistic Seismic Hazard Analysis for a New-Build Nuclear Site in South Africa
(2015)
A probabilistic seismic hazard analysis has been conducted for a potential nuclear power plant site on the coast of South Africa, a country of low-to-moderate seismicity. The hazard study was conducted as a SSHAC Level 3 process, the first application of this approach outside North America. Extensive geological investigations identified five fault sources with a non-zero probability of being seismogenic. Five area sources were defined for distributed seismicity, the least active being the host zone for which the low recurrence rates for earthquakes were substantiated through investigations of historical seismicity. Empirical ground-motion prediction equations were adjusted to a horizon within the bedrock at the site using kappa values inferred from weak-motion analyses. These adjusted models were then scaled to create new equations capturing the range of epistemic uncertainty in this region with no strong motion recordings. Surface motions were obtained by convolving the bedrock motions with site amplification functions calculated using measured shear-wave velocity profiles.
Pseudotachylitic breccia (PTB) in the form of cm-wide melt breccia veinlets locally occurs on the exposed central uplift of the 380Ma Siljan impact structure. The host rock to the PTBs is the so-called Jarna granite of quartz monzonitic to syenodioritic composition. The nearly 603m long BH-5 drill core from Hattberg, near the centre of the Siljan central uplift, contains numerous veins and pods of PTB. In particular, two major zones of 60m combined width contain extensive PTB network breccias (30% actual melt breccia component), with individual melt breccia occurrences up to >1m in length. Core logging and petrographic and geochemical analysis of the core have been performed, and the data are interpreted to suggest the following. (1) The impact event caused low to moderate (at essentially <20GPa) shock deformation in the host rock and in clasts of this lithology within the PTB. (2) Macroscopic deformation of the basement mainly comprises fracturing, with only localised cataclasis. (3) No evidence for shock melting (i.e. compression/decompression melting early in the cratering process) could be observed. (4) Optical and scanning electron microscopy showed that dark PTB contains a definite melt component. (5) Shearing has significantly affected this part of the central uplift, but its effects are limited to very short displacements and likely did not result in extensive melting. (6) A frictional heating component upon melt generation can, however, not be excluded, as many PTB samples contain clasts of a mafic (gabbroic) component, although only in one place along the entire core, a 1.2cm-wide section through such material in direct contact to host rock was observed. Consequently, we suggest that, upon uplift in the central part of the impact structure, considerable melt volumes were generated locally, especially in areas that had been affected by extensive cataclasis and where grain size comminution favoured melt formation. Rapid decompression related to central uplift formation is the preferred process for the generation of the PTB melt breccias.
The Earth's biggest magmatic events are believed to originate from massive melting when hot mantle plumes rising from the lowermost mantle reach the base of the lithosphere. Classical models predict large plume heads that cause kilometre-scale surface uplift, and narrow (100 km radius) plume tails that remain in the mantle after the plume head spreads below the lithosphere. However, in many cases, such uplifts and narrow plume tails are not observed. Here using numerical models, we show that the issue can be resolved if major mantle plumes contain up to 15-20% of recycled oceanic crust in a form of dense eclogite, which drastically decreases their buoyancy and makes it depth dependent. We demonstrate that, despite their low buoyancy, large enough thermochemical plumes can rise through the whole mantle causing only negligible surface uplift. Their tails are bulky (4200 km radius) and remain in the upper mantle for 100 millions of years.
Faulting processes in volcanic areas result from a complex interaction of pressurized fluid-filled cracks and conduits with the host rock and local and regional tectonic setting. Often, volcanic seismicity is difficult to decipher in terms of the physical processes involved, and there is a need for models relating the mechanics of volcanic sources to observations. Here we use focal mechanism data of the energetic swarm induced by the 2000 dike intrusion at Miyakejima (Izu Archipelago, Japan), to study the relation between the 3-D dike-induced stresses and the characteristics of the seismicity. We perform a clustering analysis on the focal mechanism (FM) solutions and relate them to the dike stress field and to the scaling relationships of the earthquakes. We find that the strike and rake angles of the FMs are strongly correlated and cluster on bands in a strike-rake plot. We suggest that this is consistent with optimally oriented faults according to the expected pattern of Coulomb stress changes. We calculate the frequency-size distribution of the clustered sets finding that focal mechanisms with a large strike-slip component are consistent with the Gutenberg-Richter relation with a b value of about 1. Conversely, events with large normal faulting components deviate from the Gutenberg-Richter distribution with a marked roll-off on its right-hand tail, suggesting a lack of large-magnitude events (M-w>5.5). This may result from the interplay of the limited thickness and lower rock strength of the layer of rock above the dike, where normal faulting is expected, and lower stress levels linked to the faulting style and low confining pressure.
We present a new record of eolian dust flux to the western Subarctic North Pacific (SNP) covering the past 27,000years based on a core from the Detroit Seamount. Comparing the SNP dust record to the North Greenland Ice Core Project (NGRIP) ice core record shows significant differences in the amplitude of dust changes to the two regions during the last deglaciation, while the timing of abrupt changes is synchronous. If dust deposition in the SNP faithfully records its mobilization in East Asian source regions, then the difference in the relative amplitude must reflect climate-related changes in atmospheric dust transport to Greenland. Based on the synchronicity in the timing of dust changes in the SNP and Greenland, we tie abrupt deglacial transitions in the Th-230-normalized He-4 flux record to corresponding transitions in the well-dated NGRIP dust flux record to provide a new chronostratigraphic technique for marine sediments from the SNP. Results from this technique are complemented by radiocarbon dating, which allows us to independently constrain radiocarbon paleoreservoir ages. We find paleoreservoir ages of 745140years at 11,653year B.P., 680228years at 14,630year B.P., and 790498years at 23,290year B.P. Our reconstructed paleoreservoir ages are consistent with modern surface water reservoir ages in the western SNP. Good temporal synchronicity between eolian dust records from the Subantarctic Atlantic and equatorial Pacific and the ice core record from Antarctica supports the reliability of the proposed dust tuning method to be used more widely in other global ocean regions.
We present the results of our investigations on the radiocarbon dated core sediments from the Lake Tso Moriri, NW Himalaya aimed at reconstructing palaeohydrological changes in this climatically sensitive region. Based on the detailed geochemical, mineralogical and sedimentological analysis, we recognise several short-term fluctuations superimposed upon seven major palaeohydrological stages identified in this lake since similar to 26 cal ka. Stage I (>20.2 cal ka): shallow lake characterised by input of coarse-grained detrital sediments; Stage II (20.2-16.4 cal ka): lake deepening and intensification of this trend ca. 18 cal ka; Stage III (16.4-11.2 cal ka): rising lake levels with a short term wet phase (13.1-11.7 cal ka); Stage IV (11.2-8.5 cal ka): early Holocene hydrological maxima and highest lake levels inferred to have resulted from early Holocene Indian monsoon intensification, as records from central Asia indicate weaker westerlies during this interval; Stage V (8.5-5.5 cal ka): mid-Holocene climate deterioration; Stage VI (5.5-2.7 cal ka): progressive lowering of lake level; Stage VII (2.7-0 cal ka): onset of modern conditions. The reconstructed hydrological variability in Lake Tso Moriri is governed by temperature changes (meltwater inflow) and monsoon precipitation (increased runoff). A regional comparison shows considerable differences with other palaeorecords from peninsular India during late Holocene. (C) 2014 Elsevier Ltd and INQUA. All rights reserved.
As part of ongoing research on Holocene lacustrine sediments of Lonar Crater Lake (central India), pollen assemblages in lake surface sediment and soil samples were studied to unravel pollenevegetation relationships, including pollen transport processes in tropical dry deciduous forest vegetation. Furthermore, palynological results were compared with geochemical proxies and spatial features of the lake sediments and the vegetation. The obtained data reveal strong differences in pollen assemblages and pollen concentrations between and within the studied trapping media. Local arboreal vegetation is adequately represented in the soil samples, but is less represented in the lake surface sediment samples. The composition of the lacustrine pollen assemblages is mainly influenced by patterns of transport through surface and channel runoff. Besides the relevance of our new data for reliable interpretation of fossil pollen spectra extracted from Lonar sediment cores, the results of this study are of general importance for the understanding of Quaternary pollen assemblages from tropical lacustrine archives, as well as for the implementation and selection of suitable approaches for quantitative pollen based environmental reconstructions in south Asia and beyond. (C) 2015 Elsevier Ltd and INQUA. All rights reserved.
Individual great earthquakes are posited to release the elastic strain energy that has accumulated over centuries by the gradual movement of tectonic plates(1,2). However, knowledge of plate deformation during a complete seismic cycle-two successive great earthquakes and the intervening interseismic period-remains incomplete(3). A complete seismic cycle began in south-central Chile in 1835 with an earthquake of about magnitude 8.5 (refs 4,5) and ended in 2010 with a magnitude 8.8 earthquake(6). During the first earthquake, an uplift of Isla Santa Maria by 2.4 to 3m was documented(4,5). In the second earthquake, the island was uplifted(7) by 1.8 m. Here we use nautical surveys made in 1804, after the earthquake in 1835 and in 1886, together with modern echo sounder surveys and GPS measurements made immediately before and after the 2010 earthquake, to quantify vertical deformation through the complete seismic cycle. We find that in the period between the two earthquakes, Isla Santa Maria subsided by about 1.4 m. We simulate the patterns of vertical deformation with a finite-element model and find that they agree broadly with predictions from elastic rebound theory(2). However, comparison with geomorphic and geologic records of millennial coastline emergence(8,9) reveal that 10-20% of the vertical uplift could be permanent.
In the present study, we report the results of comprehensive amino acid (AA) analyses of four Indian lakes from different climate regimes. We focus on the investigation of sediment cores retrieved from the lakes but data of modern sediment as well as vascular plant, soil, and suspended particulate matter samples from individual lakes are also presented. Commonly used degradation and organic matter source indices are tested for their applicability to the lake sediments, and we discuss potential reasons for possible limitations. A principal component analysis including the monomeric AA composition of organic matter of all analysed samples indicates that differences in organic matter sources and the environmental properties of the individual lakes are responsible for the major variability in monomeric AA distribution of the different samples. However, the PCA also gives a factor that most probably separates the samples according to their state of organic matter degradation. Using the factor loadings of the individual AA monomers, we calculate a lake sediment degradation index (LI) that might be applicable to other palaeo-lake investigations. (C) 2015 Elsevier Ltd. All rights reserved.
This study focuses on evaluating the potential of ALOS/PALSAR time-series data to analyze the activation of deep-seated landslides in the foothill zone of the high mountain Alai range in the southern Tien Shan (Kyrgyzstan). Most previous field-based landslide investigations have revealed that many landslides have indicators for ongoing slow movements in the form of migrating and newly developing cracks. L-band ALOS/PALSAR data for the period between 2007 and 2010 are available for the 484 km(2) area in this study. We analyzed these data using the Small Baseline Subset (SBAS) time-series technique to assess the surface deformation related to the activation of landslides. We observed up to +/- 17 mm/year of LOS velocity deformation rates, which were projected along the local steepest slope and resulted in velocity rates of up to -63 mm/year. The obtained rates indicate very slow movement of the deep-seated landslides during the observation time. We also compared these movements with precipitation and earthquake records. The results suggest that the deformation peaks correlate with rainfall in the 3 preceding months and with an earthquake event. Overall, the results of this study indicated the great potential of L-band InSAR time series analysis for efficient spatiotemporal identification and monitoring of slope activations in this region of high landslide activity in Southern Kyrgyzstan.
Vegetated dunes are recognized as important natural barriers that shelter inland ecosystems and coastlines suffering daily erosive impacts of the sea and extreme events, such as tsunamis. However, societal responses to erosion and shoreline retreat often result in man-made coastal defence structures that cover part of the intertidal and upper shore zones causing coastal squeeze and habitat loss, especially for upper shore biota, such as dune plants. Coseismic uplift of up to 2.0 m on the Peninsula de Arauco (South central Chile, ca. 37.5 degrees S) caused by the 2010 Maule earthquake drastically modified the coastal landscape, including major increases in the width of uplifted beaches and the immediate conversion of mid to low sandy intertidal habitat to supralittoral sandy habitat above the reach of average tides and waves. To investigate the early stage responses in species richness, cover and across-shore distribution of the hitherto absent dune plants, we surveyed two formerly intertidal armoured sites and a nearby intertidal unarmoured site on a sandy beach located on the uplifted coast of Llico (Peninsula de Arauco) over two years. Almost 2 years after the 2010 earthquake, dune plants began to recruit, then rapidly grew and produced dune hummocks in the new upper beach habitats created by uplift at the three sites. Initial vegetation responses were very similar among sites. However, over the course of the study, the emerging vegetated dunes of the armoured sites suffered a slowdown in the development of the spatial distribution process, and remained impoverished in species richness and cover compared to the unarmoured site. Our results suggest that when released from the effects of coastal squeeze, vegetated dunes can recover without restoration actions. However, subsequent human activities and management of newly created beach and dune habitats can significantly alter the trajectory of vegetated dune development. Management that integrates the effects of natural and human induced disturbances, and promotes the development of dune vegetation as natural barriers can provide societal and conservation benefits in coastal ecosystems.
Climate change will alter the forces of predation and competition in temperate ectotherm food webs. This may increase local extinction rates, change the fate of invasions and impede species reintroductions into communities. Invasion success could be modulated by traits (e.g., defenses) and adaptations to climate. We studied how different temperatures affect the time until extinction of species, using bitrophic and tritrophic planktonic food webs to evaluate the relative importance of predatory overexploitation and competitive exclusion, at 15 and 25 A degrees C. In addition, we tested how inclusion of a subtropical as opposed to a temperate strain in this model food web affects times until extinction. Further, we studied the invasion success of the temperate rotifer Brachionus calyciflorus into the planktonic food web at 15 and 25 A degrees C on five consecutive introduction dates, during which the relative forces of predation and competition differed. A higher temperature dramatically shortened times until extinction of all herbivore species due to carnivorous overexploitation in tritrophic systems. Surprisingly, warming did not increase rates of competitive exclusion among the tested herbivore species in bitrophic communities. Including a subtropical herbivore strain reduced top-down control by the carnivore at high temperature. Invasion attempts of temperate B. calyciflorus into the food web always succeeded at 15 A degrees C, but consistently failed at 25 A degrees C due to voracious overexploitation by the carnivore. Pre-induction of defenses (spines) in B. calyciflorus before the invasion attempt did not change its invasion success at the high temperature. We conclude that high temperatures may promote local extinctions in temperate ectotherms and reduce their chances of successful recovery.
This paper reports on extremely thick and massive duricrusts in soils of two basins in the Alentejo (southern Portugal). Since different types of duricrusts (calcretes, silcretes and palycretes) have been reported from other regions in the Mediterranean, the purpose of this study was to identify the cementing agents in the duricrusts and to compare their composition in the two basins. Moreover, the study aimed at identifying the processes involved in duricrust formation, and especially the role of topography and lateral water and element transport in the landscape, and drawing conclusions about environmental conditions and time of duricrust formation.
After studying an extensive number of road cuts in the field and mapping soil patterns in parts of the two basins by manual augering, ten pedons were selected for detailed description and sampling. Thin sections were analysed under a petrographic microscope, focusing on the micromorphology and optical properties of the cementing materials. Selected samples were studied by scanning electron microscopy and energy-dispersive X-ray spectroscopy to reconfirm the optical identification. The laboratory analyses included pH, carbonate contents, and X-ray diffraction analysis.
The duricrusts in the eastern Sado basin are indurated by silica. Combination of XRD and thin section analysis allowed to identify opal-CT as a major component, while opal-A is present to a lesser extent, and chalcedony is very rare. The cementing materials of the duricrusts in the Oriola basin are palygorskite and calcite, which may occur alone or in combination within a soil profile.
The thick duricrusts formed in the basins through precipitation of calcite, palygorskite and silica from lateral water flows, which ran from the Serra de Portel into the basins, during short moist seasons in a generally warm, semi-arid climate with strong evapotranspiration. Lithology of the upper catchment areas (element sources) and topography control the spatial distribution of the different duricrusts. Their formation took place mainly during the Pliocene. Palygorskite transformation to smectite in the upper parts of the palycretes indicates that palygorskite is unstable in the present (more humid, typical Mediterranean) climate. This study demonstrates the potential role of lateral water and element transport in landscapes that need to be considered in pedological studies and concepts, and the use of mineral indicators of past climates such as palygorskite and the ageing stage of silica precipitations as tools for reconstructing environmental conditions and possible time of duricrust formation.
We herein reappraise the pressure-temperature (PT) evolution of the high-pressure and low-temperature (HP-LT) Tavsanli zone (western Turkey) in order to (i) better characterize rock units exhumed along a cooling subduction interface, from birth to steady state and (ii) constrain exhumation and detachment dynamics, as well as mechanical coupling between plates. Based on PT estimates and field observations three oceanic complexes are recognized between the HP-LT continental margin and the obducted ophiolite, with PT estimates ranging from incipient metamorphism to blueschist-fades conditions. PT conditions for the continental unit are reappraised to 24 kbar and similar to 500 degrees C on the basis of pseudosection modelling and Raman spectroscopy on carbonaceous material. A tentative reconstruction of the subduction zone evolution is proposed using available radiometric and palaeogeographic data and recent thermomechanical modelling. Both PT conditions and field observations point out to the slicing of km-sized units at different preferred depths along the subduction interface, thus providing constraints on the dynamics of accretion and underplating. In particular, the comparison of PT estimates for the Tavsanli zone and for other broadly similar fossil subduction settings (i.e., Oman, Corsica, New Caledonia, Franciscan, Schistes Lustres) suggests that units are detached preferentially from the slab at specific depths of 30-40 km (i.e., downdip of the seismogenic zone) and similar to 80 km. We propose that these depths are controlled by major changes in mechanical coupling along the plate interface, whereas exhumation through time would rather be controlled by large-scale geodynamic boundary conditions. (C) 2015 Elsevier B.V. All rights reserved.
Background: GEOMAGIA50.v3 for sediments is a comprehensive online database providing access to published paleomagnetic, rock magnetic, and chronological data obtained from lake and marine sediments deposited over the past 50 ka. Its objective is to catalogue data that will improve our understanding of changes in the geomagnetic field, physical environments, and climate.
Findings: GEOMAGIA50.v3 for sediments builds upon the structure of the pre-existing GEOMAGIA50 database for magnetic data from archeological and volcanic materials. A strong emphasis has been placed on the storage of geochronological data, and it is the first magnetic archive that includes comprehensive radiocarbon age data from sediments. The database will be updated as new sediment data become available.
Conclusions: The web-based interface for the sediment database is located at http://geomagia.gfz-potsdam.de/geomagiav3/SDquery.php. This paper is a companion to Brown et al. (Earth Planets Space doi:10.1186/s40623-015-0232-0,2015) and describes the data types, structure, and functionality of the sediment database.
Hydrothermal carbonization (HTC) produces carbon-rich nano-micro size particles. In this study, magnetic hydrochar (MHC) was prepared from model compound cellulose by simply adding ferrites during HTC. The effects of ferrites on HTC were evaluated by characterizing solid MHC and corresponding process liquid. Additionally, magnetic stability of MHC was tested by magnetic susceptibility method. Finally, MHC was used as support media for anaerobic films in anaerobic digestion (AD). Ash-free mass yield was around 50% less in MHC than hydrochar produced without ferrites at any certain HTC reaction condition, where organic part of MHC is mainly carbon. In fact, amorphous hydrochar was growing on the surface of inorganic ferrites. MHC maintained magnetic susceptibility regardless of reaction time at reaction temperature 250 degrees C. Pronounced inhibitory effects of magnetic hydrochar occurred during start-up of AD but diminished with prolong AD times. Visible biofilms were observed on the MHC by laser scanning microscope after AD. (C) 2015 Elsevier Ltd. All rights reserved.
Computer-based simulation models are frequently used in hydrological research and engineering but also in other fields of environmental sciences. New case studies often require existing model concepts to be adapted. Extensions may be necessary due to the peculiarities of the studied natural system or subtleties of anthropogenic control. In other cases, simplifications must be made in response to scarce data, incomplete knowledge, or restrictions set by the spatio-temporal scale of application. This paper introduces an open-source modeling framework called ECHSE designed to cope with the above-mentioned challenges. It provides a lightweight infrastructure for the rapid development of new, reusable simulation tools and, more importantly, the safe modification of existing formulations. ECHSE-based models treat the simulated system as a collection of interacting objects. Although feedbacks are generally supported, the majority of the objects' interactions is expected to be of the feed-forward type. Therefore, the ECHSE software is particularly useful in the context of hydrological catchment modeling. Conversely, it is unsuitable, e.g., for fully hydrodynamic simulations and groundwater flow modeling. The focus of the paper is put on a comprehensible outline of the ECHSE's fundamental concepts and limitations. For the purpose of illustration, a specific, ECHSE-based solution for hydrological catchment modeling is presented which has undergone testing in a number of river basins. (C) 2015 Elsevier Ltd. All rights reserved.
The topographic signature of a mountain belt depends on the interplay of tectonic, climatic and erosional processes, whose relative importance changes over times, while quantifying these processes and their rates at specific times remains a challenge. The eastern Andes of central Bolivia offer a natural laboratory in which such interplay has been debated. Here, we investigate the Rio Grande catchment which crosses orthogonally the eastern Andes orogen from the Eastern Cordillera into the Subandean Zone, exhibiting a catchment relief of up to 5000 m. Despite an enhanced tectonic activity in the Subandes, local relief, mean and modal slopes and channel steepness indices are largely similar compared to the Eastern Cordillera and the intervening Interandean Zone. Nevertheless, a dataset of 57 new cosmogenic 10Be and 26AI catchment wide denudation rates from the Rio Grande catchment reveals up to one order of magnitude higher denudation rates in the Subandean Zone (mean 0.8 mm/yr) compared to the upstream physiographic regions. We infer that tectonic activity in the thrusting dominated Subandean belt causes higher denudation rates based on cumulative rock uplift investigations and due to the absence of a pronounced climate gradient. Furthermore, the lower rock strength of the Subandean sedimentary units correlates with mean slopes similar to the ones of the Eastern Cordillera and Interandean Zone, highlighting the fact, that lithology and rock strength can control high denudation rates at low slopes.
Low denudation rates measured at the outlet of the Rio Grande catchment (Abapo) are interpreted to be a result of a biased cosmogenic nuclide mixing that is dominated by headwater signals from the Eastern Cordillera and the Interandean zone and limited catchment sediment connectivity in the lower river reaches. Therefore, comparisons of short- (i.e., sediment yield) and millennial denudation rates require caution when postulating tectonic and/or climatic forcing without detailed studies. (C) 2015 The Authors. Published by Elsevier B.V.