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)
The mechanisms by which climate and vegetation affect erosion rates over various time scales lie at the heart of understanding landscape response to climate change. Plot-scale field experiments show that increased vegetation cover slows erosion, implying that faster erosion should occur under low to moderate vegetation cover. However, demonstrating this concept over long time scales and across landscapes has proven to be difficult, especially in settings complicated by tectonic forcing and variable slopes. We investigate this problem by measuring cosmogenic Be-10-derived catchment-mean denudation rates across a range of climate zones and hillslope gradients in the Kenya Rift, and by comparing our results with those published from the Rwenzori Mountains of Uganda. We find that denudation rates from sparsely vegetated parts of the Kenya Rift are up to 0.13 mm/yr, while those from humid and more densely vegetated parts of the Kenya Rift flanks and the Rwenzori Mountains reach a maximum of 0.08 mm/yr, despite higher median hillslope gradients. While differences in lithology and recent land-use changes likely affect the denudation rates and vegetation cover values in some of our studied catchments, hillslope gradient and vegetation cover appear to explain most of the variation in denudation rates across the study area. Our results support the idea that changing vegetation cover can contribute to complex erosional responses to climate or land-use change and that vegetation cover can play an important role in determining the steady-state slopes of mountain belts through its stabilizing effects on the land surface.
In the last decade, the number and dimensions of catastrophic flooding events in the Niger River Basin (NRB) have markedly increased. Despite the devastating impact of the floods on the population and the mainly agriculturally based economy of the riverine nations, awareness of the hazards in policy and science is still low. The urgency of this topic and the existing research deficits are the motivation for the present dissertation.
The thesis is an initial detailed assessment of the increasing flood risk in the NRB. The research strategy is based on four questions regarding (1) features of the change in flood risk, (2) reasons for the change in the flood regime, (3) expected changes of the flood regime given climate and land use changes, and (4) recommendations from previous analysis for reducing the flood risk in the NRB.
The question examining the features of change in the flood regime is answered by means of statistical analysis. Trend, correlation, changepoint, and variance analyses show that, in addition to the factors exposure and vulnerability, the hazard itself has also increased significantly in the NRB, in accordance with the decadal climate pattern of West Africa. The northern arid and semi-arid parts of the NRB are those most affected by the changes.
As potential reasons for the increase in flood magnitudes, climate and land use changes are attributed by means of a hypothesis-testing framework. Two different approaches, based on either data analysis or simulation, lead to similar results, showing that the influence of climatic changes is generally larger compared to that of land use changes. Only in the dry areas of the NRB is the influence of land use changes comparable to that of climatic alterations.
Future changes of the flood regime are evaluated using modelling results. First ensembles of statistically and dynamically downscaled climate models based on different emission scenarios are analyzed. The models agree with a distinct increase in temperature. The precipitation signal, however, is not coherent. The climate scenarios are used to drive an eco-hydrological model. The influence of climatic changes on the flood regime is uncertain due to the unclear precipitation signal. Still, in general, higher flood peaks are expected. In a next step, effects of land use changes are integrated into the model. Different scenarios show that regreening might help to reduce flood peaks. In contrast, an expansion of agriculture might enhance the flood peaks in the NRB. Similarly to the analysis of observed changes in the flood regime, the impacts of climate- and land use changes for the future scenarios are also most severe in the dry areas of the NRB.
In order to answer the final research question, the results of the above analysis are integrated into a range of recommendations for science and policy on how to reduce flood risk in the NRB. The main recommendations include a stronger consideration of the enormous natural climate variability in the NRB and a focus on so called “no-regret” adaptation strategies which account for high uncertainty, as well as a stronger consideration of regional differences. Regarding the prevention and mitigation of catastrophic flooding, the most vulnerable and sensitive areas in the basin, the arid and semi-arid Sahelian and Sudano-Sahelian regions, should be prioritized. Eventually, an active, science-based and science-guided flood policy is recommended. The enormous population growth in the NRB in connection with the expected deterioration of environmental and climatic conditions is likely to enhance the region´s vulnerability to flooding. A smart and sustainable flood policy can help mitigate these negative impacts of flooding on the development of riverine societies in West Africa.
Climate or Land Use?
(2015)
This study intends to contribute to the ongoing discussion on whether land use and land cover changes (LULC) or climate trends have the major influence on the observed increase of flood magnitudes in the Sahel. A simulation-based approach is used for attributing the observed trends to the postulated drivers. For this purpose, the ecohydrological model SWIM (Soil and Water Integrated Model) with a new, dynamic LULC module was set up for the Sahelian part of the Niger River until Niamey, including the main tributaries Sirba and Goroul. The model was driven with observed, reanalyzed climate and LULC data for the years 1950-2009. In order to quantify the shares of influence, one simulation was carried out with constant land cover as of 1950, and one including LULC. As quantitative measure, the gradients of the simulated trends were compared to the observed trend. The modeling studies showed that for the Sirba River only the simulation which included LULC was able to reproduce the observed trend. The simulation without LULC showed a positive trend for flood magnitudes, but underestimated the trend significantly. For the Goroul River and the local flood of the Niger River at Niamey, the simulations were only partly able to reproduce the observed trend. In conclusion, the new LULC module enabled some first quantitative insights into the relative influence of LULC and climatic changes. For the Sirba catchment, the results imply that LULC and climatic changes contribute in roughly equal shares to the observed increase in flooding. For the other parts of the subcatchment, the results are less clear but show, that climatic changes and LULC are drivers for the flood increase; however their shares cannot be quantified. Based on these modeling results, we argue for a two-pillar adaptation strategy to reduce current and future flood risk: Flood mitigation for reducing LULC-induced flood increase, and flood adaptation for a general reduction of flood vulnerability.
Samples of 474 forest stands in Germany were analysed for concentrations of polycyclic aromatic hydrocarbons (PAHs) in three sampling depths. Enhanced concentrations were mainly found at spots relatively close to densely industrialized and urbanized regions and at some topographically elevated areas. Average enrichment factors between mineral soil and humic layer depend on humus type i.e. decrease from mull via moder to more Based on their compound-patterns, the observed samples could be assigned to three main clusters. For some parts of our study area a uniform assignment of samples to clusters over larger regions could be identified. For instance, samples taken at vicinity to brown-coal strip-mining districts are characterized by high relative abundances of low-molecular-weight PAHs. These results suggest that PAHs are more likely originated from local and regional emitters rather than from long-range transport and that specific source-regions can be identified based on PAH fingerprints. (C) 2015 Elsevier Ltd. All rights reserved.
Central Asia is located at the confluence of large-scale atmospheric circulation systems. It is thus likely to be highly susceptible to changes in the dynamics of those systems; however, little is still known about the regional paleoclimate history. Here we present carbon and hydrogen isotopic compositions of n-alkanoic acids from a late Holocene sediment core from Lake Karakuli (eastern Pamir, Xinjiang Province, China). Instrumental evidence and isotopeenabled climate model experiments with the Laboratoire de Meteorologie Dynamique Zoom model version 4 (LMDZ4) demonstrate that delta D values of precipitation in the region are influenced by both temperature and precipitation amount. We find that these parameters are inversely correlated on an annual scale, i.e., the climate has varied between relatively cool and wet and more warm and dry over the last 50 years. Since the isotopic signals of these changes are in the same direction and therefore additive, isotopes in precipitation are sensitive recorders of climatic changes in the region. Additionally, we infer that plants use year-round precipitation (including snowmelt), and thus leaf wax delta D values must also respond to shifts in the proportion of moisture derived from westerly storms during late winter and early spring. Downcore results give evidence for a gradual shift to cooler and wetter climates between 3.5 and 2.5 cal kyr BP, interrupted by a warm and dry episode between 3.0 and 2.7 kyr BP. Further cool and wet episodes occur between 1.9 and 1.5 and between 0.6 and 0.1 kyr BP, the latter coeval with the Little Ice Age. Warm and dry episodes from 2.5 to 1.9 and 1.5 to 0.6 kyr BP coincide with the Roman Warm Period and Medieval Climate Anomaly, respectively. Finally, we find a drying tend in recent decades. Regional comparisons lead us to infer that the strength and position of the westerlies, and wider northern hemispheric climate dynamics, control climatic shifts in arid Central Asia, leading to complex local responses. Our new archive from Lake Karakuli provides a detailed record of the local signatures of these climate transitions in the eastern Pamir.
Topographic migration of 2D and 3D ground-penetrating radar data considering variable velocities
(2015)
We present a 2D/3D topographic migration scheme for ground-penetrating radar (GPR) data which is able to account for variable velocities by using the root mean square (rms) velocity approximation. We test our migration scheme using a synthetic 2D example and compare our migrated image to the results obtained using common GPR migration approaches. Furthermore, we apply it to 2D and 3D field data. These examples are recorded across common subsurface settings including surface topography and variations in the GPR subsurface velocity field caused by a shallow ground water table. In such field settings, our migration strategy provides well focused images of commonoffset GPR data without the need for a detailed interval velocity model. The synthetic and field examples demonstrate that our topographic migration scheme allows for accurate GPR imaging in the presence of variations in surface topography and subsurface velocity.
Flow phenomena in the unsaturated zone are highly variable in time and space. Thus, it is challenging to measure and monitor such processes under field conditions. Here, we present a new setup and interpretation approach for combining a dye tracer experiment with a 4D ground-penetrating radar (GPR) survey. Therefore, we designed a rainfall experiment during which we measured three surface-based 3D GPR surveys using a pair of 500 MHz antennas. Such a survey setup requires accurate acquisition and processing techniquesto extract time-lapse information supporting the interpretation of selected cross-sections photographed after excavating the site. Our results reveal patterns of traveltime changes in the measured GPR data, which are associated with soil moisture changes. As distinct horizons are present at our site, such changes can be quantified and transferred into changes in total soil moisture content. Our soil moisture estimates are similar to the amount of infiltrated water, which confirms our experimental approach and makes us confident for further developing this strategy, especially, with respect to improving the temporal and spatial resolution. (C) 2015 Elsevier B.V. All rights reserved.
Mapping urban forest leaf area index with airborne lidar using penetration metrics and allometry
(2015)
In urban areas, leaf area index (LAI) is a key ecosystem structural attribute with implications for energy and water balance, gas exchange, and anthropogenic energy use. In this study, we estimated LAI spatially using airborne lidar in downtown Santa Barbara, California, USA. We implemented two different modeling approaches. First, we directly estimated effective LAI (LAIe) using scan angle- and clump-corrected lidar laser penetration metrics (LPM). Second, we adapted existing allometric equations to estimate crown structural metrics including tree height and crown base height using lidar. The latter approach allowed for LAI estimates at the individual tree-crown scale. The LPM method, at both high and decimated point densities, resulted in good linear agreement with estimates from ground-based hemispherical photography (r(2) = 0.82, y = 0.99x) using a model that assumed a spherical leaf angle distribution. Within individual tree crown segments, the lidar estimates of crown structure closely paralleled field measurements (e.g., r(2) = 0.87 for crown length). LAI estimates based on the lidar crown measurements corresponded well with estimates from field measurements (r(2) = 0.84, y = 0.97x + 0.10). Consistency of the LPM and allometric lidar methods was also strong at 71 validation plots (r(2) = 0.88) and at 450 additional sample locations across the entire study area (r(2) = 0.72). This level of correspondence exceeded that of the canopy hemispherical photography and allometric, ground-based estimates (r(2) = 0.53). The first-order alignment of these two disparate methods may indicate that the error bounds for mapping LAI in cities are small enough to pursue large scale, spatially explicit estimation. (C) 2015 Elsevier Inc All rights reserved.
With accelerating climate cooling in the late Cenozoic, glacial and periglacial erosion became more widespread on the surface of the Earth. The resultant shift in erosion patterns significantly changed the large-scale morphology of many mountain ranges worldwide. Whereas the glacial fingerprint is easily distinguished by its characteristic fjords and U-shaped valleys, the periglacial fingerprint is more subtle but potentially prevails in some mid- to high-latitude landscapes. Previous models have advocated a frost-driven control on debris production at steep headwalls and glacial valley sides. Here we investigate the important role that periglacial processes also play in less steep parts of mountain landscapes. Understanding the influences of frost-driven processes in low-relief areas requires a focus on the consequences of an accreting soil mantle, which characterises such surfaces. We present a new model that quantifies two key physical processes: frost cracking and frost creep, as a function of both temperature and sediment thickness. Our results yield new insights into how climate and sediment transport properties combine to scale the intensity of periglacial processes. The thickness of the soil mantle strongly modulates the relation between climate and the intensity of mechanical weathering and sediment flux. Our results also point to an offset between the conditions that promote frost cracking and those that promote frost creep, indicating that a stable climate can provide optimal conditions for only one of those processes at a time. Finally, quantifying these relations also opens up the possibility of including periglacial processes in large-scale, long-term landscape evolution models, as demonstrated in a companion paper.
The Central Pontides is an accretionary-type orogenic area within the Alpine-Himalayan orogenic belt characterized by pre-collisional tectonic continental growth. The region comprises Mesozoic subduction-accretionary complexes and an accreted intra-oceanic arc that are sandwiched between the Laurasian active continental margin and Gondwana-derived the Kırşehir Block. The subduction-accretion complexes mainly consist of an Albian-Turonian accretionary wedge representing the Laurasian active continental margin. To the north, the wedge consists of slate/phyllite and metasandstone intercalation with recrystallized limestone, Na-amphibole-bearing metabasite (PT= 7–12 kbar and 400 ± 70 ºC) and tectonic slices of serpentinite representing accreted distal part of a large Lower Cretaceous submarine turbidite fan deposited on the Laurasian active continental margin that was subsequently accreted and metamorphosed. Raman spectra of carbonaceous material (RSCM) of the metapelitic rocks revealed that the metaflysch sequence consists of metamorphic packets with distinct peak metamorphic temperatures. The majority of the metapelites are low-temperature (ca. 330 °C) slates characterized by lack of differentiation of the graphite (G) and D2 defect bands. They possibly represent offscraped distal turbidites along the toe of the Albian accretionary wedge. The rest are phyllites that are characterized by slightly pronounced G band with D2 defect band occurring on its shoulder. Peak metamorphic temperatures of these phyllites are constrained to 370-385 °C. The phyllites are associated with a strip of incipient blueschist facies metabasites which are found as slivers within the offscraped distal turbidites. They possibly represent underplated continental metasediments together with oceanic crustal basalt along the basal décollement. Tectonic emplacement of the underplated rocks into the offscraped distal turbidites was possibly achieved by out-of-sequence thrusting causing tectonic thickening and uplift of the wedge. 40Ar/39Ar phengite ages from the phyllites are ca. 100 Ma, indicating Albian subduction and regional HP metamorphism.
The accreted continental metasediments are underlain by HP/LT metamorphic rocks of oceanic origin along an extensional shear zone. The oceanic metamorphic sequence mainly comprises tectonically thickened deep-seated eclogite to blueschist facies metabasites and micaschists. In the studied area, metabasites are epidote-blueschists locally with garnet (PT= 17 ± 1 kbar and 500 ± 40 °C). Lawsonite-blueschists are exposed as blocks along the extensional shear zone (PT= 14 ± 2 kbar and 370–440 °C). They are possibly associated with low shear stress regime of the initial stage of convergence. Close to the shear zone, the footwall micaschists consist of quartz, phengite, paragonite, chlorite, rutile with syn-kinematic albite porphyroblast formed by pervasive shearing during exhumation. These types of micaschists are tourmaline-bearing and their retrograde nature suggests high-fluid flux along shear zones. Peak metamorphic mineral assemblages are partly preserved in the chloritoid-micaschist farther away from the shear zone representing the zero strain domains during exhumation. Three peak metamorphic assemblages are identified and their PT conditions are constrained by pseudosections produced by Theriak-Domino and by Raman spectra of carbonaceous material: 1) garnet-chloritoid-glaucophane with lawsonite pseudomorphs (P= 17.5 ± 1 kbar, T: 390-450 °C) 2) chloritoid with glaucophane pseudomorphs (P= 16-18 kbar, T: 475 ± 40 °C) and 3) relatively high-Mg chloritoid (17%) with jadeite pseudomorphs (P= 22-25 kbar; T: 440 ± 30 °C) in addition to phengite, paragonite, quartz, chlorite, rutile and apatite. The last mineral assemblage is interpreted as transformation of the chloritoid + glaucophane assemblage to chloritoid + jadeite paragenesis with increasing pressure. Absence of tourmaline suggests that the chloritoid-micaschist did not interact with B-rich fluids during zero strain exhumation. 40Ar/39Ar phengite age of a pervasively sheared footwall micaschist is constrained to 100.6 ± 1.3 Ma and that of a chloritoid-micaschist is constrained to 91.8 ± 1.8 Ma suggesting exhumation during on-going subduction with a southward younging of the basal accretion and the regional metamorphism. To the south, accretionary wedge consists of blueschist and greenschist facies metabasite, marble and volcanogenic metasediment intercalation. 40Ar/39Ar phengite dating reveals that this part of the wedge is of Middle Jurassic age partly overprinted during the Albian. Emplacement of the Middle Jurassic subduction-accretion complexes is possibly associated with obliquity of the Albian convergence.
Peak metamorphic assemblages and PT estimates of the deep-seated oceanic metamorphic sequence suggest tectonic stacking within wedge with different depths of burial. Coupling and exhumation of the distinct metamorphic slices are controlled by decompression of the wedge possibly along a retreating slab. Structurally, decompression of the wedge is evident by an extensional shear zone and the footwall micaschists with syn-kinematic albite porphyroblasts. Post-kinematic garnets with increasing grossular content and pseudomorphing minerals within the chloritoid-micaschists also support decompression model without an extra heating.
Thickening of subduction-accretionary complexes is attributed to i) significant amount of clastic sediment supply from the overriding continental domain and ii) deep level basal underplating by propagation of the décollement along a retreating slab. Underplating by basal décollement propagation and subsequent exhumation of the deep-seated subduction-accretion complexes are connected and controlled by slab rollback creating a necessary space for progressive basal accretion along the plate interface and extension of the wedge above for exhumation of the tectonically thickened metamorphic sequences. This might be the most common mechanism of the tectonic thickening and subsequent exhumation of deep-seated HP/LT subduction-accretion complexes.
To the south, the Albian-Turonian accretionary wedge structurally overlies a low-grade volcanic arc sequence consisting of low-grade metavolcanic rocks and overlying metasedimentary succession is exposed north of the İzmir-Ankara-Erzincan suture (İAES), separating Laurasia from Gondwana-derived terranes. The metavolcanic rocks mainly consist of basaltic andesite/andesite and mafic cognate xenolith-bearing rhyolite with their pyroclastic equivalents, which are interbedded with recrystallized pelagic limestone and chert. The metavolcanic rocks are stratigraphically overlain by recrystallized micritic limestone with rare volcanogenic metaclastic rocks. Two groups can be identified based on trace and rare earth element characteristics. The first group consists of basaltic andesite/andesite (BA1) and rhyolite with abundant cognate gabbroic xenoliths. It is characterized by relative enrichment of LREE with respect to HREE. The rocks are enriched in fluid mobile LILE, and strongly depleted in Ti and P reflecting fractionation of Fe-Ti oxides and apatite, which are found in the mafic cognate xenoliths. Abundant cognate gabbroic xenoliths and identical trace and rare earth elements compositions suggest that rhyolites and basaltic andesites/andesites (BA1) are cogenetic and felsic rocks were derived from a common mafic parental magma by fractional crystallization and accumulation processes. The second group consists only of basaltic andesites (BA2) with flat REE pattern resembling island arc tholeiites. Although enriched in LILE, this group is not depleted in Ti or P.
Geochemistry of the metavolcanic rocks indicates supra-subduction volcanism evidenced by depletion of HFSE and enrichment of LILE. The arc sequence is sandwiched between an Albian-Turonian subduction-accretionary complex representing the Laurasian active margin and an ophiolitic mélange. Absence of continent derived detritus in the arc sequence and its tectonic setting in a wide Cretaceous accretionary complex suggest that the Kösdağ Arc was intra-oceanic. This is in accordance with basaltic andesites (BA2) with island arc tholeiite REE pattern.
Zircons from two metarhyolite samples give Late Cretaceous (93.8 ± 1.9 and 94.4 ± 1.9 Ma) U/Pb ages. Low-grade regional metamorphism of the intra-oceanic arc sequence is constrained 69.9 ± 0.4 Ma by 40Ar/39Ar dating on metamorphic muscovite from a metarhyolite indicating that the arc sequence became part of a wide Tethyan Cretaceous accretionary complex by the latest Cretaceous. The youngest 40Ar/39Ar phengite age from the overlying subduction-accretion complexes is 92 Ma confirming southward younging of an accretionary-type orogenic belt. Hence, the arc sequence represents an intra-oceanic paleo-arc that formed above the sinking Tethyan slab and finally accreted to Laurasian active continental margin. Abrupt non-collisional termination of arc volcanism was possibly associated with southward migration of the arc volcanism similar to the Izu-Bonin-Mariana arc system.
The intra-oceanic Kösdağ Arc is coeval with the obducted supra-subduction ophiolites in NW Turkey suggesting that it represents part of the presumed but missing incipient intra-oceanic arc associated with the generation of the regional supra-subduction ophiolites. Remnants of a Late Cretaceous intra-oceanic paleo-arc and supra-subduction ophiolites can be traced eastward within the Alp-Himalayan orogenic belt. This reveals that Late Cretaceous intra-oceanic subduction occurred as connected event above the sinking Tethyan slab. It resulted as arc accretion to Laurasian active margin and supra-subduction ophiolite obduction on Gondwana-derived terranes.
A tectonic slice of an arc sequence consisting of low-grade metavolcanic rocks and overlying metasedimentary succession is exposed in the Central Pontides north of the Izmir-Ankara-Erzincan suture 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 metasedimentary succession comprises recrystallized micritic limestone with rare volcanogenic metaclastic rocks and stratigraphically overlies the metavolcanic rocks. The geochemistry of the metavolcanic rocks indicates an arc setting evidenced by depletion of HFSE (Ti, P and Nb) and enrichment of fluid mobile LILE. Identical trace and rare earth elements compositions of basaltic andesites/andesites and rhyolites suggest that they are cogenetic and derived from a common parental magma. The arc sequence crops out between an Albian-Turonian subduction-accretionary complex representing the Laurasian active margin and an ophiolitic melange. Absence of continent derived detritus in the arc sequence and its tectonic setting in a wide Cretaceous accretionary complex suggest that the Kosdag Arc was intra-oceanic. Zircons from two metarhyolite samples give Late Cretaceous (93.8 +/- 1.9 and 94.4 +/- 1.9 Ma) U/Pb ages. These ages are the same as the age of the supra-subduction ophiolites in western Turkey, which implies that that the Kosdag Arc may represent part of the incipient arc formed during the generation of the supra-subduction ophiolites. The low-grade regional metamorphism in the Kosdag Arc is constrained to 69.9 +/- 0.4 Ma by Ar-40/Ar-39 muscovite dating indicating that the arc sequence became part of a wide Tethyan Cretaceous accretionary complex by the latest Cretaceous. Non-collisional cessation of the arc volcanism is possibly associated with southward migration of the magmatism as in the Izu-Bonin-Mariana arc system. (c) 2015 Elsevier Ltd. All rights reserved.
Albian-Turonian subduction-accretionary complexes are exposed widely in the Central Pontides. A major portion of the accretionary complexes is made up of a metaflysch sequence consisting of slate/phyllite and metasandstone intercalation with blocks of marble, Na-amphibole bearing metabasite, and serpentinite. The metaflysch sequence represents distal parts of a large Lower Cretaceous submarine turbidite fan deposited on the Laurasian active continental margin that was subsequently accreted and metamorphosed during the Albian. Raman spectra of carbonaceous material of the metapelitic rocks revealed that the metaflysch consists of metamorphic packets with distinct peak metamorphic temperatures. The majority of the metapelites are low-temperature (ca. 330 degrees 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. Other phyllites are characterized by a slightly pronounced G band with a D2 defect band occurring on its shoulder. Peak metamorphic temperatures of these phyllites are constrained to 370-385 degrees C. The phyllites are associated with a strip of incipient blueschist facies metabasites and are found as a sliver within the offscraped distal turbidites. We interpret the phyllites as underplated continental sediments together with oceanic crustal basalt along the basal decollement. 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.
The idea that climatically modulated erosion may impact orogenic processes has challenged geoscientists for decades. Although modeling studies and physical calculations have provided a solid theoretical basis supporting this interaction, to date, field-based work has produced inconclusive results. The central-western Alborz Mountains in the northern sectors of the Arabia-Eurasia collision zone constitute a promising area to explore these potential feedbacks. This region is characterized by asymmetric precipitation superimposed on an orogen with a history of spatiotemporal changes in exhumation rates, deformation patterns, and prolonged, km-scale base-level changes. Our analysis suggests that despite the existence of a strong climatic gradient at least since 17.5 Ma, the early orogenic evolution (from similar to 36 to 9-6 Ma) was characterized by decoupled orographic precipitation and tectonics. In particular, faster exhumation and sedimentation along the more arid southern orogenic flank point to a north-directed accretionary flux and underthrusting of Central Iran. Conversely, from 6 to 3 Ma, erosion rates along the northern orogenic flank became higher than those in the south, where they dropped to minimum values. This change occurred during a similar to 3-Myr-long, km-scale base-level lowering event in the Caspian Sea. We speculate that mass redistribution processes along the northern flank of the Alborz and presumably across all mountain belts adjacent to the South Caspian Basin and more stable areas of the Eurasian plate increased the sediment load in the basin and ultimately led to the underthrusting of the Caspian Basin beneath the Alborz Mountains. This underthrusting in turn triggered a new phase of northward orogenic expansion, transformed the wetter northern flank into a new pro-wedge, and led to the establishment of apparent steady-state conditions along the northern orogenic flank (i.e., rock uplift equal to erosion rates). Conversely, the southern mountain front became the retro-wedge and experienced limited tectonic activity. These observations overall raise the possibility that mass-distribution processes during a pronounced erosion phase driven by base-level changes may have contributed to the inferred regional plate-tectonic reorganization of the northern Arabia-Eurasia collision during the last similar to 5 Ma. (C) 2015 Elsevier B.V. All rights reserved.
Cosmic-Ray neutron sensing (CRS) is a unique approach to measure soil moisture at field scale filling the gap of current methodologies. However, CRS signal is affected by all the hydrogen pools on the land surface and understanding their relative importance plays an important role for the application of the method e.g., validation of remote sensing products and data assimilation. In this study, a soil moisture scaling approach is proposed to estimate directly the correct CRS soil moisture based on the soil moisture profile measured at least in one position within the field. The approach has the advantage to avoid the need to introduce one correction for each hydrogen contribution and to estimate indirectly all the related time-varying hydrogen pools. Based on the data collected in three crop seasons, the scaling approach shows its ability to identify and to quantify the seasonal biomass water equivalent. Additionally, the analysis conducted at sub-daily time resolution is able to quantify the daily vertical redistribution of the water biomass and the rainfall interception, showing promising applications of the CRS method also for these types of measurements. Overall, the study underlines how not only soil moisture but all the specific hydrological processes in the soil-plant-atmosphere continuum should be considered for a proper evaluation of the CRS signal. For this scope, the scaling approach reveals to be a simple and pragmatic analysis that can be easily extended to other experimental sites. (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.
In the Sri Lankan Highlands erosion and chemical weathering rates are among the lowest for global mountain denudation. In this tropical humid setting, highly weathered deep saprolite profiles have developed from high-grade metamorphic charnockite during spheroidal weathering of the bedrock. The spheroidal weathering produces rounded corestones and spalled rindlets at the rock-saprolite interface. We used detailed textural, mineralogical and chemical analyses to reconstruct the sequence of weathering reactions and their causes. The first mineral attacked by weathering was found to be pyroxene initiated by in situ Fe oxidation. Volumetric calculations suggest that this oxidation leads to the generation of porosity due to the formation of micro-fractures allowing for fluid transport and subsequent dissolution of biotite and plagioclase. The rapid ensuing plagioclase weathering leads to formation of high secondary porosity in the corestone over a distance of only a few cm and eventually to the final disaggregation of bedrock to saprolite. The first secondary phases are oxides or amorphous precipitates from which secondary minerals (mainly gibbsite, kaolinite and goethite) form. As oxidation is the first weathering reaction, the supply of O-2 is a rate-limiting factor for chemical weathering. Hence, the supply of O-2 and its consumption at depth connects processes at the weathering front with those at the Earth's surface in a feedback mechanism. The strength of the feedback depends on the relative weight of advective versus diffusive transport of O-2 through the weathering profile. The feedback will be stronger with dominating diffusive transport. The low weathering rate is explained by the nature of this feedback that is ultimately dependent on the transport of O-2 through the whole regolith, and on lithological factors such as low bedrock porosity and the amount of Fe-bearing primary minerals. Tectonic quiescence in this region and low pre-development erosion rate (attributed to a dense vegetation cover) minimize the rejuvenation of the thick and cohesive regolith column, finally leading to low denudation rates. (C) 2015 Elsevier B.V. All rights reserved.
Much progress has been made in estimating recurrence intervals of great and giant subduction earthquakes using terrestrial, lacustrine, and marine paleoseismic archives. Recent detailed records suggest these earthquakes may have variable recurrence periods and magnitudes forming supercycles. Understanding seismic supercycles requires long paleoseismic archives that record timing and magnitude of such events. Turbidite paleoseismic archives may potentially extend past earthquake records to the Pleistocene and can thus complement commonly shorter-term terrestrial archives. However, in order to unambiguously establish recurring seismicity as a trigger mechanism for turbidity currents, synchronous deposition of turbidites in widely spaced, isolated depocenters has to be ascertained. Furthermore, characteristics that predispose a seismically active continental margin to turbidite paleoseismology and the correct sample site selection have to be taken into account.
Here we analyze 8 marine sediment cores along 950 km of the Chile margin to test for the feasibility of compiling detailed and continuous paleoseismic records based on turbidites. Our results suggest that the deposition of areally widespread, synchronous turbidites triggered by seismicity is largely controlled by sediment supply and, hence, the climatic and geomorphic conditions of the adjacent subaerial setting. The feasibility of compiling a turbidite paleoseismic record depends on the delicate balance between sufficient sediment supply providing material to fail frequently during seismic shaking and sufficiently low sedimentation rates to allow for coeval accumulation of planktonic foraminifera for high-resolution radiocarbon dating.
We conclude that offshore northern central Chile (29-32.5 degrees S) Holocene turbidite paleoseismology is not feasible, because sediment supply from the semi-arid mainland is low and almost no Holocene turbidity-current deposits are found in the cores. In contrast, in the humid region between 36 and 38 degrees S frequent Holocene turbidite deposition may generally correspond to paleoseismic events. However, high terrigenous sedimentation rates prevent high-resolution radiocarbon dating. The climatic transition region between 32.5 and 36 degrees S appears to be best suited for turbidite paleoseismology. (C) 2015 Elsevier Ltd. All rights reserved.
Newly acquired high-resolution bathymetric data (with 5 m and 2 m grid sizes) from the continental shelf off Concepcion (Chile), in combination with seismic reflection profiles, reveal a distinctly different evolution for the Biobio submarine canyon compared to that of one of its tributaries. Both canyons are incised into the shelf of the active margin. Whereas the inner shelf appears to be mantled with unconsolidated sediment, the outer shelf shows the influence of strong bottom currents that form drifts of loose sediment and transport -material into the Biobio submarine canyon and onto the continental slope.
The main stem of the Biobio Canyon is connected to the mouth of the Biobio River and currently provides a conduit for terrestrial sediment from the continental shelf to the deep seafloor. In contrast, the head of its tributary closest to the coast is located similar to 24 km offshore of the present-day coastline at 120 m water depth, and it is subject to passive sedimentation. However, canyon activity within the study area is interpreted to be controlled not only by the direct input of fluvial sediments into the canyon head facilitated by the river-mouth to canyon-head connection, but also by input from southward-directed bottom currents and possibly longshore drift. In addition, about 24 km offshore of the present-day coastline, the main stem of the Biobio Canyon has steep canyon walls next to sites of active tectonic deformation that are prone to wall failure. Mass-failure events may also foster turbidity currents and contribute to canyon feeding. In contrast, the tributary has less steep canyon walls with limited evidence of canyon-wall failure and is located down-system of bottom currents from the Biobio Canyon. It consequently receives neither fluvial nor longshore sediments. Therefore, the canyon's connectivity to fluvial or longshore sediment delivery pathways is affected by the distance of the canyon head from the coastline and the orientation of the canyon axis relative to the direction of bottom currents.
The ability of a submarine canyon to act as an active conduit for large quantities of terrestrial sediment toward the deep sea during sea-level highstands may be controlled by several different conditions simultaneously. These include bottom current direction, structural deformation of the seafloor affecting canyon location and orientation as well as canyon-wall failure, shelf gradient and associated distance from the canyon head to the coast, and fluvial networks. The complex interplay between these factors may vary even within an individual canyon system, resulting in distinct levels of canyon activity on a regional scale.
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.
Based on high-spatiotemporal-resolution data, the authors perform a climatological study of strong rainfall events propagating from southeastern South America to the eastern slopes of the central Andes during the monsoon season. These events account for up to 70% of total seasonal rainfall in these areas. They are of societal relevance because of associated natural hazards in the form of floods and landslides, and they form an intriguing climatic phenomenon, because they propagate against the direction of the low-level moisture flow from the tropics. The responsible synoptic mechanism is analyzed using suitable composites of the relevant atmospheric variables with high temporal resolution. The results suggest that the low-level inflow from the tropics, while important for maintaining sufficient moisture in the area of rainfall, does not initiate the formation of rainfall clusters. Instead, alternating low and high pressure anomalies in midlatitudes, which are associated with an eastward-moving Rossby wave train, in combination with the northwestern Argentinean low, create favorable pressure and wind conditions for frontogenesis and subsequent precipitation events propagating from southeastern South America toward the Bolivian Andes.
Extreme Rainfall of the South American Monsoon System: A Dataset Comparison Using Complex Networks
(2015)
In this study, the authors compare six different rainfall datasets for South America with a focus on their representation of extreme rainfall during the monsoon season (December February): the gauge-calibrated TRMM 3B42 V7 satellite product; the near-real-time TRMM 3B42 V7 RT, the GPCP 1 degrees daily (1DD) V1.2 satellite gauge combination product, the Interim ECMWF Re-Analysis (ERA-Interim) product; output of a high-spatial-resolution run of the ECHAM6 global circulation model; and output of the regional climate model Eta. For the latter three, this study can be understood as a model evaluation. In addition to statistical values of local rainfall distributions, the authors focus on the spatial characteristics of extreme rainfall covariability. Since traditional approaches based on principal component analysis are not applicable in the context of extreme events, they apply and further develop methods based on complex network theory. This way, the authors uncover substantial differences in extreme rainfall patterns between the different datasets: (i) The three model-derived datasets yield very different results than the satellite gauge combinations regarding the main climatological propagation pathways of extreme events as well as the main convergence zones of the monsoon system. (ii) Large discrepancies are found for the development of mesoscale convective systems in southeastern South America. (iii) Both TRMM datasets and ECHAM6 indicate a linkage of extreme rainfall events between the central Amazon basin and the eastern slopes of the central Andes, but this pattern is not reproduced by the remaining datasets. The authors' study suggests that none of the three model-derived datasets adequately captures extreme rainfall patterns in South America.
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.
Adjustment of empirically derived ground motion prediction equations (GMPEs), from a data- rich region/site where they have been derived to a data-poor region/site, is one of the major challenges associated with the current practice of seismic hazard analysis. Due to the fre- quent use in engineering design practices the GMPEs are often derived for response spectral ordinates (e.g., spectral acceleration) of a single degree of freedom (SDOF) oscillator. The functional forms of such GMPEs are based upon the concepts borrowed from the Fourier spectral representation of ground motion. This assumption regarding the validity of Fourier spectral concepts in the response spectral domain can lead to consequences which cannot be explained physically.
In this thesis, firstly results from an investigation that explores the relationship between Fourier and response spectra, and implications of this relationship on the adjustment issues of GMPEs, are presented. The relationship between the Fourier and response spectra is explored by using random vibration theory (RVT), a framework that has been extensively used in earthquake engineering, for instance within the stochastic simulation framework and in the site response analysis. For a 5% damped SDOF oscillator the RVT perspective of response spectra reveals that no one-to-one correspondence exists between Fourier and response spectral ordinates except in a limited range (i.e., below the peak of the response spectra) of oscillator frequencies. The high oscillator frequency response spectral ordinates are dominated by the contributions from the Fourier spectral ordinates that correspond to the frequencies well below a selected oscillator frequency. The peak ground acceleration (PGA) is found to be related with the integral over the entire Fourier spectrum of ground motion which is in contrast to the popularly held perception that PGA is a high-frequency phenomenon of ground motion.
This thesis presents a new perspective for developing a response spectral GMPE that takes the relationship between Fourier and response spectra into account. Essentially, this frame- work involves a two-step method for deriving a response spectral GMPE: in the first step two empirical models for the FAS and for a predetermined estimate of duration of ground motion are derived, in the next step, predictions from the two models are combined within the same RVT framework to obtain the response spectral ordinates. In addition to that, a stochastic model based scheme for extrapolating the individual acceleration spectra beyond the useable frequency limits is also presented. To that end, recorded acceleration traces were inverted to obtain the stochastic model parameters that allow making consistent extrapola- tion in individual (acceleration) Fourier spectra. Moreover an empirical model, for a dura- tion measure that is consistent within the RVT framework, is derived. As a next step, an oscillator-frequency-dependent empirical duration model is derived that allows obtaining the most reliable estimates of response spectral ordinates. The framework of deriving the response spectral GMPE presented herein becomes a self-adjusting model with the inclusion of stress parameter (∆σ) and kappa (κ0) as the predictor variables in the two empirical models. The entire analysis of developing the response spectral GMPE is performed on recently compiled RESORCE-2012 database that contains recordings made from Europe, the Mediterranean and the Middle East. The presented GMPE for response spectral ordinates should be considered valid in the magnitude range of 4 ≤ MW ≤ 7.6 at distances ≤ 200 km.
Empirical ground-motion prediction equations (GMPEs) require adjustment to make them appropriate for site-specific scenarios. However, the process of making such adjustments remains a challenge. This article presents a holistic framework for the development of a response spectral GMPE that is easily adjustable to different seismological conditions and does not suffer from the practical problems associated with adjustments in the response spectral domain. The approach for developing a response spectral GMPE is unique, because it combines the predictions of empirical models for the two model components that characterize the spectral and temporal behavior of the ground motion. Essentially, as described in its initial form by Bora et al. (2014), the approach consists of an empirical model for the Fourier amplitude spectrum (FAS) and a model for the ground-motion duration. These two components are combined within the random vibration theory framework to obtain predictions of response spectral ordinates. In addition, FAS corresponding to individual acceleration records are extrapolated beyond the useable frequencies using the stochastic FAS model, obtained by inversion as described in Edwards and Fah (2013a). To that end, a (oscillator) frequency-dependent duration model, consistent with the empirical FAS model, is also derived. This makes it possible to generate a response spectral model that is easily adjustable to different sets of seismological parameters, such as the stress parameter Delta sigma, quality factor Q, and kappa kappa(0). The dataset used in Bora et al. (2014), a subset of the RESORCE-2012 database, is considered for the present analysis. Based upon the range of the predictor variables in the selected dataset, the present response spectral GMPE should be considered applicable over the magnitude range of 4 <= M-w <= 7.6 at distances <= 200 km.
Spectral fingerprinting
(2015)
Current research on runoff and erosion processes, as well as an increasing demand for sustainable watershed management emphasize the need for an improved understanding of sediment dynamics. This involves the accurate assessment of erosion rates and sediment transfer, yield and origin. A variety of methods exist to capture these processes at the catchment scale. Among these, sediment fingerprinting, a technique to trace back the origin of sediment, has attracted increasing attention by the scientific community in recent years. It is a two-step procedure, based on the fundamental assumptions that potential sources of sediment can be reliably discriminated based on a set of characteristic ‘fingerprint’ properties, and that a comparison of source and sediment fingerprints allows to quantify the relative contribution of each source.
This thesis aims at further assessing the potential of spectroscopy to assist and improve the sediment fingerprinting technique. Specifically, this work focuses on (1) whether potential sediment sources can be reliably identified based on spectral features (‘fingerprints’), whether (2) these spectral fingerprints permit the quantification of relative source contribution, and whether (3) in situ derived source information is sufficient for this purpose. Furthermore, sediment fingerprinting using spectral information is applied in a study catchment to (4) identify major sources and observe how relative source contributions change between and within individual flood events. And finally, (5) spectral fingerprinting results are compared and combined with simultaneous sediment flux measurements to study sediment origin, transport and storage behaviour.
For the sediment fingerprinting approach, soil samples were collected from potential sediment sources within the Isábena catchment, a meso-scale basin in the central Spanish Pyrenees. Undisturbed samples of the upper soil layer were measured in situ using an ASD spectroradiometer and subsequently sampled for measurements in the laboratory. Suspended sediment was sampled automatically by means of ISCO samplers at the catchment as well as at the five major subcatchment outlets during flood events, and stored fine sediment from the channel bed was collected from 14 cross-sections along the main river. Artificial mixtures of known contributions were produced from source soil samples. Then, all source, sediment and mixture samples were dried and spectrally measured in the laboratory. Subsequently, colour coefficients and physically based features with relation to organic carbon, iron oxide, clay content and carbonate, were calculated from all in situ and laboratory spectra. Spectral parameters passing a number of prerequisite tests were submitted to principal component analyses to study natural clustering of samples, discriminant function analyses to observe source differentiation accuracy, and a mixing model for source contribution assessment. In addition, annual as well as flood event based suspended sediment fluxes from the catchment and its subcatchments were calculated from rainfall, water discharge and suspended sediment concentration measurements using rating curves and Quantile Regression Forests. Results of sediment flux monitoring were interpreted individually with respect to storage behaviour, compared to fingerprinting source ascriptions and combined with fingerprinting to assess their joint explanatory potential.
In response to the key questions of this work, (1) three source types (land use) and five spatial sources (subcatchments) could be reliably discriminated based on spectral fingerprints. The artificial mixture experiment revealed that while (2) laboratory parameters permitted source contribution assessment, (3) the use of in situ derived information was insufficient. Apparently, high discrimination accuracy does not necessarily imply good quantification results. When applied to suspended sediment samples of the catchment outlet, the spectral fingerprinting approach was able to (4) quantify the major sediment sources: badlands and the Villacarli subcatchment, respectively, were identified as main contributors, which is consistent with field observations and previous studies. Thereby, source contribution was found to vary both, within and between individual flood events. Also sediment flux was found to vary considerably, annually as well as seasonally and on flood event base. Storage was confirmed to play an important role in the sediment dynamics of the studied catchment, whereas floods with lower total sediment yield tend to deposit and floods with higher yield rather remove material from the channel bed. Finally, a comparison of flux measurements with fingerprinting results highlighted the fact that (5) immediate transport from sources to the catchment outlet cannot be assumed. A combination of the two methods revealed different aspects of sediment dynamics that none of the techniques could have uncovered individually.
In summary, spectral properties provide a fast, non-destructive, and cost-efficient means to discriminate and quantify sediment sources, whereas, unfortunately, straight-forward in situ collected source information is insufficient for the approach. Mixture modelling using artificial mixtures permits valuable insights into the capabilities and limitations of the method and similar experiments are strongly recommended to be performed in the future. Furthermore, a combination of techniques such as e.g. (spectral) sediment fingerprinting and sediment flux monitoring can provide comprehensive understanding of sediment dynamics.
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.
Injection of nanoscale zero-valent iron (nZVI) is an innovative technology for in situ installation of a permeable reactive barrier in the subsurface. Zerovalent iron (ZVI) is highly reactive with chlorinated hydrocarbons (CHCs) and renders them into less harmful substances. Application of nZVI instead of granular ZVI can increase rates of dechlorination of CHCs by orders of magnitude, due to its higher surface area. This approach is still difficult to apply due to fast agglomeration and sedimentation of colloidal suspensions of nZVI, which leads to very short transport distances. To overcome this issue of limited mobility, polyanionic stabilisers are added to increase surface charge and stability of suspensions. In field experiments maximum transport distances of a few metres were achieved. A new approach, which is investigated in this thesis, is enhanced mobility of nZVI by a more mobile carrier colloid. The investigated composite material consists of activated carbon, which is loaded with nZVI.
In this cumulative thesis, transport characteristics of carbon-colloid supported nZVI (c-nZVI) are investigated. Investigations started with column experiments in 40 cm columns filled with various porous media to investigate on physicochemical influences on transport characteristics. The experimental setup was enlarged to a transport experiment in a 1.2-m-sized two-dimensional aquifer tank experiment, which was filled with granular porous media. Further, a field experiment was performed in a natural aquifer system with a targeted transport distance of 5.3 m. Parallel to these investigations, alternative methods for transport observations were investigated by using noninvasive tomographic methods. Experiments using synchrotron radiation and magnetic resonance (MRI) were performed to investigate in situ transport characteristics in a non-destructive way.
Results from column experiments show potentially high mobility under environmental relevant conditions. Addition of mono-and bivalent salts, e.g. more than 0.5 mM/L CaCl2, might decrease mobility. Changes in pH to values below 6 can inhibit mobility at all. Measurements of colloid size show changes in the mean particle size by a factor of ten. Measurements of zeta potential revealed an increase of –62 mV to –82 mV. Results from the 2D-aquifer test system suggest strong particle deposition in the first centimetres and only weak straining in the further travel path and no gravitational influence on particle transport. Straining at the beginning of the travel path in the porous medium was observed with tomographic investigations of transport. MRI experiments revealed similar results to the previous experiments, and observations using synchrotron radiation suggest straining of colloids at pore throats. The potential for high transport distances, which was suggested from laboratory experiments, was confirmed in the field experiment, where the transport distance of 5.3 m was reached by at least 10% of injected nZVI. Altogether, transport distances of the investigated carbon-colloid supported nZVI are higher than published results of traditional nZVI.
The application of nanoscale zero-valent iron (nZVI) for subsurface remediation of groundwater contaminants is a promising new technology, which can be understood as alternative to the permeable reactive barrier technique using granular iron. Dechlorination of organic contaminants by zero-valent iron seems promising. Currently, one limitation to widespread deployment is the fast agglomeration and sedimentation of nZVI in colloidal suspensions, even more so when in soils and sediments, which limits the applicability for the treatment of sources and plumes of contamination. Colloid-supported nZVI shows promising characteristics to overcome these limitations. Mobility of Carbo-Iron Colloids (CIC) - a newly developed composite material based on finely ground activated carbon as a carrier for nZVI - was tested in a field application: In this study, a horizontal dipole flow field was established between two wells separated by 53 m in a confined, natural aquifer. The injection/extraction rate was 500 L/h. Approximately 12 kg of CIC was suspended with the polyanionic stabilizer carboxymethyl cellulose. The suspension was introduced into the aquifer at the injection well. Breakthrough of CIC was observed visually and based on total particle and iron concentrations detected in samples from the extraction well. Filtration of water samples revealed a particle breakthrough of about 12% of the amount introduced. This demonstrates high mobility of CIC particles and we suggest that nZVI carried on CIC can be used for contaminant plume remediation by in-situ formation of reactive barriers. (C) 2015 Elsevier B.V. All rights reserved.
Endosulfan - an agricultural insecticide and banned by Stockholm Convention - is produced as a 2:1 to 7:3 mixture of isomers endosulfan I (ESI) and endosulfan II (ESII). Endosulfan is transformed under aerobic conditions into endosulfan sulfate (ESS). The study shows for 76 sampling locations in German forests that endosulfan is abundant in all samples with an opposite ratio between the ESI and ESII than the technical product, where the main metabolite ESS is found with even higher abundance. The ratio between ESI/ESII and ESS show clear dependence on the type of stands (coniferous vs. deciduous) and humus type and increases from deciduous via mixed to coniferous forest stands. The study argues for a systematic monitoring of ESI, ESII, and ESS and underlines the need for further research, specifically on the fate of endosulfan including biomagnifications and bioaccumulation in soil. (C) 2015 Elsevier Ltd. All rights reserved.
The continuously increasing demand for rare earth elements in technical components of modern technologies, brings the detection of new deposits closer into the focus of global exploration. One promising method to globally map important deposits might be remote sensing, since it has been used for a wide range of mineral mapping in the past. This doctoral thesis investigates the capacity of hyperspectral remote sensing for the detection of rare earth element deposits. The definition and the realization of a fundamental database on the spectral characteristics of rare earth oxides, rare earth metals and rare earth element bearing materials formed the basis of this thesis. To investigate these characteristics in the field, hyperspectral images of four outcrops in Fen Complex, Norway, were collected in the near-field. A new methodology (named REEMAP) was developed to delineate rare earth element enriched zones. The main steps of REEMAP are: 1) multitemporal weighted averaging of multiple images covering the sample area; 2) sharpening the rare earth related signals using a Gaussian high pass deconvolution technique that is calibrated on the standard deviation of a Gaussian-bell shaped curve that represents by the full width of half maxima of the target absorption band; 3) mathematical modeling of the target absorption band and highlighting of rare earth elements. REEMAP was further adapted to different hyperspectral sensors (EO-1 Hyperion and EnMAP) and a new test site (Lofdal, Namibia). Additionally, the hyperspectral signatures of associated minerals were investigated to serve as proxy for the host rocks. Finally, the capacity and limitations of spectroscopic rare earth element detection approaches in general and of the REEMAP approach specifically were investigated and discussed. One result of this doctoral thesis is that eight rare earth oxides show robust absorption bands and, therefore, can be used for hyperspectral detection methods. Additionally, the spectral signatures of iron oxides, iron-bearing sulfates, calcite and kaolinite can be used to detect metasomatic alteration zones and highlight the ore zone. One of the key results of this doctoral work is the developed REEMAP approach, which can be applied from near-field to space. The REEMAP approach enables rare earth element mapping especially for noisy images. Limiting factors are a low signal to noise ratio, a reduced spectral resolution, overlaying materials, atmospheric absorption residuals and non-optimal illumination conditions. Another key result of this doctoral thesis is the finding that the future hyperspectral EnMAP satellite (with its currently published specifications, June 2015) will be theoretically capable to detect absorption bands of erbium, dysprosium, holmium, neodymium and europium, thulium and samarium. This thesis presents a new methodology REEMAP that enables a spatially wide and rapid hyperspectral detection of rare earth elements in order to meet the demand for fast, extensive and efficient rare earth exploration (from near-field to space).
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.
This study investigates the spatial and temporal distributions of 14 key arboreal taxa and their driving forces during the last 22,000 calendar years before ad 1950 (kyr BP) using a taxonomically harmonized and temporally standardized fossil pollen dataset with a 500-year resolution from the eastern part of continental Asia. Logistic regression was used to estimate pollen abundance thresholds for vegetation occurrence (presence or dominance), based on modern pollen data and present ranges of 14 taxa in China. Our investigation reveals marked changes in spatial and temporal distributions of the major arboreal taxa. The thermophilous (Castanea, Castanopsis, Cyclobalanopsis, Fagus, Pterocarya) and eurythermal (Juglans, Quercus, Tilia, Ulmus) broadleaved tree taxa were restricted to the current tropical or subtropical areas of China during the Last Glacial Maximum (LGM) and spread northward since c. 14.5kyr BP. Betula and conifer taxa (Abies, Picea, Pinus), in contrast, retained a wider distribution during the LGM and showed no distinct expansion direction during the Late Glacial. Since the late mid-Holocene, the abundance but not the spatial extent of most trees decreased. The changes in spatial and temporal distributions for the 14 taxa are a reflection of climate changes, in particular monsoonal moisture, and, in the late Holocene, human impact. The post-LGM expansion patterns in eastern continental China seem to be different from those reported for Europe and North America, for example, the westward spread for eurythermal broadleaved taxa.
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.
With less than two decades of activity, research on melt inclusions (MI) in crystals from rocks that have undergone crustal anatexis - migmatites and granulites - is a recent addition to crustal petrology and geochemistry. Studies on this subject started with glassy inclusions in anatectic crustal enclaves in lavas, and then progressed to regionally metamorphosed and partially melted crustal rocks, where melt inclusions are normally crystallized into a cryptocrystalline aggregate (nanogranitoid).
Since the first paper on melt inclusions in the granulites of the Kerala Khondalite Belt in 2009, reported and studied occurrences are already a few tens. Melt inclusions in migmatites and granulites show many analogies with their more common and long studied counterparts in igneous rocks, but also display very important differences and peculiarities, which are the subject of this review. Microstructurally, melt inclusions in anatectic rocks are small, commonly 10 mu m in diameter, and their main mineral host is peritectic garnet, although several other hosts have been observed. Inclusion contents vary from glass in enclaves that were cooled very rapidly from supersolidus temperatures, to completely crystallized material in slowly cooled regional migmatites. The chemical composition of the inclusions can be analyzed combining several techniques (SEM, EMP, NanoSIMS, LA-ICP-MS), but in the case of crystallized inclusions the experimental remelting under confining pressure in a piston cylinder is a prerequisite. The melt is generally granitic and peraluminous, although granodioritic to trondhjemitic compositions have also been found.
Being mostly primary in origin, inclusions attest for the growth of their peritectic host in the presence of melt. As a consequence, the inclusions have the unique ability of preserving information on the composition of primary anatectic crustal melts, before they undergo any of the common following changes in their way to produce crustal magmas. For these peculiar features, melt inclusions in migmatites and granulites, largely overlooked so far, have the potential to become a fundamental tool for the study of crustal melting, crustal differentiation, and even the generation of the continental crust. (C) 2015 The Authors. Published by Elsevier B.V.
In this study we use the anisotropy of magnetic susceptibility (AMS) and paleomagnetic data for deciphering the origin of magnetic lineation in weakly deformed sedimentary rocks and for evaluating oroclinal processes within the Arabia-Eurasia collision zone. In particular, we have analyzed the Miocene Upper Red Formation (URF) from the outer curved front of the southern Central Alborz Mountains of north Iran, to test for the first time with paleomagnetic data the origin (primary versus secondary) of this orogenic arc. AMS data document the existence of a magnetic lineation parallel to the orientation of the major tectonic structures, which vary along strike from WNW to ENE. These directions are highly oblique to the paleoflow directions and hence suggest that the magnetic lineation in the URF was produced by compressional deformation during layer-parallel shortening. In addition, our paleomagnetic data document clockwise and anticlockwise rotations along vertical axis for the western and eastern sectors of the Central Alborz Mountains, respectively. Combined, our results suggest that the orogen represents an orocline, which formed not earlier than circa 7.6Ma most likely through bending processes caused by the relative motion between the rigid crustal blocks of the collision zone. Moreover, our study provides new insights into the Iran-Southern Caspian Basin kinematic evolution suggesting that the present-day SW motion of the South Caspian Basin with respect to Central Iran postdates oroclinal bending and hence cannot be as old as late Miocene to early Pliocene but a rather recent configuration (i.e., 3 to <1Ma).
Modelling the transfer of supraglacial meltwater to the bed of Leverett Glacier, Southwest Greenland
(2015)
Meltwater delivered to the bed of the Greenland Ice Sheet is a driver of variable ice-motion through changes in effective pressure and enhanced basal lubrication. Ice surface velocities have been shown to respond rapidly both to meltwater production at the surface and to drainage of supraglacial lakes, suggesting efficient transfer of meltwater from the supraglacial to subglacial hydrological systems. Although considerable effort is currently being directed towards improved modelling of the controlling surface and basal processes, modelling the temporal and spatial evolution of the transfer of melt to the bed has received less attention. Here we present the results of spatially distributed modelling for prediction of moulins and lake drainages on the Leverett Glacier in Southwest Greenland. The model is run for the 2009 and 2010 ablation seasons, and for future increased melt scenarios. The temporal pattern of modelled lake drainages are qualitatively comparable with those documented from analyses of repeat satellite imagery. The modelled timings and locations of delivery of meltwater to the bed also match well with observed temporal and spatial patterns of ice surface speed-ups. This is particularly true for the lower catchment (< 1000 m a.s.l.) where both the model and observations indicate that the development of moulins is the main mechanism for the transfer of surface meltwater to the bed. At higher elevations (e.g. 1250-1500 m a.s.l.) the development and drainage of supraglacial lakes becomes increasingly important. At these higher elevations, the delay between modelled melt generation and subsequent delivery of melt to the bed matches the observed delay between the peak air temperatures and subsequent velocity speed-ups, while the instantaneous transfer of melt to the bed in a control simulation does not. Although both moulins and lake drainages are predicted to increase in number for future warmer climate scenarios, the lake drainages play an increasingly important role in both expanding the area over which melt accesses the bed and in enabling a greater proportion of surface melt to reach the bed.
A detailed architectural analysis has been performed on the porcelaneous species Taberina bingistani Henson, 1948 from diverse localities, including its type locality in Iran. Characterized by a subcylindrical to slightly conical shape, with a cribbed apertural face, the Cenomanian Middle-East species bingistani is removed from the American early-Paleogene genus Taberina Keijzer, 1945 (type species: T. cubana) and ascribed to the new genus Praetaberina. In the early stages of growth, the chambers are planispiral-involute, becoming uncoiled in the later stages. The interior of the chambers shows numerous septula and pillars. The septula are radial and distributed in the marginal area of the chamber lumen, while the pillars occupy the central area and are regularly distributed between the intercameral foramina. Pillars are in alignment from one chamber to the next. A new species, from southeast Italy, Praetaberina apula, is also erected. It differs from P. bingistani in its smaller size and significantly reduced structural elements. Both species characterize upper Cenomanian shallow-water deposits.
Earthquakes occurring close to hydrocarbon fields under production are often under critical view of being induced or triggered. However, clear and testable rules to discriminate the different events have rarely been developed and tested. The unresolved scientific problem may lead to lengthy public disputes with unpredictable impact on the local acceptance of the exploitation and field operations. We propose a quantitative approach to discriminate induced, triggered, and natural earthquakes, which is based on testable input parameters. Maxima of occurrence probabilities are compared for the cases under question, and a single probability of being triggered or induced is reported. The uncertainties of earthquake location and other input parameters are considered in terms of the integration over probability density functions. The probability that events have been human triggered/induced is derived from the modeling of Coulomb stress changes and a rate and state-dependent seismicity model. In our case a 3-D boundary element method has been adapted for the nuclei of strain approach to estimate the stress changes outside the reservoir, which are related to pore pressure changes in the field formation. The predicted rate of natural earthquakes is either derived from the background seismicity or, in case of rare events, from an estimate of the tectonic stress rate. Instrumentally derived seismological information on the event location, source mechanism, and the size of the rupture plane is of advantage for the method. If the rupture plane has been estimated, the discrimination between induced or only triggered events is theoretically possible if probability functions are convolved with a rupture fault filter. We apply the approach to three recent main shock events: (1) the M-w 4.3 Ekofisk 2001, North Sea, earthquake close to the Ekofisk oil field; (2) the M-w 4.4 Rotenburg 2004, Northern Germany, earthquake in the vicinity of the Sohlingen gas field; and (3) the M-w 6.1 Emilia 2012, Northern Italy, earthquake in the vicinity of a hydrocarbon reservoir. The three test cases cover the complete range of possible causes: clearly human induced, not even human triggered, and a third case in between both extremes.
The recently proposed global monsoon hypothesis interprets monsoon systems as part of one global-scale atmospheric overturning circulation, implying a connection between the regional monsoon systems and an in-phase behaviour of all northern hemispheric monsoons on annual timescales (Trenberth et al., 2000). Whether this concept can be applied to past climates and variability on longer timescales is still under debate, because the monsoon systems exhibit different regional characteristics such as different seasonality (i. e. onset, peak and withdrawal). To investigate the interconnection of different monsoon systems during the pre-industrial Holocene, five transient global climate model simulations have been analysed with respect to the rainfall trend and variability in different sub-domains of the Afro-Asian monsoon region. Our analysis suggests that on millennial timescales with varying orbital forcing, the monsoons do not behave as a tightly connected global system. According to the models, the Indian and North African monsoons are coupled, showing similar rainfall trend and moderate correlation in centennial rainfall variability in all models. The East Asian monsoon changes independently during the Holocene. The dissimilarities in the seasonality of the monsoon sub-systems lead to a stronger response of the North African and Indian monsoon systems to the Holocene insolation forcing than of the East Asian monsoon and affect the seasonal distribution of Holocene rainfall variations. Within the Indian and North African monsoon domain, precipitation solely changes during the summer months, showing a decreasing Holocene precipitation trend. In the East Asian monsoon region, the precipitation signal is determined by an increasing precipitation trend during spring and a decreasing precipitation change during summer, partly balancing each other. A synthesis of reconstructions and the model results do not reveal an impact of the different seasonality on the timing of the Holocene rainfall optimum in the different sub-monsoon systems. Rather they indicate locally inhomogeneous rainfall changes and show that single palaeo-records should not be used to characterise the rainfall change and monsoon evolution for entire monsoon sub-systems.
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.
We have studied the focal mechanism of an M-w 4.0 earthquake that occurred on 8 October 2011 in the southwest of El Hierro (Canary Islands), the largest shock of the swarm that preceded the submarine eruption of El Hierro 2011-2012. The joint focal mechanism solution of 34 foreshocks has also been obtained. The results for the mainshock are consistent with a pure double-couple mechanism of a strike-slip motion with a reverse component and a focal depth of 12-13 km. The stress pattern obtained from the focal mechanism indicates horizontal compression in a north-northwest-south-southeast direction, parallel to the southern ridge of the island, and a quasi-horizontal extension in an east-west direction. Similar stress pattern is derived from the joint solution of the foreshocks. The occurrence of this family of earthquakes at the moment of the maximum strain rate of the pre-eruptive swarm suggests that their rupture process is related to tectonic stress, which led to the eruption only two days later, 5 km away from the mainshock epicenter.
On 11 August 2012 an earthquake doublet (M-w 6.4 and 6.2) occurred near the city of Ahar, northwest Iran. Both events were only 6 km and 11 minutes apart, producing a surface rupture of about 12 km in length. Historical and modern seismicity has so far been sparse in this area. Spatially, the region represents a transitional zone between different tectonic domains, including compression in Iran, westward extrusion of the Anatolian plate, and thrusting beneath the Caucasus.
In this study, we inverted the surface waveforms of the two mainshocks and 11 aftershocks (M-w >= 4.3) to obtain regional seismic moment tensors. The earthquakes analyzed can be grouped into pure strike slip (including the first mainshock) and oblique reverse mechanisms (including the second mainshock). The sequence provides information about faulting mechanisms at the spatial scale of the entire rock volume affected by the earthquake doublet, including coinciding deformation on minor faults (sub) parallel to the main fault and Riedel shears. It occurred on a so far unknown fault structure, which we call the Ahar fault.
Alongside the seismological data, we used geological maps, satellite images, and digital elevation data to analyze the geomorphology of the region. Our analysis suggests that the adjacent North Tabriz fault, which accomodates up to 7 mm/yr of right-lateral strike-slip faulting, does not compensate the entire lateral shear strain, and that part of it is compensated farther north. Combined, our results suggest a temporally and spatially complex style of deformation (reverse and strike slip) overprinting older reverse deformation.
Observed streamflow of headwater catchments of the Tarim River (Central Asia) increased by about 30% over the period 1957-2004. This study aims at assessing to which extent these streamflow trends can be attributed to changes in air temperature or precipitation. The analysis includes a data-based approach using multiple linear regression and a simulation-based approach using a hydrological model. The hydrological model considers changes in both glacier area and surface elevation. It was calibrated using a multiobjective optimization algorithm with calibration criteria based on glacier mass balance and daily and interannual variations of discharge. The individual contributions to the overall streamflow trends from changes in glacier geometry, temperature, and precipitation were assessed using simulation experiments with a constant glacier geometry and with detrended temperature and precipitation time series. The results showed that the observed changes in streamflow were consistent with the changes in temperature and precipitation. In the Sari-Djaz catchment, increasing temperatures and related increase of glacier melt were identified as the dominant driver, while in the Kakshaal catchment, both increasing temperatures and increasing precipitation played a major role. Comparing the two approaches, an advantage of the simulation-based approach is the fact that it is based on process-based relationships implemented in the hydrological model instead of statistical links in the regression model. However, data-based approaches are less affected by model parameter and structural uncertainties and typically fast to apply. A complementary application of both approaches is recommended.
Semi-automated detection of annual laminae (varves) in lake sediments using a fuzzy logic algorithm
(2015)
Annual laminae (varves) in lake sediments are typically visually identified, measured and counted, although numerous attempts have been made to automate this process. The reason for the failure of most of these automated algorithms for varve counting is the complexity of the seasonal laminations, typically rich in lateral fades variations and internal heterogeneities. In the manual counting of varves, the investigator acquired and interpreted flexible numbers of complex decision criteria to understand whether a particular simple lamination is a varve or not. Fuzzy systems simulate the flexible decision making process in a computer by introducing a smooth transition between true varve and false varve. In our investigation, we use an adaptive neuro fuzzy inference system (ANFIS) to detect varves on the basis of a digital image of the sediment. The results of the application of the ANFIS to laminated sediments from the Meerfelder Maar (Eifel, Germany) and from a landslide-dammed lake in the Quebrada de Cafayate of Argentina are compared with manual varve counts and possible reasons for the differences are discussed. (C) 2015 Elsevier B.V. All rights reserved.
The periglacial engine of mountain erosion - Part 2: Modelling large-scale landscape evolution
(2015)
There is growing recognition of strong periglacial control on bedrock erosion in mountain landscapes, including the shaping of low-relief surfaces at high elevations (summit flats). But, as yet, the hypothesis that frost action was crucial to the assumed Late Cenozoic rise in erosion rates remains compelling and untested. Here we present a landscape evolution model incorporating two key periglacial processes - regolith production via frost cracking and sediment transport via frost creep - which together are harnessed to variations in temperature and the evolving thickness of sediment cover. Our computational experiments time-integrate the contribution of frost action to shaping mountain topography over million-year timescales, with the primary and highly reproducible outcome being the development of flattish or gently convex summit flats. A simple scaling of temperature to marine delta O-18 records spanning the past 14 Myr indicates that the highest summit flats in mid-to high-latitude mountains may have formed via frost action prior to the Quaternary. We suggest that deep cooling in the Quaternary accelerated mechanical weathering globally by significantly expanding the area subject to frost. Further, the inclusion of subglacial erosion alongside periglacial processes in our computational experiments points to alpine glaciers increasing the long-term efficiency of frost-driven erosion by steepening hillslopes.
The Upper Devonian Munster Basin of southern Ireland has traditionally been viewed as a post-orogenic molasse deposit that was sourced from the Caledonides of central Ireland and subsequently deformed by the end Carboniferous Variscan orogenic event. The basin fill is composed of super-mature quartz arenite sandstone that clearly represents a second cycle of deposition. The source of this detritus is now recognized as Lower Devonian Dingle Basin red bed sequences to the north. This genetic link is based on the degree of similarity in the detrital mica chemistry in both of these units; micas plot in identical fields and define the same trends. In addition, the two sequences show increased textural and chemical maturity up-sequence and define indistinguishable Ar-40/Ar-39 age ranges for the detrital mica grains. Partial resetting of the Ar ages can be attributed to elevated heat flow in the region caused by Munster Basin extension and subsequent Variscan deformation. The combined evidence from southwest Ireland therefore points to a Caledonian or possibly Taconian primary source area that initially shed detritus into the Lower Devonian Dingle Basin which was subsequently recycled into the Upper Devonian Munster Basin following mid-Devonian Acadian basin inversion. (C) 2014 International Association for Gondwana Research. Published by Elsevier B.V. All rights reserved.
Hydrology is rich in methods that use information theory to evaluate monitoring networks. Yet in most existing studies, only the available data set as a whole is used, which neglects the intraannual variability of the hydrological system. In this paper, we demonstrate how this variability can be considered by extending monitoring evaluation to subsets of the available data. Therefore, we separately evaluated time windows of fixed length, which were shifted through the data set, and successively extended time windows. We used basic information theory measures and a greedy ranking algorithm based on the criterion of maximum information/minimum redundancy. The network investigated monitored surface and groundwater levels at quarter-hourly intervals and was located at an artificially drained lowland site in the Spreewald region in north-east Germany. The results revealed that some of the monitoring stations were of value permanently while others were needed only temporally. The prevailing meteorological conditions, particularly the amount of precipitation, affected the degree of similarity between the water levels measured. The hydrological system tended to act more individually during periods of no or little rainfall. The optimal monitoring setup, its stability, and the monitoring effort necessary were influenced by the meteorological forcing. Altogether, the methodology presented can help achieve a monitoring network design that has a more even performance or covers the conditions of interest (e.g., floods or droughts) best.
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.