Refine
Year of publication
- 2016 (284) (remove)
Document Type
- Article (216)
- Doctoral Thesis (35)
- Other (14)
- Postprint (13)
- Review (5)
- Habilitation Thesis (1)
Language
- English (284) (remove)
Keywords
- erosion (5)
- Pollen (4)
- climate change (4)
- August 2002 flood (3)
- Body waves (3)
- Central Europe (3)
- Floods Directive (3)
- Himalaya (3)
- Himalayas (3)
- Inversion (3)
Institute
- Institut für Geowissenschaften (284) (remove)
The 11 July 1889 Chilik earthquake (M-w 8.0-8.3) forms part of a remarkable sequence of large earthquakes in the late nineteenth and early twentieth centuries in the northern Tien Shan. Despite its importance, the source of the 1889 earthquake remains unknown, though the macroseismic epicenter is sited in the Chilik valley, similar to 100 km southeast of Almaty, Kazakhstan (similar to 2 million population). Several short fault segments that have been inferred to have ruptured in 1889 are too short on their own to account for the estimated magnitude. In this paper we perform detailed surveying and trenching of the similar to 30 km long Saty fault, one of the previously inferred sources, and find that it was formed in a single earthquake within the last 700 years, involving surface slip of up to 10 m. The scarp-forming event, likely to be the 1889 earthquake, was the only surface-rupturing event for at least 5000 years and potentially for much longer. From satellite imagery we extend the mapped length of fresh scarps within the 1889 epicentral zone to a total of similar to 175 km, which we also suggest as candidate ruptures from the 1889 earthquake. The 175 km of rupture involves conjugate oblique left-lateral and right-lateral slip on three separate faults, with step overs of several kilometers between them. All three faults were essentially invisible in the Holocene geomorphology prior to the last slip. The recurrence interval between large earthquakes on any of these faults, and presumably on other faults of the Tien Shan, may be longer than the timescale over which the landscape is reset, providing a challenge for delineating sources of future hazard.
This case study evaluates the suitability of radar-based quantitative precipitation estimates (QPEs) for the simulation of streamflow in the Marikina River Basin (MRB), the Philippines. Hourly radar-based QPEs were produced from reflectivity that had been observed by an S-band radar located about 90 km from the MRB. Radar data processing and precipitation estimation were carried out using the open source library wradlib. To assess the added value of the radar-based QPE, we used spatially interpolated rain gauge observations (gauge-only (GO) product) as a benchmark. Rain gauge observations were also used to quantify rainfall estimation errors at the point scale. At the point scale, the radar-based QPE outperformed the GO product in 2012, while for 2013, the performance was similar. For both periods, estimation errors substantially increased from daily to the hourly accumulation intervals. Despite this fact, both rainfall estimation methods allowed for a good representation of observed streamflow when used to force a hydrological simulation model of the MRB. Furthermore, the results of the hydrological simulation were consistent with rainfall verification at the point scale: the radar-based QPE performed better than the GO product in 2012, and equivalently in 2013. Altogether, we could demonstrate that, in terms of streamflow simulation, the radar-based QPE can perform as good as or even better than the GO product - even for a basin such as the MRB which has a comparatively dense rain gauge network. This suggests good prospects for using radar-based QPE to simulate and forecast streamflow in other parts of the Philippines where rain gauge networks are not as dense.
Subsurface microbial communities undertake many terminal electron-accepting processes, often simultaneously. Using a tritium-based assay, we measured the potential hydrogen oxidation catalyzed by hydrogenase enzymes in several subsurface sedimentary environments (Lake Van, Barents Sea, Equatorial Pacific, and Gulf of Mexico) with different predominant electron-acceptors. Hydrogenases constitute a diverse family of enzymes expressed by microorganisms that utilize molecular hydrogen as a metabolic substrate, product, or intermediate. The assay reveals the potential for utilizing molecular hydrogen and allows qualitative detection of microbial activity irrespective of the predominant electron-accepting process. Because the method only requires samples frozen immediately after recovery, the assay can be used for identifying microbial activity in subsurface ecosystems without the need to preserve live material. We measured potential hydrogen oxidation rates in all samples from multiple depths at several sites that collectively span a wide range of environmental conditions and biogeochemical zones. Potential activity normalized to total cell abundance ranges over five orders of magnitude and varies, dependent upon the predominant terminal electron acceptor. Lowest per-cell potential rates characterize the zone of nitrate reduction and highest per-cell potential rates occur in the methanogenic zone. Possible reasons for this relationship to predominant electron acceptor include (i) increasing importance of fermentation in successively deeper biogeochemical zones and (ii) adaptation of H(2)ases to successively higher concentrations of H-2 in successively deeper zones.
Subsurface microbial communities undertake many terminal electron-accepting processes, often simultaneously. Using a tritium-based assay, we measured the potential hydrogen oxidation catalyzed by hydrogenase enzymes in several subsurface sedimentary environments (Lake Van, Barents Sea, Equatorial Pacific, and Gulf of Mexico) with different predominant electron-acceptors. Hydrogenases constitute a diverse family of enzymes expressed by microorganisms that utilize molecular hydrogen as a metabolic substrate, product, or intermediate. The assay reveals the potential for utilizing molecular hydrogen and allows qualitative detection of microbial activity irrespective of the predominant electron-accepting process. Because the method only requires samples frozen immediately after recovery, the assay can be used for identifying microbial activity in subsurface ecosystems without the need to preserve live material. We measured potential hydrogen oxidation rates in all samples from multiple depths at several sites that collectively span a wide range of environmental conditions and biogeochemical zones. Potential activity normalized to total cell abundance ranges over five orders of magnitude and varies, dependent upon the predominant terminal electron acceptor. Lowest per-cell potential rates characterize the zone of nitrate reduction and highest per-cell potential rates occur in the methanogenic zone. Possible reasons for this relationship to predominant electron acceptor include (i) increasing importance of fermentation in successively deeper biogeochemical zones and (ii) adaptation of H(2)ases to successively higher concentrations of H-2 in successively deeper zones.
Hydrologic regionalization deals with the investigation of homogeneity in watersheds and provides a classification of watersheds for regional analysis. The classification thus obtained can be used as a basis for mapping data from gauged to ungauged sites and can improve extreme event prediction. This paper proposes a wavelet power spectrum (WPS) coupled with the self-organizing map method for clustering hydrologic catchments. The application of this technique is implemented for gauged catchments. As a test case study, monthly streamflow records observed at 117 selected catchments throughout the western United States from 1951 through 2002. Further, based on WPS of each station, catchments are classified into homogeneous clusters, which provides a representative WPS pattern for the streamflow stations in each cluster.
Analysis of time-lapse ground-penetrating radar (GPR) data can provide information regarding subsurface hydrological processes, such as preferential flow. However, the analysis of time-lapse data is often limited by data quality; for example, for noisy input data, the interpretation of difference images is often difficult. Motivated by modern image-processing tools, we have developed two robust GPR attributes, which allow us to distinguish amplitude (contrast similarity) and time-shift (structural similarity) variations related to differences between individual time-lapse GPR data sets. We tested and evaluated our attributes using synthetic data of different complexity. Afterward, we applied them to a field data example, in which subsurface flow was induced by an artificial rainfall event. For all examples, we identified our structural similarity attribute to be a robust measure for highlighting time-lapse changes also in data with low signal-to-noise ratios. We determined that our new attribute-based workflow is a promising tool to analyze time-lapse GPR data, especially for imaging subsurface hydrological processes.
Here, we study the 3-D subduction initiation process induced by the interaction between a hot thermochemical mantle plume and oceanic lithosphere using thermo-mechanical viscoplastic finite difference marker-in-cell models. Our numerical modeling results show that self-sustaining subduction is induced by plume-lithosphere interaction when the plume is sufficiently buoyant, the oceanic lithosphere is sufficiently old and the plate is weak enough to allow the buoyant plume to. pass through it. Subduction initiation occurs following penetration of the lithosphere by the hot plume and the downward displacement of broken, nearly circular segments of lithosphere (proto-slabs) as a result of partially molten plume rocks overriding the proto-slabs. Our experiments show four different deformation regimes in response to plume-lithosphere interaction: a) self-sustaining subduction initiation, in which subduction becomes self-sustaining; b) frozen subduction initiation, in which subduction stops at shallow depths; c) slab break-off, in which the subducting circular slab breaks off soon after formation; and d) plume underplating, in which the plume does not pass through the lithosphere and instead spreads beneath it (i.e., failed subduction initiation). These regimes depend on several parameters, such as the size, composition, and temperature of the plume, the brittle/plastic strength and age of the oceanic lithosphere, and the presence/absence of lithospheric heterogeneities. The results show that subduction initiates and becomes self-sustaining when the lithosphere is older than 10 Myr and the non dimensional ratio of the plume buoyancy force and lithospheric strength above the plume is higher than approximately 2. The outcomes of our numerical experiments are applicable for subduction initiation in the modern and Precambrian Earth and for the origin of plume-related corona structures on Venus. (C) 2016 Elsevier B.V. All rights reserved.
Intracontinental deformation usually is a result of tectonic forces associated with distant plate collisions. In general, the evolution of mountain ranges and basins in this environment is strongly controlled by the distribution and geometries of preexisting structures. Thus, predictive models usually fail in forecasting the deformation evolution in these kinds of settings. Detailed information on each range and basin-fill is vital to comprehend the evolution of intracontinental mountain belts and basins. In this dissertation, I have investigated the complex Cenozoic tectonic evolution of the western Tien Shan in Central Asia, which is one of the most active intracontinental ranges in the world. The work presented here combines a broad array of datasets, including thermo- and geochronology, paleoenvironmental interpretations, sediment provenance and subsurface interpretations in order to track changes in tectonic deformation. Most of the identified changes are connected and can be related to regional-scale processes that governed the evolution of the western Tien Shan.
The NW-SE trending Talas-Fergana fault (TFF) separates the western from the central Tien Shan and constitutes a world-class example of the influence of preexisting anisotropies on the subsequent structural development of a contractile orogen. While to the east most of ranges and basins have a sub-parallel E-W trend, the triangular-shaped Fergana basin forms a substantial feature in the western Tien Shan morphology with ranges on all three sides. In this thesis, I present 55 new thermochronologic ages (apatite fission track and zircon (U-Th)/He)) used to constrain exhumation histories of several mountain ranges in the western Tien Shan. At the same time, I analyzed the Fergana basin-fill looking for progressive changes in sedimentary paleoenvironments, source areas and stratal geometrical configurations in the subsurface and outcrops.
The data presented in this thesis suggests that low cooling rates (<1°C Myr-1), calm depositional environments, and low depositional rates (<10 m Myr-1) were widely distributed across the western Tien Shan, describing a quiescent tectonic period throughout the Paleogene. Increased cooling rates in the late Cenozoic occurred diachronously and with variable magnitudes in different ranges. This rapid cooling stage is interpreted to represent increased erosion caused by active deformation and constrains the onset of Cenozoic deformation in the western Tien Shan. Time-temperature histories derived from the northwestern Tien Shan samples show an increase in cooling rates by ~25 Ma. This event is correlated with a synchronous pulse
iv
in the South Tien Shan. I suggest that strike-slip motion along the TFF commenced at the Oligo-Miocene boundary, facilitating CCW rotation of the Fergana basin and enabling exhumation of the linked horsetail splays. Higher depositional rates (~150 m Myr-1) in the Oligo-Miocene section (Massaget Fm.) of the Fergana basin suggest synchronous deformation in the surrounding ranges. The central Alai Range also experienced rapid cooling around this time, suggesting that the onset of intramontane basin fragmentation and isolation is coeval. These results point to deformation starting simultaneously in the late Oligocene – early Miocene in geographically distant mountain ranges. I suggest that these early uplifts are controlled by reactivated structures (like the TFF), which are probably the frictionally weakest and most-suitably oriented for accommodating and transferring N-S horizontal shortening along the western Tien Shan.
Afterwards, in the late Miocene (~10 Ma), a period of renewed rapid cooling affected the Tien Shan and most mountain ranges and inherited structures started to actively deform. This episode is widely distributed and an increase in exhumation is interpreted in most of the sampled ranges. Moreover, the Pliocene section in the basin subsurface shows the higher depositional rates (>180 m Myr-1) and higher energy facies. The deformation and exhumation increase further contributed to intramontane basin partitioning. Overall, the interpretation is that the Tien Shan and much of Central Asia suffered a global increase in the rate of horizontal crustal shortening. Previously, stress transfer along the rigid Tarim block or Pamir indentation has been proposed to account for Himalayan hinterland deformation. However, the extent of the episode requires a different and broader geodynamic driver.
This review presents a compositional database of primary anatectic granitoid magmas, entirely based on melt inclusions (MI) in high-grade metamorphic rocks. Although MI are well known to igneous petrologists and have been extensively studied in intrusive and extrusive rocks, MI in crustal rocks that have undergone anatexis (migmatites and granulites) are a novel subject of research. They are generally trapped along the heating path by peritectic phases produced by incongruent melting reactions. Primary MI in high-grade metamorphic rocks are small, commonly 5-10 pm in diameter, and their most common mineral host is peritectic garnet. In most cases inclusions have crystallized into a cryptocrystalline aggregate and contain a granitoid phase assemblage (nanogranitoid inclusions) with quartz, K-feldspar, plagioclase, and one or two mica depending on the particular circumstances. After their experimental remelting under high-confining pressure, nanogranitoid MI can be analyzed combining several techniques (EMP, LA-ICP-MS, NanoSIMS, Raman). The trapped melt is granitic and metaluminous to peraluminous, and sometimes granodioritic, tonalitic, and trondhjemitic in composition, in agreement with the different P-T-a(H2o) conditions of melting and protolith composition, and overlap the composition of experimental glasses produced at similar conditions. Being trapped along the up-temperature trajectory as opposed to classic MI in igneous rocks formed during down-temperature magma crystallization fundamental information provided by nanogranitoid MI is the pristine composition of the natural primary anatectic melt for the specific rock under investigation. So far similar to 600 nanogranitoid MI, coming from several occurrences from different geologic and geodynamic settings and ages, have been characterized. Although the compiled MI database should be expanded to other potential sources of crustal magmas, MI data collected so far can be already used as natural "starting-point" compositions to track the processes involved in formation and evolution of granitoid magmas.
New age data have been obtained to time constrain the recent Quaternary volcanism of El Hierro (Canary Islands) and to estimate its recurrence rate. We have carried out Ar-40/Ar-39 geochronology on samples spanning the entire volcanostratigraphic sequence of the island and C-14 geochronology on the most recent eruption on the northeast rift of the island: 2280 +/- 30 yr BP. We combine the new absolute data with a revision of published ages onshore, some of which were identified through geomorphological criteria (relative data). We present a revised and updated chronology of volcanism for the last 33 ka that we use to estimate the maximum eruptive recurrence of the island. The number of events per year determined is 9.7 x 10(-4) for the emerged part of the island, which means that, as a minimum, one eruption has occurred approximately every 1000 years. This highlights the need of more geochronological data to better constrain the eruptive recurrence of El Hierro. (C) 2015 Elsevier Ltd. All rights reserved.
Terrigenous sediment supply, marine transport, and depositional processes along tectonically active margins are key to decoding turbidite successions as potential archives of climatic and seismic forcings. Sequence stratigraphic models predict coarse-grained sediment delivery to deep-marine sites mainly during sea-level fall and lowstand. Marine siliciclastic deposition during transgressions and highstands has been attributed to sustained connectivity between terrigenous sources and marine sinks facilitated by narrow shelves. To decipher the controls on Holocene highstand turbidite deposition, we analyzed 12 sediment cores from spatially discrete, coeval turbidite systems along the Chile margin (29 degrees-40 degrees S) with changing climatic and geomorphic characteristics but uniform changes in sea level. Sediment cores from intraslope basins in north-central Chile (29 degrees-33 degrees S) offshore a narrow to absent shelf record a shut-off of turbidite deposition during the Holocene due to postglacial aridification. In contrast, core sites in south-central Chile (36 degrees-40 degrees S) offshore a wide shelf record frequent turbidite deposition during highstand conditions. Two core sites are linked to the Biobio river-canyon system and receive sediment directly from the river mouth. However, intraslope basins are not connected via canyons to fluvial systems but yield even higher turbidite frequencies. High sediment supply combined with a wide shelf and an undercurrent moving sediment toward the shelf edge appear to control Holocene turbidite sedimentation and distribution. Shelf undercurrents may play an important role in lateral sediment transport and supply to the deep sea and need to be accounted for in sediment-mass balances.
Although the climate development over the Holocene in the Northern Hemisphere is well known, palaeolimnological climate reconstructions reveal spatiotemporal variability in northern Eurasia. Here we present a multi-proxy study from north-eastern Siberia combining sediment geochemistry, and diatom and pollen data from lake-sediment cores covering the last 38,000 cal. years. Our results show major changes in pyrite content and fragilarioid diatom species distributions, indicating prolonged seasonal lake-ice cover between similar to 13,500 and similar to 8900 cal. years BP and possibly during the 8200 cal. years BP cold event. A pollen-based climate reconstruction generated a mean July temperature of 17.8 degrees C during the Holocene Thermal Maximum (HTM) between similar to 8900 and similar to 4500 cal. years BP. Naviculoid diatoms appear in the late Holocene indicating a shortening of the seasonal ice cover that continues today. Our results reveal a strong correlation between the applied terrestrial and aquatic indicators and natural seasonal climate dynamics in the Holocene. Planktonic diatoms show a strong response to changes in the lake ecosystem due to recent climate warming in the Anthropocene. We assess other palaeolimnological studies to infer the spatiotemporal pattern of the HTM and affirm that the timing of its onset, a difference of up to 3000 years from north to south, can be well explained by climatic teleconnections. The westerlies brought cold air to this part of Siberia until the Laurentide ice sheet vanished 7000 years ago. The apparent delayed ending of the HTM in the central Siberian record can be ascribed to the exceedance of ecological thresholds trailing behind increases in winter temperatures and decreases in contrast in insolation between seasons during the mid to late Holocene as well as lacking differentiation between summer and winter trends in paleolimnological reconstructions. (C) 2015 Elsevier Ltd. All rights reserved.
The Pamirs represent the indented westward continuation of the northern margin of the Tibetan Plateau, dividing the Tarim and Tajik basins. Their evolution may be a key factor influencing aridification of the Asian interior, yet the tectonics of the Pamir Salient are poorly understood. We present a provenance study of the Aertashi section, a Paleogene to late Neogene clastic succession deposited in the Tarim basin to the north of the NW margin of Tibet (the West Kunlun) and to the east of the Pamirs. Our detrital zircon U-Pb ages coupled with zircon fission track, bulk rock Sm-Nd, and petrography data document changes in contributing source terranes during the Oligocene to Miocene, which can be correlated to regional tectonics. We propose a model for the evolution of the Pamir and West Kunlun (WKL), in which the WKL formed topography since at least similar to 200 Ma. By similar to 25 Ma, movement along the Pamir-bounding faults such as the Kashgar-Yecheng Transfer System had commenced, marking the onset of Pamir indentation into the Tarim-Tajik basin. This is coincident with basinward expansion of the northern WKL margin, which changed the palaeodrainage pattern within the Kunlun, progressively cutting off the more southerly WKL sources from the Tarim basin. An abrupt change in the provenance and facies of sediments at Aertashi has a maximum age of 14 Ma; this change records when the Pamir indenter had propagated sufficiently far north that the North Pamir was now located proximal to the Aertashi region.
The South American Andes are frequently exposed to intense rainfall events with varying moisture sources and precipitation-forming processes. In this study, we assess the spatiotemporal characteristics and geographical origins of rainfall over the South American continent. Using high-spatiotemporal resolution satellite data (TRMM 3B42 V7), we define four different types of rainfall events based on their (1) high magnitude, (2) long temporal extent, (3) large spatial extent, and (4) high magnitude, long temporal and large spatial extent combined. In a first step, we analyze the spatiotemporal characteristics of these events over the entire South American continent and integrate their impact for the main Andean hydrologic catchments. Our results indicate that events of type 1 make the overall highest contributions to total seasonal rainfall (up to 50%). However, each consecutive episode of the infrequent events of type 4 still accounts for up to 20% of total seasonal rainfall in the subtropical Argentinean plains. In a second step, we employ complex network theory to unravel possibly non-linear and long-ranged climatic linkages for these four event types on the high-elevation Altiplano-Puna Plateau as well as in the main river catchments along the foothills of the Andes. Our results suggest that one to two particularly large squall lines per season, originating from northern Brazil, indirectly trigger large, long-lasting thunderstorms on the Altiplano Plateau. In general, we observe that extreme rainfall in the catchments north of approximately 20 degrees S typically originates from the Amazon Basin, while extreme rainfall at the eastern Andean foothills south of 20 degrees S and the Puna Plateau originates from southeastern South America.
In contrast to recent advances in projecting sea levels, estimations about the economic impact of sea level rise are vague. Nonetheless, they are of great importance for policy making with regard to adaptation and greenhouse-gas mitigation. Since the damage is mainly caused by extreme events, we propose a stochastic framework to estimate the monetary losses from coastal floods in a confined region. For this purpose, we follow a Peak-over-Threshold approach employing a Poisson point process and the Generalised Pareto Distribution. By considering the effect of sea level rise as well as potential adaptation scenarios on the involved parameters, we are able to study the development of the annual damage. An application to the city of Copenhagen shows that a doubling of losses can be expected from a mean sea level increase of only 11 cm. In general, we find that for varying parameters the expected losses can be well approximated by one of three analytical expressions depending on the extreme value parameters. These findings reveal the complex interplay of the involved parameters and allow conclusions of fundamental relevance. For instance, we show that the damage typically increases faster than the sea level rise itself. This in turn can be of great importance for the assessment of sea level rise impacts on the global scale. Our results are accompanied by an assessment of uncertainty, which reflects the stochastic nature of extreme events. While the absolute value of uncertainty about the flood damage increases with rising mean sea levels, we find that it decreases in relation to the expected damage.
The functional form of empirical response spectral ground-motion prediction equations (GMPEs) is often derived using concepts borrowed from Fourier spectral modeling of ground motion. As these GMPEs are subsequently calibrated with empirical observations, this may not appear to pose any major problems in the prediction of ground motion for a particular earthquake scenario. However, the assumption that Fourier spectral concepts persist for response spectra can lead to undesirable consequences when it comes to the adjustment of response spectral GMPEs to represent conditions not covered in the original empirical data set. In this context, a couple of important questions arise, for example, what are the distinctions and/or similarities between Fourier and response spectra of ground motions? And, if they are different, then what is the mechanism responsible for such differences and how do adjustments that are made to Fourier amplitude spectrum (FAS) manifest in response spectra? The present article explores the relationship between the Fourier and response spectrum of ground motion by using random vibration theory (RVT). With a simple Brune (1970, 1971) source model, RVT-generated acceleration spectra for a fixed magnitude and distance scenario are used. The RVT analyses reveal that the scaling of low oscillator-frequency response spectral ordinates can be treated as being equivalent to the scaling of the corresponding Fourier spectral ordinates. However, the high oscillator-frequency response spectral ordinates are controlled by a rather wide band of Fourier spectral ordinates. In fact, the peak ground acceleration, counter to the popular perception that it is a reflection of the high-frequency characteristics of ground motion, is controlled by the entire Fourier spectrum of ground motion. Additionally, this article demonstrates how an adjustment made to FAS is similar or different to the same adjustment made to response spectral ordinates. For this purpose, two cases: adjustments to the stress parameter (Delta sigma) (source term), and adjustments to the attributes reflecting site response (V-S - kappa(0)) are considered.
Fossil oyster shells are well-suited to provide palaeotemperature proxies from geologic to seasonal timescales due to their ubiquitous occurrence from Triassic to Quaternary sediments, the seasonal nature of their shell growth and their relative strong resistance to post-mortem alteration. However, the common use to translate calcitic oxygen isotopes into palaeotemperatures is challenged by uncertainties in accounting for past seawater delta O-18, especially in shallow coastal environment where oysters calcify. In principle, the Mg/Ca ratio in oyster shells can provide an alternative palaeothermometer. Several studies provided temperature calibrations for this potential proxy based on modem species, nevertheless their application to palaeo-studies remains hitherto unexplored. Here, we show that past temperature variability in seawater can be obtained from Mg/Ca analyses from selected fossil oyster species and specimens. High-resolution Mg/Ca profiles, combined with delta O-18, were obtained along 41 fossil oyster shells of seven different species from the Palaeogene Proto-Paratethys sea (Central Asia) found in similar as well as different depositional age and environments providing comparison. Suitable Mg/Ca profiles, defined by continuous cyclicity and reproducibility within one shell, are found to be consistent for specimens of the same species but differ systematically between species, implying a dominant species-specific effect on the Mg/Ca signal. Two species studied here (Ostrea (Turkostrea) strictiplicata and Sokolowia buhsii) provide an excellent proxy for palaeoclimate reconstruction from China to Europe in Palaeogene marine sediments. More generally, the protocol developed here can be applied to identify other fossil oyster species suitable for palaeoclimate reconstructions. (C) 2015 Elsevier B.V. All rights reserved.
Synergistic applications based on integrated hyperspectral and lidar data are receiving a growing interest from the remote-sensing community. A prerequisite for the optimum sensor fusion of hyperspectral and lidar data is an accurate geometric coalignment. The simple unadjusted integration of lidar elevation and hyperspectral reflectance causes a substantial loss of information and does not exploit the full potential of both sensors. This paper presents a novel approach for the geometric coalignment of hyperspectral and lidar airborne data, based on their respective adopted return intensity information. The complete approach incorporates ray tracing and subpixel procedures in order to overcome grid inherent discretization. It aims at the correction of extrinsic and intrinsic (camera resectioning) parameters of the hyperspectral sensor. In additional to a tie-point-based coregistration, we introduce a ray-tracing-based back projection of the lidar intensities for area-based cost aggregation. The approach consists of three processing steps. First is a coarse automatic tie-point-based boresight alignment. The second step coregisters the hyperspectral data to the lidar intensities. Third is a parametric coalignment refinement with an area-based cost aggregation. This hybrid approach of combining tie-point features and area-based cost aggregation methods for the parametric coregistration of hyperspectral intensity values to their corresponding lidar intensities results in a root-mean-square error of 1/3 pixel. It indicates that a highly integrated and stringent combination of different coalignment methods leads to an improvement of the multisensor coregistration.
The humid tropics are the region with the highest rate of land-cover change worldwide. Especially prevalent is the deforestation of old-growth tropical forests to create space for cattle pastures and soybean fields.
The regional water cycle is influenced by vegetation cover in various ways. Especially evapotranspiration considerably contributes to water vapor content in the lower atmosphere. Besides active transpiration by plants, evaporation from wetted plant surfaces further known as interception loss is an important supply of water vapor. Changes in interception loss due to change in land cover and the related consequences on the regional water cycle in the humid tropics of Latin America are the research focus of my thesis. (1) In an experimental setup I assess differences in interception loss between an old-growth tropical forest and a soybean plantation. (2) In a modeling study, I examine interception losses of these two vegetation types compared to a younger secondary forest with the use of the Gash interception model, including an uncertainty analysis for the estimation of the necessary model parameters. (3) Studying the water balance of a 192-km² catchment I disentangle the influences of changes in land cover and climatic factors on interception loss.
The three different research sites in my thesis represent a currently typical spectrum for land-cover changes in Latin America. In the first example I study the consequences of deforestation of transitional forest, which forms the transition from the Brazilian tree savanna (cerrado) to tropical rain forest, for the establishment of soybean fields in the southern Amazon basin. The second study site is a young secondary forest within the “Agua Salud” project area in Panama as an example of reforestation of former pastures. The third study site is the Cirí Grande river catchment which comprises a mixture of young and old forests as well as pastures, which is typical for the southern sub-catchments of the Panama Canal.
The experimental approach consists of the indirect estimation of interception loss by measuring throughfall and stem flow. For the first experimental study I measured throughfall as well as stem flow manually. Measurements of the leaf area index of the two land covers do not show distinct differences; hence it could not serve as an explanation for the differences in the measured interception loss. The considerably higher interception loss at the soybean field is attributed to a possible underestimation of stemflow but also to the stronger ventilation within the well-structured plant rows causing higher evaporation rates. This situation is valid only for two months of the rainy season, when soybean plants are fully developed. In the annual balance evapotranspiration at the soybean site is clearly less than at the forest site, accelerating the development of fast runoff components and consequently discharge. In the medium term, a reduction of water availability in the study area can be expected.
For the modeling study, throughfall in a young secondary forest is sampled automatically. The resulting temporally high-resolution dataset allows the distinction between different precipitation and interception events. The core of this study is the sensitivity and uncertainty analysis of the Gash interception model parameters and the consequences for its results. Canopy storage capacity plays a key role for the model and parameter uncertainty. With increasing storage capacity uncertainty in parameter delineation also increases. Evaporation rate as the driving component of the interception process incorporates in this context the largest parameter uncertainty. Depending on the selected method for parameter estimation, parameter values may vary tremendously.
In the third study, I analyze the influence of interception loss on the water balance of the Cirí Grande catchment, incorporating the interlinked effects of temperature, precipitation and changes of the land use mosaic using the SWAT (soil water assessment tool) model. Constructing several land-cover scenarios I assess their influence on the catchment’s discharge. The results show that land-cover change exerts only a small influence on annual discharge in the Cirí Grande catchment whereas an increase in temperature markedly influences evapotranspiration. The temperature-induced larger transpiration and interception loss balances the simultaneous increase in annual precipitation, such that the resulting changes in annual discharge are negligible.
The results of the three studies show the considerable effect of land cover on interception. However, the magnitude of this effect can be masked by changes in local conditions, especially by an increase in temperature. Hence, the results cannot be transferred easily between the different study sites. For modeling purposes, this means that measurements of vegetation characteristics as well as interception loss at the respective sites are indispensable.
Carboniferous metagranites with U-Pb zircon crystallization ages of 331-315 Ma crop out in the Afyon zone in the northern margin of the Anatolide-Tauride Block, which is commonly regarded as part of Gondwana during the Late Palaeozoic. They are peraluminous, calc-alkaline and are characterized by increase in Rb and Ba, decrease in Nb-Ta, and enrichment in Sr and high LILE/HFSE ratios compatible with a continental arc setting. The metagranites intrude a metasedimentary sequence of phyllite, metaquartzite and marble; both the Carboniferous metagranites and metasedimentary rocks are overlain unconformably by Lower Triassic metaconglomerates, metavolcanics and Upper Triassic to Cretaceous recrystallized limestones. The low-grade metamorphism and deformation occurred at the Cretaceous-Tertiary boundary. There is no evidence for Carboniferous deformation and metamorphism in the region. Carboniferous arc-type granites and previously described Carboniferous subduction-accretion complexes on the northern margin of the Anatolide-Tauride Block suggest southward subduction of Paleotethys under Gondwana during the Carboniferous. Considering the Variscan-related arc granites in Pelagonian and Sakarya zones on the active southern margin of Laurasia, a dual subduction of Paleotethys can be envisaged between Early Carboniferous and Late Permian. However, the southward subduction was short-lived and by the Late Permian the Gondwana margin became passive. (C) 2016 Elsevier B.V. All rights reserved.
Numerous (meta-)gabbroic dikes or stocks occur within the latest Neoproterozoic-early Cambrian series of the Menderes Massif (Anatolide-Tauride Block, western Turkey). These well-preserved rocks were locally converted into eclogitic metagabbros and garnet amphibolites along the contacts or shear zones. Both bulk-rock composition and compositions of igneous clinopyroxenes suggest continental tholeiitic affinity. U-Pb dating of igneous zircons from gabbroic rocks yielded a mean age of 563 +/- 1 Ma (2 sigma), indicating emplacement during the latest Neoproterozoic (Ediacaran). On the other hand, rims of zircons from eclogitic metagabbro gave 535 +/- 3 Ma (2 sigma) (early Cambrian), in addition to 558 +/- 3 Ma (2 sigma) obtained from the igneous core of zircons. These ages are interpreted as the time of high-P metamorphism and crystallization age of gabbroic protolith, respectively. Given the estimated paleogeographic position of the Anatolide-Tauride Block during the late Neoproterozoic and early Cambrian, this orogenic event can be spatially and temporally related to the northward continuity of 600-500 Ma orogenic event (Malagasy/Kuunga orogeny) extending from western margin of India, Madagascar, via Arabia up to northern margin of Gondwana beneath thick Phanerozoic cover series in Arabian Peninsula. Therefore, the high-P evolution of the basement of the Menderes Massif and associated basic intrusions can be interpreted to mark the latest stages of consumption of the basin/oceanic branches and final amalgamation of the Gondwana during the late Neoproterozoic-early Cambrian around the Arabian region. (C) 2015 International Association for Gondwana Research. Published by Elsevier B.V. All rights reserved.
By reducing current velocity, tidal marsh vegetation can diminish storm surges and storm waves. Conversely, currents often exert high mechanical stresses onto the plants and hence affect vegetation structure and plant characteristics. In our study, we aim at analysing this interaction from both angles. On the one hand, we quantify the reduction of current velocity by Bolboschoenus maritimus, and on the other hand, we identify functional traits of B. maritimus’ ramets along environmental gradients. Our results show that tidal marsh vegetation is able to buffer a large proportion of the flow velocity at currents under normal conditions. Cross-shore current velocity decreased with distance from the marsh edge and was reduced by more than 50% after 15 m of vegetation. We were furthermore able to show that plants growing at the marsh edge had a significantly larger diameter than plants from inside the vegetation. We found a positive correlation between plant thickness and cross-shore current which could provide an adaptive value in habitats with high mechanical stress. With the adapted morphology of plants growing at the highly exposed marsh edge, the entire vegetation belt is able to better resist the mechanical stress of high current velocities. This self-adaptive effect thus increases the ability of B. maritimus to grow and persist in the pioneer zone and may hence better contribute to ecosystem-based coastal protection by reducing current velocity.
Extreme hydro-meteorological events, such as severe droughts or heavy rainstorms, constitute primary manifestations of climate variability and exert a critical impact on the natural environment and human society. This is particularly true for high-mountain areas, such as the eastern flank of the southern Central Andes of NW Argentina, a region impacted by deep convection processes that form the basis of extreme events, often resulting in floods, a variety of mass movements, and hillslope processes. This region is characterized by pronounced E-W gradients in topography, precipitation, and vegetation cover, spanning low to medium-elevation, humid and densely vegetated areas to high-elevation, arid and sparsely vegetated environments. This strong E-W gradient is mirrored by differences in the efficiency of surface processes, which mobilize and transport large amounts of sediment through the fluvial system, from the steep hillslopes to the intermontane basins and further to the foreland. In a highly sensitive high-mountain environment like this, even small changes in the spatiotemporal distribution, magnitude and rates of extreme events may strongly impact environmental conditions, anthropogenic activity, and the well-being of mountain communities and beyond. However, although the NW Argentine Andes comprise the catchments for the La Plata river that traverses one of the most populated and economically relevant areas of South America, there are only few detailed investigations of climate variability and extreme hydro-meteorological events.
In this thesis, I focus on deciphering the spatiotemporal variability of rainfall and river discharge, with particular emphasis on extreme hydro-meteorological events in the subtropical southern Central Andes of NW Argentina during the past seven decades. I employ various methods to assess and quantify statistically significant trend patterns of rainfall and river discharge, integrating high-quality daily time series from gauging stations (40 rainfall and 8 river discharge stations) with gridded datasets (CPC-uni and TRMM 3B42 V7), for the period between 1940 and 2015. Evidence for a general intensification of the hydrological cycle at intermediate elevations (~ 0.5 – 3 km asl) at the eastern flank of the southern Central Andes is found both from rainfall and river-discharge time-series analysis during the period from 1940 to 2015. This intensification is associated with the increase of the annual total amount of rainfall and the mean annual discharge. However, most pronounced trends are found at high percentiles, i.e. extreme hydro-meteorological events, particularly during the wet season from December to February.An important outcome of my studies is the recognition of a rapid increase in the amount of river discharge during the period between 1971 and 1977, most likely linked to the 1976-77 global climate shift, which is associated with the North Pacific Ocean sea surface temperature variability. Interestingly, after this rapid increase, both rainfall and river discharge decreased at low and intermediate elevations along the eastern flank of the Andes. In contrast, during the same time interval, at high elevations, extensive areas on the arid Puna de Atacama plateau have recorded increasing annual rainfall totals. This has been associated with more intense extreme hydro-meteorological events from 1979 to 2014. This part of the study reveals that low-, intermediate, and high-elevation sectors in the Andes of NW Argentina respond differently to changing climate conditions.
Possible forcing mechanisms of the pronounced hydro-meteorological variability observed in the study area are also investigated. For the period between 1940 and 2015, I analyzed modes of oscillation of river discharge from small to medium drainage basins (102 to 104 km2), located on the eastern flank of the orogen. First, I decomposed the relevant monthly time series using the Hilbert-Huang Transform, which is particularly appropriate for non-stationary time series that result from non-linear natural processes. I observed that in the study region discharge variability can be described by five quasi-periodic oscillatory modes on timescales varying from 1 to ~20 years. Secondly, I tested the link between river-discharge variations and large-scale climate modes of variability, using different climate indices, such as the BEST ENSO (Bivariate El Niño-Southern Oscillation Time-series) index. This analysis reveals that, although most of the variance on the annual timescale is associated with the South American Monsoon System, a relatively large part of river-discharge variability is linked to Pacific Ocean variability (PDO phases) at multi-decadal timescales (~20 years). To a lesser degree, river discharge variability is also linked to the Tropical South Atlantic (TSA) sea surface temperature anomaly at multi-annual timescales (~2-5 years).
Taken together, these findings exemplify the high degree of sensitivity of high-mountain environments with respect to climatic variability and change. This is particularly true for the topographic transitions between the humid, low-moderate elevations and the semi-arid to arid highlands of the southern Central Andes. Even subtle changes in the hydro-meteorological regime of these areas of the mountain belt react with major impacts on erosional hillslope processes and generate mass movements that fundamentally impact the transport capacity of mountain streams. Despite more severe storms in these areas, the fluvial system is characterized by pronounced variability of the stream power on different timescales, leading to cycles of sediment aggradation, the loss of agriculturally used land and severe impacts on infrastructure.
Recent studies have shown that the 1976-77 global climate shift strongly affected the South American climate. In our study, we observed a link between this climate shift and river-discharge variability in the subtropical Southern Central Andes. We analyzed the daily river-discharge time series between 1940 and 1999 from small to medium mountain drainage basins (10(2)-10(4) km(2) ) across a steep climatic and topographic gradient. We document that the discharge frequency distribution changed significantly, with higher percentiles exhibiting more pronounced trends. A change point between 1971 and 1977 marked an intensification of the hydrological cycle, which resulted in increased river discharge. In the upper Rio Bermejo basin of the northernmost Argentine Andes, the mean annual discharge increased by 40% over 7 years. Our findings are important for flood risk management in areas impacted by the 1976-77 climate shift; discharge frequency distribution analysis provides important insights into the variability of the hydrological cycle in the Andean realm.
The uppermost Cretaceous (upper Campanian-Maastrichtian) marine deposits of the central south Pyrenees host a rich larger benthic foraminiferal fauna and several rudist-rich levels. These marine deposits are directly overlain by the continental facies of the Aren and Tremp Formations, which are famous for their fossil dinosaur remains. Larger benthic foraminiferal distribution documents an important faunal turnover in all the carbonate platform environments within the photic zone, from open marine to littoral areas. Biostratigraphy indicates that this turnover occurred close to the Campanian-Maastrichtian boundary. This is also confirmed by strontium isotope stratigraphy which indicates an earliest Maastrichtian age for the appearance of the larger benthic foraminiferal assemblage constituted by Lepidorbitoides socialis, Clypeorbis mammillata, Wannierina cataluniensis, Orbitoides gruenbachensis, Siderolites aff. calcitrapoides, Fascispira colomi, Omphalocyclus macroporus and Laffiteina mengaudi. In particular, a numerical age of 71 Ma is obtained for the Hippurites radiosus level, just a few meters below the first continental deposits of the Aren sensu stricto Formation. The youngest marine sediments of the central south Pyrenees are early Maastrichtian in age. This is also an important constraint for the age of the end-Cretaceous dinosaur fossil localities of the Tremp basin. (C) 2015 Elsevier Ltd. All rights reserved.
The 2014 April 1, M-w 8.1 Iquique (Chile) earthquake struck in the Northern Chile seismic gap. With a rupture length of less than 200 km, it left unbroken large segments of the former gap. Early studies were able to model the main rupture features but results are ambiguous with respect to the role of aseismic slip and left open questions on the remaining hazard at the Northern Chile gap. A striking observation of the 2014 earthquake has been its extensive preparation phase, with more than 1300 events with magnitude above M-L 3, occurring during the 15 months preceding the main shock. Increasing seismicity rates and observed peak magnitudes accompanied the last three weeks before the main shock. Thanks to the large data sets of regional recordings, we assess the precursor activity, compare foreshocks and aftershocks and model rupture preparation and rupture effects. To tackle inversion challenges for moderate events with an asymmetric network geometry, we use full waveforms techniques to locate events, map the seismicity rate and derive source parameters, obtaining moment tensors for more than 300 events (magnitudes M-w 4.0-8.1) in the period 2013 January 1-2014 April 30. This unique data set of fore- and aftershocks is investigated to distinguish rupture process models and models of strain and stress rotation during an earthquake. Results indicate that the spatial distributions of foreshocks delineated the shallower part of the rupture areas of the main shock and its largest aftershock, well matching the spatial extension of the aftershocks cloud. Most moment tensors correspond to almost pure double couple thrust mechanisms, consistent with the slab orientation. Whereas no significant differences are observed among thrust mechanisms in different areas, nor among thrust foreshocks and aftershocks, the early aftershock sequence is characterized by the presence of normal fault mechanisms, striking parallel to the trench but dipping westward. These events likely occurred in the shallow wedge structure close to the slab interface and are consequence of the increased extensional stress in this region after the largest events. The overall stress inversion result suggests a minor stress rotation after the main shock, but a significant release of the deviatoric stress. The temporal change in the distribution of focal mechanisms can also be explained in terms of the spatial heterogeneity of the stress field: under such interpretation, the potential of a large megathrust earthquake breaking a larger segment offshore Northern Chile remains high.
The Alanya Massif, which is located to the south of central Taurides in Turkey, presents a typical nappe pile consisting of thrust sheets with contrasting metamorphic histories. In two thrust sheets, Sugozu and GundogmuAY nappes, HP metamorphism under eclogite (550-567 A degrees C/14-18 kbar) and blueschist facies (435-480 A degrees C/11-13 kbar) conditions have been recognized, respectively. Whereas the rest of the Massif underwent MP metamorphism under greenschist to amphibolite facies (525-555 A degrees C/6.5-7.5 kbar) conditions. Eclogite facies metamorphism in Sugozu nappe, which consists of homogeneous garnet-glaucophane-phengite schists with eclogite lenses is dated at 84.8 +/- A 0.8, 84.7 +/- A 1.5 and 82 +/- A 3 Ma (Santonian-Campanian) by Ar-40/Ar-39 phengite, U/Pb zircon and rutile dating methods, respectively. Similarly, phengites in GundogmuAY nappe representing an accretionary complex yield 82-80 Ma (Campanian) ages for blueschist facies metamorphism. During the exhumation, the retrograde overprint of the HP units under greenschist-amphibolite facies conditions and tectonic juxtaposition with the Barrovian units occurred during Campanian (75-78 Ma). Petrological and geochronological data clearly indicate a similar Late Cretaceous tectonometamorphic evolution for both Alanya (84-75 Ma) and Bitlis (84-72 Ma) Massifs. They form part of a single continental sliver (Alanya-Bitlis microcontinent), which was rifted from the southern part of the Anatolide-Tauride platform. The P-T-t coherence between two Massifs suggests that both Massifs have been derived from the closure of the same ocean (Alanya-Bitlis Ocean) located to the south of the Anatolide-Tauride block by a northward subduction. The boundary separating the autochthonous Tauride platform to the north from both the Alanya and Bitlis Massifs to the south represents a suture zone, the Pamphylian-Alanya-Bitlis suture.
Over the past decades, rapid and constant advances have motivated GNSS technology to approach the ability to monitor transient ground motions with mm to cm accuracy in real-time. As a result, the potential of using real-time GNSS for natural hazards prediction and early warning has been exploited intensively in recent years, e.g., landslides and volcanic eruptions monitoring. Of particular note, compared with traditional seismic instruments, GNSS does not saturate or tilt in terms of co-seismic displacement retrieving, which makes it especially valuable for earthquake and earthquake induced tsunami early warning. In this thesis, we focus on the application of real-time GNSS to fast seismic source inversion and tsunami early warning.
Firstly, we present a new approach to get precise co-seismic displacements using cost effective single-frequency receivers. As is well known, with regard to high precision positioning, the main obstacle for single-frequency GPS receiver is ionospheric delay. Considering that over a few minutes, the change of ionospheric delay is almost linear, we constructed a linear model for each satellite to predict ionospheric delay. The effectiveness of this method has been validated by an out-door experiment and 2011 Tohoku event, which confirms feasibility of using dense GPS networks for geo-hazard early warning at an affordable cost.
Secondly, we extended temporal point positioning from GPS-only to GPS/GLONASS and assessed the potential benefits of multi-GNSS for co-seismic displacement determination. Out-door experiments reveal that when observations are conducted in an adversary environment, adding a couple of GLONASS satellites could provide more reliable results. The case study of 2015 Illapel Mw 8.3 earthquake shows that the biases between co-seismic displacements derived from GPS-only and GPS/GLONASS vary from station to station, and could be up to 2 cm in horizontal direction and almost 3 cm in vertical direction. Furthermore, slips inverted from GPS/GLONASS co-seismic displacements using a layered crust structure on a curved plane are shallower and larger for the Illapel event.
Thirdly, we tested different inversion tools and discussed the uncertainties of using real-time GNSS for tsunami early warning. To be exact, centroid moment tensor inversion, uniform slip inversion using a single Okada fault and distributed slip inversion in layered crust on a curved plane were conducted using co-seismic displacements recorded during 2014 Pisagua earthquake. While the inversion results give similar magnitude and the rupture center, there are significant differences in depth, strike, dip and rake angles, which lead to different tsunami propagation scenarios. Even though, resulting tsunami forecasting along the Chilean coast is close to each other for all three models.
Finally, based on the fact that the positioning performance of BDS is now equivalent to GPS in Asia-Pacific area and Manila subduction zone has been identified as a zone of potential tsunami hazard, we suggested a conceptual BDS/GPS network for tsunami early warning in South China Sea. Numerical simulations with two earthquakes (Mw 8.0 and Mw 7.5) and induced tsunamis demonstrate the viability of this network. In addition, the advantage of BDS/GPS over a single GNSS system by source inversion grows with decreasing earthquake magnitudes.
The role that climate and environmental history may have played in influencing human evolution has been the focus of considerable interest and controversy among paleoanthropologists for decades. Prior attempts to understand the environmental history side of this equation have centered around the study of outcrop sediments and fossils adjacent to where fossil hominins (ancestors or close relatives of modern humans) are found, or from the study of deep sea drill cores. However, outcrop sediments are often highly weathered and thus are unsuitable for some types of paleoclimatic records, and deep sea core records come from long distances away from the actual fossil and stone tool remains. The Hominin Sites and Paleolakes Drilling Project (HSPDP) was developed to address these issues. The project has focused its efforts on the eastern African Rift Valley, where much of the evidence for early hominins has been recovered. We have collected about 2 km of sediment drill core from six basins in Kenya and Ethiopia, in lake deposits immediately adjacent to important fossil hominin and archaeological sites. Collectively these cores cover in time many of the key transitions and critical intervals in human evolutionary history over the last 4 Ma, such as the earliest stone tools, the origin of our own genus Homo, and the earliest anatomically modern Homo sapiens. Here we document the initial field, physical property, and core description results of the 2012-2014 HSPDP coring campaign.
Volunteered geographical information (VGI) and citizen science have become important sources data for much scientific research. In the domain of land cover, crowdsourcing can provide a high temporal resolution data to support different analyses of landscape processes. However, the scientists may have little control over what gets recorded by the crowd, providing a potential source of error and uncertainty. This study compared analyses of crowdsourced land cover data that were contributed by different groups, based on nationality (labelled Gondor and Non-Gondor) and on domain experience (labelled Expert and Non-Expert). The analyses used a geographically weighted model to generate maps of land cover and compared the maps generated by the different groups. The results highlight the differences between the maps how specific land cover classes were under-and over-estimated. As crowdsourced data and citizen science are increasingly used to replace data collected under the designed experiment, this paper highlights the importance of considering between group variations and their impacts on the results of analyses. Critically, differences in the way that landscape features are conceptualised by different groups of contributors need to be considered when using crowdsourced data in formal scientific analyses. The discussion considers the potential for variation in crowdsourced data, the relativist nature of land cover and suggests a number of areas for future research. The key finding is that the veracity of citizen science data is not the critical issue per se. Rather, it is important to consider the impacts of differences in the semantics, affordances and functions associated with landscape features held by different groups of crowdsourced data contributors.
P-wave polarization at the Grafenberg array (GRF) in southern Germany is analysed in terms of azimuthal deviations and deviations in the vertical polarization using 20 yr of broad-band recordings. An automated procedure for estimating P-wave polarization parameters is suggested, based on the definition of a characteristic function, which evaluates the polarization angles and their time variability as well as the amplitude, linearity and the signal-to-noise ratio of the P wave. P-wave polarization at the GRF array is shown to depend mainly on frequency and backazimuth and only slightly on epicentral distance indicating depth-dependent local anisotropy and lateral heterogeneity. A harmonic analysis is applied to the azimuthal anomalies to analyse their periodicity as a function of backazimuth. The dominant periods are 180A degrees A and 360A degrees. At low frequencies, between 0.03 and 0.1 Hz, the observed fast directions of azimuthal anisotropy inferred from the 180A degrees A periodicity are similar across the array. The average fast direction of azimuthal anisotropy at these frequencies is N20A degrees E with an uncertainty of about 8A degrees A and is consistent with fast directions of Pn-wave propagation. Lateral velocity gradients determined for the low-frequency band are compatible with the Moho topography of the area. A more complex pattern in the horizontal fast axis orientation beneath the GRF array is observed in the high-frequency band between 0.1 and 0.5 Hz, and is attributed to anisotropy in the upper crust. A remarkable rotation of the horizontal fast axis orientation across the suture between the geological units Moldanubicum and Saxothuringicum is observed. In contrast, the 360A degrees A periodicity at high frequencies is rather consistent across the array and may either point to lower velocities in the upper crust towards the Bohemian Massif and/or to anisotropy dipping predominantly in the NE-SW direction. Altogether, P-wave polarization analysis indicates the presence of layered lithospheric anisotropy in the area of the GRF array. Seismic anisotropy is more variable in the brittle upper crust compared to lower crustal and subcrustal depths.
Dynamics of mantle plumes
(2016)
Mantle plumes are a link between different scales in the Earth’s mantle: They are an important part of large-scale mantle convection, transporting material and heat from the core-mantle boundary to the surface, but also affect processes on a smaller scale, such as melt generation and transport and surface magmatism. When they reach the base of the lithosphere, they cause massive magmatism associated with the generation of large igneous provinces, and they can be related to mass extinction events (Wignall, 2001) and continental breakup (White and McKenzie, 1989).
Thus, mantle plumes have been the subject of many previous numerical modelling studies (e.g. Farnetani and Richards, 1995; d’Acremont et al., 2003; Lin and van Keken, 2005; Sobolev et al., 2011; Ballmer et al., 2013). However, complex mechanisms, such as the development and implications of chemical heterogeneities in plumes, their interaction with mid-ocean ridges and global mantle flow, and melt ascent from the source region to the surface are still not very well understood; and disagreements between observations and the predictions of classical plume models have led to a challenge of the plume concept in general (Czamanske et al., 1998; Anderson, 2000; Foulger, 2011). Hence, there is a need for more sophisticated models that can explain the underlying physics, assess which properties and processes are important, explain how they cause the observations visible at the Earth’s surface and provide a link between the different scales.
In this work, integrated plume models are developed that investigate the effect of dense recycled oceanic crust on the development of mantle plumes, plume–ridge interaction under the influence of global mantle flow and melting and melt migration in form of two-phase flow.
The presented analysis of these models leads to a new, updated picture of mantle plumes: Models considering a realistic depth-dependent density of recycled oceanic crust and peridotitic mantle material show that plumes with excess temperatures of up to 300 K can transport up to 15% of recycled oceanic crust through the whole mantle. However, due to the high density of recycled crust, plumes can only advance to the base of the lithosphere directly if they have high excess temperatures, high plume volumes and the lowermost mantle is subadiabatic, or plumes rise from the top or edges of thermo-chemical piles. They might only cause minor surface uplift, and instead of the classical head–tail structure, these low-buoyancy plumes are predicted to be broad features in the lower mantle with much less pronounced plume heads. They can form a variety of shapes and regimes, including primary plumes directly advancing to the base of the lithosphere, stagnating plumes, secondary plumes rising from the core–mantle boundary or a pool of eclogitic material in the upper mantle and failing plumes. In the upper mantle, plumes are tilted and deflected by global mantle flow, and the shape, size and stability of the melting region is influenced by the distance from nearby plate boundaries, the speed of the overlying plate and the movement of the plume tail arriving from the lower mantle. Furthermore, the structure of the lithosphere controls where hot material is accumulated and melt is generated. In addition to melting in the plume tail at the plume arrival position, hot plume material flows upwards towards opening rifts, towards mid-ocean ridges and towards other regions of thinner lithosphere, where it produces additional melt due to decompression. This leads to the generation of either broad ridges of thickened magmatic crust or the separation into multiple thinner lines of sea mount chains at the surface. Once melt is generated within the plume, it influences its dynamics, lowering the viscosity and density, and while it rises the melt volume is increased up to 20% due to decompression. Melt has the tendency to accumulate at the top of the plume head, forming diapirs and initiating small-scale convection when the plume reaches the base of the lithosphere. Together with the introduced unstable, high-density material produced by freezing of melt, this provides an efficient mechanism to thin the lithosphere above plume heads.
In summary, this thesis shows that mantle plumes are more complex than previously considered, and linking the scales and coupling the physics of different processes occurring in mantle plumes can provide insights into how mantle plumes are influenced by chemical heterogeneities, interact with the lithosphere and global mantle flow, and are affected by melting and melt migration. Including these complexities in geodynamic models shows that plumes can also have broad plume tails, might produce only negligible surface uplift, can generate one or several volcanic island chains in interaction with a mid–ocean ridge, and can magmatically thin the lithosphere.
This manuscript proposes a method to assess hydrological drought in semi-arid environments under high impoundment rate and applies it to the semi-arid Jaguaribe River basin in Brazil. It analyzes droughts (1) in the largest reservoir systems; (2) in the Upper Basin, considering 4744 reservoirs, 800 wells and almost 18,000 cisterns; and (3) in reservoirs of different sizes during multiyear droughts. Results show that the water demand is constrained in the basin; hydrological and meteorological droughts are often out of phase; there is a negative correlation between storage level and drought severity; and the small systems cannot cope with long-term droughts.
This study analyses some hydrological driving forces and their interrelation with surface-flow initiation in a semiarid Caatinga basin (12km(2)), Northeastern Brazil. During the analysis period (2005 - 2014), 118 events with precipitation higher than 10mm were monitored, providing 45 events with runoff, 25 with negligible runoff and 49 without runoff. To verify the dominant processes, 179 on-site measurements of saturated hydraulic conductivity (Ksat) were conducted. The results showed that annual runoff coefficient lay below 0.5% and discharge at the outlet has only occurred four days per annum on average, providing an insight to the surface-water scarcity of the Caatinga biome. The most relevant variables to explain runoff initiation were total precipitation and maximum 60-min rainfall intensity (I-60). Runoff always occurred when rainfall surpassed 31mm, but it never occurred for rainfall below 14mm or for I-60 below 12mmh(-1). The fact that the duration of the critical intensity is similar to the basin concentration time (65min) and that the infiltration threshold value approaches the river-bank saturated hydraulic conductivity support the assumption that Hortonian runoff prevails. However, none of the analysed variables (total or precedent precipitation, soil moisture content, rainfall intensities or rainfall duration) has been able to explain the runoff initiation in all monitored events: the best criteria, e.g. failed to explain 27% of the events. It is possible that surface-flow initiation in the Caatinga biome is strongly influenced by the root-system dynamics, which changes macro-porosity status and, therefore, initial abstraction. Copyright (c) 2016 John Wiley & Sons, Ltd.
We compare the ability of various site-condition proxies (SCPs) to reduce the aleatory variability of ground motion prediction equations (GMPEs). Three SCPs (measured V-S30, inferred V-S30, local topographic slope) and two accelerometric databases (RESORCE and NGA-West2) are considered. An artificial neural network (ANN) approach including a random-effect procedure is used to derive GMPEs setting the relationship between peak ground acceleration (PGA), peak ground velocity (PGV), pseudo-spectral acceleration [PSA(T)], and explanatory variables (M-w, R-JB, and V-S30 or Slope). The analysis is performed using both discrete site classes and continuous proxy values. All "non-measured" SCPs exhibit a rather poor performance in reducing aleatory variability, compared to the better performance of measured V-S30. A new, fully data-driven GMPE based on the NGA-West2 is then derived, with an aleatory variability value depending on the quality of the SCP. It proves very consistent with previous GMPEs built on the same data set. Measuring V-S30 allows for benefit from an aleatory variability reduction up to 15%.
Variations in the distribution of mass within an orogen may lead to transient sediment storage, which in turn might affect the state of stress and the level of fault activity. Distinguishing between different forcing mechanisms causing variations of sediment flux and tectonic activity, is therefore one of the most challenging tasks in understanding the spatiotemporal evolution of active mountain belts.
The Himalayan mountain belt is one of the most significant Cenozoic collisional mountain belt, formed due to collision between northward-bound Indian Plate and the Eurasian Plate during the last 55-50 Ma. Ongoing convergence of these two tectonic plates is accommodated by faulting and folding within the Himalayan arc-shaped orogen and the continued lateral and vertical growth of the Tibetan Plateau and mountain belts adjacent to the plateau as well as regions farther north. Growth of the Himalayan orogen is manifested by the development of successive south-vergent thrust systems. These thrust systems divide the orogen into different morphotectonic domains. From north to south these thrusts are the Main Central Thrust (MCT), the Main Boundary Thrust (MBT) and the Main Frontal Thrust (MFT). The growing topography interacts with moisture-bearing monsoonal winds, which results in pronounced gradients in rainfall, weathering, erosion and sediment transport toward the foreland and beyond. However, a fraction of this sediment is trapped and transiently stored within the intermontane valleys or ‘dun’s within the lower-elevation foothills of the range. Improved understanding of the spatiotemporal evolution of these sediment archives could provide a unique opportunity to decipher the triggers of variations in sediment production, delivery and storage in an actively deforming mountain belt and support efforts to test linkages between sediment volumes in intermontane basins and changes in the shallow crustal stress field. As sediment redistribution in mountain belts on timescales of 102-104 years can effect cultural characteristics and infrastructure in the intermontane valleys and may even impact the seismotectonics of a mountain belt, there is a heightened interest in understanding sediment-routing processes and causal relationships between tectonism, climate and topography. It is here at the intersection between tectonic processes and superposed climatic and sedimentary processes in the Himalayan orogenic wedge, where my investigation is focused on. The study area is the intermontane Kangra Basin in the northwestern Sub-Himalaya, because the characteristics of the different Himalayan morphotectonic provinces are well developed, the area is part of a region strongly influenced by monsoonal forcing, and the existence of numerous fluvial terraces provides excellent strain markers to assess deformation processes within the Himalayan orogenic wedge. In addition, being located in front of the Dhauladhar Range the region is characterized by pronounced gradients in past and present-day erosion and sediment processes associated with repeatedly changing climatic conditions. In light of these conditions I analysed climate-driven late Pleistocene-Holocene sediment cycles in this tectonically active region, which may be responsible for triggering the tectonic re-organization within the Himalayan orogenic wedge, leading to out-of-sequence thrusting, at least since early Holocene.
The Kangra Basin is bounded by the MBT and the Sub-Himalayan Jwalamukhi Thrust (JMT) in the north and south, respectively and transiently stores sediments derived from the Dhauladhar Range. The Basin contains ~200-m-thick conglomerates reflecting two distinct aggradation phases; following aggradation, several fluvial terraces were sculpted into these fan deposits. 10Be CRN surface exposure dating of these terrace levels provides an age of 53.4±3.2 ka for the highest-preserved terrace (AF1); subsequently, this surface was incised until ~15 ka, when the second fan (AF2) began to form. AF2 fan aggradation was superseded by episodic Holocene incision, creating at least four terrace levels. We find a correlation between variations in sediment transport and ∂18O records from regions affected by the Indian Summer Monsoon (ISM). During strengthened ISMs sand post-LGM glacial retreat, aggradation occurred in the Kangra Basin, likely due to high sediment flux, whereas periods of a weakened ISM coupled with lower sediment supply coincided with renewed re-incision.
However, the evolution of fluvial terraces along Sub-Himalayan streams in the Kangra sector is also forced by tectonic processes. Back-tilted, folded terraces clearly document tectonic activity of the JMT. Offset of one of the terrace levels indicates a shortening rate of 5.6±0.8 to 7.5±1.0 mm.a-1 over the last ~10 ka. Importantly, my study reveals that late Pleistocene/Holocene out-of-sequence thrusting accommodates 40-60% of the total 14±2 mm.a-1 shortening partitioned throughout the Sub-Himalaya. Importantly, the JMT records shortening at a lower rate over longer timescales hints towards out-of-sequence activity within the Sub-Himalaya. Re-activation of the JMT could be related to changes in the tectonic stress field caused by large-scale sediment removal from the basin. I speculate that the deformation processes of the Sub-Himalaya behave according to the predictions of critical wedge model and assume the following: While >200m of sediment aggradation would trigger foreland-ward propagation of the deformation front, re-incision and removal of most of the stored sediments (nearly 80-85% of the optimum basin-fill) would again create a sub-critical condition of the wedge taper and trigger the retreat of the deformation front.
While tectonism is responsible for the longer-term processes of erosion associated with steepening hillslopes, sediment cycles in this environment are mainly the result of climatic forcing. My new 10Be cosmogenic nuclide exposure dates and a synopsis of previous studies show the late Pleistocene to Holocene alluvial fills and fluvial terraces studied here record periodic fluctuations of sediment supply and transport capacity on timescales of 1000-100000 years. To further evaluate the potential influence of climate change on these fluctuations, I compared the timing of aggradation and incision phases recorded within remnant alluvial fans and terraces with continental climate archives such as speleothems in neighboring regions affected by monsoonal precipitation. Together with previously published OSL ages yielding the timing of aggradation, I find a correlation between variations in sediment transport with oxygen-isotope records from regions affected by the Indian Summer Monsoon (ISM). Accordingly, during periods of increased monsoon intensity (transitions from dry and cold to wet and warm periods – MIS4 to MIS3 and MIS2 to MIS1) (MIS=marine isotope stage) and post-Last Glacial Maximum glacial retreat, aggradation occurred in the Kangra Basin, likely due to high sediment flux. Conversely, periods of weakened monsoon intensity or lower sediment supply coincide with re-incision of the existing basin-fill.
Finally, my study entails part of a low-temperature thermochronology study to assess the youngest exhumation history of the Dhauladhar Range. Zircon helium (ZHe) ages and existing low-temperature data sets (ZHe, apatite fission track (AFT)) across this range, together with 3D thermokinematic modeling (PECUBE) reveals constraints on exhumation and activity of the range-bounding Main Boundary Thrust (MBT) since at least mid-Miocene time. The modeling results indicate mean slip rates on the MBT-fault ramp of ~2 – 3 mm.a-1 since its activation. This has lead to the growth of the >5-km-high frontal Dhauladhar Range and continuous deep-seated exhumation and erosion. The obtained results also provide interesting constraints of deformation patterns and their variation along strike. The results point towards the absence of the time-transient ‘mid-crustal ramp’ in the basal decollement and
duplexing of the Lesser Himalayan sequence, unlike the nearby regions or even the central Nepal domain. A fraction of convergence (~10-15%) is accommodated along the deep-seated MBT-ramp, most likely merging into the MHT. This finding is crucial for a rigorous assessment of the overall level of tectonic activity in the Himalayan morphotectonic provinces as it contradicts recently-published geodetic shortening estimates. In these studies, it has been proposed that the total Himalayan shortening in the NW Himalaya is accommodated within the Sub-Himalaya whereas no tectonic activity is assigned to the MBT.
Deciphering the response of sediment routing systems to climatic forcing is fundamental for understanding the impacts of climate change on landscape evolution. In the Kangra Basin (northwest Sub-Himalaya, India), upper Pleistocene to Holocene alluvial fills and fluvial terraces record periodic fluctuations of sediment supply and transport capacity on timescales of 10(3) to 10(5) yr. To evaluate the potential influence of climate change on these fluctuations, we compare the timing of aggradation and incision phases recorded within remnant alluvial fans and terraces with climate archives. New surface-exposure dating of six terrace levels with in-situ cosmogenic Be-10 indicates the onset of incision phases. Two terrace surfaces from the highest level (T1) sculpted into the oldest preserved alluvial fan (AF1) date back to 53.4 +/- 3.2 ka and 43.0 +/- 2.7 ka (1 sigma). T2 surfaces sculpted into the remnants of AF1 have exposure ages of 18.6 +/- 1.2 ka and 15.3 +/- 0.9 ka, while terraces sculpted into the upper Pleistocene-Holocene fan (AF2) provide ages of 9.3 +/- 0.4 ka (T3), 7.1 +/- 0.4 ka (T4), 5.2 +/- 0.4 ka (T5) and 3.6 +/- 0.2 ka (T6). Together with previously published OSL ages yielding the timing of aggradation, we find a correlation between variations in sediment transport with oxygen-isotope records from regions affected by the Indian Summer Monsoon. During periods of increased monsoon intensity and post-Last Glacial Maximum glacial retreat, aggradation occurred in the Kangra Basin, likely due to high sediment flux, whereas periods of weakened monsoon intensity or lower sediment supply coincide with incision. (C) 2016 Elsevier B.V. All rights reserved.
The southernmost thrust of the Himalayan orogenic wedge that separates the foreland from the orogen, the Main Frontal Thrust, is thought to accommodate most of the ongoing crustal shortening in the Sub-Himalaya. Steepened longitudinal river profile segments, terrace offsets, and back-tilted fluvial terraces within the Kangra reentrant of the NW Sub-Himalaya suggest Holocene activity of the Jwalamukhi Thrust (JMT) and other thrust faults that may be associated with strain partitioning along the toe of the Himalayan wedge. To assess the shortening accommodated by the JMT, we combine morphometric terrain analyses with in situ Be-10-based surface-exposure dating of the deformed terraces. Incision into upper Pleistocene sediments within the Kangra Basin created two late Pleistocene terrace levels (T1 and T2). Subsequent early Holocene aggradation shortly before similar to 10ka was followed by episodic reincision, which created four cut-and-fill terrace levels, the oldest of which (T3) was formed at 10.10.9ka. A vertical offset of 445m of terrace T3 across the JMT indicates a shortening rate of 5.60.8 to 7.51.1mma(-1) over the last similar to 10ka. This result suggests that thrusting along the JMT accommodates 40-60% of the total Sub-Himalayan shortening in the Kangra reentrant over the Holocene. We speculate that this out-of-sequence shortening may have been triggered or at least enhanced by late Pleistocene and Holocene erosion of sediments from the Kangra Basin.
Strong waves in the mid-latitude circulation have been linked to extreme surface weather and thus changes in waviness could have serious consequences for society. Several theories have been proposed which could alter waviness, including tropical sea surface temperature anomalies or rapid climate change in the Arctic. However, so far it remains unclear whether any changes in waviness have actually occurred. Here we propose a novel meandering index which captures the maximum waviness in geopotential height contours at any given day, using all information of the full spatial position of each contour. Data are analysed on different time scale (from daily to 11 day running means) and both on hemispheric and regional scales. Using quantile regressions, we analyse how seasonal distributions of this index have changed over 1979-2015. The most robust changes are detected for autumn which has seen a pronounced increase in strongly meandering patterns at the hemispheric level as well as over the Eurasian sector. In summer for both the hemisphere and the Eurasian sector, significant downward trends in meandering are detected on daily timescales which is consistent with the recently reported decrease in summer storm track activity. The American sector shows the strongest increase in meandering in the warm season: in particular for 11 day running mean data, indicating enhanced amplitudes of quasi-stationary waves. Our findings have implications for both the occurrence of recent cold spells and persistent heat waves in the mid-latitudes.
The water cycle of sites with shallow groundwater tables is characterized by complex interactions of hydrological and ecological processes. The water balance components, which are subject to diurnal fluctuations, are best measured with groundwater lysimeters. However, the lower boundary condition of such lysimeters affects most of the hydrological variables, particularly when considering short time scales, and has to be defined in such a way as to facilitate realistic simulations. In this paper, different means of controlling the lower boundary condition of groundwater lysimeters were compared with respect to their ability to simulate the behavior of the water balance components properly. Measurements of rain-free periods from a lysimeter station installed in the Spreewald wetland in north-east Germany were evaluated. The most common groundwater lysimeter type is controlled using a Mariotte bottle and sets the groundwater level in the soil monolith to a constant level, which here caused an alteration of the inflow to the lysimeter, with respect to both its value and diurnal behavior. Still, daily evapotranspiration values were realistic and this simple and robust approach may be used for time intervals not shorter than one day. High-resolution measurements can be gained from lysimeters that automatically adjust the groundwater level by a system of pumps and valves on an hourly basis. Still, reliable results were only obtained when the conditions in the lysimeter and the surrounding field, where the target groundwater level was measured, were in accordance. Otherwise (e.g., when the groundwater level differed) an unrealistic inflow behavior evolved. Reasonable results, even for slightly diverging conditions, were gained with a new approach that defined the lower boundary conditions by controlling the inflows and outflows of the lysimeter. This approach further enabled the groundwater level itself to be the study subject, thereby enlarging the field of possible applications of groundwater lysimeters. (C) 2015 Elsevier B.V. All rights reserved.
Effects of data and model simplification on the results of a wetland water resource management model
(2016)
This paper presents the development of a wetland water balance model for use in a large river basin with many different wetlands. The basic model was primarily developed for a single wetland with a complex water management system involving large amounts of specialized input data and water management details. The aim was to simplify the model structure and to use only commonly available data as input for the model, with the least possible loss of accuracy. Results from different variants of the model and data adaptation were tested against results from a detailed model. This shows that using commonly available data and unifying and simplifying the input data is tolerable up to a certain level. The simplification of the model has greater effects on the evaluated water balance components than the data adaptation. Because this simplification was necessary for large-scale use, we suggest that, for reasons of comparability, simpler models should always be applied with uniform data bases for large regions, though these should only be moderately simplified. Further, we recommend using these simplified models only for large-scale comparisons and using more specific, detailed models for investigations on smaller scales.
Effects of Data and Model Simplification on the Results of a Wetland Water Resource Management Model
(2016)
This paper presents the development of a wetland water balance model for use in a large river basin with many different wetlands. The basic model was primarily developed for a single wetland with a complex water management system involving large amounts of specialized input data and water management details. The aim was to simplify the model structure and to use only commonly available data as input for the model, with the least possible loss of accuracy. Results from different variants of the model and data adaptation were tested against results from a detailed model. This shows that using commonly available data and unifying and simplifying the input data is tolerable up to a certain level. The simplification of the model has greater effects on the evaluated water balance components than the data adaptation. Because this simplification was necessary for large-scale use, we suggest that, for reasons of comparability, simpler models should always be applied with uniform data bases for large regions, though these should only be moderately simplified. Further, we recommend using these simplified models only for large-scale comparisons and using more specific, detailed models for investigations on smaller scales.
Many German lakes experienced significant water level declines in recent decades that are not fully understood due to the short observation period. At a typical northeastern German groundwater-fed lake with a complex basin morphology, an acoustic sub-bottom profile was analysed together with a transect of five sediment cores, which were correlated using multiple proxies (sediment facies, -XRF, macrofossils, subfossil Cladocera). Shifts in the boundary between sand and mud deposition were controlled by lake level changes, and hence, allowed the quantification of an absolute lake level amplitude of similar to 8m for the Holocene. This clearly exceeded observed modern fluctuations of 1.3m (AD 1973-2010). Past lake level changes were traced continuously using the calcium-record. During high lake levels, massive organic muds were deposited in the deepest lake basin, whereas lower lake levels isolated the sub-basins and allowed carbonate deposition. During the beginning of the Holocene (>9700cal. a BP), lake levels were high, probably due to final melting of permafrost and dead-ice remains. The establishment of water-use intensive Pinus forests caused generally low (3-4m below modern) but fluctuating lake levels (9700-6400cal. a BP). Afterwards, the lake showed an increasing trend and reached a short-term highstand at c.5000cal. a BP (4m above modern). At the transition towards a cooler and wetter late Holocene, forests dominated by Quercus and Fagus and initial human impact probably contributed more positively to groundwater recharge. Lake levels remained high between 3800 and 800cal. a BP, but the lake system was not sensitive enough to record short-term fluctuations during this period. Lake level changes were recorded again when humans profoundly affected the drainage system, land cover and lake trophy. Hence, local Holocene water level changes reflect feedbacks between catchment and vegetation characteristics and human impact superimposed by climate change at multiple temporal scales.
Floods frequently cause substantial economic and human losses, particularly in developing countries. For the development of sound flood risk management schemes that reduce flood consequences, detailed insights into the different components of the flood risk management cycle, such as preparedness, response, flood impact analyses and recovery, are needed. However, such detailed insights are often lacking: commonly, only (aggregated) data on direct flood damage are available. Other damage categories such as losses owing to the disruption of production processes are usually not considered, resulting in incomplete risk assessments and possibly inappropriate recommendations for risk management. In this paper, data from 858 face-to-face interviews among flood-prone households and small businesses in Can Tho city in the Vietnamese Mekong Delta are presented to gain better insights into the damage caused by the 2011 flood event and its management by households and businesses.
The role of soil erosion in terrestrial carbon (C) sequestration and release remains one of the most important uncertainties in our attempts to determine the potential of soils to mediate climate change. Despite its widely recognized importance for terrestrial C sequestration, to date, no Earth System Model (ESM) implements soil erosion effects on carbon cycling in sufficient detail. So far, available studies have mostly investigated the magnitude of erosional C transport and in-situ measurements of vertical C fluxes on the catchment or regional scale. Recognizing the need to adequately represent C erosion processes and controls in ESMs, we provide a comprehensive cross-disciplinary review on lateral C redistribution in the landscape and discuss the implications for bio-geochemical cycling of carbon. We present current knowledge on the role of erosional C distribution in controlling the stabilization and release of C in soils, taking into consideration the important geomorphic, ecological, hydrologic, pedologic and micro-climatic processes and controls that affect soil organic carbon (SOC) stock, fluxes, and persistence in dynamic landscapes. Further, we provide an overview on latest experimental and modelling approaches that are being used to investigate the role of erosion in the carbon cycle. Finally, to advance our understanding of the role of soil redistribution in biogeochemical cycles of essential elements, we discuss the most promising topics for future research in this field. (C) 2015 Elsevier B.V. All rights reserved.
Event coincidence analysis for quantifying statistical interrelationships between event time series
(2016)
Studying event time series is a powerful approach for analyzing the dynamics of complex dynamical systems in many fields of science. In this paper, we describe the method of event coincidence analysis to provide a framework for quantifying the strength, directionality and time lag of statistical interrelationships between event series. Event coincidence analysis allows to formulate and test null hypotheses on the origin of the observed interrelationships including tests based on Poisson processes or, more generally, stochastic point processes with a prescribed inter-event time distribution and other higher-order properties. Applying the framework to country-level observational data yields evidence that flood events have acted as triggers of epidemic outbreaks globally since the 1950s. Facing projected future changes in the statistics of climatic extreme events, statistical techniques such as event coincidence analysis will be relevant for investigating the impacts of anthropogenic climate change on human societies and ecosystems worldwide.
Understanding the role of natural climate variability under the pressure of human induced changes of climate and landscapes, is crucial to improve future projections and adaption strategies. This doctoral thesis aims to reconstruct Holocene climate and environmental changes in NE Germany based on annually laminated lake sediments. The work contributes to the ICLEA project (Integrated CLimate and Landscape Evolution Analyses). ICLEA intends to compare multiple high-resolution proxy records with independent chronologies from the N central European lowlands, in order to disentangle the impact of climate change and human land use on landscape development during the Lateglacial and Holocene. In this respect, two study sites in NE Germany are investigated in this doctoral project, Lake Tiefer See and palaeolake Wukenfurche. While both sediment records are studied with a combination of high-resolution sediment microfacies and geochemical analyses (e.g. µ-XRF, carbon geochemistry and stable isotopes), detailed proxy understanding mainly focused on the continuous 7.7 m long sediment core from Lake Tiefer See covering the last ~6000 years. Three main objectives are pursued at Lake Tiefer See: (1) to perform a reliable and independent chronology, (2) to establish microfacies and geochemical proxies as indicators for climate and environmental changes, and (3) to trace the effects of climate variability and human activity on sediment deposition.
Addressing the first aim, a reliable chronology of Lake Tiefer See is compiled by using a multiple-dating concept. Varve counting and tephra findings form the chronological framework for the last ~6000 years. The good agreement with independent radiocarbon dates of terrestrial plant remains verifies the robustness of the age model. The resulting reliable and independent chronology of Lake Tiefer See and, additionally, the identification of nine tephras provide a valuable base for detailed comparison and synchronization of the Lake Tiefer See data set with other climate records. The sediment profile of Lake Tiefer See exhibits striking alternations between well-varved and non-varved sediment intervals. The combination of microfacies, geochemical and microfossil (i.e. Cladocera and diatom) analyses indicates that these changes of varve preservation are caused by variations of lake circulation in Lake Tiefer See. An exception is the well-varved sediment deposited since AD 1924, which is mainly influenced by human-induced lake eutrophication. Well-varved intervals before the 20th century are considered to reflect phases of reduced lake circulation and, consequently, stronger anoxic conditions. Instead, non-varved intervals indicate increased lake circulation in Lake Tiefer See, leading to more oxygenated conditions at the lake ground. Furthermore, lake circulation is not only influencing sediment deposition, but also geochemical processes in the lake. As, for example, the proxy meaning of δ13COM varies in time in response to changes of the oxygen regime in the lake hypolinion. During reduced lake circulation and stronger anoxic conditions δ13COM is influenced by microbial carbon cycling. In contrast, organic matter degradation controls δ13COM during phases of intensified lake circulation and more oxygenated conditions. The varve preservation indicates an increasing trend of lake circulation at Lake Tiefer See after ~4000 cal a BP. This trend is superimposed by decadal to centennial scale variability of lake circulation intensity. Comparison to other records in Central Europe suggests that the long-term trend is probably related to gradual changes in Northern Hemisphere orbital forcing, which induced colder and windier conditions in Central Europe and, therefore, reinforced lake circulation. Decadal to centennial scale periods of increased lake circulation coincide with settlement phases at Lake Tiefer See, as inferred from pollen data of the same sediment record. Deforestation reduced the wind shelter of the lake, which probably increased the sensitivity of lake circulation to wind stress. However, results of this thesis also suggest that several of these phases of increased lake circulation are additionally reinforced by climate changes. A first indication is provided by the comparison to the Baltic Sea record, which shows striking correspondence between major non-varved intervals at Lake Tiefer See and bioturbated sediments in the Baltic Sea. Furthermore, a preliminary comparison to the ICLEA study site Lake Czechowskie (N central Poland) shows a coincidence of at least three phases of increased lake circulation in both lakes, which concur with periods of known climate changes (2.8 ka event, ’Migration Period’ and ’Little Ice Age’). These results suggest an additional over-regional climate forcing also on short term increased of lake circulation in Lake Tiefer See.
In summary, the results of this thesis suggest that lake circulation at Lake Tiefer See is driven by a combination of long-term and short-term climate changes as well as of anthropogenic deforestation phases. Furthermore, the lake circulation drives geochemical cycles in the lake affecting the meaning of proxy data. Therefore, the work presented here expands the knowledge of climate and environmental variability in NE Germany. Furthermore, the integration of the Lake Tiefer See multi-proxy record in a regional comparison with another ICLEA side, Lake Czechowskie, enabled to better decipher climate changes and human impact on the lake system. These first results suggest a huge potential for further detailed regional comparisons to better understand palaeoclimate dynamics in N central Europe.
The Paleocene was a time of transition for the Arctic, with magmatic activity of the High Arctic Large Igneous Province (HALIP) giving way to magmatism of the North Atlantic Large Igneous Province in connection to plate tectonic changes in the Arctic and North Atlantic. In this study we investigate the Paleocene magmatic record and sediment pathways of the Basilika Formation exposed in the Central Tertiary Basin of Svalbard. By means of geochemistry, SmNd isotopic signatures, and zircon UPb geochronology we investigate the characteristics of several bentonite layers contained in the Basilika Formation, as well as the provenance of the intercalated clastic sediments. Our data show that the volcanic ash layers of the Basilika Formation, which were diagenetically altered to bentonites, originate from alkaline continental-rift magmatism such as the last, explosive stages of the HALIP in North Greenland and the Canadian Arctic. The volcanic ash layers were deposited on Svalbard in a flat shelf environment with dominant sediment supply from the east. Dating of detrital zircons suggests that the detritus was derived from Siberian sources, primarily from the Verkhoyansk Fold-and-Thrust Belt, which would require transport over similar to 3000 km across the Arctic.
Bedrock landslides, by excavating deep below saprolite-rock interfaces, create conditions for weathering in which all mineral phases in a lithology are initially unweathered within landslide deposits. As a result, the most labile phases dominate the weathering immediately after mobilisation and during a transient period of depletion. This mode of dissolution can strongly alter the overall output of solutes from catchments and their contribution to global chemical cycles if landslide-derived material is retained in catchments for extended periods after mass wasting.
A link between chemical weathering and physical erosion exists at the catchment scale over a wide range of erosion rates(1,2). However, in mountain environments, where erosion rates are highest, weathering may be kinetically limited(3-5) and therefore decoupled from erosion. In active mountain belts, erosion is driven by bedrock landsliding(6) at rates that depend strongly on the occurrence of extreme rainfall or seismicity(7). Although landslides affect only a small proportion of the landscape, bedrock landsliding can promote the collection and slow percolation of surface runoff in highly fragmented rock debris and create favourable conditions for weathering. Here we show from analysis of surface water chemistry in the Southern Alps of New Zealand that weathering in bedrock landslides controls the variability in solute load of these mountain rivers. We find that systematic patterns in surface water chemistry are strongly associated with landslide occurrence at scales from a single hillslope to an entire mountain belt, and that landslides boost weathering rates and river solute loads over decades. We conclude that landslides couple erosion and weathering in fast-eroding uplands and, thus, mountain weathering is a stochastic process that is sensitive to climatic and tectonic controls on mass wasting processes.
Lake Meerfelder Maar (MFM) is the northernmost Western European sediment record with annual laminations across the Younger Dryas (YD), and the onset of the YD in the record of MFM has previously been defined as an increase in non-arboreal pollen abundance at ca. 12 680 varve a BP. Here we present a palynological record at unprecedented subdecadal resolution for MFM, covering the Allerod-YD transition. Our results show a fluctuation in pollen accumulation rates (PARs) before the onset of the YD, with lower rates between ca. 12 725 and 12 685 varve a BP. The fluctuation in PARs occurs simultaneous with a previously undescribed short fluctuation in sediment composition and varve thickness, as well as with changes in biogeochemical proxies. The combined evidence indicates signs of climatic instability ca. 45 years before the onset of the YD. The PAR records of Betula and Pinus furthermore show earlier and more abrupt changes at the onset of the YD than the percentage-records do. Finally, heliophilous herbaceous taxa show a delayed increase following the onset of the YD of ca. 145 years. This paper illustrates the potential to identify previously unrecognized climate variability and vegetation change when using subdecadal-resolution analyses. Copyright (C) 2016 John Wiley & Sons, Ltd.
Over short and long timescales, wildfires can be an important cause of hydrological and geomorphological change. Mediterranean rivers are part of a fire-prone and high-energy environment in which the timing of major storms in relation to fire influences the impact on fluvial systems; accordingly, the identification of major sources, stores and fluxes of sediments is essential for providing more effective post-fire management strategies. In this study, caesium-137 and excess lead-210 were used as tracers to quantify the proportional contributions of fine sediment from hillslope surface and channel bank sources to suspended sediment and channel bed deposits before the impact of a forest wildfire in Na Borges, a Mediterranean groundwater-dominated river. It also compared burnt and unburnt spatial sources of sediment within a single catchment and the extent to which burnt material was transported downstream. The study focused on two small and steep sub-catchments, where just one of the catchments was partially affected by a wildfire. The pre-fire dynamics indicated that surface soils were the main sediment source in these ephemeral creeks. Post-fire sediment dynamics were characterised by a single flood event with a short recurrence interval (i.e. return period ca. <1 year). Sediment generated from the burnt area contributed 12% on average to bed-stored sediments within the burnt catchment, which reduced downstream to 5% along the main channel of the Na Borges River. The findings demonstrate the potential for using fallout radionuclide tracers to understand the wider impacts of wildfires on fluvial environments located outside of the burn area. Copyright (C) 2015 John Wiley & Sons, Ltd.
Despite a large number of dated glacial landforms in the Himalaya, the ice extent during the global Last Glacial Maximum (LGM) from 19 to 23 ka is only known to first order. New cosmogenic Be-10 exposure ages from well-preserved glacially polished surfaces, combined with published data, and an improved production rate scaling model allow reconstruction of the LGM ice extent and subsequent deglaciation in the Chandra Valley of NW India. We show that a >1000 m thick valley glacier retreated >150 km within a few thousand years after the onset of LGM deglaciation. By comparing the recession of the Chandra Valley Glacier and other Himalayan glaciers with those of Northern and Southern Hemisphere glaciers, we demonstrate that post-LGM deglaciation was similar and nearly finished prior to the Bolling/Allerod interstadial. Our study supports the view that many Himalayan glaciers advanced during the LGM, likely in response to global variations in temperature.
In the past, floods were basically managed by flood control mechanisms. The focus was set on the reduction of flood hazard. The potential consequences were of minor interest. Nowadays river flooding is increasingly seen from the risk perspective, including possible consequences. Moreover, the large-scale picture of flood risk became increasingly important for disaster management planning, national risk developments and the (re-) insurance industry. Therefore, it is widely accepted that risk-orientated flood management ap-proaches at the basin-scale are needed. However, large-scale flood risk assessment methods for areas of several 10,000 km² are still in early stages. Traditional flood risk assessments are performed reach wise, assuming constant probabilities for the entire reach or basin. This might be helpful on a local basis, but where large-scale patterns are important this approach is of limited use. Assuming a T-year flood (e.g. 100 years) for the entire river network is unrealistic and would lead to an overestimation of flood risk at the large scale. Due to the lack of damage data, additionally, the probability of peak discharge or rainfall is usually used as proxy for damage probability to derive flood risk. With a continuous and long term simulation of the entire flood risk chain, the spatial variability of probabilities could be consider and flood risk could be directly derived from damage data in a consistent way.
The objective of this study is the development and application of a full flood risk chain, appropriate for the large scale and based on long term and continuous simulation. The novel approach of ‘derived flood risk based on continuous simulations’ is introduced, where the synthetic discharge time series is used as input into flood impact models and flood risk is directly derived from the resulting synthetic damage time series.
The bottleneck at this scale is the hydrodynamic simu-lation. To find suitable hydrodynamic approaches for the large-scale a benchmark study with simplified 2D hydrodynamic models was performed. A raster-based approach with inertia formulation and a relatively high resolution of 100 m in combination with a fast 1D channel routing model was chosen.
To investigate the suitability of the continuous simulation of a full flood risk chain for the large scale, all model parts were integrated into a new framework, the Regional Flood Model (RFM). RFM consists of the hydrological model SWIM, a 1D hydrodynamic river network model, a 2D raster based inundation model and the flood loss model FELMOps+r. Subsequently, the model chain was applied to the Elbe catchment, one of the largest catchments in Germany. For the proof-of-concept, a continuous simulation was per-formed for the period of 1990-2003. Results were evaluated / validated as far as possible with available observed data in this period. Although each model part introduced its own uncertainties, results and runtime were generally found to be adequate for the purpose of continuous simulation at the large catchment scale.
Finally, RFM was applied to a meso-scale catchment in the east of Germany to firstly perform a flood risk assessment with the novel approach of ‘derived flood risk assessment based on continuous simulations’. Therefore, RFM was driven by long term synthetic meteorological input data generated by a weather generator. Thereby, a virtual time series of climate data of 100 x 100 years was generated and served as input to RFM providing subsequent 100 x 100 years of spatially consistent river discharge series, inundation patterns and damage values. On this basis, flood risk curves and expected annual damage could be derived directly from damage data, providing a large-scale picture of flood risk. In contrast to traditional flood risk analysis, where homogenous return periods are assumed for the entire basin, the presented approach provides a coherent large-scale picture of flood risk. The spatial variability of occurrence probability is respected. Additionally, data and methods are consistent. Catchment and floodplain processes are repre-sented in a holistic way. Antecedent catchment conditions are implicitly taken into account, as well as physical processes like storage effects, flood attenuation or channel–floodplain interactions and related damage influencing effects. Finally, the simulation of a virtual period of 100 x 100 years and consequently large data set on flood loss events enabled the calculation of flood risk directly from damage distributions. Problems associated with the transfer of probabilities in rainfall or peak runoff to probabilities in damage, as often used in traditional approaches, are bypassed.
RFM and the ‘derived flood risk approach based on continuous simulations’ has the potential to provide flood risk statements for national planning, re-insurance aspects or other questions where spatially consistent, large-scale assessments are required.
Understanding how compression affects the distribution of liquid water and gaseous oxygen in the polymer electrolyte membrane fuel cell gas diffusion layer (GDL) is vital for informing the design of improved porous materials for effective water management strategies. Pore networks extracted from synchrotron-based micro-computed tomography images of compressed GDLs were employed to simulate liquid water transport in GDL materials over a range of compression pressures. The oxygen transport resistance was predicted for each sample under dry and partially saturated conditions. A favorable GDL compression value for a preferred liquid water distribution and oxygen diffusion was found for Toray TGP-H-090 (10%), yet an optimum compression value was not recognized for SGL Sigracet 25BC. SGL Sigracet 25BC exhibited lower transport resistance values compared to Toray TGP-H-090, and this is attributed to the additional diffusion pathways provided by the microporous layer (MPL), an effect that is particularly significant under partially saturated conditions. (C) 2016 Elsevier B.V. All rights reserved.
Carbonatites are peculiar magmatic rocks with mantle-related genesis, commonly interpreted as the products of melting of CO2-bearing peridotites, or resulting from the chemical evolution of mantle derived magmas, either through extreme "differentiation or secondary immiscibility. Here we report the first finding of anatectic carbonatites of crustal origin, preserved as calcite-rich polycrystalline inclusions in garnet from low-to-medium pressure migmatites of the Oberpfalz area, SW Bohemian Massif (Central Europe). These inclusions originally trapped a melt of calciocarbonatitic composition with a characteristic enrichment in Ba, Sr and LREE. This interpretation is supported by the results of a detailed microstructural and microchemical investigation, as well as re-melting experiments using a piston cylinder apparatus. Carbonatitic inclusions coexist in the same cluster with crystallized silicate melt inclusions (nanogranites) and COH fluid inclusions, suggesting conditions of primary immiscibility between two melts and a fluid during anatexis. The production of both carbonatitic and granitic melts during the same anatectic event requires a suitable heterogeneous protolith. This may be represented by a sedimentary sequence containing marble lenses of limited extension, similar to the one still visible in the adjacent central Moldanubian Zone. The presence of CO2-rich fluid inclusions suggests furthermore that high CO2 activity during anatexis may be required to stabilize a carbonate-rich melt in a silica-dominated, system. This natural occurrence displays a remarkable similarity with experiments on carbonate-silicate melt immiscibility, where CO2 saturation is a condition commonly imposed. In conclusion, this study shows how the investigation of partial melting through melt inclusion studies may unveil unexpected processes whose evidence, while preserved in stiff minerals such as garnet, is completely obliterated in the rest of the rock due to metamorphic re-equilibration. Our results thus provide invaluable new insights into the processes which shape the geochemical evolution of our planet, such as the redistribution of carbon and strategic metals during orogenesis. (C) 2016 Elsevier B.V. All rights reserved.
A unique assemblage including kumdykolite and kokchetavite, polymorphs of albite and K-feldspar, respectively, together with cristobalite, micas, and calcite has been identified in high-pressure granulites of the Orlica-Snieznik dome (Bohemian Massif) as the product of partial melt crystallization in preserved nanogranites. Previous reports of both kumdykolite and kokchetavite in natural rocks are mainly from samples that passed through the diamond stability field. However, because the maximum pressure recorded in these host rocks is <3 GPa, our observations indicate that high pressure is not required for the formation of kumdykolite and kokchetavite, and their presence is not therefore an indicator of ultrahigh-pressure conditions. Detailed microstructural and microchemical investigation of these inclusions indicates that such phases should instead be regarded as (1) a direct mineralogical criteria to identify former melt inclusions with preserved original compositions, including H2O and CO2 contents and (2) indicators of rapid cooling of the host rocks. Thus, the present study provides novel criteria for the interpretation of melt inclusions in natural rocks and allows a more rigorous characterization of partial melts during deep subduction to mantle depth as well as their behavior on exhumation.
The common carp (Cyprinus carpio) is a highly invasive species and an ecological engineer. It has been repeatedly shown to increase nutrient concentrations and phytoplankton biomass while destroying submerged macrophytes, although there are few studies from the Mediterranean region. We studied its impact on the zooplankton community in Laguna de Medina lake, a shallow lake in Jerez de la Frontera, south-west Spain. Carp were removed with rotenone in 2007 but returned in 2010-2011. We compared zooplankton sampled monthly from 8 points from May to December in 2008 (without carp) and 2012 (with carp). Extensive macrophyte beds present in 2008 were absent in 2012. As expected, chlorophyll-a concentrations, turbidity, total suspended solids and total phosphorus were much higher in 2012. Zooplankton richness decreased from 21 taxa in 2008 to 8 taxa in 2012, accompanied by a decrease in Shannon-Wiener diversity, an increase in Evenness and a change in size distribution with loss of larger taxa. In 2008, the crustaceans were dominated by the macrocladocerans Daphnia magna and Moina brachiata and the large calanoid copepod Arctodiaptomus salinus. In 2012, these three taxa were completely absent and the zooplankton was dominated by the alien cyclopoid Acanthocyclops americanus and the rotifers Brachionus plicatilis and Keratella quadrata. Our results confirm the disappearance of macrocladocera reported by others in mesocosm experiments with carp, and suggest that alien carp facilitate the spread of the alien copepod A. americanus.
Human development has far-reaching impacts on the surface of the globe. The transformation of natural land cover occurs in different forms, and urban growth is one of the most eminent transformative processes. We analyze global land cover data and extract cities as defined by maximally connected urban clusters. The analysis of the city size distribution for all cities on the globe confirms Zipf’s law. Moreover, by investigating the percolation properties of the clustering of urban areas we assess the closeness to criticality for various countries. At the critical thresholds, the urban land cover of the countries undergoes a transition from separated clusters to a gigantic component on the country scale. We study the Zipf-exponents as a function of the closeness to percolation and find a systematic dependence, which could be the reason for deviating exponents reported in the literature. Moreover, we investigate the average size of the clusters as a function of the proximity to percolation and find country specific behavior. By relating the standard deviation and the average of cluster sizes—analogous to Taylor’s law—we suggest an alternative way to identify the percolation transition. We calculate spatial correlations of the urban land cover and find long-range correlations. Finally, by relating the areas of cities with population figures we address the global aspect of the allometry of cities, finding an exponent δ ≈ 0.85, i.e., large cities have lower densities.
The Greater and Lesser Caucasus mountains and their associated foreland basins contain similar rock types, experience a similar two-fold, along-strike variation in mean annual precipitation, and were affected by extreme base-level drops of the neighboring Caspian Sea. However, the two Caucasus ranges are characterized by decidedly different tectonic regimes and rates of deformation that are subject to moderate (less than an order of magnitude) gradients in climate, and thus allow for a unique opportunity to isolate the effects of climate and tectonics in the evolution of topography within active orogens. There is an apparent disconnect between modern climate, shortening rates, and topography of both the Greater Caucasus and Lesser Caucasus which exhibit remarkably similar topography along-strike despite the gradients in forcing. By combining multiple datasets, we examine plausible causes for this disconnect by presenting a detailed analysis of the topography of both ranges utilizing established relationships between catchment-mean erosion rates and topography (local relief, hillslope gradients, and channel steepness) and combining it with a synthesis of previously published low-temperature thermochronologic data. Modern climate of the Caucasus region is assessed through an analysis of remotely-sensed data (TRMM and MODIS) and historical streamflow data. Because along-strike variation in either erosional efficiency or thickness of accreted material fail to explain our observations, we suggest that the topography of both the western Lesser and Greater Caucasus are partially supported by different geodynamic forces. In the western Lesser Caucasus, high relief portions of the landscape likely reflect uplift related to ongoing mantle lithosphere delamination beneath the neighboring East Anatolian Plateau. In the Greater Caucasus, maintenance of high topography in the western portion of the range despite extremely low (<2-4 mm/y) modern convergence rates may be related to dynamic topography from detachment of the north-directed Greater Caucasus slab or to a recent slowing of convergence rates. Large-scale spatial gradients in climate are not reflected in the topography of the Caucasus and do not seem to exert any significant control on the tectonics or structure of either range. (C) 2016 Elsevier B.V. All rights reserved.
A Paleoarchean impact spherule-bearing interval of the 763 m long International Continental Scientific Drilling Program (ICDP) drill core BARB5 from the lower Mapepe Formation of the Fig Tree Group, Barberton Mountain Land (South Africa) was investigated using nondestructive analytical techniques. The results of visual observation, infrared (IR) spectroscopic imaging, and micro-X-ray fluorescence (lXRF) of drill cores are presented. Petrographic and sedimentary features, as well as major and trace element compositions of lithologies from the micrometer to kilometer-scale, assisted in the localization and characterization of eight spherule-bearing intervals between 512.6 and 510.5 m depth. The spherule layers occur in a strongly deformed section between 517 and 503 m, and the rocks in the core above and below are clearly less disturbed. The lXRF element maps show that spherule layers have similar petrographic and geochemical characteristics but differences in (1) sorting of two types of spherules and (2) occurrence of primary minerals (Ni-Cr spinel and zircon). We favor a single impact scenario followed by postimpact reworking, and subsequent alteration. The spherule layers are Al2O3-rich and can be distinguished from the Al2O3-poor marine sediments by distinct Al-OH absorption features in the short wave infrared (SWIR) region of the electromagnetic spectrum. Infrared images can cover tens to hundreds of square meters of lithologies and, thus, may be used to search for Al-OH-rich spherule layers in Al2O3-poor sediments, such as Eoarchean metasediments, where the textural characteristics of the spherule layers are obscured by metamorphism.
Ice-wedge polygon (IWP) peatlands in the Arctic and Subarctic are extremely vulnerable to climatic and environmental change. We present the results of a multidisciplinary paleoenvironmental study on IWPs in the northern Yukon, Canada. High-resolution laboratory analyses were carried out on a permafrost core and the overlying seasonally thawed (active) layer, from an IWP located in a drained lake basin on Herschel Island. In relation to 14 Accelerator Mass Spectrometry (AMS) radiocarbon dates spanning the last 5000 years, we report sedimentary data including grain size distribution and biogeochemical parameters (organic carbon, nitrogen, C/N ratio, delta C-13), stable water isotopes (delta O-18, delta D), as well as fossil pollen, plant macrofossil and diatom assemblages. Three sediment units (SUS) correspond to the main stages of deposition (1) in a thermokarst lake (SW : 4950 to 3950 cal yrs BP), (2) during transition from lacustrine to palustrine conditions after lake drainage (SU2: 3950 to 3120 cal yrs BP), and (3) in palustrine conditions of the IWP field that developed after drainage (SU3: 3120 cal yrs BP to 2012 CE). The lacustrine phase (pre 3950 cal yrs BP) is characterized by planktonic-benthic and pioneer diatom species indicating circumneutral waters, and very few plant macrofossils. The pollen record has captured a regional signal of relatively stable vegetation composition and climate for the lacustrine stage of the record until 3950 cal yrs BP. Palustrine conditions with benthic and acidophilic diatom species characterize the peaty shallow-water environments of the low-centered IWP. The transition from lacustrine to palustrine conditions was accompanied by acidification and rapid revegetation of the lake bottom within about 100 years. Since the palustrine phase we consider the pollen record as a local vegetation proxy dominated by the plant communities growing in the IWP. Ice-wedge cracking in water-saturated sediments started immediately after lake drainage at about 3950 cal yrs BP and led to the formation of an IWP mire. Permafrost aggradation through downward closed-system freezing of the lake talik is indicated by the stable water isotope record. The originally submerged IWP center underwent gradual drying during the past 2000 years. This study highlights the sensitivity of permafrost landscapes to climate and environmental change throughout the Holocene. (C) 2016 Elsevier Ltd. All rights reserved.
Leaf wax n-alkane H-2 values carry important information about environmental and ecophysiological processes in plants. However, the physiological and biochemical drivers that shape leaf wax n-alkane H-2 values are not completely understood. It is particularly unclear why n-alkanes in grasses are typically H-2-depleted compared with plants from other taxonomic groups such as dicotyledonous plants and why C3 grasses are H-2-depleted compared with C4 grasses. To resolve these uncertainties, we quantified the effects of leaf water evaporative H-2-enrichment and biosynthetic hydrogen isotope fractionation on n-alkane H-2 values for a range of C3 and C4 grasses grown in climate-controlled chambers. We found that only a fraction of leaf water evaporative H-2-enrichment is imprinted on the leaf wax n-alkane H-2 values in grasses. This is interesting, as previous studies have shown in dicotyledonous plants a nearly complete transfer of this H-2-enrichment to the n-alkane H-2 values. We thus infer that the typically observed H-2-depletion of n-alkanes in grasses (as opposed to dicots) is because only a fraction of the leaf water evaporative H-2-enrichment is imprinted on the H-2 values. Our experiments also show that differences in n-alkane H-2 values between C3 and C4 grasses are largely the result of systematic differences in biosynthetic fractionation between these two plant groups, which was on average -198 and-159 parts per thousand for C3 and C4 grasses, respectively. We present novel and exciting data on how leaf wax n-alkane 2H values from grasses are affected by plant physiological (leaf water evaporative 2H-enrichment) or biochemical processes (biosynthetic hydrogen isotope fractionation). These results are very interesting because they shed new light on how naturally observed differences between 2H values from C3 and C4 grasses and dicots can be explained by systematic differences in the biosynthesis of n-alkanes between these plant group (i.e. largely driven by NADPH origins).
We applied the geomorphic indices (hypsometry and stream length gradient) to evaluate the differential uplift of the central and southern Longmenshan, a mountain range characterized by rapid erosion, strong tectonic uplift, and devastating seismic hazards. The results of the geomorphic analysis indicate that the Beichuan-Yingxiu fault and the Shuangshi-Dachuan fault act as major tectonic boundaries separating areas experiencing rapid uplift from slow uplift. The results of the geomorphic analysis also suggest that the Beichuan-Yingxiu fault is the most active fault with the largest relative uplift rates compared to the rest of the faults in the Longmenshan fault system. We compared reflected relative uplift rates based on the hypsometry and stream length gradient indices with geological/geodetic absolute rates. Along-strike and across-strike variations in the hypsometry and stream length gradient correlate with the spatial patterns derived from the apatite fission track exhumation rates, the leveling-derived uplift rate, and coseismic vertical displacements during the 2008 Wenchuan earthquake. These data defined multiple fault relationships in a complex thrust zone and provided geomorphic evidence to evaluate the potential seismic hazards of the southern Longmenshan range.
The main scope of the InterPACIFIC (Intercomparison of methods for site parameter and velocity profile characterization) project is to assess the reliability of in-hole and surface-wave methods, used for estimating shear wave velocity. Three test-sites with different subsurface conditions were chosen: a soft soil, a stiff soil and a rock outcrop. This paper reports the surface-wave methods results. Specifically 14 teams of expert users analysed the same experimental surface-wave datasets, consisting of both passive and active data. Each team adopted their own strategy to retrieve the dispersion curve and the shear-wave velocity profile at each site. Despite different approaches, the dispersion curves are quite in agreement with each other. Conversely, the shear-wave velocity profiles show a certain variability that increases in correspondence of major stratigraphic interfaces. This larger variability is mainly due to non-uniqueness of the solution and lateral variability. As expected, the observed variability in V-s,V-30 estimatesis small, as solution non-uniqueness plays a limited role. (C) 2015 Elsevier Ltd. All rights reserved.
The collision of bathymetric anomalies, such as oceanic spreading centers, at convergent plate margins can profoundly affect subduction dynamics, magmatism, and the structural and geomorphic evolution of the overriding plate. The Southern Patagonian Andes of South America are a prime example for sustained oceanic ridge collision and the successive formation and widening of an extensive asthenospheric slab window since the Middle Miocene. Several of the predicted upper-plate geologic manifestations of such deep-seated geodynamic processes have been studied in this region, but many topics remain highly debated. One of the main controversial topics is the interpretation of the regional low-temperature thermochronology exhumational record and its relationship with tectonic and/or climate-driven processes, ultimately manifested and recorded in the landscape evolution of the Patagonian Andes. The prominent along-strike variance in the topographic characteristics of the Andes, combined with coupled trends in low-temperature thermochronometer cooling ages have been interpreted in very contrasting ways, considering either purely climatic (i.e. glacial erosion) or geodynamic (slab-window related) controlling factors.
This thesis focuses on two main aspects of these controversial topics. First, based on field observations and bedrock low-temperature thermochronology data, the thesis addresses an existing research gap with respect to the neotectonic activity of the upper plate in response to ridge collision - a mechanism that has been shown to affect the upper plate topography and exhumational patterns in similar tectonic settings. Secondly, the qualitative interpretation of my new and existing thermochronological data from this region is extended by inverse thermal modelling to define thermal histories recorded in the data and evaluate the relative importance of surface vs. geodynamic factors and their possible relationship with the regional cooling record.
My research is centered on the Northern Patagonian Icefield (NPI) region of the Southern Patagonian Andes. This site is located inboard of the present-day location of the Chile Triple Junction - the juncture between the colliding Chile Rise spreading center and the Nazca and Antarctic Plates along the South American convergent margin. As such this study area represents the region of most recent oceanic-ridge collision and associated slab window formation. Importantly, this location also coincides with the abrupt rise in summit elevations and relief characteristics in the Southern Patagonian Andes. Field observations, based on geological, structural and geomorphic mapping, are combined with bedrock apatite (U-Th)/He and apatite fission track (AHe and AFT) cooling ages sampled along elevation transects across the orogen. This new data reveals the existence of hitherto unrecognized neotectonic deformation along the flanks of the range capped by the NPI.
This deformation is associated with the closely spaced oblique collision of successive oceanic-ridge segments in this region over the past 6 Ma. I interpret that this has caused a crustal-scale partitioning of deformation and the decoupling, margin-parallel migration, and localized uplift of a large crustal sliver (the NPI block) along the subduction margin. The location of this uplift coincides with a major increase of summit elevations and relief at the northern edge of the NPI massif. This mechanism is compatible with possible extensional processes along the topographically subdued trailing edge of the NPI block as documented by very recent and possibly still active normal faulting. Taken together, these findings suggest a major structural control on short-wavelength variations in topography in the Southern Patagonian Andes - the region affected by ridge collision and slab window formation.
The second research topic addressed here focuses on using my new and existing bedrock low-temperature cooling ages in forward and inverse thermal modeling. The data was implemented in the HeFTy and QTQt modeling platforms to constrain the late Cenozoic thermal history of the Southern Patagonian Andes in the region of the most recent upper-plate sectors of ridge collision. The data set combines AHe and AFT data from three elevation transects in the region of the Northern Patagonian Icefield. Previous similar studies claimed far-reaching thermal effects of the approaching ridge collision and slab window to affect patterns of Late Miocene reheating in the modelled thermal histories. In contrast, my results show that the currently available data can be explained with a simpler thermal history than previously proposed. Accordingly, a reheating event is not needed to reproduce the observations. Instead, the analyzed ensemble of modelled thermal histories defines a Late Miocene protracted cooling and Pliocene-to-recent stepwise exhumation. These findings agree with the geological record of this region. Specifically, this record indicates an Early Miocene phase of active mountain building associated with surface uplift and an active fold-and-thrust belt, followed by a period of stagnating deformation, peneplanation, and lack of synorogenic deposition in the Patagonian foreland. The subsequent period of stepwise exhumation likely resulted from a combination of pulsed glacial erosion and coeval neotectonic activity. The differences between the present and previously published interpretation of the cooling record can be reconciled with important inconsistencies of previously used model setup. These include mainly the insufficient convergence of the models and improper assumptions regarding the geothermal conditions in the region. This analysis puts a methodological emphasis on the prime importance of the model setup and the need for its thorough examination to evaluate the robustness of the final outcome.
The subduction of bathymetric anomalies at convergent margins can profoundly affect subduction dynamics, magmatism, and the structural and geomorphic evolution of the overriding plate. The Northern Patagonian Icefield (NPI) is located east of the Chile Triple Junction at similar to 47 degrees S, where the Chile Rise spreading center collides with South America. This region is characterized by an abrupt increase in summit elevations and relief that has been controversially debated in the context of geodynamic versus glacial erosion effects on topography. Here we present geomorphic, thermochronological, and structural data that document neotectonic activity along hitherto unrecognized faults along the flanks of the NPI. New apatite (U-Th)/He bedrock cooling ages suggest faulting since 2-3 Ma. We infer the northward translation of an similar to 140 km long fore-arc sliver-the NPI block-results from enhanced partitioning of oblique plate convergence due to the closely spaced collision of three successive segments of the Chile Rise. In this model, greater uplift occurs in the hanging wall of the Exploradores thrust at the northern leading edge of the NPI block, whereas the Cachet and Liquine-Ofqui dextral faults decouple the NPI block along its eastern and western flanks, respectively. Localized extension possibly occurs at its southern trailing edge along normal faults associated with margin-parallel extension, tectonic subsidence, and lower elevations along the Andean crest line. Our neotectonic model provides a novel explanation for the abrupt topographic variations inland of the Chile Triple Junction and emphasizes the fundamental effects of local tectonics on exhumation and topographic patterns in this glaciated landscape.
Hummocky soil landscapes are characterized by 3D spatial patterns of soil types that result from erosion-affected pedogenesis. Due to tillage and water erosion, truncated profiles have been formed at steep and mid slopes and colluvial soils at hollows, while intact profiles remained at plateau positions. Pedogenetic variations in soil horizons lead to spatial differences in the soil water balance at hillslope positions. Here, possible interactions between erosion affected soil properties, the water balances, and the crop growth and feedback effects of erosion on the leaching rates were assumed. The hypothesis was tested by water balance simulations comparing uniform with hillslope position-specific crop and root growths for soils at plateau, flat mid slope, steep slope, and hollow using the Hydrus-1D program. The boundary condition data were monitored at the CarboZALF-D experimental field site, which was cropped with perennial lucerne (Medicago sativa L.) in 2013 and 2014. Crop and root growth at the four hillslope positions was assumed proportional to observed leaf area index (LAI). Fluxes of dissolved organic and inorganic carbon (DOC, DIC) were obtained from simulated water fluxes and measured DOC and DIC concentrations. For the colluvic soil at hollow, the crop growth was initially highest and later limited by an increasing water table; here the predominately upward flow led to a net input in DIC and DOC. For the truncated soils at steep slopes, simulations support the hypothesis that reduced crop growth caused an increase in percolation and DIC leaching from the subsoil horizons, which in turn led to accelerated soil development and more soil variations along eroding hillslopes in arable soil landscapes. (C) 2015 Elsevier B.V. All rights reserved.
Variational methods are employed in situations where exact Bayesian inference becomes intractable due to the difficulty in performing certain integrals. Typically, variational methods postulate a tractable posterior and formulate a lower bound on the desired integral to be approximated, e.g. marginal likelihood. The lower bound is then optimised with respect to its free parameters, the so-called variational parameters. However, this is not always possible as for certain integrals it is very challenging (or tedious) to come up with a suitable lower bound. Here, we propose a simple scheme that overcomes some of the awkward cases where the usual variational treatment becomes difficult. The scheme relies on a rewriting of the lower bound on the model log-likelihood. We demonstrate the proposed scheme on a number of synthetic and real examples, as well as on a real geophysical model for which the standard variational approaches are inapplicable.
Adequate irrigation inputs are essential for the application of hydrological models in irrigated catchments, but reliable data on both the amount and the frequency of irrigation applications are often missing at an appropriate spatial scale. In this paper, we demonstrate and test approaches to estimate irrigation inputs for distributed hydrological modelling. In this context, the Soil and Water Assessment Tool was applied to simulate water balances for an irrigated catchment in southeast Australia during the period 2008–2010. Two methods for estimating irrigation inputs were tested. One method was based on a fixed irrigation application rate, whereas the other one had variable irrigation rates depending on season and the irrigated crop. These two approaches were also compared with the ‘auto-irrigation’ method within the Soil and Water Assessment Tool model. The method with variable irrigation rates resulted in the most reasonable interpretation of the readily available irrigation data, consistent estimates of irrigation runoff coefficients throughout the year and the best fit to observed data on both drain flows at the catchment outlet and spatial evapotranspiration patterns. We also found that the different irrigation inputs significantly affected simulated water balances, in particular deep percolation under relatively dry climatic conditions. All these results suggest that it is possible to infer irrigation inputs from readily available data and local knowledge, adequate for hydrological modelling in irrigated catchments. Our study also demonstrates that, in order to predict reliable water balances in irrigated catchments, an accurate knowledge of irrigation scheduling and irrigation runoff is required. Copyright © 2015 John Wiley & Sons, Ltd.
Decreasing groundwater levels in many parts of Germany and decreasing low flows in Central Europe have created a need for adaptation measures to stabilize the water balance and to increase low flows. The objective of our study was to estimate the impact of ditch water level management on stream-aquifer interactions in small lowland catchments of the mid-latitudes. The water balance of a ditch-irrigated area and fluxes between the subsurface and the adjacent stream were modeled for three runoff recession periods using the Hydrus-2D software package. The results showed that the subsurface flow to the stream was closely related to the difference between the water level in the ditch system and the stream. Evapotranspiration during the growing season additionally reduced base flow. It was crucial to stop irrigation during a recession period to decrease water withdrawal from the stream and enhance the base flow by draining the irrigated area. Mean fluxes to the stream were between 0.04 and 0.64 ls(-1) for the first 20 days of the low-flow periods. This only slightly increased the flow in the stream, whose mean was 57 ls(-1) during the period with the lowest flows. Larger areas would be necessary to effectively increase flows in mesoscale catchments.
Decreasing groundwater levels in many parts of Germany and decreasing low flows in Central Europe have created a need for adaptation measures to stabilize the water balance and to increase low flows. The objective of our study was to estimate the impact of ditch water level management on stream-aquifer interactions in small lowland catchments of the mid-latitudes. The water balance of a ditch-irrigated area and fluxes between the subsurface and the adjacent stream were modeled for three runoff recession periods using the Hydrus-2D software package. The results showed that the subsurface flow to the stream was closely related to the difference between the water level in the ditch system and the stream. Evapotranspiration during the growing season additionally reduced base flow. It was crucial to stop irrigation during a recession period to decrease water withdrawal from the stream and enhance the base flow by draining the irrigated area. Mean fluxes to the stream were between 0.04 and 0.64 ls(-1) for the first 20 days of the low-flow periods. This only slightly increased the flow in the stream, whose mean was 57 ls(-1) during the period with the lowest flows. Larger areas would be necessary to effectively increase flows in mesoscale catchments.
The lakes in the Kenyan Rift Valley offer the unique opportunity to study a wide range of hydrochemical environmental conditions, ranging from freshwater to highly saline and alkaline lakes. Because little is known about the hydro- and biogeochemical conditions in the underlying lake sediments, it was the aim of this study to extend the already existing data sets with data from porewater and biomarker analyses. Additionally, reduced sulphur compounds and sulphate reduction rates in the sediment were determined. The new data was used to examine the anthropogenic and microbial influence on the lakes sediments as well as the influence of the water chemistry on the degradation and preservation of organic matter in the sediment column. The lakes discussed in this study are: Logipi, Eight (a small crater lake in the region of Kangirinyang), Baringo, Bogoria, Naivasha, Oloiden, and Sonachi.
The biomarker compositions were similar in all studied lake sediments; nevertheless, there were some differences between the saline and freshwater lakes. One of those differences is the occurrence of a molecule related to β-carotene, which was only found in the saline lakes. This molecule most likely originates from cyanobacteria, single-celled organisms which are commonly found in saline lakes. In the two freshwater lakes, stigmasterol, a sterol characteristic for freshwater algae, was found. In this study, it was shown that Lakes Bogoria and Sonachi can be used for environmental reconstructions with biomarkers, because the absence of oxygen at the lake bottoms slowed the degradation process. Other lakes, like for example Lake Naivasha, cannot be used for such reconstructions, because of the large anthropogenic influence. But the biomarkers proved to be a useful tool to study those anthropogenic influences. Additionally, it was observed that horizons with a high concentration of elemental sulphur can be used as temporal markers. Those horizons were deposited during times when the lake levels were very low. The sulphur was deposited by microorganisms which are capable of anoxygenic photosynthesis or sulphide oxidation.
Zircon (ZrSiO4), hafnon (HfSiO4) and five intermediate compositions were synthesized from a Pb silicate melt. The resulting crystals were 20-300 mu m in size and displayed sector and growth zoning. Raman spectra were acquired at locations in the sample for which preceding electron microprobe (EMP) analyses revealed sufficient compositional homogeneity. The dataset documents shifts of Raman bands with changing composition. In this study, bands that have previously not been reported were found for the intermediate compositions and for pure hafnon, in particular at wavenumbers less than 200 cm(-1). For these external modes, the dataset provides new insight into the compositional dependence of their frequencies. Density-functional theory calculations support the observations and are used for a detailed interpretation of the spectra. The pitfalls of the EMP analysis along the zircon-hafnon join are highlighted.
Cities play a vital role in the global climate change mitigation agenda. City population density is one of the key factors that influence urban energy consumption and the subsequent GHG emissions. However, previous research on the relationship between population density and GHG emissions led to contradictory results due to urban/rural definition conundrum and the varying methodologies for estimating GHG emissions. This work addresses these ambiguities by employing the City Clustering Algorithm (CCA) and utilizing the gridded CO2 emissions data. Our results, derived from the analysis of all inhabited areas in the US, show a sub-linear relationship between population density and the total emissions (i.e. the sum of on-road and building emissions) on a per capita basis. Accordingly, we find that doubling the population density would entail a reduction in the total CO2 emissions in buildings and on-road sectors typically by at least 42%. Moreover, we find that population density exerts a higher influence on on-road emissions than buildings emissions. From an energy consumption point of view, our results suggest that on-going urban sprawl will lead to an increase in on-road energy consumption in cities and therefore stresses the importance of developing adequate local policy measures to limit urban sprawl. (C) 2016 Elsevier Ltd. All rights reserved.
We present an algorithm that performs sequentially one-dimensional inversion of subsurface magnetic permeability and electrical conductivity by using multi-configuration electromagnetic induction sensor data. The presented method is based on the conversion of the in-phase and out-of-phase data into effective magnetic permeability and electrical conductivity of the equivalent homogeneous half-space. In the case of small-offset systems, such as portable electromagnetic induction sensors, for which in-phase and out-of-phase data are moderately coupled, the effective half-space magnetic permeability and electrical conductivity can be inverted sequentially within an iterative scheme. We test and evaluate the proposed inversion strategy using synthetic and field examples. First, we apply it to synthetic data for some highly magnetic environments. Then, the method is tested on real field data acquired in a basaltic environment to image a formation of archaeological interest. These examples demonstrate that a joint interpretation of in-phase and out-of-phase data leads to a better characterisation of the subsurface in magnetic environments such as volcanic areas.
Steppe vegetation represents a key marker of past Asian aridification and is associated with monsoonal intensification. Little is, however, known about the origin of this pre-Oligocene vegetation, its specific composition and how it changed over time and responded to climatic variations. Here, we describe the morphological characters of Ephedraceae pollen in Eocene strata of the Xining Basin and compare the pollen composition with the palynological composition of Late Cretaceous and Paleocene deposits of the Xining Basin and the Quaternary deposits of the Qaidam Basin. We find that the Late Cretaceous steppe was dominated by Gnetaceaepollenites; in the transition from the Cretaceous to the Paleocene, Gnetaceaepollenites became extinct and Ephedripites subgenus Ephedripites dominated the flora with rare occurrences of Ephedripites subgen. Distachyapites; the middle to late Eocene presents a strong increase of Ephedripites subgen. Distachyapites; and the Quaternary/Recent is marked by a significantly lower diversity of Ephedraceae (and Nitrariaceae) compared to the Eocene. In the modern landscape of China, only a fraction of the Paleogene species diversity of Ephedraceae remains and we propose that these alterations in Ephedreaceae composition occurred in response to the climatic changes at least since the Eocene. In particular, the strong Eocene monsoons that enhanced the continental aridification may have played an important role in the evolution of Ephedripites subgen. Distachyapites triggering an evolutionary shift to wind-pollination in this group. Conceivably, the Ephedraceae/Nitrariaceae dominated steppe ended during the Eocene/Oligocene climatic cooling and aridification, which favoured other plant taxa.
Oceanic lithospheric S-wave velocities from the analysis of P-wave polarization at the ocean floor
(2016)
Our knowledge of the absolute S-wave velocities of the oceanic lithosphere is mainly based on global surface wave tomography, local active seismic or compliance measurements using oceanic infragravity waves. The results of tomography give a rather smooth picture of the actual S-wave velocity structure and local measurements have limitations regarding the range of elastic parameters or the geometry of the measurement. Here, we use the P-wave polarization (apparent P-wave incidence angle) of teleseismic events to investigate the S-wave velocity structure of the oceanic crust and the upper tens of kilometres of the mantle beneath single stations. In this study, we present an up to our knowledge new relation of the apparent P-wave incidence angle at the ocean bottom dependent on the half-space S-wave velocity. We analyse the angle in different period ranges at ocean bottom stations (OBSs) to derive apparent S-wave velocity profiles. These profiles are dependent on the S-wave velocity as well as on the thickness of the layers in the subsurface. Consequently, their interpretation results in a set of equally valid models. We analyse the apparent P-wave incidence angles of an OBS data set which was collected in the Eastern Mid Atlantic. We are able to determine reasonable S-wave-velocity-depth models by a three-step quantitative modelling after a manual data quality control, although layer resonance sometimes influences the estimated apparent S-wave velocities. The apparent S-wave velocity profiles are well explained by an oceanic PREM model in which the upper part is replaced by four layers consisting of a water column, a sediment, a crust and a layer representing the uppermost mantle. The obtained sediment has a thickness between 0.3 and 0.9 km with S-wave velocities between 0.7 and 1.4 km s(-1). The estimated total crustal thickness varies between 4 and 10 km with S-wave velocities between 3.5 and 4.3 km s(-1). We find a slight increase of the total crustal thickness from similar to 5 to similar to 8 km towards the South in the direction of a major plate boundary, the Gloria Fault. The observed crustal thickening can be related with the known dominant compression in the vicinity of the fault. Furthermore, the resulting mantle S-wave velocities decrease from values around 5.5 to 4.5 km s(-1) towards the fault. This decrease is probably caused by serpentinization and indicates that the oceanic transform fault affects a broad region in the uppermost mantle. Conclusively, the presented method is useful for the estimation of the local S-wave velocity structure beneath ocean bottom seismic stations. It is easy to implement and consists of two main steps: (1) measurement of apparent P-wave incidence angles in different period ranges for real and synthetic data, and (2) comparison of the determined apparent S-wave velocities for real and synthetic data to estimate S-wave velocity-depth models.
Restoring degraded lands in rural environments that are heavily managed to meet subsistence needs is a challenge due to high rates of disturbance and resource extraction. This study investigates the efficacy of erosion control structures (ECSs) as restoration tools in the context of a watershed rehabilitation and wet meadow (bofedal) restoration program in the Bolivian Andes. In an effort to enhance water security and increase grazing stability, Aymara indigenous communities built over 15,000 check dams, 9,100 terraces, 5,300 infiltration ditches, and 35 pasture improvement trials. Communities built ECSs at different rates, and we compared vegetation change in the highest restoration management intensity, lowest restoration management intensity, and nonproject control communities. We used line transects to measure changes in vegetation cover and standing water in gullies with check dams and without check dams, and related these ground measurements to a time series (1986-2009) of normalized difference vegetation index derived from Landsat TM5 images. Evidence suggests that check dams increase bofedal vegetation and standing water at a local scale, and lead to increased greenness at a basin scale when combined with other ECSs. Watershed rehabilitation enhances ecosystem services significant to local communities (grazing stability, water security), which creates important synergies when conducting land restoration in rural development settings.
Glacial legacies on interglacial vegetation at the Pliocene-Pleistocene transition in NE Asia
(2016)
Broad-scale climate control of vegetation is widely assumed. Vegetation-climate lags are generally thought to have lasted no more than a few centuries. Here our palaeoecological study challenges this concept over glacial–interglacial timescales. Through multivariate analyses of pollen assemblages from Lake El’gygytgyn, Russian Far East and other data we show that interglacial vegetation during the Plio-Pleistocene transition mainly reflects conditions of the preceding glacial instead of contemporary interglacial climate. Vegetation–climate disequilibrium may persist for several millennia, related to the combined effects of permafrost persistence, distant glacial refugia and fire. In contrast, no effects from the preceding interglacial on glacial vegetation are detected. We propose that disequilibrium was stronger during the Plio-Pleistocene transition than during the Mid-Pliocene Warm Period when, in addition to climate, herbivory was important. By analogy to the past, we suggest today’s widespread larch ecosystem on permafrost is not in climate equilibrium. Vegetation-based reconstructions of interglacial climates used to assess atmospheric CO2–temperature relationships may thus yield misleading simulations of past global climate sensitivity.
The global carbon cycle is closely linked to Earth’s climate. In the context of continuously unchecked anthropogenic CO₂ emissions, the importance of natural CO₂ bond and carbon storage is increasing. An important biogenic mechanism of natural atmospheric CO₂ drawdown is the photosynthetic carbon fixation in plants and the subsequent longterm deposition of plant detritus in sediments.
The main objective of this thesis is to identify factors that control mobilization and transport of plant organic matter (pOM) through rivers towards sedimentation basins. I investigated this aspect in the eastern Nepalese Arun Valley. The trans-Himalayan Arun River is characterized by a strong elevation gradient (205 − 8848 m asl) that is accompanied by strong changes in ecology and climate ranging from wet tropical conditions in the Himalayan forelad to high alpine tundra on the Tibetan Plateau. Therefore, the Arun is an excellent natural laboratory, allowing the investigation of the effect of vegetation cover, climate, and topography on plant organic matter mobilization and export in tributaries along the gradient.
Based on hydrogen isotope measurements of plant waxes sampled along the Arun River and its tributaries, I first developed a model that allows for an indirect quantification of pOM contributed to the mainsetm by the Arun’s tributaries. In order to determine the role of climatic and topographic parameters of sampled tributary catchments, I looked for significant statistical relations between the amount of tributary pOM export and tributary characteristics (e.g. catchment size, plant cover, annual precipitation or runoff, topographic measures). On one hand, I demonstrated that pOMsourced from the Arun is not uniformly derived from its entire catchment area. On the other, I showed that dense vegetation is a necessary, but not sufficient, criterion for high tributary pOM export. Instead, I identified erosion and rainfall and runoff as key factors controlling pOM sourcing in the Arun Valley. This finding is supported by terrestrial cosmogenic nuclide concentrations measured on river sands along the Arun and its tributaries in order to quantify catchment wide denudation rates. Highest denudation rates corresponded well with maximum pOM mobilization and export also suggesting the link between erosion and pOM sourcing.
The second part of this thesis focusses on the applicability of stable isotope records such as plant wax n-alkanes in sediment archives as qualitative and quantitative proxy for the variability of past Indian Summer Monsoon (ISM) strength. First, I determined how ISM strength affects the hydrogen and oxygen stable isotopic composition (reported as δD and δ18O values vs. Vienna Standard Mean Ocean Water) of precipitation in the Arun Valley and if this amount effect (Dansgaard, 1964) is strong enough to be recorded in potential paleo-ISM isotope proxies. Second, I investigated if potential isotope records across the Arun catchment reflect ISM strength dependent precipitation δD values only, or if the ISM isotope signal is superimposed by winter precipitation or glacial melt. Furthermore, I tested if δD values of plant waxes in fluvial deposits reflect δD values of environmental waters in the respective catchments.
I showed that surface water δD values in the Arun Valley and precipitation δD from south of the Himalaya both changed similarly during two consecutive years (2011 & 2012) with distinct ISM rainfall amounts (~20% less in 2012). In order to evaluate the effect of other water sources (Winter-Westerly precipitation, glacial melt) and evapotranspiration in the Arun Valley, I analysed satellite remote sensing data of rainfall distribution (TRMM 3B42V7), snow cover (MODIS MOD10C1), glacial coverage (GLIMSdatabase, Global Land Ice Measurements from Space), and evapotranspiration (MODIS MOD16A2). In addition to the predominant ISM in the entire catchment I found through stable isotope analysis of surface waters indications for a considerable amount of glacial melt derived from high altitude tributaries and the Tibetan Plateau. Remotely sensed snow cover data revealed that the upper portion of the Arun also receives considerable winter precipitation, but the effect of snow melt on the Arun Valley hydrology could not be evaluated as it takes place in early summer, several months prior to our sampling campaigns. However, I infer that plant wax records and other potential stable isotope proxy archives below the snowline are well-suited for qualitative, and potentially quantitative, reconstructions of past changes of ISM strength.
It is commonly recognized that soil moisture exhibits spatial heterogeneities occurring in a wide range of scales. These heterogeneities are caused by different factors ranging from soil structure at the plot scale to land use at the landscape scale. There is an urgent need for effi-cient approaches to deal with soil moisture heterogeneity at large scales, where manage-ment decisions are usually made. The aim of this dissertation was to test innovative ap-proaches for making efficient use of standard soil hydrological data in order to assess seep-age rates and main controls on observed hydrological behavior, including the role of soil het-erogeneities.
As a first step, the applicability of a simplified Buckingham-Darcy method to estimate deep seepage fluxes from point information of soil moisture dynamics was assessed. This was done in a numerical experiment considering a broad range of soil textures and textural het-erogeneities. The method performed well for most soil texture classes. However, in pure sand where seepage fluxes were dominated by heterogeneous flow fields it turned out to be not applicable, because it simply neglects the effect of water flow heterogeneity. In this study a need for new efficient approaches to handle heterogeneities in one-dimensional water flux models was identified.
As a further step, an approach to turn the problem of soil moisture heterogeneity into a solu-tion was presented: Principal component analysis was applied to make use of the variability among soil moisture time series for analyzing apparently complex soil hydrological systems. It can be used for identifying the main controls on the hydrological behavior, quantifying their relevance, and describing their particular effects by functional averaged time series. The ap-proach was firstly tested with soil moisture time series simulated for different texture classes in homogeneous and heterogeneous model domains. Afterwards, it was applied to 57 mois-ture time series measured in a multifactorial long term field experiment in Northeast Germa-ny.
The dimensionality of both data sets was rather low, because more than 85 % of the total moisture variance could already be explained by the hydrological input signal and by signal transformation with soil depth. The perspective of signal transformation, i.e. analyzing how hydrological input signals (e.g., rainfall, snow melt) propagate through the vadose zone, turned out to be a valuable supplement to the common mass flux considerations. Neither different textures nor spatial heterogeneities affected the general kind of signal transfor-mation showing that complex spatial structures do not necessarily evoke a complex hydro-logical behavior. In case of the field measured data another 3.6% of the total variance was unambiguously explained by different cropping systems. Additionally, it was shown that dif-ferent soil tillage practices did not affect the soil moisture dynamics at all.
The presented approach does not require a priori assumptions about the nature of physical processes, and it is not restricted to specific scales. Thus, it opens various possibilities to in-corporate the key information from monitoring data sets into the modeling exercise and thereby reduce model uncertainties.
Classification of clouds, cirrus, snow, shadows and clear sky areas is a crucial step in the pre-processing of optical remote sensing images and is a valuable input for their atmospheric correction. The Multi-Spectral Imager on board the Sentinel-2's of the Copernicus program offers optimized bands for this task and delivers unprecedented amounts of data regarding spatial sampling, global coverage, spectral coverage, and repetition rate. Efficient algorithms are needed to process, or possibly reprocess, those big amounts of data. Techniques based on top-of-atmosphere reflectance spectra for single-pixels without exploitation of external data or spatial context offer the largest potential for parallel data processing and highly optimized processing throughput. Such algorithms can be seen as a baseline for possible trade-offs in processing performance when the application of more sophisticated methods is discussed. We present several ready-to-use classification algorithms which are all based on a publicly available database of manually classified Sentinel-2A images. These algorithms are based on commonly used and newly developed machine learning techniques which drastically reduce the amount of time needed to update the algorithms when new images are added to the database. Several ready-to-use decision trees are presented which allow to correctly label about 91% of the spectra within a validation dataset. While decision trees are simple to implement and easy to understand, they offer only limited classification skill. It improves to 98% when the presented algorithm based on the classical Bayesian method is applied. This method has only recently been used for this task and shows excellent performance concerning classification skill and processing performance. A comparison of the presented algorithms with other commonly used techniques such as random forests, stochastic gradient descent, or support vector machines is also given. Especially random forests and support vector machines show similar classification skill as the classical Bayesian method.
The effect of methods of statistical downscaling of daily precipitation on changes in extreme flow indices under a plausible future climate change scenario was investigated in 11 catchments selected from 9 countries in different parts of Europe. The catchments vary from 67 to 6171 km(2) in size and cover different climate zones. 15 regional climate model outputs and 8 different statistical downscaling methods, which are broadly categorized as change factor and bias correction based methods, were used for the comparative analyses. Different hydrological models were implemented in different catchments to simulate daily runoff. A set of flood indices were derived from daily flows and their changes have been evaluated by comparing their values derived from simulations corresponding to the current and future climate. Most of the implemented downscaling methods project an increase in the extreme flow indices in most of the catchments. The catchments where the extremes are expected to increase have a rainfall dominated flood regime. In these catchments, the downscaling methods also project an increase in the extreme precipitation in the seasons when the extreme flows occur. In catchments where the flooding is mainly caused by spring/summer snowmelt, the downscaling methods project a decrease in the extreme flows in three of the four catchments considered. A major portion of the variability in the projected changes in the extreme flow indices is attributable to the variability of the climate model ensemble, although the statistical downscaling methods contribute 35-60% of the total variance. (C) 2016 Elsevier B.V. All rights reserved.
Knowledge of the quality factor of near-surface materials is of fundamental interest in various applications. Attenuation can be very strong close to the surface and thus needs to be properly assessed. In recent years, several researchers have studied the retrieval of attenuation coefficients from the cross correlation of ambient seismic noise. Yet, the determination of exact amplitude information from noise-correlation functions is, in contrast to the extraction of traveltimes, not trivial. Most of the studies estimated attenuation coefficients on the regional scale and within the microseism band. In this paper, we investigate the possibility to derive attenuation coefficients from seismic noise at much shallower depths and higher frequencies (> 1 Hz). The Euroseistest area in northern Greece offers ideal conditions to study quality factor retrieval from ambient noise for different rock types. Correlations are computed between the stations of a small scale array experiment (station spacings < 2 km) that was carried out in the Euroseistest area in 2011. We employ the correlation of the coda of the correlation (C-3) method instead of simple cross correlations to mitigate the effect of uneven noise source distributions on the correlation amplitude. Transient removal and temporal flattening are applied instead of 1-bit normalization in order to retain relative amplitudes. The C-3 method leads to improved correlation results (higher signal-to-noise ratio and improved time symmetry) compared to simple cross correlations. The C-3 functions are rotated from the ZNE to the ZRT system and we focus on Love wave arrivals on the transverse component and on Love wave quality factors Q(L). The analysis is performed for selected stations being either situated on soft soil or on weathered rock. Phase slowness is extracted using a slant-stack method. Attenuation parameters are inferred by inspecting the relative amplitude decay of Love waves with increasing interstation distance. We observe that the attenuation coefficient gamma and Q(L) can be reliably extracted for stations situated on soft soil whereas the derivation of attenuation parameters is more problematic for stations that are located on weathered rock. The results are in acceptable conformance with theoretical Love wave attenuation curves that were computed using 1-D shear wave velocity and quality factor profiles from the Euroseistest area.
Rainfall-induced attenuation is a major source of underestimation for radar-based precipitation estimation at C-band. Unconstrained gate-by-gate correction procedures are known to be inherently unstable and thus not suited for unsupervised attenuation correction. In this study, we evaluate three different procedures to constrain gate-by-gate attenuation correction using reflectivity as the only input. These procedures are benchmarked against rainfall estimates from uncorrected radar data, using six years of radar observations from the single-polarized C-band radar in South-West Germany. The precipitation estimation error is obtained by comparing the radar-based estimates to rain gauge observations. All attenuation correction procedures benchmarked in this study lead to an effective improvement of precipitation estimation. The first method caps the corrections if the rain intensity increase exceeds a factor of two. The second method decreases the parameters of the attenuation correction iteratively for every radar beam calculation until attaining a stability criterion. The second method outperforms the first method and leads to a consistent distribution of path-integrated attenuation along the radar beam. As a third method, we propose a slight modification of Kraemer's approach which allows users to exert better control over attenuation correction by introducing an additional constraint that prevents unplausible corrections in cases of dramatic signal losses.
The molecular biomarker composition of two sediment cores from Sanabria Lake (NW Iberian Peninsula) and a survey of modern plants in the watershed provide a reconstruction of past vegetation and landscape dynamics since deglaciation. During a proglacial stage in Lake Sanabria (prior to 14.7 cal ka BP), very low biomarker concentration and carbon preference index (CPI) values similar to 1 suggest that the n-alkanes could have derived from eroded ancient sediment sources or older organic matter with high degree of maturity. During the Late glacial (14.7-11.7 cal ka BP) and the Holocene (last 11.7 cal ka BP) intervals with higher biomarker and triterpenoid concentrations (high %nC(29) , nC(31) alkanes), higher CPI and average carbon length (ACL), and lower P-aq (proportion of aquatic plants) are indicative of major contribution of vascular land plants from a more forested watershed (e.g. Mid Holocene period 7.0-4.0 cal ka BP). Lower biomarker concentrations (high %nC(27) alkanes), CPI and ACL values responded to short phases with decreased allochthonous contribution into the lake that correspond to centennial-scale periods of regional forest decline (e.g. 4-3 ka BP, Roman deforestation after 2.0 ka, and some phases of the LIA, seventeenth-nineteenth centuries). Human activities in the watershed were significant during early medieval times (1.3-1.0 cal ka BP) and since 1960 CE, in both cases associated with relatively higher productivity stages in the lake (lower biomarker and triterpenoid concentrations, high %nC(23) and %nC(31) respectively, lower ACL and CPI values and higher P-aq). The lipid composition of Sanabria Lake sediments indicates a major allochthonous (watershed-derived) contribution to the organic matter budget since deglaciation, and a dominant oligotrophic status during the lake history. The study constrains the climate and anthropogenic forcings and watershed versus lake sources in organic matter accumulation processes and helps to design conservation and management policies in mountain, oligotrophic lakes.
Rapidly uplifting coastlines are frequently associated with convergent tectonic boundaries, like subduction zones, which are repeatedly breached by giant megathrust earthquakes. The coastal relief along tectonically active realms is shaped by the effect of sea-level variations and heterogeneous patterns of permanent tectonic deformation, which are accumulated through several cycles of megathrust earthquakes. However, the correlation between earthquake deformation patterns and the sustained long-term segmentation of forearcs, particularly in Chile, remains poorly understood. Furthermore, the methods used to estimate permanent deformation from geomorphic markers, like marine terraces, have remained qualitative and are based on unrepeatable methods. This contrasts with the increasing resolution of digital elevation models, such as Light Detection and Ranging (LiDAR) and high-resolution bathymetric surveys.
Throughout this thesis I study permanent deformation in a holistic manner: from the methods to assess deformation rates, to the processes involved in its accumulation. My research focuses particularly on two aspects: Developing methodologies to assess permanent deformation using marine terraces, and comparing permanent deformation with seismic cycle deformation patterns under different spatial scales along the M8.8 Maule earthquake (2010) rupture zone. Two methods are developed to determine deformation rates from wave-built and wave-cut terraces respectively. I selected an archetypal example of a wave-built terrace at Santa Maria Island studying its stratigraphy and recognizing sequences of reoccupation events tied with eleven radiocarbon sample ages (14C ages). I developed a method to link patterns of reoccupation with sea-level proxies by iterating relative sea level curves for a range of uplift rates. I find the best fit between relative sea-level and the stratigraphic patterns for an uplift rate of 1.5 +- 0.3 m/ka.
A Graphical User Interface named TerraceM® was developed in Matlab®. This novel software tool determines shoreline angles in wave-cut terraces under different geomorphic scenarios. To validate the methods, I select test sites in areas of available high-resolution LiDAR topography along the Maule earthquake rupture zone and in California, USA. The software allows determining the 3D location of the shoreline angle, which is a proxy for the estimation of permanent deformation rates. The method is based on linear interpolations to define the paleo platform and cliff on swath profiles. The shoreline angle is then located by intersecting these interpolations. The
accuracy and precision of TerraceM® was tested by comparing its results with previous assessments, and through an experiment with students in a computer lab setting at the University
of Potsdam.
I combined the methods developed to analyze wave-built and wave-cut terraces to assess regional patterns of permanent deformation along the (2010) Maule earthquake rupture. Wave-built terraces are tied using 12 Infra Red Stimulated luminescence ages (IRSL ages) and shoreline angles in wave-cut terraces are estimated from 170 aligned swath profiles. The comparison of coseismic slip, interseismic coupling, and permanent deformation, leads to three areas of high permanent uplift, terrace warping, and sharp fault offsets. These three areas correlate with regions of high slip and low coupling, as well as with the spatial limit of at least eight historical megathrust ruptures (M8-9.5). I propose that the zones of upwarping at Arauco and Topocalma reflect changes in frictional properties of the megathrust, which result in discrete boundaries for the propagation of mega earthquakes.
To explore the application of geomorphic markers and quantitative morphology in offshore areas I performed a local study of patterns of permanent deformation inferred from hitherto unrecognized drowned shorelines at the Arauco Bay, at the southern part of the (2010) Maule earthquake rupture zone. A multidisciplinary approach, including morphometry, sedimentology, paleontology, 3D morphoscopy, and a landscape Evolution Model is used to recognize, map, and assess local rates and patterns of permanent deformation in submarine environments. Permanent deformation patterns are then reproduced using elastic models to assess deformation rates of an active submarine splay fault defined as Santa Maria Fault System. The best fit suggests a reverse structure with a slip rate of 3.7 m/ka for the last 30 ka. The register of land level changes during the earthquake cycle at Santa Maria Island suggest that most of the deformation may be accrued through splay fault reactivation during mega earthquakes, like the (2010) Maule event. Considering a recurrence time of 150 to 200 years, as determined from historical and geological observations, slip between 0.3 and 0.7 m per event would be required to account for the 3.7 m/ka millennial slip rate. However, if the SMFS slips only every ~1000 years, representing a few megathrust earthquakes, then a slip of ~3.5 m per event would be required to account for the long- term rate. Such event would be equivalent to a magnitude ~6.7 earthquake capable to generate a local tsunami.
The results of this thesis provide novel and fundamental information regarding the amount of permanent deformation accrued in the crust, and the mechanisms responsible for this accumulation at millennial time-scales along the M8.8 Maule earthquake (2010) rupture zone. Furthermore, the results of this thesis highlight the application of quantitative geomorphology and the use of repeatable methods to determine permanent deformation, improve the accuracy of marine terrace assessments, and estimates of vertical deformation rates in tectonically active coastal areas. This is vital information for adequate coastal-hazard assessments and to anticipate realistic earthquake and tsunami scenarios.
High-resolution topographic data greatly facilitate the remote identification of geomorphic features, furnishing valuable information concerning surface processes and characterization of reference markers for quantifying tectonic deformation. Marine terraces have been used as long baseline geodetic markers of relative past sea-level positions, reflecting the interplay between vertical crustal movements and sea-level oscillations. Uplift rates may be determined from the terrace age and the elevation of its shoreline angle, a geomorphic feature that can be correlated with past sea-levels positions. A precise definition of the shoreline angle in time and space is essential to obtain reliable uplift rates with coherent spatial correlation. To improve our ability to rapidly assess and map shoreline angles at regional and local scales, we have developed TerraceM, a MATLAB (R) graphical user interface that allows the shoreline angle and its associated error to be estimated using high-resolution topography. TerraceM uses topographic swath profiles oriented orthogonally to the terrace riser. Four functions are included to analyze the swath profiles and extract the shoreline angle, from both staircase sequences of multiple terraces and rough coasts characterized by eroded remnants of emerged terrace surfaces. The former are measured by outlining the paleocliffs and paieo-platforms and finding their intersection by extrapolating linear regressions, whereas the latter are assessed by automatically detecting peaks of sea-stack tops and back-projecting them to the modern sea cliff. In the absence of rigorous absolute age determinations of marine terraces, their geomorphic age may be estimated using previously published diffusion models. Postprocessing functions are included to obtain first-order statistics of shoreline-angle elevations and their spatial distribution. TerraceM has the ability to process series of profiles from several sites in an efficient and structured workflow. Results may be exported in Google Earth and ESRI shapefile formats. The precision and accuracy of the method have been estimated from a case study at Santa Cruz, California, by comparing TerraceM results with published field measurements. The repeatability was evaluated using multiple measurements made by inexperienced users. TerraceM will improve the efficiency and precision of estimating shoreline-angle elevations in wave-cut terraces in both marine and lacustrine environments.
Mineral topsoils possess large organic carbon (OC) contents but there is only limited knowledge on the mechanisms controlling the preservation of organic matter (OM) against microbial decay. Samples were taken from the uppermost mineral topsoil horizon (0 to 5 cm) of seven sites under mature deciduous forest showing OC contents between 69 and 164 g kg(-1) and a wide range in mineral characteristics. At first, organic particles and the water-extractable OM were removed from the soil samples. Thereafter, Na-pyrophosphate extractable organic matter (OM(PY)), assumed to be indicative for OM bound via cation mediated interactions, and the OM remaining in the extraction residue (OM(ER)), supposed to be indicative for OM occluded in mechanically highly stable micro-aggregates, were sequentially separated and quantified. The composition of OM(PY) and OM(ER) was analyzed by FTIR and their stability by C-14 measurements. The OC remaining in the extraction residues accounted for 38 to 59% of the bulk soil OC (SOC) suggesting a much larger relevance of OM(ER) for the OM dynamic in the analyzed soils as compared with OM(PY) that accounted for 1.6 to 7.5% of the SOC. The FUR analyses revealed a lower relative proportion of C=O groups in OM(ER) compared to OM(PY) indicating differences in the degree of microbial processing between these fractions. Correlation analyses suggest an increase in the stability of OM(PY) with the soil pH and contents of Na-pyrophosphate soluble Fe, Al, and Mg and an increase in the stability of OM(ER) with the soil pH and the contents of clay and oxalate-soluble Fe and Al. Despite the detected influence of soil mineral characteristics on the turnover of OM(PY) and OM(ER), the Delta C-14 signatures indicated mean residence times less than 100 years. The presence of less stabilized OM in these fractions can be derived from methodological uncertainties and/or the fast cycling compartment of mineral-associated OM. (C) 2015 Elsevier B.V. All rights reserved.