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.
Our understanding of the impact of climate-driven environmental change on prehistoric human populations is hampered by the scarcity of continuous paleoenvironmental records in the vicinity of archaeological sites. Here we compare a continuous paleoclimatic record of the last 20 ka before present from the Chew Bahir basin, southwest Ethiopia, with the available archaeological record of human presence in the region. The correlation of this record with orbitally-driven insolation variations suggests a complex nonlinear response of the environment to climate forcing, reflected in several long-term and short-term transitions between wet and dry conditions, resulting in abrupt changes between favorable and unfavorable living conditions for humans. Correlating the archaeological record in the surrounding region of the Chew Bahir basin, presumably including montane and lake-marginal refugia for human populations, with our climate record suggests a complex interplay between humans and their environment during the last 20 ka. The result may contribute to our understanding of how a dynamic environment may have impacted the adaptation and dispersal of early humans in eastern Africa. (C) 2015 Elsevier Ltd. All rights reserved.
Extreme droughts, heat waves, frosts, precipitation, wind storms and other climate extremes may impact the structure, composition and functioning of terrestrial ecosystems, and thus carbon cycling and its feedbacks to the climate system. Yet, the interconnected avenues through which climate extremes drive ecological and physiological processes and alter the carbon balance are poorly understood. Here, we review the literature on carbon cycle relevant responses of ecosystems to extreme climatic events. Given that impacts of climate extremes are considered disturbances, we assume the respective general disturbance-induced mechanisms and processes to also operate in an extreme context. The paucity of well-defined studies currently renders a quantitative meta-analysis impossible, but permits us to develop a deductive framework for identifying the main mechanisms (and coupling thereof) through which climate extremes may act on the carbon cycle. We find that ecosystem responses can exceed the duration of the climate impacts via lagged effects on the carbon cycle. The expected regional impacts of future climate extremes will depend on changes in the probability and severity of their occurrence, on the compound effects and timing of different climate extremes, and on the vulnerability of each land-cover type modulated by management. Although processes and sensitivities differ among biomes, based on expert opinion, we expect forests to exhibit the largest net effect of extremes due to their large carbon pools and fluxes, potentially large indirect and lagged impacts, and long recovery time to regain previous stocks. At the global scale, we presume that droughts have the strongest and most widespread effects on terrestrial carbon cycling. Comparing impacts of climate extremes identified via remote sensing vs. ground-based observational case studies reveals that many regions in the (sub-)tropics are understudied. Hence, regional investigations are needed to allow a global upscaling of the impacts of climate extremes on global carbon-climate feedbacks.
Shallow-water carbonates are invaluable archives of past global change. They hold the record of how neritic biologic communities reacted to palaeoenvironmental changes. However, attempts to decipher these geological archives are often severely hampered by the low stratigraphic resolution attained by biostratigraphy. This is particularly the case for the Upper Cretaceous carbonate platforms of the central Tethyan realm: their biostratigraphy suffers from very low resolution and poor correlation with the standard biochronologic scales based on ammonites, planktic foraminifers and calcareous nannoplankton.
In this paper we show how this problem can be tackled by integrating biostratigraphy with isotope stratigraphy. We present a detailed record of the benthic foraminiferal biostratigraphy and carbon and strontium isotope stratigraphy of three upper Cenomanian-middle Campanian sections belonging to the Apennine Carbonate Platform of southern Italy. For the upper Cenomanian-Turonian interval, the carbon isotope curves of the studied sections are easily correlated to the reference curve of the English Chalk. The correlation is facilitated by the matching of the prominent positive excursion corresponding to the Oceanic Anoxic Event 2. For the Coniacian-middle Campanian interval, the correlation is mainly based on strontium isotope stratigraphy. We use the Sr-87/Sr-86 ratios of the low-Mg calcite of well preserved rudist shells to obtain accurate chronostratigraphic ages for many levels of the three studied sections. The ages obtained by Sr isotope stratigraphy are then used to better constrain the matching of the carbon isotope curves.
From the high-resolution chronostratigraphic age-model stablished by isotope stratigraphy, we derive the chronostratigraphic calibration of benthic foraminiferal biostratigraphic events. For the first time the benthic foraminiferal biozones of the Apennine Carbonate Platform can be accurately correlated to the standard ammonite biozonation. This result is of great relevance because the biostratigraphic schemes of other carbonate platforms in the central and southern Tethyan realm are largely based on the same biostratigraphic events. (C) 2014 Elsevier Ltd. All rights reserved.
Accelerating climate change and increased economic and environmental interests in permafrost-affected regions have resulted in an acute need for more directed permafrost research. In June 2014, 88 early career researchers convened to identify future priorities for permafrost research. This multidisciplinary forum concluded that five research topics deserve greatest attention: permafrost landscape dynamics, permafrost thermal modeling, integration of traditional knowledge, spatial distribution of ground ice, and engineering issues. These topics underline the need for integrated research across a spectrum of permafrost-related domains and constitute a contribution to the Third International Conference on Arctic Research Planning (ICARP III).
Thermal permafrost degradation and coastal erosion in the Arctic remobilize substantial amounts of organic carbon (OC) and nutrients which have accumulated in late Pleistocene and Holocene unconsolidated deposits. Permafrost vulnerability to thaw subsidence, collapsing coastlines and irreversible landscape change are largely due to the presence of large amounts of massive ground ice such as ice wedges. However, ground ice has not, until now, been considered to be a source of dissolved organic carbon (DOC), dissolved inorganic carbon (DIC) and other elements which are important for ecosystems and carbon cycling. Here we show, using biogeochemical data from a large number of different ice bodies throughout the Arctic, that ice wedges have the greatest potential for DOC storage, with a maximum of 28.6 mg L-1 (mean: 9.6 mg L-1). Variation in DOC concentration is positively correlated with and explained by the concentrations and relative amounts of typically terrestrial cations such as Mg2+ and K+. DOC sequestration into ground ice was more effective during the late Pleistocene than during the Holocene, which can be explained by rapid sediment and OC accumulation, the prevalence of more easily degradable vegetation and immediate incorporation into permafrost. We assume that pristine snowmelt is able to leach considerable amounts of well-preserved and highly bioavailable DOC as well as other elements from surface sediments, which are rapidly frozen and stored in ground ice, especially in ice wedges, even before further degradation. We found that ice wedges in the Yedoma region represent a significant DOC (45.2 Tg) and DIC (33.6 Tg) pool in permafrost areas and a freshwater reservoir of 4200 km(2). This study underlines the need to discriminate between particulate OC and DOC to assess the availability and vulnerability of the permafrost car-bon pool for ecosystems and climate feedback upon mobilization.
Upwelling hot mantle plumes are thought to disintegrate continental lithosphere and are considered to be drivers of active continental breakup. The formation of the Walvis Ridge during the opening of the South Atlantic is related to a putative plume-induced breakup. We investigated the crustal structure of the Walvis Ridge (southeast Atlantic Ocean) at its intersection with the continental margin and searched for anomalies related to the possible plume head. The overall structure we identify suggests that no broad plume head existed during opening of the South Atlantic and anomalous mantle melting occurred only locally. We therefore question the importance of a plume head as a driver of continental breakup and further speculate that the hotspot was present before the rifting, leaving a track of kimberlites in the African craton.
The upper Middle to lower Upper Miocene (Serravallian to Tortonian) sedimentary succession in Salento (southern Italy) includes glauconite- and phosphate-rich deposits, which are associated with pelagic micrite. In Baia del Ciolo and Marittima (southern Salento), the succession is composed of shallow-water platform carbonates of Late Oligocene age (Chattian; Porto Badisco Formation), which are overlain by a 20- to 30-cm-thick level of glauconite-rich micrite with abundant reworked particles and fossils of the underlying Porto Badisco Formation. This interval is in turn covered by an up to 15 cm thick phosphatic crust ("Livello ad Aturia"), which itself is overlain either by a hemipelagic chalk-like carbonate of Middle to Late Miocene age ("Pietra Leccese"; Marittima) or directly by a micrite of Late Miocene age (Messinian; Novaglie Formation; Baia del Ciolo), which shallows upwards into a shallow-water platform carbonate. A large hiatus is present in this succession, which likely includes the Lower and lower Middle Miocene. In the region of Lecce, two discrete levels enriched in glauconite and phosphate-each associated with a major discontinuity-occur within the Pietra Leccese. The strontium-isotope ages derived on phosphate nodules and phosphatized and non-phosphatized fossils and calcareous nannofossil ages indicate a time interval of phosphogenesis between 13.5 and 7.5 Ma, with two clusters at 12 and 10.5 Ma. The glauconite and phosphate-rich sediments resulted from a current-dominated regime, which was characterized by low overall sedimentation rates, erosion and sediment reworking, and authigenesis. This regime was likely related to a generally westward-directed bottom current, which was forced to upwell once it arrived at the western border of the eastern Mediterranean basin. The timing of the principal phosphogenic phases can only partly be correlated to those of other occurrences in this part of the Mediterranean (Malta, Gozo, southern Sicily, Matese, Latium-Abruzzi platform). The partial diachrony in phosphogenesis may express the effect of lateral switching in and/or focusing of upwelling zones. (C) 2015 Elsevier B.V. All rights reserved.
Probing the core-mantle boundary beneath Europe and Western Eurasia: A detailed study using PcP
(2015)
We use PcP (the core reflected P phase) recordings of deep earthquakes and nuclear explosions from the Grafenberg (Germany) and NORSAR (Norway) arrays to investigate the core-mantle boundary region beneath Europe and western Eurasia. We find evidence for a previously unknown ultra-low velocity zone 600 km south-east of Moscow, located at the edge of a middle-size low shear- velocity region imaged in seismic tomography that is located beneath the Volga river region. The observed amplitude variations of PcP can be modelled by velocity reductions of P and S-waves of -5% and -15%, respectively, with a density increase of +15%. Travel time delays of pre-and postcursors are indicating a thickness of about 13 km for this ultra-low velocity region (ULVZ). However, our modelling also reveals highly ambiguous amplitude variations of PcP and a reflection off the top of the anomaly for various ULVZs and topography models. Accordingly, large velocity contrasts of up to -10% in V-P and -20% in Vs cannot be excluded. In general, the whole Volga river region shows a complex pattern of PcP amplitudes caused most likely by CMB undulations. Further PcP probes beneath Paris, Kiev and northern Italy indicate likely normal CMB conditions, whereas the samples below Finland and the Hungary-Slovakia border yield strongly amplified PcP signals suggesting strong CMB topography effects.
We evaluate the amplitude behaviour of PcP as a function of distance and several ULVZ models using the 1D reflectivity and the 2D Gauss beam method. The influence of the velocity and density perturbations is analysed as well as the anomaly thickness, the dominant period of the source wavelet and interface topographies. Strong variation of the PcP amplitude are obtained as a function of distance and of the impedance contrast. We also consider two types of topographies: undulations atop the CMB in the presence of flat ULVZs and vice versa. Where a broad range of CMB topography dimensions lead to large PcP amplitude variations, only large ULVZ undulations generate significant amplitude scattering. Consequently, this indicates that topography effects of anomalies may mask the true medium parameters as well as the ULVZ thickness. Moreover, there might be a possibility of misinterpreting the precursor as PcP, in particular for thin ULVZs. (C) 2015 Elsevier B.V. All rights reserved.
Scientific theories of how subduction and plate tectonics began on Earth-and what the tectonic structure of Earth was before this-remain enigmatic and contentious(1). Understanding viable scenarios for the onset of subduction and plate tectonics(2,3) is hampered by the fact that subduction initiation processes must have been markedly different before the onset of global plate tectonics because most present-day subduction initiation mechanisms require acting plate forces and existing zones of lithospheric weakness, which are both consequences of plate tectonics(4). However, plume-induced subduction initiation(5-9) could have started the first subduction zone without the help of plate tectonics. Here, we test this mechanism using high-resolution three-dimensional numerical thermomechanical modelling. We demonstrate that three key physical factors combine to trigger self-sustained subduction: (1) a strong, negatively buoyant oceanic lithosphere; (2) focused magmatic weakening and thinning of lithosphere above the plume; and (3) lubrication of the slab interface by hydrated crust. We also show that plume-induced subduction could only have been feasible in the hotter early Earth for old oceanic plates. In contrast, younger plates favoured episodic lithospheric drips rather than self-sustained subduction and global plate tectonics.
On 2012 August 11, a pair of large, damaging earthquakes struck the Varzaghan-Ahar region in northwest Iran, in a region where there was no major mapped fault or any well-documented historical seismicity. To investigate the active tectonics of the source region we applied a combination of seismological methods (local aftershock network, calibrated multiple event relocation and focal mechanism studies), field observations (structural geology and geomorphological) and inversions for the regional stress field. The epicentral region is north of the North Tabriz Fault. The first main shock is characterized by right-lateral strike-slip motion on an almost E-W fault plane of about 23 km length extending from the surface to a depth of about 14 km. The second main shock occurred on an ENE-striking fault that dips at 60-70A degrees to the NW. Independent inversions of focal mechanisms and geologically determined fault kinematic data for the active stress state yield a transpressional tectonic regime with sigma(1) oriented N132E. For the region northeast of the North Tabriz Fault, the presence of rigid lithosphere of the South Caspian Basin implies the kinematic adjustment by northward transferring of the contracted masses through both distributed deformation and structural deflections. Our results suggest that the kinematic adjustment inside a contracting wedge may occur along interacting crosswise or conjugate faults to accommodate low rates of internal deformation. At a global scale, our results indicate that despite the basic assumption of 'rigid blocks' in geodetic plate modelling, internal deformation of block-like regions could control the kinematics of deformation and the level of seismic hazard within and around such regions of low deformation rate.
We propose a RAndom Interacting Network (RAIN) model to study the interactions between a pair of complex networks. The model involves two major steps: (i) the selection of a pair of nodes, one from each network, based on intra-network node-based characteristics, and (ii) the placement of a link between selected nodes based on the similarity of their relative importance in their respective networks. Node selection is based on a selection fitness function and node linkage is based on a linkage probability defined on the linkage scores of nodes. The model allows us to relate within-network characteristics to between-network structure. We apply the model to the interaction between the USA and Schengen airline transportation networks (ATNs). Our results indicate that two mechanisms: degree-based preferential node selection and degree-assortative link placement are necessary to replicate the observed inter-network degree distributions as well as the observed inter-network assortativity. The RAIN model offers the possibility to test multiple hypotheses regarding the mechanisms underlying network interactions. It can also incorporate complex interaction topologies. Furthermore, the framework of the RAIN model is general and can be potentially adapted to various real-world complex systems.
Most airborne electromagnetic (EM) processing programs assume a flat ground surface. However, in mountainous areas, the system can be at an angle with regard to the ground. As the system is no longer parallel to the ground surface, the measured magnetic field has to be corrected and the ground induced eddy current has to be modelled in a better way when performing a very fine interpretation of the data. We first recall the theoretical background for the modelling of a magnetic dipole source and study it in regard to the case of an arbitrarily oriented magnetic dipole. We show in particular how transient central loop helicopter borne data are influenced by this inclination. The result shows that the effect of topography on airborne EM is more important at early time windows and for systems using a short cut-off source. In this paper, we suggest that an estimate be made off the locally averaged inclination of the system to the ground and then to correct the data for this before inverting it (whether the inversion assumes a flat 1D, 2D or 3D sub-surface). Both 1D and 2D inversions are applied to synthetic and real data sets with such a correction. The consequence on the ground imaging is small for slopes with an angle less than 25 degrees but the correction factor can be useful for improving the estimation of depths in mountainous areas.
Ground loop-loop electromagnetic surveys are often conducted to fulfill the low-induction-number condition. To image the distribution of electric conductivity inside the ground, it is then necessary to collect a multioffset data set. We considered that less time-consuming constant offset measurements can also reach this objective. This can be achieved by performing multifrequency soundings, which are commonly performed for the airborne electromagnetic method. Ground multifrequency soundings have to be interpreted carefully because they contain high-induction-number data. These data are interpreted in two steps. First, the in-phase and out-of-phase data are converted into robust apparent conductivities valid for all the induction numbers. Second, the apparent conductivity data are inverted in 1D and 2D to obtain the true distribution of the ground conductivity. For the inversion, we used a general half-space Jacobian for the apparent conductivity valid for all the induction numbers. This method was applied and validated on synthetic data computed with the full Maxwell theory. The method was then applied on field data acquired in the test site of Provins, in the Parisian basin, France. The result revealed good agreement with borehole and geologic information, demonstrating the applicability of our method.
For near surface geophysical surveys, small-fixed offset loop-loop electromagnetic induction (EMI) sensors are usually placed parallel to the ground surface (i.e., both loops are at the same height above ground). In this study, we evaluate the potential of making measurements with a system that is not parallel to the ground; i.e., by positioning the system at different inclinations with respect to ground surface. First, we present the Maxwell theory for inclined magnetic dipoles over a homogeneous half space. By analyzing the sensitivities of such configurations, we,show that varying the angle of the system would result in improved imaging capabilities. For example, we show that acquiring data with a vertical system allows detection of a conductive body with a better lateral resolution compared to data acquired using standard horizontal configurations. The synthetic responses are presented for a heterogeneous medium and compared to field data acquired in the historical Park Sanssouci in Potsdam, Germany. After presenting a detailed sensitivity analysis and synthetic examples of such ground conductivity measurements, we suggest a new strategy of acquisition that allows to better estimate the true distribution of electrical conductivity using instruments with a fixed, small offset between the loops. This strategy is evaluated using field data collected at a well-constrained test-site in Horstwalde (Germany). Here, the target buried utility pipes are best imaged using vertical system configurations demonstrating the potential of our approach for typical applications. (C) 2015 Elsevier B.V. Pill rights reserved.
Interlayered blueschists and greenschists of the Coastal Cordillera (Chile) are part of a Late Palaeozoic accretionary complex. They represent metavolcanic rocks with oceanic affinities based on predominantly 01B-type REE patterns and immobile trace element ratios. Both rock types have similar mineralogies, albeit with different mineral modal abundances. Amphibole is the major mafic mineral and varies compositionally from glaucophane to actinolite. The presence of glaucophane relicts as cores in zoned amphiboles in both blueschists and greenschists is evidence for a pervasive high-pressure metamorphic stage, indicating that tectonic juxtaposition is an unlikely explanation for the cm-dm scale interlayering. During exhumation, a retrograde greenschist-facies overprint stabilized chlorite + albite + winchitic/actinolitic amphibole + phengitic white mica +/- epidote +/- K-feldspar at 0.4 +/- 0.1 GPa. Geochemical variability can be partly ascribed to primary magmatic and partly to secondary metasomatic processes that occurred under greenschist-facies conditions. Isocon diagrams of several adjacent blueschist-greenschist pairs with similar protolith geochemistry were used to evaluate metasomatic changes due to retrograde fluid-rock interaction. The most important geochemical changes are depletion of Si and Na and addition of water in the greenschists compared to the blueschists. Transition metals and LILE are mobilized to varying degrees. The unsystematic deviations from magmatic fractionation trends suggest open system conditions and influx of an external fluid. Pseudosection and water isopleth calculations show that the rocks were dehydrating during most of their exhumation history and remained at water-saturated conditions. The mineralogical changes, in particular breakdown of blue amphibole and replacement by chlorite, albite and calcic/sodic-calcic amphibole, are the prime cause for the distinct coloring. Pseudo-binary phase diagrams were used as a means to link bulk rock geochemical variability to modal and chemical changes in the mineralogy. The geochemical changes induced by fluid-rock interaction are important in two ways: First, the bulk rock chemistry is altered, leading to the stabilization of higher modal proportions of chlorite in the greenschists. Second, the retrograde overprint is a selective, layer-parallel fluid infiltration process, causing more intense greenschist-facies recrystallization in greenschist layers and therefore preferential preservation of blue amphibole in blueschist layers. Hence, the distinct colors were acquired by a combination of compositional variability, both primary magmatic and secondary metasomatic, and the different intensity of retrograde fluid infiltration. (C) 2014 Elsevier B.V. All rights reserved.
Geodetic and seismologic studies support a tectonic model for the central Himalaya wherein similar to 2 cm/yr of Indo-Asian convergence is accommodated along the primary decollement under the range, the Main Himalayan thrust. A steeper midcrustal ramp in the Main Himalayan thrust is commonly invoked as driving rapid rock uplift along a range-parallel band in the Greater Himalaya. This tectonic model, developed primarily from studies in central Nepal, is commonly assumed to project along strike with little lateral variation in Main Himalayan thrust geometry or associated rock uplift patterns. Here, we synthesize multiple lines of evidence for a major discontinuity in the Main Himalayan thrust in western Nepal. Analysis of topography and seismicity indicates that west of similar to 82.5 degrees E, the single band of steep topography and seismicity along the Main Himalayan thrust ramp in central Nepal bifurcates around a high-elevation, low-relief landscape, resulting in a two-step topographic front along an similar to 150 km segment of the central Himalaya. Although multiple models could explain this bifurcation, the full suite of data appears to be most consistent with a northward bend to the Main Himalayan thrust ramp and activation of a young duplex horse to the south. This poorly documented segmentation of the Main Himalayan thrust has important implications for the seismogenic potential of the western Nepal seismic gap and for models of the ongoing evolution of the orogen.
Applications of graph theory have proliferated across the academic spectrum in recent years. Whereas geosciences and landscape ecology have made rich use of graph theory, its use seems limited in physical geography, and particularly in geomorphology. Common applications of graph theory analyses of connectivity, path or transport efficiencies, subnetworks, network structure, system behaviour and dynamics, and network optimization or engineering all have uses or potential uses in geomorphology and closely related fields. In this paper, we give a short introduction to graph theory and review previous geomorphological applications or works in related fields that have been particularly influential. Network-like geomorphic systems can be classified into nonspatial or spatially implicit system components linked by statistical/causal relationships and spatial units linked by some spatial relationship, for example by fluxes of matter and/or energy. We argue that, if geomorphic system properties and behaviour (e.g., complexity, sensitivity, synchronisability, historical contingency, connectivity etc.) depend on system structure and if graph theory is able to quantitatively describe the configuration of system components, then graph theory should provide us with tools that help in quantifying system properties and in inferring system behaviour. (C) 2015 Elsevier B.V. All rights reserved.
Weather radar analysis has become increasingly sophisticated over the past 50 years, and efforts to keep software up to date have generally lagged behind the needs of the users. We argue that progress has been impeded by the fact that software has not been developed and shared as a community.
Recently, the situation has been changing. In this paper, the developers of a number of open-source software (OSS) projects highlight the potential of OSS to advance radar-related research. We argue that the community-based development of OSS holds the potential to reduce duplication of efforts and to create transparency in implemented algorithms while improving the quality and scope of the software. We also conclude that there is sufficiently mature technology to support collaboration across different software projects. This could allow for consolidation toward a set of interoperable software platforms, each designed to accommodate very specific user requirements.
In a recent BAMS article, it is argued that community-based Open Source Software (OSS) could foster scientific progress in weather radar research, and make weather radar software more affordable, flexible, transparent, sustainable, and interoperable.
Nevertheless, it can be challenging for potential developers and users to realize these benefits: tools are often cumbersome to install; different operating systems may have particular issues, or may not be supported at all; and many tools have steep learning curves.
To overcome some of these barriers, we present an open, community-based virtual machine (VM). This VM can be run on any operating system, and guarantees reproducibility of results across platforms. It contains a suite of independent OSS weather radar tools (BALTRAD, Py-ART, wradlib, RSL, and Radx), and a scientific Python stack. Furthermore, it features a suite of recipes that work out of the box and provide guidance on how to use the different OSS tools alone and together. The code to build the VM from source is hosted on GitHub, which allows the VM to grow with its community.
We argue that the VM presents another step toward Open (Weather Radar) Science. It can be used as a quick way to get started, for teaching, or for benchmarking and combining different tools. It can foster the idea of reproducible research in scientific publishing. Being scalable and extendable, it might even allow for real-time data processing.
We expect the VM to catalyze progress toward interoperability, and to lower the barrier for new users and developers, thus extending the weather radar community and user base.
In a recent BAMS article, it is argued that community-based Open Source Software (OSS) could foster scientific progress in weather radar research, and make weather radar software more affordable, flexible, transparent, sustainable, and interoperable.
Nevertheless, it can be challenging for potential developers and users to realize these benefits: tools are often cumbersome to install; different operating systems may have particular issues, or may not be supported at all; and many tools have steep learning curves.
To overcome some of these barriers, we present an open, community-based virtual machine (VM). This VM can be run on any operating system, and guarantees reproducibility of results across platforms. It contains a suite of independent OSS weather radar tools (BALTRAD, Py-ART, wradlib, RSL, and Radx), and a scientific Python stack. Furthermore, it features a suite of recipes that work out of the box and provide guidance on how to use the different OSS tools alone and together. The code to build the VM from source is hosted on GitHub, which allows the VM to grow with its community.
We argue that the VM presents another step toward Open (Weather Radar) Science. It can be used as a quick way to get started, for teaching, or for benchmarking and combining different tools. It can foster the idea of reproducible research in scientific publishing. Being scalable and extendable, it might even allow for real-time data processing.
We expect the VM to catalyze progress toward interoperability, and to lower the barrier for new users and developers, thus extending the weather radar community and user base.
In a recent BAMS article, it is argued that community-based Open Source Software (OSS) could foster scientific progress in weather radar research, and make weather radar software more affordable, flexible, transparent, sustainable, and interoperable.Nevertheless, it can be challenging for potential developers and users to realize these benefits: tools are often cumbersome to install; different operating systems may have particular issues, or may not be supported at all; and many tools have steep learning curves.To overcome some of these barriers, we present an open, community-based virtual machine (VM). This VM can be run on any operating system, and guarantees reproducibility of results across platforms. It contains a suite of independent OSS weather radar tools (BALTRAD, Py-ART, wradlib, RSL, and Radx), and a scientific Python stack. Furthermore, it features a suite of recipes that work out of the box and provide guidance on how to use the different OSS tools alone and together. The code to build the VM from source is hosted on GitHub, which allows the VM to grow with its community.We argue that the VM presents another step toward Open (Weather Radar) Science. It can be used as a quick way to get started, for teaching, or for benchmarking and combining different tools. It can foster the idea of reproducible research in scientific publishing. Being scalable and extendable, it might even allow for real-time data processing.We expect the VM to catalyze progress toward interoperability, and to lower the barrier for new users and developers, thus extending the weather radar community and user base.
A seismological network was operated at the junction of the aseismic Walvis Ridge with the northwestern Namibian coast. We mapped crustal thickness and bulk V-p/V-s ratio by the H-k analysis of receiver functions. In the Damara Belt, the crustal thickness is similar to 35km with a V-p/V-s ratio of <1.75. The crust is similar to 30km thick at the coast in the Kaoko Belt. Strong variations in crustal thickness and V-p/V-s ratios are found at the landfall of the Walvis Ridge. Here and at similar to 150km northeast of the coast, the crustal thickness increases dramatically reaching 44km and the V-p/V-s ratios are extremely high (similar to 1.89). These anomalies are interpreted as magmatic underplating produced by the mantle plume during the breakup of Gondwana. The area affected by the plume is smaller than 300km in diameter, possibly ruling out the existence of a large plume head under the continent during the breakup.
This thesis contains three experimental studies addressing the interplay between deformation and the mineral reaction between natural calcite and magnesite. The solid-solid mineral reaction between the two carbonates causes the formation of a magnesio-calcite precursor layer and a dolomite reaction rim in every experiment at isostatic annealing and deformation conditions.
CHAPTER 1 briefly introduces general aspects concerning mineral reactions in nature and diffusion pathways for mass transport. Moreover, results of previous laboratory studies on the influence of deformation on mineral reactions are summarized. In addition, the main goals of this study are pointed out.
In CHAPTER 2, the reaction between calcite and magnesite single crystals is examined at isostatic annealing conditions. Time series performed at a fixed temperature revealed a diffusion-controlled dolomite rim growth. Two microstructural domains could be identified characterized by palisade-shaped dolomite grains growing into the magnesite and granular dolomite growing towards calcite. A model was provided for the dolomite rim growth based on the counter-diffusion of CaO and MgO. All reaction products exhibited a characteristic crystallographic relationship with respect to the calcite reactant. Moreover, kinetic parameters of the mineral reaction were determined out of a temperature series at a fixed time. The main goal of the isostatic test series was to gain information about the microstructure evolution, kinetic parameters, chemical composition and texture development of the reaction products. The results were used as a reference to quantify the influence of deformation on the mineral reaction.
CHAPTER 3 deals with the influence of non-isostatic deformation on dolomite and magnesio-calcite layer production between calcite and magnesite single crystals. Deformation was achieved by triaxial compression and by torsion. Triaxial compression up to 38 MPa axial stress at a fixed time showed no significant influence of stress and strain on dolomite formation. Time series conducted at a fixed stress yield no change in growth rates for dolomite and magnesio-calcite at low strains. Slightly larger magnesio-calcite growth rates were observed at strains above >0.1. High strains at similar stresses were caused by the activation of additional glide systems in the calcite single crystal and more mobile dislocations in the magnesio-calcite grains, providing fast diffusion pathways. In torsion experiments a gradual decrease in dolomite and magnesio-calcite layer thickness was observed at a critical shear strain. During deformation, crystallographic orientations of reaction products rearranged with respect to the external framework. A direct effect of the mineral reaction on deformation could not be recognized due to the relatively small reaction product widths.
In CHAPTER 4, the influence of starting material microfabrics and the presence of water on the reaction kinetics was evaluated. In these experimental series polycrystalline material was in contact with single crystals or two polycrystalline materials were used as reactants. Isostatic annealing resulted in different dolomite and magnesio-calcite layer thicknesses, depending on starting material microfabrics. The reaction progress at the magnesite interface was faster with smaller magnesite grain size, because grain boundaries provided fast pathways for diffusion and multiple nucleation sites for dolomite formation. Deformation by triaxial compression and torsion yield lower dolomite rim thicknesses compared to annealed samples for the same time. This was caused by grain coarsening of polycrystalline magnesite during deformation. In contrast, magnesio-calcite layers tended to be larger during deformation, which triggered enhanced diffusion along grain boundaries. The presence of excess water had no significant influence on the reaction kinetics, at least if the reactants were single crystals.
In CHAPTER 5 general conclusions about the interplay between deformation and the mineral reaction in the carbonate system are presented.
Finally, CHAPTER 6 highlights possible future work in the carbonate system based on the results of this study.
Contrasting response of two shallow eutrophic cold temperate lakes to a partial winterkill of fish
(2015)
Food-web effects of winterkill are difficult to predict as the enhanced mortality of planktivorous fish may be counterbalanced by an even higher mortality of piscivores. We hypothesised that a winterkill in a clear and a turbid shallow lake would equalise their fish community composition, but seasonal plankton successions would differ between lakes. After a partial winterkill, we observed a reduction of fish biomass by 16 and 43% in a clear-water and a turbid small temperate lake, respectively. Fish biomass and piscivore shares (5% of fish biomass) were similar in both lakes after this winterkill, but young-of-the-year (YOY) abundances were higher in the turbid lake. Top-down control by crustaceans was only partly responsible for low phytoplankton biomass at the end of May following the winterkill in both lakes. Summer phytoplankton biomass remained low in the clear-water lake despite high abundances of YOY fish (mainly roach). In contrast, the crustacean biomass of the turbid lake was reduced in summer by a high YOY abundance (sunbleak and roach), leading to a strong increase in phytoplankton biomass. The YOY abundance of fish in shallow eutrophic lakes may thus be more important for their summer phytoplankton development after winterkill than the relative abundance of piscivores.
Within the scope of Russian German palaeoenvironmental research, Two-Yurts Lake (TYL, Dvuh-Yurtochnoe in Russian) was chosen as the main scientific target area to decipher Holocene climate variability on Kamchatka. The 5 x 2 km large and 26 m deep lake is of proglacial origin and situated on the eastern flank of Sredinny Ridge at the northwestern end of the Central Kamchatka Valley, outside the direct influence of active volcanism. Here, we present results of a multi-proxy study on sediment cores, spanning about the last 7000 years. The general tenor of the TYL record is an increase in continentality and winter snow cover in conjunction with a decrease in temperature, humidity, and biological productivity after 5000-4500 cal yrs BP, inferred from pollen and diatom data and the isotopic composition of organic carbon. The TYL proxy data also show that the late Holocene was punctuated by two colder spells, roughly between 4500 and 3500 cal yrs BP and between 1000 and 200 cal yrs BP, as local expressions of the Neoglacial and Little Ice Age, respectively. These environmental changes can be regarded as direct and indirect responses to climate change, as also demonstrated by other records in the regional terrestrial and marine realm. Long-term climate deterioration was driven by decreasing insolation, while the short-term climate excursions are best explained by local climatic processes. The latter affect the configuration of atmospheric pressure systems that control the sources as well as the temperature and moisture of air masses reaching Kamchatka. (C) 2015 Elsevier B.V. All rights reserved.
Does textural heterogeneity matter? Quantifying transformation of hydrological signals in soils
(2015)
Textural heterogeneity causes complex water flow patterns and soil moisture dynamics in soils that hamper monitoring and modeling soil hydrological processes. These patterns can be generated by process based models considering soil texture heterogeneities. However, there is urgent need for tools for the inverse approach, that is, to analyze observed dynamics in a quantitative way independent from any model approach in order to identify effects of soil texture heterogeneity. Here, studying the transformation of hydrological input signals (e.g., rainfall, snow melt) propagating through the vadose zone is a promising supplement to the common perspective of mass flux considerations. In this study we applied a recently developed new approach for quantitative analysis of hydrological time series (i) to investigate the effect of soil texture on the signal transformation behavior and (ii) to analyze to what degree soil moisture dynamics from a heterogeneous profile can be reproduced by a corresponding homogenous substrate. We used simulation models to generate three data sets of soil moisture time series considering homogeneous substrates (HOM), homogeneous substrates with noise added (NOISE), and heterogeneous substrates (HET). The soil texture classes sand, loamy sand, clay loam and silt were considered. We applied a principal component analysis (also called empirical orthogonal functions) to identify predominant functional patterns and to measure the degree of signal transformation of single time series. For the HOM case 86.7% of the soil moisture dynamics were reproduced by the first two principal components. Based on these results a quantitative measure for the degree of transformation of the input signal was derived. The general nature of signal transformation was nearly identical in all textures, but the intensity of signal damping per depth interval decreased from fine to coarse textures. The same functional patterns occurred in the HET data set. However, here the signal damping of time series did not increase monotonically with soil depth. The analysis succeeded in extracting the same signal transformation behavior from the NOISE data set compared to that of the HOM case in spite of being blurred by random noise. Thus, principal component analysis proved to be a very robust tool to disentangle between independent effects and to measure the degree of transformation of the input signal. The suggested approach can be used for (i) data processing, including subtracting measurement noise (ii) identification of factors controlling soil water dynamics, (iii) assessing the mean signal transformation in heterogeneous soils based on observed soil moisture time series, and (iv) model building, calibration and evaluation. (C) 2015 Elsevier B.V. All rights reserved.
Paleomagnetic dating of the India-Asia collision hinges on determining the Paleogene latitude of the Lhasa terrane (southern Tibet). Reported latitudes range from 5 degrees N to 30 degrees N, however, leading to contrasting paleogeographic interpretations. Here we report new data from the Eocene Linzizong volcanic rocks in the Nanmulin Basin, which previously yielded data suggesting a low paleolatitude (similar to 10 degrees N). New zircon U-Pb dates indicate an age of similar to 52Ma. Negative fold tests, however, demonstrate that the isolated characteristic remanent magnetizations, with notably varying inclinations, are not primary. Rock magnetic analyses, end-member modeling of isothermal remanent magnetization acquisition curves, and petrographic observations are consistent with variable degrees of posttilting remagnetization due to low-temperature alteration of primary magmatic titanomagnetite and the formation of secondary pigmentary hematite that unblock simultaneously. Previously reported paleomagnetic data from the Nanmulin Basin implying low paleolatitude should thus not be used to estimate the time and latitude of the India-Asia collision. We show that the paleomagnetic inclinations vary linearly with the contribution of secondary hematite to saturation isothermal remanent magnetization. We tentatively propose a new method to recover a primary remanence with inclination of 38.1 degrees (35.7 degrees, 40.5 degrees) (95% significance) and a secondary remanence with inclination of 42.9 degrees (41.5 degrees,44.4 degrees) (95% significance). The paleolatitude defined by the modeled primary remanence21 degrees N (19.8 degrees N, 23.1 degrees N)is consistent with the regional compilation of published results from pristine volcanic rocks and sedimentary rocks of the upper Linzizong Group corrected for inclination shallowing. The start of the Tibetan Himalaya-Asia collision was situated at similar to 20 degrees N and took place by similar to 50Ma.
The Paleogene latitude of the Lhasa terrane (southern Tibet) can constrain the age of the onset of the India-Asia collision. Estimates for this latitude, however, vary from 5 degrees N to 30 degrees N, and thus, here, we reassess the geochronology and paleomagnetism of Paleogene volcanic rocks from the Linzizong Group in the Linzhou basin. The lower and upper parts of the section previously yielded particularly conflicting ages and paleolatitudes. We report consistent Ar-40/Ar-39 and U-Pb zircon dates of similar to 52Ma for the upper Linzizong, and Ar-40/Ar-39 dates (similar to 51Ma) from the lower Linzizong are significantly younger than U-Pb zircon dates (64-63Ma), suggesting that the lower Linzizong was thermally and/or chemically reset. Paleomagnetic results from 24 sites in lower Linzizong confirm a low apparent paleolatitude of similar to 5 degrees N, compared to the upper part (similar to 20 degrees N) and to underlying Cretaceous strata (similar to 20 degrees N). Detailed rock magnetic analyses, end-member modeling of magnetic components, and petrography from the lower and upper Linzizong indicate widespread secondary hematite in the lower Linzizong, whereas hematite is rare in upper Linzizong. Volcanic rocks of the lower Linzizong have been hydrothermally chemically remagnetized, whereas the upper Linzizong retains a primary remanence. We suggest that remagnetization was induced by acquisition of chemical and thermoviscous remanent magnetizations such that the shallow inclinations are an artifact of a tilt correction applied to a secondary remanence in lower Linzizong. We estimate that the Paleogene latitude of Lhasa terrane was 204 degrees N, consistent with previous results suggesting that India-Asia collision likely took place by similar to 52Ma at similar to 20 degrees N.
The Tibetan Himalaya represents the northernmost continental unit of the Indian plate that collided with Asia in the Cenozoic. Paleomagnetic studies on the Tibetan Himalaya can help constrain the dimension and paleogeography of "Greater India,' the Indian plate lithosphere that subducted and underthrusted below Asia after initial collision. Here we present a paleomagnetic investigation of a Jurassic (limestones) and Lower Cretaceous (volcaniclastic sandstones) section of the Tibetan Himalaya. The limestones yielded positive fold test, showing a prefolding origin of the isolated remanent magnetizations. Detailed paleomagnetic analyses, rock magnetic tests, end-member modeling of acquisition curves of isothermal remanent magnetization, and petrographic investigation reveal that the magnetic carrier of the Jurassic limestones is authigenic magnetite, whereas the dominant magnetic carrier of the Lower Cretaceous volcaniclastic sandstones is detrital magnetite. Our observations lead us to conclude that the Jurassic limestones record a prefolding remagnetization, whereas the Lower Cretaceous volcaniclastic sandstones retain a primary remanence. The volcaniclastic sandstones yield an Early Cretaceous paleolatitude of 55.5 degrees S [52.5 degrees S, 58.6 degrees S] for the Tibetan Himalaya, suggesting it was part of the Indian continent at that time. The size of "Greater India' during Jurassic time cannot be estimated from these limestones. Instead, a paleolatitude of the Tibetan Himalaya of 23.8 degrees S [21.8 degrees S, 26.1 degrees S] during the remagnetization process is suggested. It is likely that the remagnetization, caused by the oxidation of early diagenetic pyrite to magnetite, was induced during 103-83 or 77-67 Ma. The inferred paleolatitudes at these two time intervals imply very different tectonic consequences for the Tibetan Himalaya.
Several solutions have been proposed to explain the long-standing kinematic observation that postcollisional upper crustal shortening within the Himalaya and Asia is much less than the magnitude of India-Asia convergence. Here we implement these hypotheses in global plate reconstructions and test paleolatitudes predicted by the global apparent polar wander path against independent, and the most robust paleomagnetic data. Our tests demonstrate that (1) reconstructed 600-750km postcollisional intra-Asian shortening is a minimum value; (2) a 52Ma collision age is only consistent with paleomagnetic data if intra-Asian shortening was 900km; a 56-58Ma collision age requires greater intra-Asian shortening; (3) collision ages of 34 or 65Ma incorrectly predict Late Cretaceous and Paleogene paleolatitudes of the Tibetan Himalaya (TH); and (4) Cretaceous counterclockwise rotation of India cannot explain the paleolatitudinal divergence between the TH and India. All hypotheses, regardless of collision age, require major Cretaceous extension within Greater India.
The India-Asia suture zone of southern Tibet exposes Lower Cretaceous Xigaze ophiolites and radiolarian cherts, and time-equivalent Asian-derived clastic forearc sedimentary rocks (Xigaze Group). These ophiolites have been interpreted to have formed in the forearc of the north-dipping subduction zone below Tibet that produced the Gangdese magmatic arc around 15-20 degrees N, or in the forearc of a subequatorial intra-oceanic subduction zone. To better constrain the latitude of the ophiolites, we carried out an integrated paleomagnetic, geochronologic and stratigraphical study on epi-ophiolitic radiolarites (Chongdui and Bainang sections), and Xigaze Group turbiditic sandstones unconformably overlying the ophiolite's mantle units (Sangsang section). Detrital zircon U-Pb geochronology of tuffaceous layers from the Chongdui section and sandstones of the Xigaze Group at the Sangsang section provides maximum depositional ages of 116.5 +/- 3.1 Ma and 128.8 +/- 3.4 Ma, respectively, for the Chongdui section and an Asian provenance signature for the Xigaze Group. Paleomagnetic analyses, integrated with rock magnetic experiments, indicate significant compaction-related inclination 'shallowing' of the remanence within the studied rocks. Two independent methods are applied for the inclination shallowing correction of the paleomagnetic directions from the Sangsang section, yielding consistent mean paleolatitudes of 16.2 degrees N 113 degrees N, 20.9 degrees N] and 16.8 degrees N [11.1 degrees N, 23.3 degrees N], respectively. These results are indistinguishable from recent paleolatitude estimates for the Gangdese arc in southern Tibet. Radiolarites from the Chongdui and Bainang sections yield low paleomagnetic inclinations that would suggest a sub-equatorial paleolatitude, but the distribution of the paleomagnetic directions in these rocks strongly suggests a low inclination bias by compaction. Our data indicate that spreading of the Xigaze ophiolite occurred in the Gangdese forearc, and formed the basement of the forearc strata. (C) 2015 Elsevier B.V. All rights reserved.
The Brazilian Cerrado is recognized as one of the most threatened biomes in the world, as the region has experienced a striking change from natural Cerrado vegetation to intense cash crop production. This paper reviews the history of land conversion in the Cerrado and the development of soil properties and water resources under past and ongoing land use. We compared soil and water quality parameters from different land uses considering 80 soil and 18 water studies conducted in different regions across the Cerrado to provide quantitative evidence of soil and water alterations from land use change.
Following the conversion of native Cerrado, significant effects on soil pH, bulk density and available P and K for croplands and less-pronounced effects on pastures were evident. Soil total N did not differ between land uses because most of the sites classified as croplands were nitrogen-fixing soybeans, which are not artificially fertilized with N. In contrast, water quality studies showed nitrogen enrichment in agricultural catchments, indicating fertilizer impacts and potential susceptibility to eutrophication. Regardless of the land use, P is widely absent because of the high-fixing capacities of deeply weathered soils and the filtering capacity of riparian vegetation. Pesticides, however, were consistently detected throughout the entire aquatic system. In several case studies, extremely high-peak concentrations exceeded Brazilian and European Union (EU) water quality limits, which were potentially accompanied by serious health implications. Land use intensification is likely to continue, particularly in regions where less annual rainfall and severe droughts are projected in the northeastern and western Cerrado. Thus, the leaching risk and displacement of agrochemicals are expected to increase, particularly because the current legislation has caused a reduction in riparian vegetation. We conclude that land use intensification is likely to seriously limit the Cerrado's future regarding both agricultural productivity and ecosystem stability. Because only limited data are available, we recommend further field studies to understand the interaction between terrestrial and aquatic systems. This study may serve as a valuable database for integrated modelling to investigate the impact of land use and climate change on soil and water resources and to test and develop mitigation measures for the Cerrado. Copyright (C) 2014 John Wiley & Sons, Ltd.
The Brazilian Cerrado is recognised as one of the most threatened biomes in the world, as the region has experienced a striking change from natural vegetation to intense cash crop production. The impacts of rapid agricultural expansion on soil and water resources are still poorly understood in the region. Therefore, the overall aim of the thesis is to improve our understanding of the ecohydrological processes causing water and soil degradation in the Brazilian Cerrado.
I first present a metaanalysis to provide quantitative evidence and identifying the main impacts of soil and water alterations resulting from land use change. Second, field studies were conducted to (i) examine the effects of land use change on soils of natural cerrado transformed to common croplands and pasture and (ii) indicate how agricultural production affects water quality across a meso-scale catchment. Third, the ecohydrological process-based model SWAT was tested with simple scenario analyses to gain insight into the impacts of land use and climate change on the water cycling in the upper São Lourenço catchment which experienced decreasing discharges in the last 40 years.
Soil and water quality parameters from different land uses were extracted from 89 soil and 18 water studies in different regions across the Cerrado. Significant effects on pH, bulk density and available P and K for croplands and less-pronounced effects on pastures were evident. Soil total N did not differ between land uses because most of the cropland sites were N-fixing soybean cultivations, which are not artificially fertilized with N. By contrast, water quality studies showed N enrichment in agricultural catchments, indicating fertilizer impacts and potential susceptibility to eutrophication. Regardless of the land use, P is widely absent because of the high-fixing capacities of deeply weathered soils and the filtering capacity of riparian vegetation. Pesticides, however, were consistently detected throughout the entire aquatic system. In several case studies, extremely high-peak concentrations exceeded Brazilian and EU water quality limits, which pose serious health risks.
My field study revealed that land conversion caused a significant reduction in infiltration rates near the soil surface of pasture (–96 %) and croplands (–90 % to –93 %). Soil aggregate stability was significantly reduced in croplands than in cerrado and pasture. Soybean crops had extremely high extractable P (80 mg kg–1), whereas pasture N levels declined. A snapshot water sampling showed strong seasonality in water quality parameters. Higher temperature, oxi-reduction potential (ORP), NO2–, and very low oxygen concentrations (<5 mg•l–1) and saturation (<60 %) were recorded during the rainy season. By contrast, remarkably high PO43– concentrations (up to 0.8 mg•l–1) were measured during the dry season. Water quality parameters were affected by agricultural activities at all sampled sub-catchments across the catchment, regardless of stream characteristic. Direct NO3– leaching appeared to play a minor role; however, water quality is affected by topsoil fertiliser inputs with impact on small low order streams and larger rivers. Land conversion leaving cropland soils more susceptible to surface erosion by increased overland flow events.
In a third study, the field data were used to parameterise SWAT. The model was tested with different input data and calibrated in SWAT-CUP using the SUFI-2 algorithm. The model was judged reliable to simulate the water balance in the Cerrado. A complete cerrado, pasture and cropland cover was used to analyse the impact of land use on water cycling as well as climate change projections (2039–2058) according to the projections of the RCP 8.5 scenario. The actual evapotranspiration (ET) for the cropland scenario was higher compared to the cerrado cover (+100 mm a–1). Land use change scenarios confirmed that deforestation caused higher annual ET rates explaining partly the trend of decreased streamflow. Taking all climate change scenarios into account, the most likely effect is a prolongation of the dry season (by about one month), with higher peak flows in the rainy season. Consequently, potential threats for crop production with lower soil moisture and increased erosion and sediment transport during the rainy season are likely and should be considered in adaption plans.
From the three studies of the thesis I conclude that land use intensification is likely to seriously limit the Cerrado’s future regarding both agricultural productivity and ecosystem stability. Because only limited data are available for the vast biome, we recommend further field studies to understand the interaction between terrestrial and aquatic systems. This thesis may serve as a valuable database for integrated modelling to investigate the impact of land use and climate change on soil and water resources and to test and develop mitigation measures for the Cerrado in the future.
In probabilistic seismic hazard analysis, different ground-motion prediction equations (GMPEs) are commonly combined within a logic tree framework. The selection of appropriate GMPEs, however, is a non-trivial task, especially for regions where strong motion data are sparse and where no indigenous GMPE exists because the set of models needs to capture the whole range of ground-motion uncertainty. In this study we investigate the aggregation of GMPEs into a mixture model with the aim to infer a backbone model that is able to represent the center of the ground-motion distribution in a logic tree analysis. This central model can be scaled up and down to obtain the full range of ground-motion uncertainty. The combination of models into a mixture is inferred from observed ground-motion data. We tested the new approach for Northern Chile, a region for which no indigenous GMPE exists. Mixture models were calculated for interface and intraslab type events individually. For each source type we aggregated eight subduction zone GMPEs using mainly new strong-motion data that were recorded within the Plate Boundary Observatory Chile project and that were processed within this study. We can show that the mixture performs better than any of its component GMPEs, and that it performs comparable to a regression model that was derived for the same dataset. The mixture model seems to represent the median ground motions in that region fairly well. It is thus able to serve as a backbone model for the logic tree.
Exploring, exploiting and evolving diversity of aquatic ecosystem models: a community perspective
(2015)
Here, we present a community perspective on how to explore, exploit and evolve the diversity in aquatic ecosystem models. These models play an important role in understanding the functioning of aquatic ecosystems, filling in observation gaps and developing effective strategies for water quality management. In this spirit, numerous models have been developed since the 1970s. We set off to explore model diversity by making an inventory among 42 aquatic ecosystem modellers, by categorizing the resulting set of models and by analysing them for diversity. We then focus on how to exploit model diversity by comparing and combining different aspects of existing models. Finally, we discuss how model diversity came about in the past and could evolve in the future. Throughout our study, we use analogies from biodiversity research to analyse and interpret model diversity. We recommend to make models publicly available through open-source policies, to standardize documentation and technical implementation of models, and to compare models through ensemble modelling and interdisciplinary approaches. We end with our perspective on how the field of aquatic ecosystem modelling might develop in the next 5-10 years. To strive for clarity and to improve readability for non-modellers, we include a glossary.
The architecture of coastal sequences in tectonically-active regions results mostly from a combination of sea-level and land-level changes. The objective of this study is to unravel these signals by combining sequence stratigraphy and sedimentology of near-shore sedimentary sequences in wave-built terraces. We focus on Santa Maria Island at the south-central Chile margin, which hosts excellent exposures of coastal sediments from Marine Isotope Stage 3. A novel method based on statistical analysis of grain-size distributions coupled with fades descriptions provided a detailed account of transgressive-regressive cycles. Radiocarbon ages from paleosols constrain the chronology between >53 and similar to 31 cal ka BP. Because the influence of glaciations can be neglected, we calculated relative sea-level curves by tying the onset of deposition on a bedrock abrasion platform to a global sea-level curve. The observed depositional cycles match those predicted for uplift rates between 1.2 and 1.8 m/ka. The studied sedimentary units represent depositional cycles that resulted in reoccupation events of an existing marine terrace. Our study demonstrates wave-built marine terrace deposits along clastic shorelines in temperate regions can be used to distinguish between tectonic uplift and climate-induced sea-level changes.
The segmentation of major fault systems in subduction zones controls earthquake magnitude and location, but the causes for the existence of segment boundaries and the relationships between long-term deformation and the extent of earthquake rupture, are poorly understood. We compare permanent and seismic-cycle deformation patterns along the rupture zone of the 2010 Maule earthquake (M8.8), which ruptured 500 km of the Chile subduction margin. We analyzed the morphology of MIS-5 marine terraces using LiDAR topography and established their chronology and coeval origin with twelve luminescence ages, stratigraphy and geomorphic correlation, obtaining a virtually continuous distribution of uplift rates along the entire rupture zone. The mean uplift rate for these terraces is 0.5 m/ka. This value is exceeded in three areas, which have experienced rapid emergence of up to 1.6 m/ka; they are located at the northern, central, and southern sectors of the rupture zone, referred to as Topocalma, Carranza and Arauco, respectively. The three sectors correlate with boundaries of eight great earthquakes dating back to 1730. The Topocalma and Arauco sectors, located at the boundaries of the 2010 rupture, consist of broad zones of crustal warping with wavelengths of 60 and 90 km, respectively. These two regions coincide with the axes of oroclinal bending of the entire Andean margin and correlate with changes in curvature of the plate interface. Rapid uplift at Carranza, in turn, is of shorter wavelength and associated with footwall flexure of three crustal-scale normal faults. The uplift rate at Carranza is inversely correlated with plate coupling as well as with coseismic slip, suggesting permanent deformation may accumulate interseismically. We propose that the zones of upwarping at Arauco and Topocalma reflect changes in frictional properties of the megathrust resulting in barriers to the propagation of great earthquakes. Slip during the 1960 (M9.5) and 2010 events overlapped with the similar to 90-km-long zone of rapid uplift at Arauco; similarly, slip in 2010 and 1906 extended across the similar to 60-km-long section of the megathrust at Topocalma, but this area was completely breached by the 1730 (M similar to 9) event, which propagated southward until Carranza. Both Arauco and Topocalma show evidence of sustained rapid uplift since at least the middle Pleistocene. These two sectors might thus constitute discrete seismotectonic boundaries restraining most, but not all great earthquake ruptures. Based on our observations, such barriers might be breached during multi-segment super-cycle events. (C) 2015 Elsevier Ltd. All rights reserved.
Plan oblique relief shows terrain with a side view on a two-dimensional map, resulting in visualizations where the third dimension of the terrain is more explicit than on traditional two-dimensional maps. Existing plan oblique maps are static: the angle of terrain inclination is not adjustable and the orientation of plan oblique inclination does not change with the orientation of the map. This article introduces two complementary methods that address these issues by using the 3D graphics pipeline to render plan oblique relief for tile-based web maps. The goal is to allow users to adjust the terrain inclination and map rotation angles to better visualize the third dimension of the terrain. The first method pre-renders plan oblique tiles with a server-side application. The tiles are visualized with a standard web mapping framework. The second method renders plan oblique relief on-the-fly in a web browser using WebGL and a customized version of OpenLayers 3, which enables users to select arbitrary terrain inclination and map rotation angles. The second method uses a tiled digital terrain model that is loaded by the web browser. The browser applies the plan oblique transformation, computes a shaded relief, and texturizes the terrain with tiled map layers.
We present ostracod data from the Middle Palaeolithic open air site of Nahal Mahanayeem Outlet (NMO) at the southern edge of the Hula Basin, northern Israel.
Sediments of two Pleistocene water bodies are exposed at the site. The first one is an archaeologically sterile, light-colored limnic carbonate with an Early Pleistocene age. It contains an ostracod fauna assemblage dominated by Candona neglecta, Candonopsis kingsleii, and Pseudocandona sp., and, in minor abundances, Cypria ophtalmica, Cyprideis sp., Humphcypris sp., Fabaeformiscandona cf. fabaeformis and Ilyocypris sp. These sediments were deposited in a shallow, freshwater to oligohaline lake under stable conditions.
Sediments of the second water body are silty and dark-colored with a depositional age of 65 ka, belonging to the Late Pleistocene Ashmura Formation. The unit covers a geologically complex topography of tectonically uplifted limnic deposits and a hill-like gravel bar at the site. The most important archaeological layer is situated at its base, containing a lithic assemblage ascribed primarily to the Middle Palaeolithic Mousterian tradition and very well preserved flora and fauna. In the sediments from the archaeological layers, the brackish water ostracod Cyprideis torosa and the foraminifer Ammonia tepida could be identified. In sediments of the Ashmura Formation taken near the site, C. neglecta, Candona angulata, Ilyogpris sp., C kingsleii, Pseudocandona sp., C ophtalmica, Darwinula stevensoni, Trajancypris sp. and Potamogpris smaragdina were found, indicating a freshwater to slightly oligohaline stagnant water body. The ostracod fauna of the NMO site, together with geochemical data, allow us to reconstruct a depositional environment of the margin at a shallow lake with brackish or saline springs nearby during the site's occupation by Middle Palaeolithic humans.
Additionally, Candona weltneri, Candona cf. meerfeldiana, C kingsleii, Cyclocypris laevis, C ophtalmica, Cyprideis sp., Fabaeformiscandona cf. fabaeformis, P. smaragdina, Pseudocandona depressa, Trajancypris sp, Zonocypris cf. costata and A. tepida could be recorded for the first time for the Pleistocene limnic strata of the Hula Basin. (C) 2014 Elsevier Ltd. All rights reserved.
The ICDP "PaleoVan" drilling campaign at Lake Van, Turkey, provided a long (> 100 m) record of lacustrine subsurface sedimentary microbial cell abundance. After the ICDP campaign at Potrok Aike, Argentina, this is only the second time deep lacustrine cell counts have been documented. Two sites were cored and revealed a strikingly similar cell distribution despite differences in organic matter content and microbial activity. Although shifted towards higher values, cell counts from Lake Potrok Aike, Argentina, reveal very similar distribution patterns with depth. The lacustrine cell count data are significantly different from published marine records; the most probable cause is differences in sedimentary organic matter composition with marine sediments containing a higher fraction of labile organic matter. Previous studies showed that microbial activity and abundance increase centimetres to metres around geologic interfaces. The finely laminated Lake Van sediment allowed studying this phenomenon on the microscale. We sampled at the scale of individual laminae, and in some depth intervals, we found large differences in microbial abundance between the different laminae. This small-scale heterogeneity is normally overlooked due to much larger sampling intervals that integrate over several centimetres. However, not all laminated intervals exhibit such large differences in microbial abundance, and some non-laminated horizons show large variability on the millimetre scale as well. The reasons for such contrasting observations remain elusive, but indicate that heterogeneity of microbial abundance in subsurface sediments has not been taken into account sufficiently. These findings have implications not just for microbiological studies but for geochemistry as well, as the large differences in microbial abundance clearly show that there are distinct microhabitats that deviate considerably from the surrounding layers.
The ICDP "PaleoVan" drilling campaign at Lake Van, Turkey, provided a long (> 100 m) record of lacustrine subsurface sedimentary microbial cell abundance. After the ICDP campaign at Potrok Aike, Argentina, this is only the second time deep lacustrine cell counts have been documented. Two sites were cored and revealed a strikingly similar cell distribution despite differences in organic matter content and microbial activity. Although shifted towards higher values, cell counts from Lake Potrok Aike, Argentina, reveal very similar distribution patterns with depth. The lacustrine cell count data are significantly different from published marine records; the most probable cause is differences in sedimentary organic matter composition with marine sediments containing a higher fraction of labile organic matter. Previous studies showed that microbial activity and abundance increase centimetres to metres around geologic interfaces. The finely laminated Lake Van sediment allowed studying this phenomenon on the microscale. We sampled at the scale of individual laminae, and in some depth intervals, we found large differences in microbial abundance between the different laminae. This small-scale heterogeneity is normally overlooked due to much larger sampling intervals that integrate over several centimetres. However, not all laminated intervals exhibit such large differences in microbial abundance, and some non-laminated horizons show large variability on the millimetre scale as well. The reasons for such contrasting observations remain elusive, but indicate that heterogeneity of microbial abundance in subsurface sediments has not been taken into account sufficiently. These findings have implications not just for microbiological studies but for geochemistry as well, as the large differences in microbial abundance clearly show that there are distinct microhabitats that deviate considerably from the surrounding layers.
Earth's life-support systems are in rapid decline, yet we have few metrics or indicators with which to track these changes. The world's governments are calling for biodiversity and ecosystem-service monitoring to guide and evaluate international conservation policy as well as to incorporate natural capital into their national accounts. The Group on Earth Observations Biodiversity Observation Network (GEO BON) has been tasked with setting up this monitoring system. Here we explore the immediate feasibility of creating a global ecosystem-service monitoring platform under the GEO BON framework through combining data from national statistics, global vegetation models, and production function models. We found that nine ecosystem services could be annually reported at a national scale in the short term: carbon sequestration, water supply for hydropower, and non-fisheries marine products, crop, livestock, game meat, fisheries, mariculture, and timber production. Reported changes in service delivery over time reflected ecological shocks (e.g., droughts and disease outbreaks), highlighting the immediate utility of this monitoring system. Our work also identified three opportunities for creating a more comprehensive monitoring system. First, investing in input data for ecological process models (e.g., global land-use maps) would allow many more regulating services to be monitored. Currently, only 1 of 9 services that can be reported is a regulating service. Second, household surveys and censuses could help evaluate how nature affects people and provides non-monetary benefits. Finally, to forecast the sustainability of service delivery, research efforts could focus on calculating the total remaining biophysical stocks of provisioning services. Regardless, we demonstrated that a preliminary ecosystem-service monitoring platform is immediately feasible. With sufficient international investment, the platform could evolve further into a much-needed system to track changes in our planet's life-support systems. (C) 2015 Elsevier Ltd. All rights reserved.
Models for estimating flood losses to infrastructure are rare and their reliability is seldom investigated although infrastructure losses might contribute considerably to the overall flood losses. In this paper, an empirical modelling approach for estimating direct structural flood damage to railway infrastructure and associated financial losses is presented. Via a combination of event data, i.e. photo-documented damage on the Northern Railway in Lower Austria caused by the March River flood in 2006, and simulated flood characteristics, i.e. water levels, flow velocities and combinations thereof, the correlations between physical flood impact parameters and damage occurred to the railway track were investigated and subsequently rendered into a damage model. After calibrating the loss estimation using recorded repair costs of the Austrian Federal Railways, the model was applied to three synthetic scenarios with return periods of 30, 100 and 300 years of March River flooding. Finally, the model results are compared to depth-damage-curve-based approaches for the infrastructure sector obtained from the Rhine Atlas damage model and the Damage Scanner model. The results of this case study indicate a good performance of our two-stage model approach. However, due to a lack of independent event and damage data, the model could not yet be validated. Future research in natural risk should focus on the development of event and damage documentation procedures to overcome this significant hurdle in flood damage modelling.
In many procedures of seismic risk mitigation, ground motion simulations are needed to test systems or improve their effectiveness. For example they may be used to estimate the level of ground shaking caused by future earthquakes. Good physical models for ground motion simulation are also thought to be important for hazard assessment, as they could close gaps in the existing datasets. Since the observed ground motion in nature shows a certain variability, part of which cannot be explained by macroscopic parameters such as magnitude or position of an earthquake, it would be desirable that a good physical model is not only able to produce one single seismogram, but also to reveal this natural variability.
In this thesis, I develop a method to model realistic ground motions in a way that is computationally simple to handle, permitting multiple scenario simulations. I focus on two aspects of ground motion modelling. First, I use deterministic wave propagation for the whole frequency range – from static deformation to approximately 10 Hz – but account for source variability by implementing self-similar slip distributions and rough fault interfaces. Second, I scale the source spectrum so that the modelled waveforms represent the correct radiated seismic energy. With this scaling I verify whether the energy magnitude is suitable as an explanatory variable, which characterises the amount of energy radiated at high frequencies – the advantage of the energy magnitude being that it can be deduced from observations, even in real-time.
Applications of the developed method for the 2008 Wenchuan (China) earthquake, the 2003 Tokachi-Oki (Japan) earthquake and the 1994 Northridge (California, USA) earthquake show that the fine source discretisations combined with the small scale source variability ensure that high frequencies are satisfactorily introduced, justifying the deterministic wave propagation approach even at high frequencies. I demonstrate that the energy magnitude can be used to calibrate the high-frequency content in ground motion simulations.
Because deterministic wave propagation is applied to the whole frequency range, the simulation method permits the quantification of the variability in ground motion due to parametric uncertainties in the source description. A large number of scenario simulations for an M=6 earthquake show that the roughness of the source as well as the distribution of fault dislocations have a minor effect on the simulated variability by diminishing directivity effects, while hypocenter location and rupture velocity more strongly influence the variability. The uncertainty in energy magnitude, however, leads to the largest differences of ground motion amplitude between different events, resulting in a variability which is larger than the one observed.
For the presented approach, this dissertation shows (i) the verification of the computational correctness of the code, (ii) the ability to reproduce observed ground motions and (iii) the validation of the simulated ground motion variability. Those three steps are essential to evaluate the suitability of the method for means of seismic risk mitigation.
After the flood in 2002, the level of private precautions taken increased considerably. One contributing factor is the fact that, in general, a larger proportion of people knew that they were at risk of flooding. The best level of precaution was found before the flood events in 2006 and 2011. The main reason for this might be that residents had more experience with flooding than residents affected in 2005 or 2010. Yet, overall, flood experience and knowledge did not necessarily result in building retrofitting or flood-proofing measures, which are considered as mitigating damages most effectively. Hence, investments still need to be stimulated in order to reduce future damage more efficiently.
Computer-based simulation models are frequently used in hydrological research and engineering but also in other fields of environmental sciences. New case studies often require existing model concepts to be adapted. Extensions may be necessary due to the peculiarities of the studied natural system or subtleties of anthropogenic control. In other cases, simplifications must be made in response to scarce data, incomplete knowledge, or restrictions set by the spatio-temporal scale of application. This paper introduces an open-source modeling framework called ECHSE designed to cope with the above-mentioned challenges. It provides a lightweight infrastructure for the rapid development of new, reusable simulation tools and, more importantly, the safe modification of existing formulations. ECHSE-based models treat the simulated system as a collection of interacting objects. Although feedbacks are generally supported, the majority of the objects' interactions is expected to be of the feed-forward type. Therefore, the ECHSE software is particularly useful in the context of hydrological catchment modeling. Conversely, it is unsuitable, e.g., for fully hydrodynamic simulations and groundwater flow modeling. The focus of the paper is put on a comprehensible outline of the ECHSE's fundamental concepts and limitations. For the purpose of illustration, a specific, ECHSE-based solution for hydrological catchment modeling is presented which has undergone testing in a number of river basins. (C) 2015 Elsevier Ltd. All rights reserved.