Refine
Year of publication
- 2018 (314) (remove)
Document Type
- Article (237)
- Doctoral Thesis (25)
- Postprint (20)
- Other (17)
- Review (10)
- Habilitation Thesis (2)
- Conference Proceeding (1)
- Master's Thesis (1)
- Report (1)
Keywords
- climate change (8)
- Himalaya (7)
- permafrost (7)
- Germany (5)
- landslides (4)
- Climate change (3)
- Kalahari (3)
- Naturgefahren (3)
- Palynology (3)
- Permafrost (3)
Institute
- Institut für Geowissenschaften (314) (remove)
The behaviour of individuals, businesses, and government entities before, during, and immediately after a disaster can dramatically affect the impact and recovery time. However, existing risk-assessment methods rarely include this critical factor. In this Perspective, we show why this is a concern, and demonstrate that although initial efforts have inevitably represented human behaviour in limited terms, innovations in flood-risk assessment that integrate societal behaviour and behavioural adaptation dynamics into such quantifications may lead to more accurate characterization of risks and improved assessment of the effectiveness of risk-management strategies and investments. Such multidisciplinary approaches can inform flood-risk management policy development.
The climate is a complex dynamical system involving interactions and feedbacks among different processes at multiple temporal and spatial scales. Although numerous studies have attempted to understand the climate system, nonetheless, the studies investigating the multiscale characteristics of the climate are scarce. Further, the present set of techniques are limited in their ability to unravel the multi-scale variability of the climate system. It is completely plausible that extreme events and abrupt transitions, which are of great interest to climate community, are resultant of interactions among processes operating at multi-scale. For instance, storms, weather patterns, seasonal irregularities such as El Niño, floods and droughts, and decades-long climate variations can be better understood and even predicted by quantifying their multi-scale dynamics. This makes a strong argument to unravel the interaction and patterns of climatic processes at different scales. With this background, the thesis aims at developing measures to understand and quantify multi-scale interactions within the climate system.
In the first part of the thesis, I proposed two new methods, viz, multi-scale event synchronization (MSES) and wavelet multi-scale correlation (WMC) to capture the scale-specific features present in the climatic processes. The proposed methods were tested on various synthetic and real-world time series in order to check their applicability and replicability. The results indicate that both methods (WMC and MSES) are able to capture scale-specific associations that exist between processes at different time scales in a more detailed manner as compared to the traditional single scale counterparts.
In the second part of the thesis, the proposed multi-scale similarity measures were used in constructing climate networks to investigate the evolution of spatial connections within climatic processes at multiple timescales. The proposed methods WMC and MSES, together with complex network were applied to two different datasets.
In the first application, climate networks based on WMC were constructed for the univariate global sea surface temperature (SST) data to identify and visualize the SSTs patterns that develop very similarly over time and distinguish them from those that have long-range teleconnections to other ocean regions. Further investigations of climate networks on different timescales revealed (i) various high variability and co-variability regions, and (ii) short and long-range teleconnection regions with varying spatial distance. The outcomes of the study not only re-confirmed the existing knowledge on the link between SST patterns like El Niño Southern Oscillation and the Pacific Decadal Oscillation, but also suggested new insights into the characteristics and origins of long-range teleconnections.
In the second application, I used the developed non-linear MSES similarity measure to quantify the multivariate teleconnections between extreme Indian precipitation and climatic patterns with the highest relevance for Indian sub-continent. The results confirmed significant non-linear influences that were not well captured by the traditional methods. Further, there was a substantial variation in the strength and nature of teleconnection across India, and across time scales.
Overall, the results from investigations conducted in the thesis strongly highlight the need for considering the multi-scale aspects in climatic processes, and the proposed methods provide robust framework for quantifying the multi-scale characteristics.
Quantifying the roles of single stations within homogeneous regions using complex network analysis
(2018)
Regionalization and pooling stations to form homogeneous regions or communities are essential for reliable parameter transfer, prediction in ungauged basins, and estimation of missing information. Over the years, several clustering methods have been proposed for regional analysis. Most of these methods are able to quantify the study region in terms of homogeneity but fail to provide microscopic information about the interaction between communities, as well as about each station within the communities. We propose a complex network-based approach to extract this valuable information and demonstrate the potential of our approach using a rainfall network constructed from the Indian gridded daily precipitation data. The communities were identified using the network-theoretical community detection algorithm for maximizing the modularity. Further, the grid points (nodes) were classified into universal roles according to their pattern of within- and between-community connections. The method thus yields zoomed-in details of individual rainfall grids within each community.
Python is at the forefront of scientific computation for seismologists and therefore should be introduced to students interested in becoming seismologists. On its own, Python is open source and well designed with extensive libraries. However, Python code can also be executed, visualized, and communicated to others with "Jupyter Notebooks". Thus, Jupyter Notebooks are ideal for teaching students Python and scientific computation. In this article, we designed an openly available Python library and collection of Jupyter Notebooks based on defined scientific computation learning goals for seismology students. The Notebooks cover topics from an introduction to Python to organizing data, earthquake catalog statistics, linear regression, and making maps. Our Python library and collection of Jupyter Notebooks are meant to be used as course materials for an upper-division data analysis course in an Earth Science Department, and the materials were tested in a Probabilistic Seismic Hazard course. However, seismologists or anyone else who is interested in Python for data analysis and map making can use these materials.
Mechanical and/or chemical removal of material from the subsurface may generate large subsurface cavities, the destabilisation of which can lead to ground collapse and the formation of sinkholes. Numerical simulation of the interaction of cavity growth, host material deformation and overburden collapse is desirable to better understand the sinkhole hazard but is a challenging task due to the involved high strains and material discontinuities. Here, we present 2-D distinct element method numerical simulations of cavity growth and sinkhole development. Firstly, we simulate cavity formation by quasi-static, stepwise removal of material in a single growing zone of an arbitrary geometry and depth. We benchmark this approach against analytical and boundary element method models of a deep void space in a linear elastic material. Secondly, we explore the effects of properties of different uniform materials on cavity stability and sinkhole development. We perform simulated biaxial tests to calibrate macroscopic geotechnical parameters of three model materials representative of those in which sinkholes develop at the Dead Sea shoreline: mud, alluvium and salt. We show that weak materials do not support large cavities, leading to gradual sagging or suffusion-style subsidence. Strong materials support quasi-stable to stable cavities, the overburdens of which may fail suddenly in a caprock or bedrock collapse style. Thirdly, we examine the consequences of layered arrangements of weak and strong materials. We find that these are more susceptible to sinkhole collapse than uniform materials not only due to a lower integrated strength of the overburden but also due to an inhibition of stabilising stress arching. Finally, we compare our model sinkhole geometries to observations at the Ghor Al-Haditha sinkhole site in Jordan. Sinkhole depth / diameter ratios of 0.15 in mud, 0.37 in alluvium and 0.33 in salt are reproduced successfully in the calibrated model materials. The model results suggest that the observed distribution of sinkhole depth / diameter values in each material type may partly reflect sinkhole growth trends.
Mechanical and/or chemical removal of material from the subsurface may generate large subsurface cavities, the destabilisation of which can lead to ground collapse and the formation of sinkholes. Numerical simulation of the interaction of cavity growth, host material deformation and overburden collapse is desirable to better understand the sinkhole hazard but is a challenging task due to the involved high strains and material discontinuities. Here, we present 2-D distinct element method numerical simulations of cavity growth and sinkhole development. Firstly, we simulate cavity formation by quasi-static, stepwise removal of material in a single growing zone of an arbitrary geometry and depth. We benchmark this approach against analytical and boundary element method models of a deep void space in a linear elastic material. Secondly, we explore the effects of properties of different uniform materials on cavity stability and sinkhole development. We perform simulated biaxial tests to calibrate macroscopic geotechnical parameters of three model materials representative of those in which sinkholes develop at the Dead Sea shoreline: mud, alluvium and salt. We show that weak materials do not support large cavities, leading to gradual sagging or suffusion-style subsidence. Strong materials support quasi-stable to stable cavities, the overburdens of which may fail suddenly in a caprock or bedrock collapse style. Thirdly, we examine the consequences of layered arrangements of weak and strong materials. We find that these are more susceptible to sinkhole collapse than uniform materials not only due to a lower integrated strength of the overburden but also due to an inhibition of stabilising stress arching. Finally, we compare our model sinkhole geometries to observations at the Ghor Al-Haditha sinkhole site in Jordan. Sinkhole depth ∕ diameter ratios of 0.15 in mud, 0.37 in alluvium and 0.33 in salt are reproduced successfully in the calibrated model materials. The model results suggest that the observed distribution of sinkhole depth ∕ diameter values in each material type may partly reflect sinkhole growth trends.
The concept of hydrologic connectivity summarizes all flow processes that link separate regions of a landscape. As such, it is a central theme in the field of catchment hydrology, with influence on neighboring disciplines such as ecology and geomorphology. It is widely acknowledged to be an important key in understanding the response behavior of a catchment and has at the same time inspired research on internal processes over a broad range of scales. From this process-hydrological point of view, hydrological connectivity is the conceptual framework to link local observations across space and scales.
This is the context in which the four studies this thesis comprises of were conducted. The focus was on structures and their spatial organization as important control on preferential subsurface flow. Each experiment covered a part of the conceptualized flow path from hillslopes to the stream: soil profile, hillslope, riparian zone, and stream.
For each study site, the most characteristic structures of the investigated domain and scale, such as slope deposits and peat layers were identified based on preliminary or previous investigations or literature reviews. Additionally, further structural data was collected and topographical analyses were carried out. Flow processes were observed either based on response observations (soil moisture changes or discharge patterns) or direct measurement (advective heat transport). Based on these data, the flow-relevance of the characteristic structures was evaluated, especially with regard to hillslope to stream connectivity.
Results of the four studies revealed a clear relationship between characteristic spatial structures and the hydrological behavior of the catchment. Especially the spatial distribution of structures throughout the study domain and their interconnectedness were crucial for the establishment of preferential flow paths and their relevance for large-scale processes. Plot and hillslope-scale irrigation experiments showed that the macropores of a heterogeneous, skeletal soil enabled preferential flow paths at the scale of centimeters through the otherwise unsaturated soil. These flow paths connected throughout the soil column and across the hillslope and facilitated substantial amounts of vertical and lateral flow through periglacial slope deposits.
In the riparian zone of the same headwater catchment, the connectivity between hillslopes and stream was controlled by topography and the dualism between characteristic subsurface structures and the geomorphological heterogeneity of the stream channel. At the small scale (1 m to 10 m) highest gains always occurred at steps along the longitudinal streambed profile, which also controlled discharge patterns at the large scale (100 m) during base flow conditions (number of steps per section). During medium and high flow conditions, however, the impact of topography and parafluvial flow through riparian zone structures prevailed and dominated the large-scale response patterns.
In the streambed of a lowland river, low permeability peat layers affected the connectivity between surface water and groundwater, but also between surface water and the hyporheic zone. The crucial factor was not the permeability of the streambed itself, but rather the spatial arrangement of flow-impeding peat layers, causing increased vertical flow through narrow “windows” in contrast to predominantly lateral flow in extended areas of high hydraulic conductivity sediments.
These results show that the spatial organization of structures was an important control for hydrological processes at all scales and study areas. In a final step, the observations from different scales and catchment elements were put in relation and compared. The main focus was on the theoretical analysis of the scale hierarchies of structures and processes and the direction of causal dependencies in this context. Based on the resulting hierarchical structure, a conceptual framework was developed which is capable of representing the system’s complexity while allowing for adequate simplifications.
The resulting concept of the parabolic scale series is based on the insight that flow processes in the terrestrial part of the catchment (soil and hillslopes) converge. This means that small-scale processes assemble and form large-scale processes and responses. Processes in the riparian zone and the streambed, however, are not well represented by the idea of convergence. Here, the large-scale catchment signal arrives and is modified by structures in the riparian zone, stream morphology, and the small-scale interactions between surface water and groundwater. Flow paths diverge and processes can better be represented by proceeding from large scales to smaller ones. The catchment-scale representation of processes and structures is thus the conceptual link between terrestrial hillslope processes and processes in the riparian corridor.
Here we present high-resolution delta O-18 records obtained from speleothems collected in the eastern Bolivian Andes. The stable isotope records are related to the regional- to large-scale atmospheric circulation over South America and allow interpreting changes in delta O-18 during the last 1400 yr as a function of changes in precipitation regimes over the southern tropical Andes. Two distinct phases with more negative delta O-18 values, interpreted as periods of increased convective activity over the eastern Andean Cordillera in Bolivia are observed concomitantly with periods of global climate anomalies during the last millennium, such as the Medieval Climate Anomaly (MCA) and the Little Ice Age (LIA) respectively. Changes in the Bolivian delta O-18 record during the LIA are apparently related to a southward displacement of the Intertropical Convergence Zone (ITCZ), which acts as a main moisture driver to intensify convection over the tropical continent. During the MCA, however, the increased convective activity observed in the Bolivian record is likely the result of a different mechanism, which implies moisture sourced mainly from the southern tropical Atlantic. This interpretation is consistent with paleoclimate records further to the north in the tropical Andes that show progressively drier conditions during this time period, indicating a more northerly position of the ITCZ. The transition period between the MCA and the LIA shows a slight tendency toward increased delta O-18 values, indicating weakened convective activity. Our results also reveal a non-stationary anti-phased behavior between the delta O-18 reconstructions from Bolivia and northeastern Brazil that confirms a continental-scale east-west teleconnection across South America during the LIA.
TDS-1 GNSS Reflectometry
(2018)
This study presents the development and a systematic evaluation study of GNSS reflectometry wind speeds. After establishing a wind speed retrieval algorithm, UK TechDemoSat-1 (TDS-1) derived winds, from May 2015 to July 2017, are compared to the Advanced Scatterometer (ASCAT). ERA-Interim wind fields of the European Centre for Medium-range Weather Forecasts (ECMWF) and in situ observation from Tropical Atmosphere Ocean buoy array in the Pacific are taken as reference. One-year averaged TDS-1 global winds demonstrate small differences with ECMWF in a majority of areas as well as discuss under- and overestimations. The pioneering TDS-1 winds demonstrate a root-mean-squared error (RMSE) and bias of 2.77 and -0.33 m/s, which are comparable to the RMSE and bias derived by ASCAT winds, as large as 2.31 and 0.25 m/s, respectively. Using buoys measurements as reference, RMSE and bias of 2.23 and -0.03 m/s for TDS-1 as well as 1.40 and -0.68 m/s for ASCAT are obtained. Utilizing rain microwave-infrared estimates of the Tropical Rainfall Measuring Mission, rain-affected observation of both ASCAT and TDS-1 are collected and evaluated. Although ASCAT winds show a significant performance degradation resulting in an RMSE and bias of 3.16 and 1.03 m/s, respectively, during rain condition, TDS-1 shows a more reliable performance with an RMSE and bias of 2.94 and -0.21 m/s, respectively, which indicates the promising capability of GNSS forward scattering for wind retrievals during rain. A decrease in TDS-1-derived bistatic radar cross sections during rain events, at weak winds, is also demonstrated.
For the first time, a rain signature in Global Navigation Satellite System Reflectometry (GNSS-R) observations is demonstrated. Based on the argument that the forward quasi-specular scattering relies upon surface gravity waves with lengths larger than several wavelengths of the reflected signal, a commonly made conclusion is that the scatterometric GNSS-R measurements are not sensitive to the surface small-scale roughness generated by raindrops impinging on the ocean surface. On the contrary, this study presents an evidence that the bistatic radar cross section sigma(0) derived from TechDemoSat-1 data is reduced due to rain at weak winds, lower than approximate to 6 m/s. The decrease is as large as approximate to 0.7 dB at the wind speed of 3 m/s due to a precipitation of 0-2 mm/hr. The simulations based on the recently published scattering theory provide a plausible explanation for this phenomenon which potentially enables the GNSS-R technique to detect precipitation over oceans at low winds.
Woodlands and steppes
(2018)
Based on fossil organism remains including plant macrofossils, charcoal, pollen, and invertebrates preserved in syngenetic deposits of the Batagay permafrost sequence in the Siberian Yana Highlands, we reconstructed the environmental history during marine isotope stages (MIS) 6 to 2. Two fossil assemblages, exceptionally rich in plant remains, allowed for a detailed description of the palaeo-vegetation during two climate extremes of the Late Pleistocene, the onset of the last glacial maximum (LGM) and the last interglacial. In addition, altogether 41 assemblages were used to outline the vegetation history since the penultimate cold stage of MIS 6. Accordingly, meadow steppes analogue to modern communities of the phytosociological order Festucetalia lenensis formed the primary vegetation during the Saalian and Weichselian cold stages. Cold-resistant tundra-steppe communities (Carici rupestris-Kobresietea bellardii) as they occur above the treeline today were, in contrast to more northern locations, mostly lacking. During the last interglacial, open coniferous woodland similar to modern larch taiga was the primary vegetation at the site. Abundant charcoal indicates wildfire events during the last interglacial. Zoogenic disturbances of the local vegetation were indicated by the presence of ruderal plants, especially by abundant Urtica dioica, suggesting that the area was an interglacial refugium for large herbivores. Meadow steppes, which formed the primary vegetation during cold stages and provided potentially suitable pastures for herbivores, were a significant constituent of the plant cover in the Yana Highlands also under the full warm stage conditions of the last interglacial. Consequently, meadow steppes occurred in the Yana Highlands during the entire investigated timespan from MIS 6 to MIS 2 documenting a remarkable environmental stability. Thus, the proportion of meadow steppe vegetation merely shifted in response to the respectively prevailing climatic conditions. Their persistence indicates low precipitation and a relatively warm growing season throughout and beyond the late Pleistocene. The studied fossil record also proves that modern steppe occurrences in the Yana Highlands did not establish as late as in the Holocene but instead are relicts of a formerly continuous steppe belt extending from Central Siberia to Northeast Yakutia during the Pleistocene. The persistence of plants and invertebrates characteristic of meadow steppe vegetation in interior Yakutia throughout the late Quaternary indicates climatic continuity and documents the suitability of this region as a refugium also for other organisms of the Pleistocene mammoth steppe including the iconic large herbivores. (C)2018 Elsevier Ltd. All rights reserved.
The SISAL database
(2018)
Stable isotope records from speleothems provide information on past climate changes, most particularly information that can be used to reconstruct past changes in precipitation and atmospheric circulation. These records are increasingly being used to provide "out-of-sample" evaluations of isotope-enabled climate models. SISAL (Speleothem Isotope Synthesis and Analysis) is an international working group of the Past Global Changes (PAGES) project. The working group aims to provide a comprehensive compilation of speleothem isotope records for climate reconstruction and model evaluation. The SISAL database contains data for individual speleothems, grouped by cave system. Stable isotopes of oxygen and carbon (delta O-18, delta C-13) measurements are referenced by distance from the top or bottom of the speleothem. Additional tables provide information on dating, including information on the dates used to construct the original age model and sufficient information to assess the quality of each data set and to erect a standardized chronology across different speleothems. The metadata table provides location information, information on the full range of measurements carried out on each speleothem and information on the cave system that is relevant to the interpretation of the records, as well as citations for both publications and archived data.
Tomato (Solanum lycopersicum L.) being an important vegetable is cultivated and used throughout the world. It not only contributes in fulfilling the basic nutritional requirements of the human body but also has many health benefits due to its rich biochemical composition. However, its production at large scale is hampered by many limiting factors such as biotic and abiotic stresses. Among the different abiotic stresses, drought poses drastic impact on tomato yield. Drought stress is genetically regulated by many transcription factors that not only regulate the stress responsive mechanism but also facilitate the growth and development of tomato plants. NAC is an important stress related transcription factor genes family, and the ATAF1 gene, a member of this family, is involved in ABA signaling and stress response. In this study, tomato variety Rio Drande was transformed with drought resistant ATAF1 gene via Agrobacterium mediated gene transformation method. The ATAF1 gene was first cloned in the pK7WFG2 vector having kanamycin selectable marker and then it was introduced in the Agrobacterium tumefaciens strain GV3101 through heat shock method. The tomato cotyledon and hypocotyl ex-plants of variety "Rio Ggrande" were cultured on callus induction medium (MS + 2.5 mg/L IAA + 2 mg/L BAP). The calli were then infected with Agrobacterium tumefaciens strain GV3101 containing ATAF1 gene and selection was carried out on the kanamycin selectable medium (MS + 100 mg/L Kan), and were regenerated on MS medium with 1 mg/L IAA + 1 mg/L BAP. Out of 216 putative transformed calli, 13 calli were able to regenerate on the selection medium. Of the 13 calli, three transgenic tomato plantlets were recovered, and these were confirmed through PCR analysis for the presence of 432 bp fragment of ATAF1 gene. The transformation protocol reported here can be used to generate drought resistant tomato plants in future.
Terrestrial environmental systems are characterised by numerous feedback links between their different compartments. However, scientific research is organized into disciplines that focus on processes within the respective compartments rather than on interdisciplinary links. Major feedback mechanisms between compartments might therefore have been systematically overlooked so far. Without identifying these gaps, initiatives on future comprehensive environmental monitoring schemes and experimental platforms might fail. We performed a comprehensive overview of feedbacks between compartments currently represented in environmental sciences and explores to what degree missing links have already been acknowledged in the literature. We focused on process models as they can be regarded as repositories of scientific knowledge that compile findings of numerous single studies. In total, 118 simulation models from 23 model types were analysed. Missing processes linking different environmental compartments were identified based on a meta-review of 346 published reviews, model inter-comparison studies, and model descriptions. Eight disciplines of environmental sciences were considered and 396 linking processes were identified and ascribed to the physical, chemical or biological domain. There were significant differences between model types and scientific disciplines regarding implemented interdisciplinary links. The most wide-spread interdisciplinary links were between physical processes in meteorology, hydrology and soil science that drive or set the boundary conditions for other processes (e.g., ecological processes). In contrast, most chemical and biological processes were restricted to links within the same compartment. Integration of multiple environmental compartments and interdisciplinary knowledge was scarce in most model types. There was a strong bias of suggested future research foci and model extensions towards reinforcing existing interdisciplinary knowledge rather than to open up new interdisciplinary pathways. No clear pattern across disciplines exists with respect to suggested future research efforts. There is no evidence that environmental research would clearly converge towards more integrated approaches or towards an overarching environmental systems theory. (c) 2017 Elsevier B.V. All rights reserved.
The Aral Sea desiccation and related changes in hydroclimatic conditions on a regional level is a hot topic for past decades. The key problem of scientific research projects devoted to an investigation of modern Aral Sea basin hydrological regime is its discontinuous nature - the only limited amount of papers takes into account the complex runoff formation system entirely. Addressing this challenge we have developed a continuous prediction system for assessing freshwater inflow into the Small Aral Sea based on coupling stack of hydrological and data-driven models. Results show a good prediction skill and approve the possibility to develop a valuable water assessment tool which utilizes the power of classical physically based and modern machine learning models both for territories with complex water management system and strong water-related data scarcity. The source code and data of the proposed system is available on a Github page (https://github.com/SMASHIproject/IWRM2018).
Runoff predictions in ungauged arctic basins using conceptual models forced by reanalysis data
(2018)
Due to global warming, the problem of assessing water resources and their vulnerability to climate drivers in the Arctic region has become a focus in the recent years. This study is aimed at investigating three lumped hydrological models to predict daily runoff of large-scale Arctic basins in the case of substantial data scarcity. All models were driven only by meteorological forcing reanalysis dataset without any additional information about landscape, soil, or vegetation cover properties of the studied basins. Model parameter regionalization based on transferring the whole parameter set showed good efficiency for predictions in ungauged basins. We run a blind test of the proposed methodology for ensemble runoff predictions on five sub-basins, for which only monthly observations were available, and obtained promising results for current water resources assessment for a broad domain of ungauged basins in the Russian Arctic.
Pre-existing weakness zones in the lithosphere such as transform faults/fracture zones and extinct mid-oceanic ridges have been suggested to facilitate subduction initiation in an intra-oceanic environment. Here, we propose that the additional forcing coming from the mantle suction flow is required to trigger the conversion of a fracture zone/transform fault into a converging plate boundary. This suction flow can be induced either from the slab remnants of former converging plate boundaries or/and from slabs of neighbouring active subduction zones. Using 2-D coupled thermo-mechanical models, we show that a sufficiently strong mantle flow is able to convert a fracture zone/transform fault into a subduction zone. However, this process is feasible only if the fracture zone/transform fault is very close to the mid-oceanic ridge. Our numerical model results indicate that time of subduction initiation depends on the velocity, domain size and location of mantle suction flow and age of the oceanic plate.
Linking deep seismic profiles with regional-scale gravity inversion is a powerful tool to deduce the architecture of rifted margins and their structural evolution. Here we map upper and lower crustal thicknesses of the northern South China Sea (SCS) margin in order to investigate the occurrence of depth-dependent crustal extension from the proximal to the distal margin. By comparing upper and lower crustal stretching factors, we find that the northern margin of the SCS is segmented in three parts: (1) sedimentary basins where upper crust is stretched more than lower crust, (2) distal margin where lower crust is stretched more than upper crust, (3) mostly proximal margin regions where the two layers have similar stretching factors. Our results suggest that sedimentary basins and distal margin prominently feature depth-dependent extension, however accommodated by different processes. While differential thinning within sedimentary basins appears to be governed by lateral pressure variations inducing lower crustal flow, we suggest the distal margin to be affected by a combination of mantle flow-induced lower crustal shearing and sequential fault activity during crustal hyper-extension.
The Central Pontides of N Turkey represents a mobile orogenic belt of the southern Eurasian margin that experienced several phases of exhumation associated with the consumption of different branches of the Neo-Tethys Ocean and the amalgamation of continental domains. Our new low-temperature thermochronology data help to constrain the timing of these episodes, providing new insights into associated geodynamic processes. In particular, our data suggest that exhumation occurred at (1) similar to 110 to 90Ma, most likely during tectonic accretion and exhumation of metamorphic rocks from the subduction zone; (2) from similar to 60 to 40Ma, during the collision of the Kirehir and Anatolide-Tauride microcontinental domains with the Eurasian margin; (3) from similar to 0 to 25Ma, either during the early stages of the Arabia-Eurasia collision (soft collision) when the Arabian passive margin reached the trench, implying 70 to 530km of subduction of the Arabian passive margin, or during a phase of trench advance predating hard collision at similar to 20Ma; and (4) similar to 11Ma to the present, during transpression associated with the westward motion of Anatolia. Our findings document the punctuated nature of fault-related exhumation, with episodes of fast cooling followed by periods of slow cooling or subsidence, the role of inverted normal faults in controlling the Paleogene exhumation pattern, and of the North Anatolian Fault in dictating the most recent pattern of exhumation.
The lack of process-based classification procedures may lead to unrealistic hyetograph design due to complex oscillation of rainfall depths when assimilated at high temporal resolutions. Four consecutive years of sub-hourly rainfall data were assimilated in three study areas (Guaraira, GEB, Sao Joao do Cariri, CEB, and Aiuaba, AEB) under distinct climates (very hot semi-arid and tropical wet). This study aimed to define rainfall events (for Minimum Inter-event Time, MIT, and Minimum Rainfall Depth, MRD, equal to 30 min and 1.016 mm, respectively), classify their hyetograph types (rectangular, R, unimodal with left-skewed, UL, right-skewed, UR, and centred peaks, UC, bimodal, B, and shapeless, SL), and compare their key rainfall properties (frequency, duration, depth, rate and peak). A rain pulse aggregation process allowed for reshaping SL-events for six different time spans varying from 2 to 30 min. The results revealed that the coastal area held predominantly R-events (64% events and 49% rainfall depth), in western semi-arid prevailed UL-events (57% events and 63% rainfall depth), whereas in eastern semi-arid mostly were R-events (61% events and 30% rainfall depth) similar to coastal area. It is concluded that each cloud formation type had important effects on hyetograph properties, differentiating them even within the same climate.
Cosmic-ray neutron sensing (CRNS) is a promising proximal soil sensing technique to estimate soil moisture at intermediate scale and high temporal resolution. However, the signal shows complex and non-unique response to all hydrogen pools near the land surface, providing some challenges for soil moisture estimation in practical applications. Aims of the study were 1) to assess the uncertainty of CRNS as a stand-alone approach to estimate volumetric soil moisture in cropped field 2) to identify the causes of this uncertainty 3) and possible improvements. Two experimental sites in Germany were equipped with a CRNS probe and point-scale soil moisture network. Additional monitoring activities were conducted during the crop growing season to characterize the soil-plant systems. This data is used to identify and quantify the different sources of uncertainty (factors). An uncertainty analysis, based on Monte Carlo approach, is applied to propagate these uncertainties to CRNS soil moisture estimations. In addition, a sensitivity analysis based on the Sobol’ method is performed to identify the most important factors explaining this uncertainty. Results show that CRNS soil moisture compares well to the soil moisture network when these point-scale values are weighted to account for the spatial sensitivity of the signal and other sources of hydrogen (lattice water and organic carbon) are added to the water content. However, the performance decreases when CRNS is considered as a stand-alone method to retrieve the actual (non-weighted) volumetric soil moisture. The support volume (penetration depth and radius) shows also a considerable uncertainty, especially in relatively dry soil moisture conditions. Four of the seven factors analyzed (the vertical soil moisture profile, bulk density, incoming neutron correction and the calibrated parameter N0) were found to play an important role. Among the possible improvements identified, a simple correction factor based on vertical point-scale soil moisture profiles shows to be a promising approach to account for the sensitivity of the CRNS signal to the upper soil layers.
Causes for slow weathering and erosion in the steep, warm, monsoon-subjected Highlands of Sri Lanka
(2018)
In the Highlands of Sri Lanka, erosion and chemical weathering rates are among the lowest for global mountain denudation. In this tropical humid setting, highly weathered deep saprolite profiles have developed from high-grade metamorphic charnockite during spheroidal weathering of the bedrock. The spheroidal weathering produces rounded corestones and spalled rindlets at the rock-saprolite interface. I used detailed textural, mineralogical, chemical, and electron-microscopic (SEM, FIB, TEM) analyses to identify the factors limiting the rate of weathering front advance in the profile, the sequence of weathering reactions, and the underlying mechanisms. The first mineral attacked by weathering was found to be pyroxene initiated by in situ Fe oxidation, followed by in situ biotite oxidation. Bulk dissolution of the primary minerals is best described with a dissolution – re-precipitation process, as no chemical gradients towards the mineral surface and sharp structural boundaries are observed at the nm scale. Only the local oxidation in pyroxene and biotite is better described with an ion by ion process. The first secondary phases are oxides and amorphous precipitates from which secondary minerals (mainly smectite and kaolinite) form. Only for biotite direct solid state transformation to kaolinite is likely. The initial oxidation of pyroxene and biotite takes place in locally restricted areas and is relatively fast: log J = -11 molmin/(m2 s). However, calculated corestone-scale mineral oxidation rates are comparable to corestone-scale mineral dissolution rates: log R = -13 molpx/(m2 s) and log R = -15 molbt/(m2 s). The oxidation reaction results in a volume increase. Volumetric calculations suggest that this observed oxidation leads to the generation of porosity due to the formation of micro-fractures in the minerals and the bedrock allowing for fluid transport and subsequent dissolution of plagioclase. At the scale of the corestone, this fracture reaction is responsible for the larger fractures that lead to spheroidal weathering and to the formation of rindlets. Since these fractures have their origin from the initial oxidational induced volume increase, oxidation is the rate limiting parameter for weathering to take place. The ensuing plagioclase weathering leads to formation of high secondary porosity in the corestone over a distance of only a few cm and eventually to the final disaggregation of bedrock to saprolite. As oxidation is the first weathering reaction, the supply of O2 is a rate-limiting factor for chemical weathering. Hence, the supply of O2 and its consumption at depth connects processes at the weathering front with erosion at the surface in a feedback mechanism. The strength of the feedback depends on the relative weight of advective versus diffusive transport of O2 through the weathering profile. The feedback will be stronger with dominating diffusive transport. The low weathering rate ultimately depends on the transport of O2 through the whole regolith, and on lithological factors such as low bedrock porosity and the amount of Fe-bearing primary minerals. In this regard the low-porosity charnockite with its low content of Fe(II) bearing minerals impedes fast weathering reactions. Fresh weatherable surfaces are a pre-requisite for chemical weathering. However, in the case of the charnockite found in the Sri Lankan Highlands, the only process that generates these surfaces is the fracturing induced by oxidation. Tectonic quiescence in this region and low pre-anthropogenic erosion rate (attributed to a dense vegetation cover) minimize the rejuvenation of the thick and cohesive regolith column, and lowers weathering through the feedback with erosion.
Subduction zone magmas are more oxidised on eruption than those at mid-ocean ridges. This is attributed either to oxidising components, derived from subducted lithosphere (slab) and added to the mantle wedge, or to oxidation processes occurring during magma ascent via differentiation. Here we provide direct evidence for contributions of oxidising slab agents to melts trapped in the sub-arc mantle. Measurements of sulfur (S) valence state in sub-arc mantle peridotites identify sulfate, both as crystalline anhydrite (CaSO4) and dissolved SO42- in spinel-hosted glass (formerly melt) inclusions. Copper-rich sulfide precipitates in the inclusions and increased Fe3+/Sigma Fe in spinel record a S6+-Fe2+ redox coupling during melt percolation through the sub-arc mantle. Sulfate-rich glass inclusions exhibit high U/Th, Pb/Ce, Sr/Nd and delta S-34 (+ 7 to + 11%), indicating the involvement of dehydration products of serpentinised slab rocks in their parental melt sources. These observations provide a link between liberated slab components and oxidised arc magmas.
Three spinel lherzolite xenoliths from Mt. Quincan (Queensland, northeastern Australia) were studied with special attention to their enclosed fluid inclusions. The xenoliths are deformed, have porphyroclastic textures and overall show very similar petrographic features. The only significant difference is manifested in the abundance of fluid inclusions in the samples, mostly in orthopyroxene porphyroclasts. Xenolith JMTQ11 is fluid inclusion-free, whereas xenolith JMTQ20 shows a high abundance of fluid inclusions (fluid inclusion-rich). Xenolith JMTQ45 represents a transitional state between the previous two, as it contains only a small amount of fluid inclusions (fluid inclusion-bearing). Previous studies revealed that these xenoliths and the entrapped fluid inclusions represent a former addition of a MORB-type fluid to the pre-existing lithosphere, resulting from asthenosphere upwelling. There is a progressive enrichment in LREE, Nb, Sr and Ti from the fluid inclusion-free xenolith through the fluid inclusion-bearing one to the fluid inclusion-rich lherzolite. This suggests an increase in the extent of the interaction between the fluid-rich melt and the lherzolite wallrock. In addition, the same interaction is considered to be responsible for the formation of pargasitic amphibole as well. The presence of fluid inclusions indicates fluid migration at mantle depth, and their association with exsolution lamellae in orthopyroxene suggests fluid entrapment following the continental rifting (thermal relaxation) during cooling. A series of analyses, including microthermometry coupled with Raman spectroscopy, FTIR hyperspectral imaging, and Focused Ion Beam-Scanning Electron Microscopy (FIB-SEM) was carried out on the fluid inclusions. Based on the results, the entrapped high-density fluid is composed of 7589 mol% CO2, 918 mol% H2O, 0.11.7 mol% N-2 and <= 0.5 mol% H2S with dissolved trace elements (melt component). Our findings suggest that the metasomatic fluid phase could have been either a fluid/fluid-rich silicate melt released from the deeper asthenosphere, or a coexisting incipient fluid-rich silicate melt. Further cooling, possibly due to thermal relaxation and the upward migration of the fluid phase, caused the investigated lherzolites to reach pargasite stability conditions. We conclude that pargasite, even if only present in very limited modal proportions, can be a common phase at spinel lherzolite stability in the lithospheric upper mantle in continental rift back-arc settings. Studies of fluid inclusions indicate that significant CO2 release from the asthenosphere in a continental rifting environment is resulting from asthenosphere upwelling and its addition to the lithospheric mantle together with fluid-rich melt lherzolite interaction that leaves a CO2-rich fluid behind.
The effects of climate and topography on soil physico-chemical and microbial parameters were studied along an extensive latitudinal climate gradient in the Coastal Cordillera of Chile (26 degrees-38 degrees S). The study sites encompass arid (Pan de Azucar), semiarid (Santa Gracia), mediterranean (La Campana) and humid (Nahuelbuta) climates and vegetation, ranging from arid desert, dominated by biological soil crusts (biocrusts), semiarid shrubland and mediterranean sclerophyllous forest, where biocrusts are present but do have a seasonal pattern to temperate-mixed forest, where biocrusts only occur as an early pioneering development stage after disturbance. All soils originate from granitic parent materials and show very strong differences in pedogenesis intensity and soil depth. Most of the investigated physical, chemical and microbiological soil properties showed distinct trends along the climate gradient. Further, abrupt changes between the arid northernmost study site and the other semi-arid to humid sites can be shown, which indicate non-linearity and thresholds along the climate gradient. Clay and total organic carbon contents (TOC) as well as Ah horizons and solum depths increased from arid to humid climates, whereas bulk density (BD), pH values and base saturation (BS) decreased. These properties demonstrate the accumulation of organic matter, clay formation and element leaching as key-pedogenic processes with increasing humidity. However, the soils in the northern arid climate do not follow this overall latitudinal trend, because texture and BD are largely controlled by aeolian input of dust and sea salts spray followed by the formation of secondary evaporate minerals. Total soil DNA concentrations and TOC increased from arid to humid sites, while areal coverage by biocrusts exhibited an opposite trend. Relative bacterial and archaeal abundances were lower in the arid site, but for the other sites the local variability exceeds the variability along the climate gradient. Differences in soil properties between topographic positions were most pronounced at the study sites with the mediterranean and humid climate, whereas microbial abundances were independent on topography across all study sites. In general, the regional climate is the strongest controlling factor for pedogenesis and microbial parameters in soils developed from the same parent material. Topographic position along individual slopes of limited length augmented this effect only under humid conditions, where water erosion likely relocated particles and elements downward. The change from alkaline to neutral soil pH between the arid and the semi-arid site coincided with qualitative differences in soil formation as well as microbial habitats. This also reflects non-linear relationships of pedogenic and microbial processes in soils depending on climate with a sharp threshold between arid and semi-arid conditions. Therefore, the soils on the transition between arid and semi-arid conditions are especially sensitive and may be well used as indicators of long and medium-term climate changes. Concluding, the unique latitudinal precipitation gradient in the Coastal Cordillera of Chile is predestined to investigate the effects of the main soil forming factor - climate - on pedogenic processes.
Ground‐motion prediction equations (GMPEs) are calibrated to predict the intensity of ground shaking at any given location, based on earthquake magnitude, source‐to‐site distance, local soil amplifications, and other parameters. GMPEs are generally assumed to be independent of time; however, evidence is increasing that large earthquakes modify the shallow soil conditions and those of the fault zone for months or years. These changes may affect the intensity of shaking and result in time‐dependent effects that can potentially be resolved by analyzing between‐event residuals (residuals between observed and predicted ground motion for individual earthquakes averaged over all stations). Here, we analyze a data set of about 65,000 recordings for about 1400 earthquakes in the moment magnitude range 2.5–6.5 that occurred in central Italy from 2008 to 2017 to capture the temporal variability of the ground shaking at high frequency. We first compute between‐event residuals for each earthquake in the Fourier domain with respect to a GMPE developed ad hoc for the analyzed data set. The between‐events show large changes after the occurrence of mainshocks such as the 2009 Mw 6.3 L'Aquila, the 2016 Mw 6.2 Amatrice, and Mw 6.5 Norcia earthquakes. Within the time span of a few months after the mainshocks, the between‐event contribution to the ground shaking varies by a factor 7. In particular, we find a large drop in the between‐events in the aftermath of the L'Aquila earthquake, followed by a slow positive trend that leads to a recovery interrupted by a new drop at the beginning of 2014. We also quantify the frequency‐dependent correlation between the Brune stress drop Δσ and the between‐events. We find that the temporal changes of Δσ resemble those of the between‐event residuals; in particular, during the period when the between‐events show the positive trend, the average logarithm of Δσ increases with an annual rate of 0.19 (i.e., the amplification factor for Δσ is 1.56 per year). Breakpoint analysis located a change in the linear trend coefficients of Δσ versus time in February 2014, although no large earthquakes occurred at that time. Finally, the temporal variability of Δσ mirrors the relative seismic‐velocity variations observed in previous studies for the same area and period, suggesting that both crack healing along the main fault system and healing of microcracks distributed at shallow depths throughout the surrounding region might be necessary to explain the wider observations of postearthquake recovery.
We present the results of a consistency check performed over the flatfile extracted from the engineering strong motion (ESM) database. The flatfile includes 23,014 recordings from 2179 earthquakes in the magnitude range from 3.5 to 7.8 that occurred since the 1970s in Europe and Middle East, as presented in the companion article by Lanzano et al. (Bull Earthq Eng, 2018a). The consistency check is developed by analyzing different residual distributions obtained from ad-hoc ground motion prediction equations for the absolute spectral acceleration (SA), displacement and Fourier amplitude spectra (FAS). Only recordings from earthquakes shallower than 40 km are considered in the analysis. The between-event, between-station and event-and-station corrected residuals are computed by applying a mixed-effect regression. We identified those earthquakes, stations, and recordings showing the largest deviations from the GMPE median predictions, and also evaluated the statistical uncertainty on the median model to get insights on the applicable magnitude–distance ranges and the usable period (or frequency) range. We observed that robust median predictions are obtained up to 8 s for SA and up to 20 Hz for FAS, although median predictions for Mw ≥ 7 show significantly larger uncertainties with ‘bumps’ starting above 5 s for SA and below 0.3 Hz for FAS. The between-station variance dominates over the other residual variances, and the dependence of the between-station residuals on logarithm of Vs30 is well-described by a piece-wise linear function with period-dependent slopes and hinge velocity around 580 m/s. Finally, we compared the between-event residuals obtained by considering two different sources of moment magnitude. The results show that, at long periods, the between-event terms from the two regressions have a weak correlation and the overall between-event variability is dissimilar, highlighting the importance of magnitude source in the regression results.
Natural thermoluminescence (TL) in rocks reflects a dynamic equilibrium between radiation-induced TL growth and decay via thermal and athermal pathways. When rocks exhume through Earth's crust and cool from high to low temperature, this equilibrium level increases as the temperature dependent thermal decay decreases. This phenomenon can be exploited to extract thermal histories of rocks. The main advantage of TL is that a single TL glow curve has a wide range of thermal stabilities (lifetime <ka to Ba), and hence can provide multiple constraints on thermal histories. Here we constrain the distribution of kinetic parameters of TL in feldspar using a glow curve deconvolution method and fitting infinitesimal sub-peaks using a general order kinetic model. Each peak corresponds to a different but closely located energy level E. Forward modeling is applied for different time-temperature histories to estimate the sensitivity and limitation of each signal for different cooling rates. The results show that it is possible to constrain thermal histories between ∼30 °C and ∼80 °C. The results also illustrate that shallower traps, i.e. with lower activation energies, can be exploited to constrain lower cooling histories >100 °C/Ma, whereas deeper traps, i.e. with higher activation energies, provide constraints on thermal histories for higher cooling rates (>300 °C/Ma). Finally, we show how the path of rock exhumation (i.e., depth vs. time) can be constrained using an inverse approach. The newly developed methodology is applied to rapidly cooled samples from the Namche Barwa massif, eastern Himalaya to suggest a trend in exhumation rate with time that follows an inverse correlation with global temperature and glaciers equilibrium altitude line (ELA).
Rock glaciers in semiarid mountains contain large amounts of ice and might be important water stores aside from glaciers, lakes, and rivers. Yet whether and how rock glaciers interact with river channels in mountain valleys remains largely unresolved. We examine the potential for rock glaciers to block or disrupt river channels, using a new inventory of more than 2000 intact rock glaciers that we mapped from remotely sensed imagery in the Karakoram (KR), Tien Shan (TS), and Altai (ALT) mountains. We find that between 5% and 14% of the rock glaciers partly buried, blocked, diverted or constricted at least 95 km of mountain rivers in the entire study area. We use a Bayesian robust logistic regression with multiple topographic and climatic inputs to discern those rock glaciers disrupting mountain rivers from those with no obvious impacts. We identify elevation and potential incoming solar radiation (PISR), together with the size of feeder basins, as dominant predictors, so that lower-lying and larger rock glaciers from larger basins are more likely to disrupt river channels. Given that elevation and PISR are key inputs for modelling the regional distribution of mountain permafrost from the positions of rock-glacier toes, we infer that river-blocking rock glaciers may be diagnostic of non-equilibrated permafrost. Principal component analysis adds temperature evenness and wet-season precipitation to the controls that characterise rock glaciers impacting on rivers. Depending on the choice of predictors, the accuracy of our classification is moderate to good with median posterior area-under-the-curve values of 0.71-0.89. Clarifying whether rapidly advancing rock glaciers can physically impound rivers, or fortify existing dams instead, deserves future field investigation. We suspect that rock-glacier dams are conspicuous features that have a polygenetic history and encourage more research on the geomorphic coupling between permafrost lobes, river channels, and the sediment cascades of semiarid mountain belts. (c) 2018 John Wiley & Sons, Ltd.
Soil pollutants such as hydrocarbons can induce toxic effects in plants and associated arbuscular mycorrhizal fungi (AMF). This study was conducted to evaluate if the legume Lotus corniculatus and the grass Elymus trachycaulus and arbuscular mycorrhizal fungi could grow in two oil sands processing by-products after bitumen extraction from the oil sands in northern Alberta, Canada. Substrate treatments were coarse tailings sand (CTS), a mix of dry mature fine tailings (MFT) with CTS (1: 1) and Pleistocene sandy soil (hydrocarbon free); microbial treatments were without AMF, with AMF and AMF plus soil bacteria isolated from oil sands reclamation sites. Plant biomass, root morphology, leaf water content, shoot tissue phosphorus content and mycorrhizal colonization were evaluated. Both plant species had reduced growth in CTS and tailings mix relative to sandy soil. AMF frequency and intensity in roots of E. trachycaulus was not influenced by soil hydrocarbons; however, it decreased significantly over time in roots of L. corniculatus without bacteria in CTS. Mycorrhizal inoculation alone did not significantly improve plant growth in CTS and tailings mix; however, inoculation with mycorrhizae plus bacteria led to a significantly positive response of both plant species in CTS. Thus, combined inoculation with selected mycorrhizae and bacteria led to synergistic effects. Such combinations may be used in future to improve plant growth in reclamation of CTS and tailings mix.
Granitic melt inclusions were found in layers of garnet clinopyroxenites from orogenic peridotites hosted in high-pressure felsic granulites of the Granulitgebirge, central Europe. The inclusions are both glassy and crystallized, and occur as clusters in the garnet. Microstructural features suggest that the inclusions formed while garnet was growing as a peritectic phase, likely alongside clinopyroxene. The chemistry of the melt, in particular its trace element signature, shows a crustal contribution, probably due to the involvement of phengite in the melt-producing reaction, most likely in the presence of a fluid. The presence of a granitoid melt in mantle rocks may be the result of localized melting of a phengite-bearing protolith either already present in the peridotites or, more likely, within the local deeply subducted crustal units. In the latter case, the melt would have infiltrated the peridotites and generated pyroxenite via metasomatism. In either case, the presence of granitoid inclusions in orogenic peridotite provides direct evidence for a genetic connection between a high-pressure crustal melt and garnet pyroxenites. The in situ characterization of these remnants of natural melt provides direct quantitative constraints on (one of) the agents responsible for the interaction between crust and mantle.
One of the main purposes of detrital thermochronology is to provide constraints on the regional-scale exhumation rate and its spatial variability in actively eroding mountain ranges. Procedures that use cooling age distributions coupled with hypsometry and thermal models have been developed in order to extract quantitative estimates of erosion rate and its spatial distribution, assuming steady state between tectonic uplift and erosion. This hypothesis precludes the use of these procedures to assess the likely transient response of mountain belts to changes in tectonic or climatic forcing. Other methods are based on an a priori knowledge of the in situ distribution of ages to interpret the detrital age distributions. In this paper, we describe a simple method that, using the observed detrital mineral age distributions collected along a river, allows us to extract information about the relative distribution of erosion rates in an eroding catchment without relying on a steady-state assumption, the value of thermal parameters or an a priori knowledge of in situ age distributions. The model is based on a relatively low number of parameters describing lithological variability among the various sub-catchments and their sizes and only uses the raw ages. The method we propose is tested against synthetic age distributions to demonstrate its accuracy and the optimum conditions for it use. In order to illustrate the method, we invert age distributions collected along the main trunk of the Tsangpo-Siang-Brahmaputra river system in the eastern Himalaya. From the inversion of the cooling age distributions we predict present-day erosion rates of the catchments along the Tsangpo-Siang-Brahmaputra river system, as well as some of its tributaries. We show that detrital age distributions contain dual information about present-day erosion rate, i. e., from the predicted distribution of surface ages within each catchment and from the relative contribution of any given catchment to the river distribution. The method additionally allows comparing modern erosion rates to long-term exhumation rates. We provide a simple implementation of the method in Python code within a Jupyter Notebook that includes the data used in this paper for illustration purposes.
The flash-flood in Braunsbach in the north-eastern part of Baden-Wuerttemberg/Germany was a particularly strong and concise event which took place during the floods in southern Germany at the end of May/early June 2016. This article presents a detailed analysis of the hydro-meteorological forcing and the hydrological consequences of this event. A specific approach, the "forensic hydrological analysis" was followed in order to include and combine retrospectively a variety of data from different disciplines. Such an approach investigates the origins, mechanisms and course of such natural events if possible in a "near real time" mode, in order to follow the most recent traces of the event. The results show that it was a very rare rainfall event with extreme intensities which, in combination with catchment properties, led to extreme runoff plus severe geomorphological hazards, i.e. great debris flows, which together resulted in immense damage in this small rural town Braunsbach. It was definitely a record-breaking event and greatly exceeded existing design guidelines for extreme flood discharge for this region, i.e. by a factor of about 10. Being such a rare or even unique event, it is not reliably feasible to put it into a crisp probabilistic context. However, one can conclude that a return period clearly above 100 years can be assigned for all event components: rainfall, peak discharge and sediment transport. Due to the complex and interacting processes, no single flood cause or reason for the very high damage can be identified, since only the interplay and the cascading characteristics of those led to such an event. The roles of different human activities on the origin and/or intensification of such an extreme event are finally discussed. (C) 2018 Elsevier B.V. All rights reserved.
How information about sediment transport processes is transmitted to the sedimentary record remains a complex problem for the interpretation of fluvial stratigraphy. Alluvial fan deposits represent the condensed archive of sediment transport, which is at least partly controlled by tectonics and climate. For three coupled catchment-fan systems in northern Death Valley, California, we measure grain size across 12 well-preserved Holocene and late-Pleistocene surfaces, mapped in detail from field observations and remote sensing. Our results show that fan surfaces correlated to the late Pleistocene are, on average, 30-50% coarser than active or Holocene fan surfaces. We adopt a self-similar form of grain size distribution based on the observed stability of the ratio between mean grain size and standard deviation downstream. Using statistical analysis, we show that fan surface grain size distributions are self-similar. We derive a relative mobility function using our self-similar grain size distributions, which describes the relative probability of a given grain size being transported. We show that the largest mobile grain sizes are between 20 and 35mm, a value that varies over time and is clearly lower in the Holocene than in the Pleistocene; a change we suggest is due to a drier climate in the Holocene. These results support recent findings that alluvial fan sedimentology can record past environmental change and that these landscapes are potentially sensitive to climatic change over a glacial-interglacial cycle. We demonstrate that the self-similarity methodology offers a means to explore changes in relative mobility of grain sizes from preserved fluvial deposits. Plain Language Summary A key challenge in Earth Science is understanding how landscapes respond to climate. It may be possible to observe measurable differences in certain landscapes settings such as alluvial fans in desert regions. Alluvial fans are believed to be effective recorders of climate, representing a cumulative store of material transported downstream by rainfall-sensitive river systems. In northern Death Valley, California, we measure at high resolution grain size on three alluvial fans with surfaces that date from the Holocene and the arid climate of today to the 20-40% wetter late-Pleistocene epoch. We find that older late-Pleistocene surfaces are coarser on average than surfaces deposited during the modern and Holocene dry period, suggesting a changing sediment transport regime potentially in response to precipitation. We also show that measured grain size distributions within and between surfaces can be successfully normalized based on the decay in mean grain size and variance downstream, exhibiting a self-similar pattern. Finally, we employ a grain size relative mobility model using our field data to establish which grain sizes are likely to be in transport or locked in the substrate. This model predicts that during the wetter late-Pleistocene mobile grain sizes are up to 40% larger than during the Holocene.
Continental rift systems open up unique possibilities to study the geodynamic system of our planet: geodynamic localization processes are imprinted in the morphology of the rift by governing the time-dependent activity of faults, the topographic evolution of the rift or by controlling whether a rift is symmetric or asymmetric. Since lithospheric necking localizes strain towards the rift centre, deformation structures of previous rift phases are often well preserved and passive margins, the end product of continental rifting, retain key information about the tectonic history from rift inception to continental rupture.
Current understanding of continental rift evolution is based on combining observations from active rifts with data collected at rifted margins. Connecting these isolated data sets is often accomplished in a conceptual way and leaves room for subjective interpretation. Geodynamic forward models, however, have the potential to link individual data sets in a quantitative manner, using additional constraints from rock mechanics and rheology, which allows to transcend previous conceptual models of rift evolution. By quantifying geodynamic processes within continental rifts, numerical modelling allows key insight to tectonic processes that operate also in other plate boundary settings, such as mid ocean ridges, collisional mountain chains or subduction zones.
In this thesis, I combine numerical, plate-tectonic, analytical, and analogue modelling approaches, whereas numerical thermomechanical modelling constitutes the primary tool. This method advanced rapidly during the last two decades owing to dedicated software development and the availability of massively parallel computer facilities. Nevertheless, only recently the geodynamical modelling community was able to capture 3D lithospheric-scale rift dynamics from onset of extension to final continental rupture.
The first chapter of this thesis provides a broad introduction to continental rifting, a summary of the applied rift modelling methods and a short overview of previews studies. The following chapters, which constitute the main part of this thesis feature studies on plate boundary dynamics in two and three dimension followed by global scale analyses (Fig. 1).
Chapter II focuses on 2D geodynamic modelling of rifted margin formation. It highlights the formation of wide areas of hyperextended crustal slivers via rift migration as a key process that affected many rifted margins worldwide. This chapter also contains a study of rift velocity evolution, showing that rift strength loss and extension velocity are linked through a dynamic feed-back. This process results in abrupt accelerations of the involved plates during rifting illustrating for the first time that rift dynamics plays a role in changing global-scale plate motions. Since rift velocity affects key processes like faulting, melting and lower crustal flow, this study also implies that the slow-fast velocity evolution should be imprinted in rifted margin structures.
Chapter III relies on 3D Cartesian rift models in order to investigate various aspects of rift obliquity. Oblique rifting occurs if the extension direction is not orthogonal to the rift trend. Using 3D lithospheric-scale models from rift initialisation to breakup I could isolate a characteristic evolution of dominant fault orientations. Further work in Chapter III addresses the impact of rift obliquity on the strength of the rift system. We illustrate that oblique rifting is mechanically preferred over orthogonal rifting, because the brittle yielding requires a lower tectonic force. This mechanism elucidates rift competition during South Atlantic rifting, where the more oblique Equatorial Atlantic Rift proceeded to breakup while the simultaneously active but less oblique West African rift system became a failed rift. Finally this Chapter also investigates the impact of a previous rift phase on current tectonic activity in the linkage area of the Kenyan with Ethiopian rift. We show that the along strike changes in rift style are not caused by changes in crustal rheology. Instead the rift linkage pattern in this area can be explained when accounting for the thinned crust and lithosphere of a Mesozoic rift event.
Chapter IV investigates rifting from the global perspective. A first study extends the oblique rift topic of the previous chapter to global scale by investigating the frequency of oblique rifting during the last 230 million years. We find that approximately 70% of all ocean-forming rift segments involved an oblique component of extension where obliquities exceed 20°. This highlights the relevance of 3D approaches in modelling, surveying, and interpretation of many rifted margins. In a final study, we propose a link between continental rift activity, diffuse CO2 degassing and Mesozoic/Cenozoic climate changes. We used recent CO2 flux measurements in continental rifts to estimate worldwide rift-related CO2 release, which we based on the global extent of rifts through time. The first-order correlation to paleo-atmospheric CO2 proxy data suggests that rifts constitute a major element of the global carbon cycle.
Movements of tectonic plates often induce oblique deformation at divergent plate boundaries. This is in striking contrast with traditional conceptual models of rifting and rifted margin formation, which often assume 2-D deformation where the rift velocity is oriented perpendicular to the plate boundary. Here we quantify the validity of this assumption by analysing the kinematics of major continent-scale rift systems in a global plate tectonic reconstruction from the onset of Pangea breakup until the present day. We evaluate rift obliquity by joint examination of relative extension velocity and local rift trend using the script-based plate reconstruction software pyGPlates. Our results show that the global mean rift obliquity since 230 Ma amounts to 34 degrees with a standard deviation of 24 degrees, using the convention that the angle of obliquity is spanned by extension direction and rift trend normal. We find that more than similar to 70 % of all rift segments exceeded an obliquity of 20 degrees demonstrating that oblique rifting should be considered the rule, not the exception. In many cases, rift obliquity and extension velocity increase during rift evolution (e.g. Australia-Antarctica, Gulf of California, South Atlantic, India-Antarctica), which suggests an underlying geodynamic correlation via obliquity-dependent rift strength. Oblique rifting produces 3-D stress and strain fields that cannot be accounted for in simplified 2-D plane strain analysis. We therefore highlight the importance of 3-D approaches in modelling, surveying, and interpretation of most rift segments on Earth where oblique rifting is the dominant mode of deformation.
Protection motivation theory (PMT) has become a popular theory to explain the risk-reducing behavior of residents against natural hazards. PMT captures the two main cognitive processes that individuals undergo when faced with a threat, namely, threat appraisal and coping appraisal. The latter describes the evaluation of possible response measures that may reduce or avert the perceived threat. Although the coping appraisal component of PMT was found to be a better predictor of protective intentions and behavior, little is known about the factors that influence individuals’ coping appraisals of natural hazards. More insight into flood-coping appraisals of PMT, therefore, are needed to better understand the decision-making process of individuals and to develop effective risk communication strategies. This study presents the results of two surveys among more than 1,600 flood-prone households in Germany and France. Five hypotheses were tested using multivariate statistics regarding factors related to flood-coping appraisals, which were derived from the PMT framework, related literature, and the literature on social vulnerability. We found that socioeconomic characteristics alone are not sufficient to explain flood-coping appraisals. Particularly, observational learning from the social environment, such as friends and neighbors, is positively related to flood-coping appraisals. This suggests that social norms and networks play an important role in flood-preparedness decisions. Providing risk and coping information can also have a positive effect. Given the strong positive influence of the social environment on flood-coping appraisals, future research should investigate how risk communication can be enhanced by making use of the observed social norms and network effects.
We propose a reduced dynamical system describing the coupled evolution of fluid flow and magnetic field at the top of the Earth's core between the years 1900 and 2014. The flow evolution is modeled with a first-order autoregressive process, while the magnetic field obeys the classical frozen flux equation. An ensemble Kalman filter algorithm serves to constrain the dynamics with the geomagnetic field and its secular variation given by the COV-OBS.x1 model. Using a large ensemble with 40,000 members provides meaningful statistics including reliable error estimates. The model highlights two distinct flow scales. Slowly varying large-scale elements include the already documented eccentric gyre. Localized short-lived structures include distinctly ageostophic features like the high-latitude polar jet on the Northern Hemisphere. Comparisons with independent observations of the length-of-day variations not only validate the flow estimates but also suggest an acceleration of the geostrophic flows over the last century. Hindcasting tests show that our model outperforms simpler predictions bases (linear extrapolation and stationary flow). The predictability limit, of about 2,000 years for the magnetic dipole component, is mostly determined by the random fast varying dynamics of the flow and much less by the geomagnetic data quality or lack of small-scale information.
Crop yield variations are strongly influenced by the spatial and temporal availabilities of water and nitrogen in the soil during the crop growth season. To estimate the quantities and distributions of water and nitrogen within a given soil, process-oriented soil models have often been used. These models require detailed information about the soil characteristics and profile architecture (e.g., soil depth, clay content, bulk density, field capacity and wilting point), but high resolution information about these soil properties, both vertically and laterally, is difficult to obtain through conventional approaches. However, on-the-go electrical resistivity tomography (ERT) measurements of the soil and data inversion tools have recently improved the lateral resolutions of the vertically distributed measurable information. Using these techniques, nearly 19,000 virtual soil profiles with defined layer depths were successfully created for a 30 ha silty cropped soil over loamy and sandy substrates in Central Germany, which were used to initialise the CArbon and Nitrogen DYnamics (CANDY) model. The soil clay content was derived from the electrical resistivity (ER) and the collected soil samples using a simple linear regression approach (the mean R-2 of clay = 0.39). The additional required structural and hydrological properties were derived from pedotransfer functions. The modelling results, derived soil texture distributions and original ER data were compared with the spatial winter wheat yield distribution in a relatively dry year using regression and boundary line analysis. The yield variation was best explained by the simulated soil water content (R-2 = 0.18) during the grain filling and was additionally validated by the measured soil water content with a root mean square error (RMSE) of 7.5 Vol%.
The literature contains a sizable number of publications where weather types are used to decompose climate shifts or trends into contributions of frequency and mean of those types. They are all based on the product rule, that is, a transformation of a product of sums into a sum of products, the latter providing the decomposition. While there is nothing to argue about the transformation itself, its interpretation as a climate shift or trend decomposition is bound to fail. While the case of a climate shift may be viewed as an incomplete description of a more complex behaviour, trend decomposition indeed produces bogus trends, as demonstrated by a synthetic counterexample with well-defined trends in type frequency and mean. Consequently, decompositions based on that transformation, be it for climate shifts or trends, must not be used.
Prior to the development of the modern Amazonian drainage network during the Late Miocene-Early Pliocene, large areas of western Amazonia may have been occupied by an extensive lacustrine and wetland environment known as the Pebas and Acre systems. These depositional systems are thought to have formed in response to foreland subsidence east of the uplifting Andes. Based on the occurrence at a present-day elevation of 1600 m of fossil pearly freshwater mussels of the genus Anodontites (indet. species) in intermontane Chota Basin of northern Ecuador (the westernmost South American fossil record for this genus), we discuss their potential implications for understanding of westernmost limit of the Miocene wetland ecosystem and consequently a later timing for the regional uplift of the Eastern Cordillera.
The Atlantic meridional overturning circulation (AMOC)—a system of ocean currents in the North Atlantic—has a major impact on climate, yet its evolution during the industrial era is poorly known owing to a lack of direct current measurements. Here we provide evidence for a weakening of the AMOC by about 3 ± 1 sverdrups (around 15 per cent) since the mid-twentieth century. This weakening is revealed by a characteristic spatial and seasonal sea-surface temperature ‘fingerprint’—consisting of a pattern of cooling in the subpolar Atlantic Ocean and warming in the Gulf Stream region—and is calibrated through an ensemble of model simulations from the CMIP5 project. We find this fingerprint both in a high-resolution climate model in response to increasing atmospheric carbon dioxide concentrations, and in the temperature trends observed since the late nineteenth century. The pattern can be explained by a slowdown in the AMOC and reduced northward heat transport, as well as an associated northward shift of the Gulf Stream. Comparisons with recent direct measurements from the RAPID project and several other studies provide a consistent depiction of record-low AMOC values in recent years.
The oldest ice core records are obtained from the East Antarctic Plateau. Water isotopes are key proxies to reconstructing past climatic conditions over the ice sheet and at the evaporation source. The accuracy of climate reconstructions depends on knowledge of all processes affecting water vapour, precipitation and snow isotopic compositions. Fractionation processes are well understood and can be integrated in trajectory-based Rayleigh distillation and isotope-enabled climate models. However, a quantitative understanding of processes potentially altering snow isotopic composition after deposition is still missing. In low-accumulation sites, such as those found in East Antarctica, these poorly constrained processes are likely to play a significant role and limit the interpretability of an ice core's isotopic composition.
By combining observations of isotopic composition in vapour, precipitation, surface snow and buried snow from Dome C, a deep ice core site on the East Antarctic Plateau, we found indications of a seasonal impact of metamorphism on the surface snow isotopic signal when compared to the initial precipitation. Particularly in summer, exchanges of water molecules between vapour and snow are driven by the diurnal sublimation-condensation cycles. Overall, we observe in between precipitation events modification of the surface snow isotopic composition. Using high-resolution water isotopic composition profiles from snow pits at five Antarctic sites with different accumulation rates, we identified common patterns which cannot be attributed to the seasonal variability of precipitation. These differences in the precipitation, surface snow and buried snow isotopic composition provide evidence of post-deposition processes affecting ice core records in low-accumulation areas.
Small repeating earthquakes are thought to represent rupture of isolated asperities loaded by surrounding creep. The observed scaling between recurrence interval and seismic moment, T-r approximate to M-1/6, contrasts with expectation assuming constant stress drop and no aseismic slip (T-r approximate to M-1/3). Here we demonstrate that simple crack models of velocity-weakening asperities in a velocity-strengthening fault predict the M-1/6 scaling; however, the mechanism depends on asperity radius, R. For small asperities ( , where R is the nucleation radius) numerical simulations with rate-state friction show interseismic creep penetrating inward from the edge, and earthquakes nucleate in the center and rupture the entire asperity. Creep penetration accounts for approximate to 25% of the slip budget, the nucleation phase takes up a larger fraction of slip. Stress drop increases with increasing R; the lack of self-similarity being due to the finite nucleation dimension. For 2R<R less than or similar to 6Rsimulations exhibit simple cycles with ruptures nucleating from the edge. Asperities with R6R exhibit complex cycles of partial and full ruptures. Here T-r is explained by an energy criterion: full rupture requires that the energy release rate everywhere on the asperity at least equals the fracture energy, leading to the scaling T-r approximate to M-1/6. Remarkably, in spite of the variability in behavior with source dimension, the scaling of T-r with stress drop , nucleation length and creep rate v(pl) is the same across all regimes: Tr approximate to vpl. This supports the use of repeating earthquakes as creepmeters and provides a physical interpretation for the scaling observed in nature. Plain Language Summary While most earthquake sequences have complex temporal patterns, some small earthquakes are quite predictable: they repeat periodically. The time between consecutive events (recurrence interval) grows with earthquake size: as intuitive, it takes longer to accumulate the mechanical energy for large earthquakes. However, the scaling between the recurrence interval and earthquake energy (seismic moment) is not what simple physical considerations predict. It is often assumed that faults are locked between events and seismic slip must therefore keep up with long-term plate motion. This leads to the scaling: Tr approximate to M01/3, but the observed scaling is . In fact, faults are not fully locked between earthquakes: they can slip slowly, or release part of the energy in smaller quakes between the larger ones. Here we use numerical simulations, and ideas from fracture mechanics, to understand what controls the time between repeating quakes. The main results are (1) analytical expressions of the recurrence interval as a function of earthquake size, predicting the observed scaling; (2) explanation of the differences between the cycle of small and large earthquakes (fraction of slow slip, direction of rupture propagation, and the occurrence of smaller quakes between large ones) and the quantities determining these transitions.
Kim et al. recently measured the structure factor of deeply supercooled water droplets (Reports, 22 December 2017, p. 1589). We raise several concerns about their data analysis and interpretation. In our opinion, the reported data do not lead to clear conclusions about the origins of water’s anomalies.
A transparent and data-driven global tectonic regionalization model for seismic hazard assessment
(2018)
A key concept that is common to many assumptions inherent within seismic hazard assessment is that of tectonic similarity. This recognizes that certain regions of the globe may display similar geophysical characteristics, such as in the attenuation of seismic waves, the magnitude scaling properties of seismogenic sources or the seismic coupling of the lithosphere. Previous attempts at tectonic regionalization, particularly within a seismic hazard assessment context, have often been based on expert judgements; in most of these cases, the process for delineating tectonic regions is neither reproducible nor consistent from location to location. In this work, the regionalization process is implemented in a scheme that is reproducible, comprehensible from a geophysical rationale, and revisable when new relevant data are published. A spatial classification-scheme is developed based on fuzzy logic, enabling the quantification of concepts that are approximate rather than precise. Using the proposed methodology, we obtain a transparent and data-driven global tectonic regionalization model for seismic hazard applications as well as the subjective probabilities (e.g. degree of being active/degree of being cratonic) that indicate the degree to which a site belongs in a tectonic category.
Coastal ecosystems in the Arctic are affected by climate change. As summer rainfall frequency and intensity are projected to increase in the future, more organic matter, nutrients and sediment could bemobilized and transported into the coastal nearshore zones. However, knowledge of current processes and future changes is limited. We investigated streamflow dynamics and the impacts of summer rainfall on lateral fluxes in a small coastal catchment on Herschel Island in the western Canadian Arctic. For the summer monitoring periods of 2014-2016, mean dissolved organic matter flux over 17 days amounted to 82.7 +/- 30.7 kg km(-2) and mean total dissolved solids flux to 5252 +/- 1224 kg km(-2). Flux of suspended sediment was 7245 kg km(-2) in 2015, and 369 kg km(-2) in 2016. We found that 2.0% of suspended sediment was composed of particulate organic carbon. Data and hysteresis analysis suggest a limited supply of sediments; their interannual variability is most likely caused by short-lived localized disturbances. In contrast, our results imply that dissolved organic carbon is widely available throughout the catchment and exhibits positive linear relationship with runoff. We hypothesize that increased projected rainfall in the future will result in a similar increase of dissolved organic carbon fluxes.
With the growing size and use of night light time series from the Visible Infrared Imaging Radiometer Suite Day/Night Band (DNB), it is important to understand the stability of the dataset. All satellites observe differences in pixel values during repeat observations. In the case of night light data, these changes can be due to both environmental effects and changes in light emission. Here we examine the stability of individual locations of particular large scale light sources (e.g., airports and prisons) in the monthly composites of DNB data from April 2012 to September 2017. The radiances for individual pixels of most large light emitters are approximately normally distributed, with a standard deviation of typically 15-20% of the mean. Greenhouses and flares, however, are not stable sources. We observe geospatial autocorrelation in the monthly variations for nearby sites, while the correlation for sites separated by large distances is small. This suggests that local factors contribute most to the variation in the pixel radiances and furthermore that averaging radiances over large areas will reduce the total variation. A better understanding of the causes of temporal variation would improve the sensitivity of DNB to lighting changes.
Himalayan rivers are frequently hit by catastrophic floods that are caused by the failure of glacial lake and landslide dams; however, the dynamics and long-term impacts of such floods remain poorly understood. We present a comprehensive set of observations that capture the July 2016 glacial lake outburst flood (GLOF) in the Bhotekoshi/Sunkoshi River of Nepal. Seismic records of the flood provide new insights into GLOF mechanics and their ability to mobilize large boulders that otherwise prevent channel erosion. Because of this boulder mobilization, GLOF impacts far exceed those of the annual summer monsoon, and GLOFs may dominate fluvial erosion and channel-hillslope coupling many tens of kilometers downstream of glaciated areas. Long-term valley evolution in these regions may therefore be driven by GLOF frequency and magnitude, rather than by precipitation.
Reducing uncertainties about carbon cycling is important in the Arctic where rapid environmental changes contribute to enhanced mobilization of carbon. Here we quantify soil organic carbon (SOC) contents of permafrost soils along the Yukon Coastal Plain and determine the annual fluxes from coastal erosion. Different terrain units were assessed based on surficial geology, morphology, and ground ice conditions. To account for the volume of wedge ice and massive ice in a unit, SOC contents were reduced by 19% and sediment contents by 16%. The SOC content in a 1m(2) column of soil varied according to the height of the bluff, ranging from 30 to 662kg, with a mean value of 183kg. Forty-four per cent of the SOC was within the top 1m of soil and values varied based on surficial materials, ranging from 30 to 53kg C/m(3), with a mean of 41kg. Eighty per cent of the shoreline was erosive with a mean annual rate of change of -0.7m/yr. This resulted in a SOC flux per meter of shoreline of 132kg C/m/yr, and a total flux for the entire 282km of the Yukon coast of 35.5 x 10(6) kg C/yr (0.036 Tg C/yr). The mean flux of sediment per meter of shoreline was 5.3 x 10(3) kg/m/yr, with a total flux of 1,832 x 10(6)kg/yr (1.832 Tg/yr). Sedimentation rates indicate that approximately 13% of the eroded carbon was sequestered in nearshore sediments, where the overwhelming majority of organic carbon was of terrestrial origin. Plain Language Summary The oceans help slow the buildup of carbon dioxide (CO2) because they absorb much of this greenhouse gas. However, if carbon from other sources is added to the oceans, it can affect their ability to absorb atmospheric CO2. Our study examines the organic carbon added to the Canadian Beaufort Sea from eroding permafrost along the Yukon coast, a region quite vulnerable to erosion. Understanding carbon cycling in this area is important because environmental changes in the Arctic such as longer open water seasons, rising sea levels, and warmer air, water and soil temperatures are likely to increase coastal erosion and, thus, carbon fluxes to the sea. We measured the carbon in different types of permafrost soils and applied corrections to account for the volume taken up by various types of ground ice. By determining how quickly the shoreline is eroding, we assessed how much organic carbon is being transferred to the ocean each year. Our results show that 36 x 10(6) kg of carbon is added annually from this section of the coast. If we extrapolate these results to other coastal areas along the Canadian Beaufort Sea, the flux of organic carbon is nearly 3 times what was previously thought.
Thermokarst lakes are prevalent in Arctic coastal lowland regions and sublake permafrost degradation and talik development contributes to greenhouse gas emissions by tapping the large permafrost carbon pool. Whereas lateral thermokarst lake expansion is readily apparent through remote sensing and shoreline measurements, sublake thawed sediment conditions and talik growth are difficult to measure. Here we combine transient electromagnetic surveys with thermal modeling, backed up by measured permafrost properties and radiocarbon ages, to reveal closed-talik geometry associated with a thermokarst lake in continuous permafrost. To improve access to talik geometry data, we conducted surveys along three transient electromagnetic transects perpendicular to lakeshores with different decadal-scale expansion rates of 0.16, 0.38, and 0.58m/year. We modeled thermal development of the talik using boundary conditions based on field data from the lake, surrounding permafrost and a borehole, independent of the transient electromagnetics. A talik depth of 91m was determined from analysis of the transient electromagnetic surveys. Using a lake initiation age of 1400years before present and available subsurface properties the results from thermal modeling of the lake center arrived at a best estimate talk depth of 80m, which is on the same order of magnitude as the results from the transient electromagnetic survey. Our approach has provided a noninvasive estimate of talik geometry suitable for comparable settings throughout circum-Arctic coastal lowland regions.
We explore the potential of spaceborne radar (SR) observations from the Ku-band precipitation radars onboard the Tropical Rainfall Measuring Mission (TRMM) and Global Precipitation Measurement (GPM) satellites as a reference to quantify the ground radar (GR) reflectivity bias. To this end, the 3-D volume-matching algorithm proposed by Schwaller and Morris (2011) is implemented and applied to 5 years (2012–2016) of observations. We further extend the procedure by a framework to take into account the data quality of each ground radar bin. Through these methods, we are able to assign a quality index to each matching SR–GR volume, and thus compute the GR calibration bias as a quality-weighted average of reflectivity differences in any sample of matching GR–SR volumes. We exemplify the idea of quality-weighted averaging by using the beam blockage fraction as the basis of a quality index. As a result, we can increase the consistency of SR and GR observations, and thus the precision of calibration bias estimates. The remaining scatter between GR and SR reflectivity as well as the variability of bias estimates between overpass events indicate, however, that other error sources are not yet fully addressed. Still, our study provides a framework to introduce any other quality variables that are considered relevant in a specific context. The code that implements our analysis is based on the wradlib open-source software library, and is, together with the data, publicly available to monitor radar calibration or to scrutinize long series of archived radar data back to December 1997, when TRMM became operational.
We explore the potential of spaceborne radar (SR) observations from the Ku-band precipitation radars onboard the Tropical Rainfall Measuring Mission (TRMM) and Global Precipitation Measurement (GPM) satellites as a reference to quantify the ground radar (GR) reflectivity bias. To this end, the 3-D volume-matching algorithm proposed by Schwaller and Morris (2011) is implemented and applied to 5 years (2012–2016) of observations. We further extend the procedure by a framework to take into account the data quality of each ground radar bin. Through these methods, we are able to assign a quality index to each matching SR–GR volume, and thus compute the GR calibration bias as a quality-weighted average of reflectivity differences in any sample of matching GR–SR volumes. We exemplify the idea of quality-weighted averaging by using the beam blockage fraction as the basis of a quality index. As a result, we can increase the consistency of SR and GR observations, and thus the precision of calibration bias estimates. The remaining scatter between GR and SR reflectivity as well as the variability of bias estimates between overpass events indicate, however, that other error sources are not yet fully addressed. Still, our study provides a framework to introduce any other quality variables that are considered relevant in a specific context. The code that implements our analysis is based on the wradlib open-source software library, and is, together with the data, publicly available to monitor radar calibration or to scrutinize long series of archived radar data back to December 1997, when TRMM became operational.
Due to the challenges in upscaling daily climatic forcing to geological time, physically realistic models describing how rainfall drives fluvial erosion are lacking. To bridge this gap between short-term hydrology and long-term geomorphology, we derive a theoretical framework for long-term fluvial erosion rates driven by realistic climate by integrating an established stochastic-mechanistic model of hydrology into a threshold-stochastic formulation of stream power. The hydrological theory provides equations for the daily streamflow probability distribution as a function of climatic boundary conditions. The new parameters introduced are rooted firmly in established climatic and hydrological theory. This allows us to account for how fluvial erosion rates respond to changes in rainfall intensity, frequency, evapotranspiration rates, and soil moisture dynamics in a way that is consistent with existing theories. We use this framework to demonstrate how hydroclimatic conditions and erosion threshold magnitude control the degree of nonlinearity between steepness index and erosion rate. We find that hydrological processes can have a significant influence on how erosive a particular climatic forcing will be. By accounting for the influence of hydrology on fluvial erosion, we conclude that climate is an important control on erosion rates and long-term landscape evolution.
GrassPlot is a collaborative vegetation-plot database organised by the Eurasian Dry Grassland Group (EDGG) and listed in the Global Index of Vegetation-Plot Databases (GIVD ID EU-00-003). GrassPlot collects plot records (releves) from grasslands and other open habitats of the Palaearctic biogeographic realm. It focuses on precisely delimited plots of eight standard grain sizes (0.0001; 0.001;... 1,000 m(2)) and on nested-plot series with at least four different grain sizes. The usage of GrassPlot is regulated through Bylaws that intend to balance the interests of data contributors and data users. The current version (v. 1.00) contains data for approximately 170,000 plots of different sizes and 2,800 nested-plot series. The key components are richness data and metadata. However, most included datasets also encompass compositional data. About 14,000 plots have near-complete records of terricolous bryophytes and lichens in addition to vascular plants. At present, GrassPlot contains data from 36 countries throughout the Palaearctic, spread across elevational gradients and major grassland types. GrassPlot with its multi-scale and multi-taxon focus complements the larger international vegetationplot databases, such as the European Vegetation Archive (EVA) and the global database " sPlot". Its main aim is to facilitate studies on the scale-and taxon-dependency of biodiversity patterns and drivers along macroecological gradients. GrassPlot is a dynamic database and will expand through new data collection coordinated by the elected Governing Board. We invite researchers with suitable data to join GrassPlot. Researchers with project ideas addressable with GrassPlot data are welcome to submit proposals to the Governing Board.
One common approach to cope with floods is the implementation of structural flood protection measures, such as levees or flood-control reservoirs, which substantially reduce the probability of flooding at the time of implementation. Numerous scholars have problematized this approach. They have shown that increasing the levels of flood protection can attract more settlements and high-value assets in the areas protected by the new measures. Other studies have explored how structural measures can generate a sense of complacency, which can act to reduce preparedness. These paradoxical risk changes have been described as "levee effect", "safe development paradox" or "safety dilemma". In this commentary, we briefly review this phenomenon by critically analysing the intended benefits and unintended effects of structural flood protection, and then we propose an interdisciplinary research agenda to uncover these paradoxical dynamics of risk.
Climate and environmental changes since the Last Glacial Maximum in the tropical zone of West Africa are usually inferred from marine and continental records. In this study, the potential of carbonate pedo-sedimentary geosystems, i.e. Vertisol relics, to record paleoenvironmental changes in the southwestern part of Chad Basin are investigated. A multi-dating approach was applied on different pedogenic organo-mineral constituents. Optically stimulated luminescence (OSL) dating was performed on the soil K-rich feldspars and was combined with radiocarbon dating on both the inorganic (C-14(inorg)) and organic carbon (C-14(org)) soil fractions. Three main pedo-sedimentary processes were assessed over the last 20 ka BP: 1) the soil parent material deposition, from 18 ka to 12 ka BP (OSL), 2) the soil organic matter integration, from 11 cal ka to 8 cal ka BP (C-14(org)), and 3) the pedogenic carbonate nodule precipitation, from 7 cal ka to 5 cal ka BP (C-14(inorg)). These processes correlate well with the Chad Basin stratigraphy and West African records and are shown to be related to significant changes in the soil water balance responding to the evolution of continental hydrology during the Late Quaternary. The last phase affecting the Vertisol relics is the increase of erosion, which is hypothesized to be due to a decrease of the vegetation cover triggered by (i) the onset of drier conditions, possibly strengthened by (ii) anthropogenic pressure. Archaeological data from Far North Cameroon and northern Nigeria, as well as sedimentation times in Lake Tilla (northeastern Nigeria), were used to test these relationships. The increase of erosion is suggested to possibly occur between c. 3 cal ka and 1 cal ka BP. Finally, satellite images revealed similar geosystems all along the Sudano-Sahelian belt, and initial C-14(inorg) ages of the samples collected in four sites gave similar ages to those reported in this study. Consequently, the carbonate pedo-sedimentary geosystems are valuable continental paleoenvironmental archives and soil water balance proxies of the semiarid tropics of West Africa. (C) 2018 Elsevier Ltd. All rights reserved.
Sulfur is an important component in volcanic gases at the Earth surface but also present in the deep Earth in hydrothermal or magmatic fluids. Little is known about the evolution of such fluids during ascent in the crust. A new optical cell was developed for in situ Raman spectroscopic investigations on fluids allowing abrupt or continuous changes of pressure up to 200 MPa at temperatures up to 750 degrees C. The concept is based on a flexible gold bellow, which separates the sample fluid from the pressure medium water. To avoid reactions between aggressive fluids and the pressure cell, steel components in contact with the fluid are shielded by gold foil. The cell was tested to study redox reactions in fluids using aqueous ammonium sulfate solutions as a model system. During heating at constant pressure of 130 MPa, sulfate ions transform first to HSO4- ions and then to molecular units such as H2SO4. Variation of pressure shows that the stability of sulfate species relies on fluid density, i.e., highly charged species are stable only in high-density fluids. Partial decomposition of ammonium was evident above 550 degrees C by the occurrence of a nitrogen peak in the Raman spectra. Reduced sulfur species were observed above 700 degrees C by Raman signals near 2590 cm(-1) assigned to HS- and H2S. No clear evidence for the formation of sulfur dioxide was found in contrary to previous studies on aqueous H2SO4, suggesting very reducing conditions in our experiments. Fluid-mineral interaction was studied by inserting into the cell a small, semi-open capsule filled with a mixture of pyrite and pyrrhotite. Oxidation of the sample assembly was evident by transformation of pyrite to pyrrhotite. As a consequence, sulfide species were observed in the fluid already at temperatures of similar to 600 degrees C.
Reconstructing how rivers respond to changes in runoff or sediment supply by incising or aggrading has been pivotal in gauging the role of the Indian Summer Monsoon (ISM) as a geomorphic driver in the Himalayas. Here we present new chronological data for fluvial aggradation and incision from the Donga alluvial fan and the upper Alaknanda River, as well as a compilation of previous work. In addition to conventional OSL-SAR (Single-Aliquot Regenerative-Dose) dating method, we have tested and applied pulsed OSL (POSL) dating for quartz samples that include K-rich feldspar inclusions, which is expected to improve the applicability and validity of OSL ages in the Lesser Himalaya. For previously dated deposits, our OSL ages are shown to be systematically older than previously reported ages. These results suggest periods of aggradation in the Alaknanda and Dehradun Valleys mainly between similar to 25 and 35 ka. This most likely reflects decreased stream power during periods of weakened monsoon. In addition, in-situ cosmogenic beryllium-10 was used to infer bedrock surface exposure ages, which are interpreted as episodes of active fluvial erosion. Resulting exposure ages span from 3 to 6 ka, suggesting that strath terraces were exposed relatively recently, and incision was dominant through most of the Holocene. In combination, our results support precipitation-driven fluvial dynamics, which regulates the balance between stream power and sediment supply. On a larger spatial scale, however, fluvial dynamics are probably not spatially homogeneous as aggradation could have been taking place in adjacent catchments while incision dominated in the study area. (C) 2018 Elsevier Ltd. All rights reserved.
This study analyzes the influence of local and regional climatic factors on the stable isotopic composition of rainfall in the Vietnamese Mekong Delta (VMD) as part of the Asian monsoon region. It is based on 1.5 years of weekly rainfall samples. In the first step, the isotopic composition of the samples is analyzed by local meteoric water lines (LMWLs) and single-factor linear correlations. Additionally, the contribution of several regional and local factors is quantified by multiple linear regression (MLR) of all possible factor combinations and by relative importance analysis. This approach is novel for the interpretation of isotopic records and enables an objective quantification of the explained variance in isotopic records for individual factors. In this study, the local factors are extracted from local climate records, while the regional factors are derived from atmospheric backward trajectories of water particles. The regional factors, i.e., precipitation, temperature, relative humidity and the length of backward trajectories, are combined with equivalent local climatic parameters to explain the response variables delta O-18, delta H-2, and d-excess of precipitation at the station of measurement. The results indicate that (i) MLR can better explain the isotopic variation in precipitation (R-2 = 0.8) compared to single-factor linear regression (R-2 = 0.3); (ii) the isotopic variation in precipitation is controlled dominantly by regional moisture regimes (similar to 70 %) compared to local climatic conditions (similar to 30 %); (iii) the most important climatic parameter during the rainy season is the precipitation amount along the trajectories of air mass movement; (iv) the influence of local precipitation amount and temperature is not sig-nificant during the rainy season, unlike the regional precipitation amount effect; (v) secondary fractionation processes (e.g., sub-cloud evaporation) can be identified through the d-excess and take place mainly in the dry season, either locally for delta O-18 and delta H-2, or along the air mass trajectories for d-excess. The analysis shows that regional and local factors vary in importance over the seasons and that the source regions and transport pathways, and particularly the climatic conditions along the pathways, have a large influence on the isotopic composition of rainfall. Although the general results have been reported qualitatively in previous studies (proving the validity of the approach), the proposed method provides quantitative estimates of the controlling factors, both for the whole data set and for distinct seasons. Therefore, it is argued that the approach constitutes an advancement in the statistical analysis of isotopic records in rainfall that can supplement or precede more complex studies utilizing atmospheric models. Due to its relative simplicity, the method can be easily transferred to other regions, or extended with other factors. The results illustrate that the interpretation of the isotopic composition of precipitation as a recorder of local climatic conditions, as for example performed for paleorecords of water isotopes, may not be adequate in the southern part of the Indochinese Peninsula, and likely neither in other regions affected by monsoon processes. However, the presented approach could open a pathway towards better and seasonally differentiated reconstruction of paleoclimates based on isotopic records.
In this study, we analyze interactions in lake and lake catchment systems of a continuous permafrost area. We assessed colored dissolved organic matter (CDOM) absorption at 440 nm (a(440)(CDOM)) and absorption slope (S300-500) in lakes using field sampling and optical remote sensing data for an area of 350 km(2) in Central Yamal, Siberia. Applying a CDOM algorithm (ratio of green and red band reflectance) for two high spatial resolution multispectral GeoEye-1 and Worldview-2 satellite images, we were able to extrapolate the a()(CDOM) data from 18 lakes sampled in the field to 356 lakes in the study area (model R-2 = 0.79). Values of a(440)(CDOM) in 356 lakes varied from 0.48 to 8.35 m(-1) with a median of 1.43 m(-1). This a()(CDOM) dataset was used to relate lake CDOM to 17 lake and lake catchment parameters derived from optical and radar remote sensing data and from digital elevation model analysis in order to establish the parameters controlling CDOM in lakes on the Yamal Peninsula. Regression tree model and boosted regression tree analysis showed that the activity of cryogenic processes (thermocirques) in the lake shores and lake water level were the two most important controls, explaining 48.4% and 28.4% of lake CDOM, respectively (R-2 = 0.61). Activation of thermocirques led to a large input of terrestrial organic matter and sediments from catchments and thawed permafrost to lakes (n = 15, mean a(440)(CDOM) = 5.3 m(-1)). Large lakes on the floodplain with a connection to Mordy-Yakha River received more CDOM (n = 7, mean a(440)(CDOM) = 3.8 m(-1)) compared to lakes located on higher terraces.
In this study, we analyze interactions in lake and lake catchment systems of a continuous permafrost area. We assessed colored dissolved organic matter (CDOM) absorption at 440 nm (a(440)(CDOM)) and absorption slope (S300-500) in lakes using field sampling and optical remote sensing data for an area of 350 km(2) in Central Yamal, Siberia. Applying a CDOM algorithm (ratio of green and red band reflectance) for two high spatial resolution multispectral GeoEye-1 and Worldview-2 satellite images, we were able to extrapolate the a()(CDOM) data from 18 lakes sampled in the field to 356 lakes in the study area (model R-2 = 0.79). Values of a(440)(CDOM) in 356 lakes varied from 0.48 to 8.35 m(-1) with a median of 1.43 m(-1). This a()(CDOM) dataset was used to relate lake CDOM to 17 lake and lake catchment parameters derived from optical and radar remote sensing data and from digital elevation model analysis in order to establish the parameters controlling CDOM in lakes on the Yamal Peninsula. Regression tree model and boosted regression tree analysis showed that the activity of cryogenic processes (thermocirques) in the lake shores and lake water level were the two most important controls, explaining 48.4% and 28.4% of lake CDOM, respectively (R-2 = 0.61). Activation of thermocirques led to a large input of terrestrial organic matter and sediments from catchments and thawed permafrost to lakes (n = 15, mean a(440)(CDOM) = 5.3 m(-1)). Large lakes on the floodplain with a connection to Mordy-Yakha River received more CDOM (n = 7, mean a(440)(CDOM) = 3.8 m(-1)) compared to lakes located on higher terraces.
It is well-known that prey species often face trade-offs between defense against predation and competitiveness, enabling predator-mediated coexistence. However, we lack an understanding of how the large variety of different defense traits with different competition costs affects coexistence and population dynamics. Our study focusses on two general defense mechanisms, that is, pre-attack (e.g., camouflage) and post-attack defenses (e.g., weaponry) that act at different phases of the predator—prey interaction. We consider a food web model with one predator, two prey types and one resource. One prey type is undefended, while the other one is pre- or post-attack defended paying costs either by a higher half-saturation constant for resource uptake or a lower maximum growth rate. We show that post-attack defenses promote prey coexistence and stabilize the population dynamics more strongly than pre-attack defenses by interfering with the predator's functional response: Because the predator spends time handling “noncrackable” prey, the undefended prey is indirectly facilitated. A high half-saturation constant as defense costs promotes coexistence more and stabilizes the dynamics less than a low maximum growth rate. The former imposes high costs at low resource concentrations but allows for temporally high growth rates at predator-induced resource peaks preventing the extinction of the defended prey. We evaluate the effects of the different defense mechanisms and costs on coexistence under different enrichment levels in order to vary the importance of bottom-up and top-down control of the prey community.
Bedrock landsliding in mountain belts can elevate overall chemical weathering rates through rapid dissolution of exhumed reactive mineral phases in transiently stored deposits. This link between a key process of erosion and the resultant weathering affects the sequestering of carbon dioxide through weathering of silicate minerals and broader links between erosion in active orogens and climate change. Here we address the effect on the carbon cycle of weathering induced by bedrock landsliding in Taiwan and the Western Southern Alps of New Zealand. Using solute chemistry data from samples of seepage from landslide deposits and river discharge from catchments with variable proportions of landsliding, we model the proportion of silicate and carbonate weathering and the balance of sulfuric and carbonic acids that act as weathering agents. We correct for secondary precipitation, geothermal, and cyclic input, to find a closer approximation of the weathering explicitly occurring within landslide deposits. We find highly variable proportions of sulfuric and carbonic acids driving weathering in landslides and stable hillslopes. Despite this variability, the predominance of rapid carbonate weathering within landslides and catchments where mass wasting is prevalent results at best in limited sequestration of carbon dioxide by this process of rapid erosion. In many cases where sulfuric acid is a key weathering agent, a net release of CO2 to the atmosphere occurs. This suggests that a causal link between erosion in mountain belts and climate change through the sequestration of CO2, if it exists, must operate through a process other than chemical weathering driven by landsliding. Plain Language Summary There is a long-standing debate surrounding the link between erosion and climate. It is often suggested that as temperatures increase, rainier and stormier weather could increase erosion of rock; as that rock is exposed, silicate minerals within could break down, which, on long time scales, can remove CO2 from the atmosphere, lowering global temperatures and acting as a negative feedback. Recent studies have shown that landslide deposits are key locations for the link between chemical weathering and physical erosion in some mountain belts. To test how landslides affect the erosion-climate link, we used samples of water seeping through landslides in Taiwan and New Zealand to calculate the amount of carbon dioxide that is either absorbed or released through this chemical reaction. We find that the large amount of freshly exposed rock in Taiwanese landslide deposits contains significant carbonate rock and sulfide minerals; the net result of the weathering of these minerals is a release of carbon dioxide, which inverts the traditional perspective on the role erosion plays in controlling carbon dioxide release. In some mountain belts, it seems that increased erosion and resulting landsliding may act to increase carbon dioxide in the air, opening further questions into the nature of erosional-climatic links.
The scientific drilling campaign PALEOVAN was conducted in the summer of 2010 and was part of the international continental drilling programme (ICDP). The main goal of the campaign was the recovery of a sensitive climate archive in the East of Anatolia. Lacustrine deposits underneath the lake floor of ‘Lake Van’ constitute this archive. The drilled core material was recovered from two locations: the Ahlat Ridge and the Northern Basin. A composite core was constructed from cored material of seven parallel boreholes at the Ahlat Ridge and covers an almost complete lacustrine history of Lake Van. The composite record offered sensitive climate proxies such as variations of total organic carbon, K/Ca ratios, or a relative abundance of arboreal pollen. These proxies revealed patterns that are similar to climate proxy variations from Greenland ice cores. Climate variations in Greenland ice cores have been dated by modelling the timing of orbital forces to affect the climate. Volatiles from melted ice aliquots are often taken as high-resolution proxies and provide a base for fitting the according temporal models.
The ICDP PALEOVAN scientific team fitted proxy data from the lacustrine drilling record to ice core data and constructed an age model. Embedded volcaniclastic layers had to be dated radiometrically in order to provide independent age constraints to the climate-stratigraphic age model. Solving this task by an application of the 40Ar/39Ar method was the main objective of this thesis. Earlier efforts to apply the 40Ar/39Ar dating resulted in inaccuracies that could not be explained satisfactorily.
The absence of K-rich feldspars in suitable tephra layers implied that feldspar crystals needed to be 500 μm in size minimum, in order to apply single-crystal 40Ar/39Ar dating. Some of the samples did not contain any of these grain sizes or only very few crystals of that size. In order to overcome this problem this study applied a combined single-crystal and multi-crystal approach with different crystal fractions from the same sample. The preferred method of a stepwise heating analysis of an aliquot of feldspar crystals has been applied to three samples. The Na-rich crystals and their young geological age required 20 mg of inclusion-free, non-corroded feldspars. Small sample volumes (usually 25 % aliquots of 5 cm3 of sample material – a spoon full of tephra) and the widespread presence of melt-inclusion led to the application of combined single- and multigrain total fusion analyses. 40Ar/39Ar analyses on single crystals have the advantage of being able to monitor the presence of excess 40Ar and detrital or xenocrystic contamination in the samples. Multigrain analyses may hide the effects from these obstacles. The results from the multigrain analyses are therefore discussed with respect to the findings from the respective cogenetic single crystal ages. Some of the samples in this study were dated by 40Ar/39Ar on feldspars on multigrain separates and (if available) in combination with only a few single crystals. 40Ar/39Ar ages from two of the samples deviated statistically from the age model. All other samples resulted in identical ages. The deviations displayed older ages than those obtained from the age model. t-Tests compared radiometric ages with available age control points from various proxies and from the relative paleointensity of the earth magnetic field within a stratigraphic range of ± 10 m. Concordant age control points from different relative chronometers indicated that deviations are a result of erroneous 40Ar/39Ar ages. The thesis argues two potential reasons for these ages: (1) the irregular appearance of 40Ar from rare melt- and fluid- inclusions and (2) the contamination of the samples with older crystals due to a rapid combination of assimilation and ejection.
Another aliquot of feldspar crystals that underwent separation for the application of 40Ar/39Ar dating was investigated for geochemical inhomogeneities. Magmatic zoning is ubiquitous in the volcaniclastic feldspar crystals. Four different types of magmatic zoning were detected. The zoning types are compositional zoning (C-type zoning), pseudo-oscillatory zoning of trace ele- ment concentrations (PO-type zoning), chaotic and patchy zoning of major and trace element concentrations (R-type zoning) and concentric zoning of trace elements (CC-type zoning). Sam- ples that deviated in 40Ar/39Ar ages showed C-type zoning, R-type zoning or a mix of different types of zoning (C-type and PO-type). Feldspars showing PO-type zoning typically represent the smallest grain size fractions in the samples. The constant major element compositions of these crystals are interpreted to represent the latest stages in the compositional evolution of feldspars in a peralkaline melt. PO-type crystals contain less melt- inclusions than other zoning types and are rarely corroded. This thesis concludes that feldspars that show PO-type zoning are most promising chronometers for the 40Ar/39Ar method, if samples provide mixed zoning types of Quaternary anorthoclase feldspars.
Five samples were dated by applying the 40Ar/39Ar method to volcanic glass. High fractions of atmospheric Ar (typically > 98%) significantly hampered the precision of the 40Ar/39Ar ages and resulted in rough age estimates that widely overlap the age model. Ar isotopes indicated that the glasses bore a chorine-rich Ar-end member. The chlorine-derived 38Ar indicated chlorine-rich fluid-inclusions or the hydration of the volcanic glass shards. This indication strengthened the evidence that irregularly distributed melt-inclusions and thus irregular distributed excess 40Ar influenced the problematic feldspar 40Ar/39Ar ages. Whether a connection between a corrected initial 40Ar/36Ar ratio from glasses to the 40Ar/36Ar ratios from pore waters exists remains unclear.
This thesis offers another age model, which is similarly based on the interpolation of the temporal tie points from geophysical and climate-stratigraphic data. The model used a PCHIP- interpolation (piecewise cubic hermite interpolating polynomial) whereas the older age model used a spline-interpolation. Samples that match in ages from 40Ar/39Ar dating of feldspars with the earlier published age model were additionally assigned with an age from the PCHIP- interpolation. These modelled ages allowed a recalculation of the Alder Creek sanidine mineral standard. The climate-stratigraphic calibration of an 40Ar/39Ar mineral standard proved that the age versus depth interpolations from PAELOVAN drilling cores were accurate, and that the applied chronometers recorded the temporal evolution of Lake Van synchronously.
Petrochemical discrimination of the sampled volcaniclastic material is also given in this thesis. 41 from 57 sampled volcaniclastic layers indicate Nemrut as their provenance. Criteria that served for the provenance assignment are provided and reviewed critically. Detailed correlations of selected PALEOVAN volcaniclastics to onshore samples that were described in detail by earlier studies are also discussed. The sampled volcaniclastics dominantly have a thickness of < 40 cm and have been ejected by small to medium sized eruptions. Onshore deposits from these types of eruptions are potentially eroded due to predominant strong winds on Nemrut and Süphan slopes. An exact correlation with the data presented here is therefore equivocal or not possible at all.
Deviating feldspar 40Ar/39Ar ages can possibly be explained by inherited 40Ar from feldspar xenocrysts contaminating the samples. In order to test this hypothesis diffusion couples of Ba were investigated in compositionally zoned feldspar crystals. The diffusive behaviour of Ba in feldspar is known, and gradients in the changing concentrations allowed for the calculation of the duration of the crystal’s magmatic development since the formation of the zoning interface. Durations were compared with degassing scenarios that model the Ar-loss during assimilation and subsequent ejection of the xenocrystals. Diffusive equilibration of the contrasting Ba concentrations is assumed to generate maximum durations as the gradient could have been developed in several growth and heating stages. The modelling does not show any indication of an involvement of inherited 40Ar in any of the deviating samples. However, the analytical set-up represents the lower limit of the required spatial resolution. Therefore, it cannot be excluded that the degassing modelling relies on a significant overestimation of the maximum duration of the magmatic history. Nevertheless, the modelling of xenocrystal degassing evidences that the irregular incorporation of excess 40Ar by melt- and fluid inclusions represents the most critical problem that needs to be overcome in dating volcaniclastic feldspars from the PALEOVAN drill cores. This thesis provides the complete background in generating and presenting 40Ar/39Ar ages that are compared to age data from a climate-stratigraphic model. Deviations are identified statistically and then discussed in order to find explanations from the age model and/or from 40Ar/39Ar geochronology. Most of the PALEOVAN stratigraphy provides several chronometers that have been proven for their synchronicity. Lacustrine deposits from Lake Van represent a key archive for reconstructing climate evolution in the eastern Mediterranean and in the Near East. The PALEOVAN record offers a climate-stratigraphic age model with a remarkable accuracy and resolution.
Changes in species' distributions are classically projected based on their climate envelopes. For Siberian forests, which have a tremendous significance for vegetation-climate feedbacks, this implies future shifts of each of the forest-forming larch (Larix) species to the north-east. However, in addition to abiotic factors, reliable projections must assess the role of historical biogeography and biotic interactions. Here, we use sedimentary ancient DNA and individual-based modelling to investigate the distribution of larch species and mitochondrial haplotypes through space and time across the treeline ecotone on the southern Taymyr peninsula, which at the same time presents a boundary area of two larch species. We find spatial and temporal patterns, which suggest that forest density is the most influential driver determining the precise distribution of species and mitochondrial haplotypes. This suggests a strong influence of competition on the species' range shifts. These findings imply possible climate change outcomes that are directly opposed to projections based purely on climate envelopes. Investigations of such fine-scale processes of biodiversity change through time are possible using paleoenvironmental DNA, which is available much more readily than visible fossils and can provide information at a level of resolution that is not reached in classical palaeoecology.
The Himalayan arc stretches >2500 km from east to west at the southern edge of the Tibetan Plateau, representing one of the most important Cenozoic continent-continent collisional orogens. Internal deformation processes and climatic factors, which drive weathering, denudation, and transport, influence the growth and erosion of the orogen. During glacial times wet-based glaciers sculpted the mountain range and left overdeepend and U-shaped valleys, which were backfilled during interglacial times with paraglacial sediments over several cycles. These sediments partially still remain within the valleys because of insufficient evacuation capabilities into the foreland. Climatic processes overlay long-term tectonic processes responsible for uplift and exhumation caused by convergence. Possible processes accommodating convergence within the orogenic wedge along the main Himalayan faults, which divide the range into four major lithologic units, are debated. In this context, the identification of processes shaping the Earth’s surface on short- and on long-term are crucial to understand the growth of the orogen and implications for landscape development in various sectors along the arc. This thesis focuses on both surface and tectonic processes that shape the landscape in the western Indian Himalaya since late Miocene.
In my first study, I dated well-preserved glacially polished bedrock on high-elevated ridges and valley walls in the upper of the Chandra Valley the by means of 10Be terrestrial cosmogenic radionuclides (TCN). I used these ages and mapped glacial features to reconstruct the extent and timing of Pleistocene glaciation at the southern front of the Himalaya. I was able to reconstruct an extensive valley glacier of ~200 km length and >1000 m thickness. Deglaciation of the Chandra Valley glacier started subsequently to insolation increase on the Northern Hemisphere and thus responded to temperature increase. I showed that the timing this deglaciation onset was coeval with retreat of further midlatitude glaciers on the Northern and Southern Hemispheres. These comparisons also showed that the post-LGM deglaciation very rapid, occurred within a few thousand years, and was nearly finished prior to the Bølling/Allerød interstadial.
A second study (co-authorship) investigates how glacial advances and retreats in high mountain environments impact the landscape. By 10Be TCN dating and geomorphic mapping, we obtained maximal length and height of the Siachen Glacier within the Nubra Valley. Today the Shyok and Nubra confluence is backfilled with sedimentary deposits, which are attributed to the valley blocking of the Siachen Glacier 900 m above the present day river level. A glacial dam of the Siachen Glacier blocked the Shyok River and lead to the evolution of a more than 20 km long lake. Fluvial and lacustrine deposits in the valley document alternating draining and filling cycles of the lake dammed by the Siachen Glacier. In this study, we can show that glacial incision was outpacing fluvial incision.
In the third study, which spans the million-year timescale, I focus on exhumation and erosion within the Chandra and Beas valleys. In this study the position and discussed possible reasons of rapidly exhuming rocks, several 100-km away from one of the main Himalayan faults (MFT) using Apatite Fission Track (AFT) thermochronometry. The newly gained AFT ages indicate rapid exhumation and confirm earlier studies in the Chandra Valley. I assume that the rapid exhumation is most likely related to uplift over subsurface structures. I tested this hypothesis by combining further low-temperature thermochronometers from areas east and west of my study area. By comparing two transects, each parallel to the Beas/Chandra Valley transect, I demonstrate similarities in the exhumation pattern to transects across the Sutlej region, and strong dissimilarities in the transect crossing the Dhauladar Range. I conclude that the belt of rapid exhumation terminates at the western end of the Kullu-Rampur window. Therewith, I corroborate earlier studies suggesting changes in exhumation behavior in the western Himalaya. Furthermore, I discussed several causes responsible for the pronounced change in exhumation patterns along strike: 1) the role of inherited pre-collisional features such as the Proterozoic sedimentary cover of the Indian basement, former ridges and geological structures, and 2) the variability of convergence rates along the Himalayan arc due to an increased oblique component towards the syntaxis.
The combination of field observations (geological and geomorphological mapping) and methods to constrain short- and long-term processes (10Be, AFT) help to understand the role of the individual contributors to exhumation and erosion in the western Indian Himalaya. With the results of this thesis, I emphasize the importance of glacial and tectonic processes in shaping the landscape by driving exhumation and erosion in the studied areas.
Despite remarkable tectonostratigraphic similarities along the Himalayan arc, pronounced topographic and exhumational variability exists in different morphotectonic segments. The processes responsible for this segmentation are debated. Of particular interest is a 30- to 40-km-wide orogen-parallel belt of rapid exhumation that extends from central Nepal to the western Himalaya and its possible linkage to a midcrustal ramp in the basal decollement, and the related growth of Lesser Himalayan duplex structures. Here we present 26 new apatite fission track cooling ages from the Beas-Lahul region, at the transition from the Central to the Western Himalaya (77 degrees-78 degrees E) to investigate segmentation in the Himalayan arc from a thermochronologic perspective. Together with previously published data from this part of the orogen, we document significant lateral changes in exhumation between the Dhauladar Range to the west, the Beas-Lahul region, and the Sutlej area to the east of the study area. In contrast to the Himalayan front farther east, exhumation in the far western sectors is focused at the frontal parts of the mountain range and associated with the hanging wall of the Main Boundary Thrust fault ramp. Our results allow us to spatially correlate the termination of the rapid exhumation belt with a midcrustal ramp to the west. We suggest that a plunging anticline at the northwestern edge of the Larji-Kullu-Rampur window represents the termination of the Central Himalayan segment, which is related to the evolution of the Lesser Himalayan duplex. Key Points
Micropetrology
(2018)
Inclusions in minerals, whether fluids, melts or crystalline phases, are small pieces of the large-scale puzzle of Nature, time-consuming to investigate and often of difficult interpretation. Yet they are windows into the past of their host mineral. Mineral inclusions provide the opportunity to unravel the genesis of their host, and the increasingly refined understanding of their elastic behaviour provides the basis for alternative, equilibrium-independent geobarometry. Fluid and melt inclusions reveal information about material transfer in the Earth system, from shallow mineralization to mantle re-fertilization via subduction. The study of inclusions is thus one of the most intriguing and fertile branches of micropetrology. In this contribution, we focus on two recent developments: the use of elasticity models to extract the formation conditions of the host crystal, and the discovery and investigation of melt inclusions in metamorphic rocks. We also discuss how to evaluate the information provided by inclusions, given that they are no longer at the pressure and temperature conditions of entrapment. We discuss how to understand and quantify the changes undergone during cooling and depressurization, and how metastability-related phenomena in inclusions, such as crystallization of rare polymorphs and preservation of the original content of volatiles in fluid and melt inclusions, provide direct evidence that inclusions represent closed systems. The field of study of inclusions in minerals still has a largely untapped potential. The most fruitful avenues for future research will emerge from continuous technological innovation in analytical and imaging techniques, the application of experimental petrology, and the development and application of new theoretical models for coupled mineral behaviour under changing P-T conditions.
The central European Bohemian Massif has undergone over two centuries of scientific investigation which has made it a pivotal area for the development and testing of modern geological theories. The discovery of melt inclusions in high-grade rocks, either crystallized as nanogranitoids or as glassy inclusions, prompted the re-evaluation of the area with an ‘inclusionist’ eye. Melt inclusions have been identified in a wide range of rocks, including felsic/perpotassic granulites, migmatites, eclogites and garnet clinopyroxenites, all the result of melting events albeit over a wide range of pressure/temperature conditions (800–1000°C/0.5–5 GPa). This contribution provides an overview of such inclusions and discusses the qualitative and quantitative constraints they provide for melting processes, and the nature of melts and fluids involved in these processes. In particular, data on trace-element signatures of melt inclusions trapped at mantle depths are presented and discussed. Moreover, experimental re-homogenization of nanogranitoids provided microstructural criteria allowing assessment of the conditions at which melt and host are mutually stable during melting. Overall this work aims to provide guidelines and suggestions for petrologists wishing to explore the fascinating field of melt inclusions in metamorphic terranes worldwide, based on the newest discoveries from the still-enigmatic Bohemian Massif.
In the Pan-African belt of the Dronning Maud Land, Antarctica, crystallized melt inclusions (nanogranitoids) occur in garnet from ultramafic granulites. The granulites contain the peak assemblage pargasite+garnet+clinopyroxene with rare relict orthopyroxene and biotite, and retrograde symplectites at contacts between garnet and amphibole. Garnet contains two generations of melt inclusions. Type 1 inclusions, interpreted as primary, are isolated, < 10 mu m in size, and generally have negative crystal shapes. They contain kokchetavite, kumdykolite, and phlogopite, with quartz and zoisite as minor phases, and undevitrified glass was identified in one inclusion. Type 2 inclusions are < 30 mu m in size, secondary, and contain amphibole, feldspars, and zoisite. Type 2 inclusions appear to be the crystallization products of a melt that coexisted with an immiscible CO2-rich fluid. The nanogranitoids were re-homogenized after heating in a piston-cylinder in a series of four experiments to investigate their composition. The conditions ranged between 900 and 950 degrees C at 1.5-2.4 GPa. Type 1 inclusions are trachytic and ultrapotassic, whereas type 2 melts are dacitic to rhyolitic. Thermodynamic modeling of the ultramafic composition in the MnNCKFMASHTO system shows that anatexis occurred at the end of the prograde P-T path, between the solidus (at ca. 860 degrees C-1.4 GPa) and the peak conditions (at ca. 960 degrees C-1.7 GPa). The model melt composition is felsic and similar to that of type 1 inclusions, particularly when the melting degree is low (< 1 mol%), close to the solidus. However the modeling fails to reproduce the highly potassic signature of the melt and its low H2O content. The combination of petrology, melt inclusion study, and thermodynamic modeling supports the interpretation that melt was produced by anatexis of the ultramafic boudins near peak P-T conditions, and that type 1 inclusions contain the anatectic melt that was present during garnet growth. The felsic, ultrapotassic composition of the primary anatectic melts is compatible with low melting degrees in the presence of biotite and amphibole as reactants.
The Gutenberg-Richter relation for earthquake magnitudes is the most famous empirical law in seismology. It states that the frequency of earthquake magnitudes follows an exponential distribution; this has been found to be a robust feature of seismicity above the completeness magnitude, and it is independent of whether global, regional, or local seismicity is analyzed. However, the exponent b of the distribution varies significantly in space and time, which is important for process understanding and seismic hazard assessment; this is particularly true because of the fact that the Gutenberg-Richter b-value acts as a proxy for the stress state and quantifies the ratio of large-to-small earthquakes. In our work, we focus on the automatic detection of statistically significant temporal changes of the b-value in seismicity data. In our approach, we use Bayes factors for model selection and estimate multiple change-points of the frequency-magnitude distribution in time. The method is first applied to synthetic data, showing its capability to detect change-points as function of the size of the sample and the b-value contrast. Finally, we apply this approach to examples of observational data sets for which b-value changes have previously been stated. Our analysis of foreshock and after-shock sequences related to mainshocks, as well as earthquake swarms, shows that only a portion of the b-value changes is statistically significant.
Tillage erosion on arable land is a very important process leading to a net downslope movement of soil and soil constitutes. Tillage erosion rates are commonly in the same order of magnitude as water erosion rates and can be even higher, especially under highly mechanized agricultural soil management. Despite its prevalence and magnitude, tillage erosion is still understudied compared to water erosion. The goal of this study was to bring together experts using different techniques to determine tillage erosion and use the different results to discuss and quantify uncertainties associated with tillage erosion measurements. The study was performed in northeastern Germany on a 10 m by 50 m plot with a mean slope of 8%. Tillage erosion was determined after two sequences of seven tillage operations. Two different micro-tracers (magnetic iron oxide mixed with soil and fluorescent sand) and one macro-tracer (passive radio-frequency identification transponders (RFIDs), size: 4 x 22 mm) were used to directly determine soil fluxes. Moreover, tillage induced changes in topography were measured for the entire plot with two different terrestrial laser scanners and an unmanned aerial system for structure from motion topography analysis. Based on these elevation differences, corresponding soil fluxes were calculated. The mean translocation distance of all techniques was 0.57 m per tillage pass, with a relatively wide range of mean soil translocation distances ranging from 039 to 0.72 m per pass. A benchmark technique could not be identified as all used techniques have individual error sources, which could not be quantified. However, the translocation distances of the macro-tracers used were consistently smaller than the translocation distances of the micro-tracers (mean difference = -26 +/- 12%), which questions the widely used assumption of non-selective soil transport via tillage operations. This study points out that tillage erosion measurements, carried out under almost optimal conditions, are subject to major uncertainties that are far from negligible. (C) 2018 Elsevier B.V. All rights reserved.
Subsurface lateral flow in hillslope soils depends on lower permeability or texture-contrasting soil horizons. In the arable hummocky soil landscape, erosion processes caused glacial till appearance closer to the soil surface at upslope positions. The objective of this work was to quantify the potential for subsurface lateral flow depending on the erosion-affected spatial hydropedological complexity. The eroded Haplic Luvisol profile was studied due to the presence of a relatively dense C horizon that varied in depth, thickness, and sloping angle. A two-dimensional numerical modeling and sensitivity analysis for the saturated hydraulic conductivity (K-s) of the C horizon and the depth to C horizon (i.e., soil solum thickness) was performed for rainstorms in 2011 and 2012 using HYDRUS-2D. A K-s value of <2.5 cm d(-1) for the C horizon was required for lateral flow initiation. Lateral flow was (i) increasing with decreasing solum thickness, indicating an erosion-induced feedback on subsurface lateral flow, and (ii) dependent on the soil moisture regime prior to rainstorms. The effect of lateral flow on the movement of a conservative tracer was simulated in the form of a "virtual experiment". Simulation scenarios revealed only a relatively small lateral shift of the tracer plume caused by a local decoupling of water (already lateral) from subsequent tracer movement (still vertical). Longer term simulations suggested that, for the present conditions, lateral flow was limited mostly to occasional summer storm events. Even without considering preferential flow contribution to lateral flow, highly complex hydropedologic interactions are present in erosion-affected heterogeneous soil profiles.
Whilst correlations have been made between the loess of Europe and China, deposits in Central Asia have remained largely overlooked by scientific investigation. The nature of the relationship between loess accumulation and palaeoclimate at the core of the Eurasian loess belt is particularly poorly understood. Here we reconstruct palaeoenvironmental change in Central Asia over the last 40 ky based on data from the Remizovka loess profile, in the northern foothills of the Tian Shan in southern Kazakhstan. Our interpretations are based on synthesis of chronostratigraphic, colour and magnetic susceptibility data, supported by chronostratigraphies from two additional sites nearby, Maibulak and Valikhanova. All three sites record substantially increased loess accumulation during late MIS 3 into the global last glacial maximum (gLGM). At Remizovka, increased loess flux occurred in two pulses at c. 38-25 ka and 22-18 ka, with the intervening period involving incipient pedogenesis. At Maibulak, two loess pulses at c. 40-30 ka and c. 28-22 ka are separated by a weakly developed paleosol which may date to the same time as pedogenesis at Remizovka. There is additional possible periglacial influence at Maibulak from c. 40-33.5 ka. At Valikhanova, there is some age overlap between paleosol and loess samples, but overall loess accumulation appears to have increased at c. 42-35 ka, c. 30 ka and the gLGM, with pedogenesis occurring >40 ka and c. 32 ka. At all three sites, Holocene loess accumulation is minimal; this period is characterised by pedogenesis. The chronostratigraphic variability between our sites highlights a need to interrogate climate-driven models for loess formation in piedmont environments. We interpret our data in the context of regional palaeoenvironmental archives to indicate that loess accumulation increased coeval with MIS 3 glacial advance in the Tian Shan, which was facilitated by northward expansion of the Asian monsoon and associated increase in precipitation. We hypothesise that increased ice volume impeded teleconnections with the temperate zone westerlies to the north; these were compressed against the piedmont resulting in increased wind strength and facilitating increased loess flux. Peak loess accumulation during the gLGM occurred under colder, drier climatic conditions, with reduced but sustained glacial ice volume and persistent influence of the westerlies in the arid Central Asian piedmont loess belt. In the absence of more widespread, reliably dated palaeoenvironmental records from the region, our data become of critical importance for understanding past environmental conditions in Central Asia, relative to elsewhere in Eurasia and globally.
Towards an understanding of climate proxy formation in the Chew Bahir basin, southern Ethiopian Rift
(2018)
Deciphering paleoclimate from lake sediments is a challenge due to the complex relationship between climate parameters and sediment composition. Here we show the links between potassium (K) concentrations in the sediments of the Chew Bahir basin in the Southern Ethiopian Rift and fluctuations in the catchment precipitation/evaporation balance. Our micro-X-ray fluorescence and X-ray diffraction results suggest that the most likely process linking climate with potassium concentrations is the authigenic illitization of smectites during episodes of higher alkalinity and salinity in the closed -basin lake, due to a drier climate. Whole-rock and clay size fraction analyses suggest that illitization of the Chew Bahir clay minerals with increasing evaporation is enhanced by octahedral Al-to-Mg substitution in the clay minerals, with the resulting layer charge increase facilitating potassium-fixation. Linking mineralogy with geochemistry shows the links between hydroclimatic control, process and formation of the Chew Bahir K patterns, in the context of well-known and widely documented eastern African climate fluctuations over the last 45,000 years. These results indicate characteristic mineral alteration patterns associated with orbitally controlled wet-dry cycles such as the African Humid Period (similar to 15-5 ka) or high-latitude controlled climate events such as the Younger Dryas (similar to 12.8-11.6 ka) chronozone. Determining the impact of authigenic mineral alteration on the Chew Bahir records enables the interpretation of the previously established pXRF-derived aridity proxy K and provides a better paleohydrological understanding of complex climate proxy formation.
The precise determination of the stable C and O isotope fractionation between water and calcite (CC) and water and aragonite (AR) is of special interest for climate reconstructions, e.g. paleotemperatures. Previous studies reported results from both laboratory and field experiments, but their results are only partly consistent. Here we present C and O isotope data of a stalagmite from the Swiss Alps, which shows CC-AR transitions along individual growth layers. Using detailed analyses both laterally and perpendicular to such layers we examined the difference in the C and O isotope fractionation factor of the HCO3- - CC and the HCO3- - AR system. For O this difference is similar to the water-CC and water-AR offset provided in experimental studies. The O isotope fractionation difference in the water-CC and water-AR system is comparable to those determined in laboratory studies but shows a statistically significant correlation with the CaCO3 precipitation rate. For C we found a fractionation difference, which is independent of CaCO3 precipitation rate and with slightly smaller values for the fractionation offset between HCO3- - CC and HCO3- - AR compared to literature values. However, we also found an unexpected decrease in delta C-13 along growth layers, which contradicts the widely used concept of Rayleigh fractionation during CO2 degassing and CaCO3 precipitation. The results of this study can be used e.g., to correct stable isotope time series of stalagmites showing CC-AR transitions along their growth axes. (C) 2018 Elsevier Ltd. All rights reserved.
Modelers can improve a model by addressing the causes for the model errors (data errors and structural errors). This leads to implementing model enhancements (MEs), for example, meteorological data based on more monitoring stations, improved calibration data, and/or modifications in process formulations. However, deciding on which MEs to implement remains a matter of expert knowledge. After implementing multiple MEs, any improvement in model performance is not easily attributed, especially when considering different objectives or aspects of this improvement (e.g., better dynamics vs. reduced bias). We present an approach for comparing the effect of multiple MEs based on real observations and considering multiple objectives (MMEMO). A stepwise selection approach and structured plots help to address the multidimensionality of the problem. Tailored analyses allow a differentiated view on the effect of MEs and their interactions. MMEMO is applied to a case study employing the mesoscale hydro-sedimentological model WASA-SED for the Mediterranean-mountainous Isabena catchment, northeast Spain. The investigated seven MEs show diverse effects: some MEs (e.g., rainfall data) cause improvements for most objectives, while other MEs (e.g., land use data) only affect a few objectives or even decrease model performance. Interaction of MEs was observed for roughly half of the MEs, confirming the need to address them in the analysis. Calibration and increasing the temporal resolution showed by far stronger impact than any of the other MEs. The proposed framework can be adopted in other studies to analyze the effect of MEs and, thus, facilitate the identification and implementation of the most promising MEs for comparable cases.
A comprehensive hydro-sedimentological dataset for the Isabena catchment, northeastern (NE) Spain, for the period 2010-2018 is presented to analyse water and sediment fluxes in a Mediterranean mesoscale catchment. The dataset includes rainfall data from 12 rain gauges distributed within the study area complemented by meteorological data of 12 official meteo-stations. It comprises discharge data derived from water stage measurements as well as suspended sediment concentrations (SSCs) at six gauging stations of the River Isabena and its sub-catchments. Soil spectroscopic data from 351 suspended sediment samples and 152 soil samples were collected to characterize sediment source regions and sediment properties via fingerprinting analyses. The Isabena catchment (445 km(2)) is located in the southern central Pyrenees ranging from 450 m to 2720 m a.s.l.; together with a pronounced topography, this leads to distinct temperature and precipitation gradients. The River Isabena shows marked discharge variations and high sediment yields causing severe siltation problems in the downstream Barasona Reservoir. The main sediment source is badland areas located on Eocene marls that are well connected to the river network. The dataset features a comprehensive set of variables in a high spatial and temporal resolution suitable for the advanced process understanding of water and sediment fluxes, their origin and connectivity and sediment budgeting and for the evaluation and further development of hydro-sedimentological models in Mediterranean mesoscale mountainous catchments.
Ice-rich yedoma-dominated landscapes store considerable amounts of organic carbon (C) and nitrogen (N) and are vulnerable to degradation under climate warming. We investigate the C and N pools in two thermokarst-affected yedoma landscapes - on Sobo-Sise Island and on Bykovsky Peninsula in the north of eastern Siberia. Soil cores up to 3m depth were collected along geomorphic gradients and analysed for organic C and N contents. A high vertical sampling density in the profiles allowed the calculation of C and N stocks for short soil column intervals and enhanced understanding of within-core parameter variability. Profile-level C and N stocks were scaled to the landscape level based on landform classifications from 5 m resolution, multispectral RapidEye satellite imagery. Mean landscape C and N storage in the first metre of soil for Sobo-Sise Island is estimated to be 20.2 kg C m(-2) and 1.8 kg N m(-2) and for Bykovsky Peninsula 25.9 kg C m(-2) and 2.2 kg N m(-2). Radiocarbon dating demonstrates the Holocene age of thermokarst basin deposits but also suggests the presence of thick Holoceneage cover layers which can reach up to 2 m on top of intact yedoma landforms. Reconstructed sedimentation rates of 0.10-0.57 mm yr(-1) suggest sustained mineral soil accumulation across all investigated landforms. Both yedoma and thermokarst landforms are characterized by limited accumulation of organic soil layers (peat). We further estimate that an active layer deepening of about 100 cm will increase organic C availability in a seasonally thawed state in the two study areas by similar to 5.8 Tg (13.2 kg C m(-2)). Our study demonstrates the importance of increasing the number of C and N storage inventories in ice-rich yedoma and thermokarst environments in order to account for high variability of permafrost and thermokarst environments in pan-permafrost soil C and N pool estimates.
Thermal conductivity (lambda) is an essential physical property of minerals and rocks and fundamental in constraining the thermal field of the lithosphere. In case that adequate samples to measure lambda are not available, it could be indirectly inferred from calculation. One of the most widely applied indirect methods for rocks involve modal mineralogy and porosity as parameters that are incorporated into mathematical mean or mixing models. Robust inferences from these approaches for crystalline rocks were impeded by a small number of studied samples or restriction to certain rock types. We employ this method and examine its applicability to low-porosity plutonic rocks by calculating bulk thermal conductivity lambda(b) for 45 samples covering the entire range from gabbro/diorite to granite. We show that the use of the harmonic-mean model for both rock matrix and porosity provided a good match between lambda(b.meas) and lambda(b.calc) of <10% deviation (2 sigma), with relative and absolute errors amounting to 1.49.7% and 4.44.9%, respectively. The results of our study constitute a big step forward to a robust conclusion on the overall applicability of the harmonic-mean model for inferring lambda(b) of isotropic, low-porosity, mafic to silicic plutonic and metamorphic rocks with an acceptable magnitude of error. Drill cuttings and enclaves form particularly interesting objects for application of this method, as they are poorly suited for direct measurement. Well-derived lambda values for those rocks would permit to calculate heat flow and to model more profoundly the thermal state of the deeper lithosphere.
Rainfall data from the Global Precipitation Measurement (GPM) mission provide a new source of information with high spatiotemporal resolution that overcomes the limitations of ground-based rainfall information worldwide. This study evaluates the performance of the Integrated multi-satellitE Retrievals for GPM (IMERG) Final Run product over Brazil by means of multi-temporal and -spatial analyses. The assessment of the IMERG Final Run product is based on six statistics obtained for the period between January-December 2016 (daily, monthly, and annual basis). The analysis consisted of comparing the satellite-based estimates against a ground-based gridded rainfall product created using daily records from 4911 rain gauges distributed throughout Brazil. Overall, the results show that the IMERG product can effectively capture the spatial patterns of rainfall across Brazil. However, the IMERG product presents a slight tendency in overestimating the ground-based rainfall at all timescales. Furthermore, the performance of the satellite product varies throughout the region. The higher errors and biases are found in the North and Central-West regions, but the low density of rain gauges in those regions can be a source of large deviations between IMERG estimates and observations. A large underestimation of the IMERG data is evident along the coastal zone of the North-east region, probably due to the inability of the passive microwave and infrared sensors to detect warm-rain processes over land. This study shows that the IMERG product can be a good source of rainfall data to complement the ground precipitation measurements in most of Brazil, although some uncertainties are found and need to be further studied
Annual greenhouse gas emissions have increased more than threefold between 1950 and 2014, posing a major threat to the integrity of the entire earth system and subsequently to humankind. Consequently, roadmaps towards low-carbon pathways are urgently needed. Our study contributes to a more detailed understanding of the dynamics of country based emission patterns and uses them to discuss prospective low-carbon pathways for countries. As availability of databases on sectoral emissions substantially increased, we employ machine learning techniques to classify emission features and pathways. By doing so, 18 representative emission patterns are derived. Overall emissions from seven sectors and for 167 countries covering the time span from 1950 to 2014 have been used in the analyses. The following significant trends can be observed: a) increasing per capita emissions due to growing fossil fuel use in many parts of the world, b) a decline in per capita emissions in some countries, and c) a shift in the emission shares, i.e., a reduction of agricultural and land use contributions in certain regions. Using the emission patterns, their dynamics, and best performing countries as role models, we show the possibility for gaining a decent human development without significantly increasing per capita emissions.
Like almost all fields of science, hydrology has benefited to a large extent from the tremendous improvements in scientific instruments that are able to collect long-time data series and an increase in available computational power and storage capabilities over the last decades. Many model applications and statistical analyses (e.g., extreme value analysis) are based on these time series. Consequently, the quality and the completeness of these time series are essential. Preprocessing of raw data sets by filling data gaps is thus a necessary procedure. Several interpolation techniques with different complexity are available ranging from rather simple to extremely challenging approaches. In this paper, various imputation methods available to the hydrological researchers are reviewed with regard to their suitability for filling gaps in the context of solving hydrological questions. The methodological approaches include arithmetic mean imputation, principal component analysis, regression-based methods and multiple imputation methods. In particular, autoregressive conditional heteroscedasticity (ARCH) models which originate from finance and econometrics will be discussed regarding their applicability to data series characterized by non-constant volatility and heteroscedasticity in hydrological contexts. The review shows that methodological advances driven by other fields of research bear relevance for a more intensive use of these methods in hydrology. Up to now, the hydrological community has paid little attention to the imputation ability of time series models in general and ARCH models in particular.
A potential human footprint on Western Central African rainforests before the Common Era has become the focus of an ongoing controversy. Between 3,000 y ago and 2,000 y ago, regional pollen sequences indicate a replacement of mature rainforests by a forest-savannah mosaic including pioneer trees. Although some studies suggested an anthropogenic influence on this forest fragmentation, current interpretations based on pollen data attribute the "rainforest crisis" to climate change toward a drier, more seasonal climate. A rigorous test of this hypothesis, however, requires climate proxies independent of vegetation changes. Here we resolve this controversy through a continuous 10,500-y record of both vegetation and hydrological changes from Lake Barombi in Southwest Cameroon based on changes in carbon and hydrogen isotope compositions of plant waxes. delta C-13-inferred vegetation changes confirm a prominent and abrupt appearance of C-4 plants in the Lake Barombi catchment, at 2,600 calendar years before AD 1950 (cal y BP), followed by an equally sudden return to rainforest vegetation at 2,020 cal y BP. delta D values from the same plant wax compounds, however, show no simultaneous hydrological change. Based on the combination of these data with a comprehensive regional archaeological database we provide evidence that humans triggered the rainforest fragmentation 2,600 y ago. Our findings suggest that technological developments, including agricultural practices and iron metallurgy, possibly related to the large-scale Bantu expansion, significantly impacted the ecosystems before the Common Era.
An essential, respected, and critical aspect of the modern practice of science and scientific publishing is peer review. The process of peer review facilitates best practices in scientific conduct and communication, ensuring that manuscripts published as accurate, valuable, and clearly communicated. The over 152 papers published in Tectonics in 2017 benefit from the time, effort, and expertise of our reviewers who have provided thoughtfully considered advice on each manuscript. This role is critical to advancing our understanding of the evolution of the continents and their margins, as these reviews lead to even clearer and higher-quality papers. In 2017, the over 423 papers submitted to Tectonics were the beneficiaries of more than 786 reviews provided by 562 members of the tectonics community and related disciplines. To everyone who has volunteered their time and intellect to peer reviewing, thank you for helping Tectonics and all other AGU Publications provide the best science possible.
Global climate change is one of the greatest challenges of the 21st century, with influence on the environment, societies, politics and economies. The (semi-)arid areas of Southern Africa already suffer from water scarcity. There is a great variety of ongoing research related to global climate history but important questions on regional differences still exist.
In southern African regions terrestrial climate archives are rare, which makes paleoclimate studies challenging. Based on the assumption that continental pans (sabkhas) represent a suitable geo-archive for the climate history, two different pans were studied in the southern and western Kalahari Desert. A combined approach of molecular biological and biogeochemical analyses is utilized to investigate the diversity and abundance of microorganisms and to trace temporal and spatial changes in paleoprecipitation in arid environments. The present PhD thesis demonstrates the applicability of pan sediments as a late Quaternary geo-archive based on microbial signature lipid biomarkers, such as archaeol, branched and isoprenoid glycerol dialkyl glycerol tetraethers (GDGTs) as well as phospholipid fatty acids (PLFA). The microbial signatures contained in the sediment provide information on the current or past microbial community from the Last Glacial Maximum to the recent epoch, the Holocene. The results are discussed in the context of regional climate evolution in southwestern Africa. The seasonal shift of the Innertropical Convergence Zone (ITCZ) along the equator influences the distribution of precipitation- and climate zones. The different expansion of the winter- and summer rainfall zones in southern Africa was confirmed by the frequency of certain microbial biomarkers. A period of increased precipitation in the south-western Kalahari could be described as a result of the extension of the winter rainfall zone during the last glacial maximum (21 ± 2 ka). Instead a period of increased paleoprecipitation in the western Kalahari was indicated during the Late Glacial to Holocene transition. This was possibly caused by a southwestern shift in the position of the summer rainfall zone associated to the southward movement of the ITCZ.
Furthermore, for the first time this study characterizes the bacterial and archaeal life based on 16S rRNA gene high-throughput sequencing in continental pan sediments and provides an insight into the recent microbial community structure. Near-surface processes play an important role for the modern microbial ecosystem in the pans. Water availability as well as salinity might determine the abundance and composition of the microbial communities. The microbial community of pan sediments is dominated by halophilic and dry-adapted archaea and bacteria. Frequently occurring microorganisms such as, Halobacteriaceae, Bacillus and Gemmatimonadetes are described in more detail in this study.
More than 41% of the Earth’s land area is covered by permanent or seasonally arid dryland ecosystems. Global development and human activity have led to an increase in aridity, resulting in ecosystem degradation and desertification around the world. The objective of the present work was to investigate and compare the microbial community structure and geochemical characteristics of two geographically distinct saline pan sediments in the Kalahari Desert of southern Africa. Our data suggest that these microbial communities have been shaped by geochemical drivers, including water content, salinity, and the supply of organic matter. Using Illumina 16S rRNA gene sequencing, this study provides new insights into the diversity of bacteria and archaea in semi-arid, saline, and low-carbon environments. Many of the observed taxa are halophilic and adapted to water-limiting conditions. The analysis reveals a high relative abundance of halophilic archaea (primarily Halobacteria), and the bacterial diversity is marked by an abundance of Gemmatimonadetes and spore-forming Firmicutes. In the deeper, anoxic layers, candidate division MSBL1, and acetogenic bacteria (Acetothermia) are abundant. Together, the taxonomic information and geochemical data suggest that acetogenesis could be a prevalent form of metabolism in the deep layers of a saline pan.
More than 41% of the Earth’s land area is covered by permanent or seasonally arid dryland ecosystems. Global development and human activity have led to an increase in aridity, resulting in ecosystem degradation and desertification around the world. The objective of the present work was to investigate and compare the microbial community structure and geochemical characteristics of two geographically distinct saline pan sediments in the Kalahari Desert of southern Africa. Our data suggest that these microbial communities have been shaped by geochemical drivers, including water content, salinity, and the supply of organic matter. Using Illumina 16S rRNA gene sequencing, this study provides new insights into the diversity of bacteria and archaea in semi-arid, saline, and low-carbon environments. Many of the observed taxa are halophilic and adapted to water-limiting conditions. The analysis reveals a high relative abundance of halophilic archaea (primarily Halobacteria), and the bacterial diversity is marked by an abundance of Gemmatimonadetes and spore-forming Firmicutes. In the deeper, anoxic layers, candidate division MSBL1, and acetogenic bacteria (Acetothermia) are abundant. Together, the taxonomic information and geochemical data suggest that acetogenesis could be a prevalent form of metabolism in the deep layers of a saline pan.
Due to a lack of well-preserved terrestrial climate archives, paleoclimate studies are sparse in southwestern Africa. Because there are no perennial lacustrine systems in this region, this study relies on a saline pan as an archive for climate information in the western Kalahari (Namibia). Molecular biological and biogeochemical analyses were combined to examine the response of indigenous microbial communities to modern and past climate-induced environmental conditions. The 16S rRNA gene high-throughput sequencing was applied to sediment samples from Omongwa pan to characterize the modern microbial diversity. Highest diversity of microorganisms, dominated by the extreme halophilic archaeon Halobacteria and by the bacterial phylum Gemmatimonadetes, was detected in the near-surface sediments of Omongwa pan. In deeper sections abundance and diversity significantly decreases and Bacillus, known to form spores, become dominant. Lipid biomarkers for living and past microbial life were analyzed to track the influence of climate variation on the abundance of microbial communities from the Last Glacial Maximum to Holocene time. Since water is an inevitable requirement for microbial life, in this dry region the abundance of past microbial biomarkers was evaluated to conclude on periods of increased paleoprecipitation in the past. The data point to a period of increased humidity in the western Kalahari during the Last Glacial to Holocene transition indicating a southward shift of the Intertropical Convergence Zone during this period. Comparison with results from a southwestern Kalahari pan suggests complex displacements of the regional atmospheric systems since the Last Glacial Maximum.
Abstract. The aim of this study is to investigate the shallow thermal field differences for two differently aged passive continental margins by analyzing regional variations in geothermal gradient and exploring the controlling factors for these variations. Hence, we analyzed two previously published 3-D conductive and lithospheric-scale thermal models of the Southwest African and the Norwegian passive margins. These 3-D models differentiate various sedimentary, crustal, and mantle units and integrate different geophysical data such as seismic observations and the gravity field. We extracted the temperature–depth distributions in 1 km intervals down to 6 km below the upper thermal boundary condition. The geothermal gradient was then calculated for these intervals between the upper thermal boundary condition and the respective depth levels (1, 2, 3, 4, 5, and 6 km below the upper thermal boundary condition). According to our results, the geothermal gradient decreases with increasing depth and shows varying lateral trends and values for these two different margins. We compare the 3-D geological structural models and the geothermal gradient variations for both thermal models and show how radiogenic heat production, sediment insulating effect, and thermal lithosphere–asthenosphere boundary (LAB) depth influence the shallow thermal field pattern. The results indicate an ongoing process of oceanic mantle cooling at the young Norwegian margin compared with the old SW African passive margin that seems to be thermally equilibrated in the present day.
Abstract. The aim of this study is to investigate the shallow thermal field differences for two differently aged passive continental margins by analyzing regional variations in geothermal gradient and exploring the controlling factors for these variations. Hence, we analyzed two previously published 3-D conductive and lithospheric-scale thermal models of the Southwest African and the Norwegian passive margins. These 3-D models differentiate various sedimentary, crustal, and mantle units and integrate different geophysical data such as seismic observations and the gravity field. We extracted the temperature–depth distributions in 1 km intervals down to 6 km below the upper thermal boundary condition. The geothermal gradient was then calculated for these intervals between the upper thermal boundary condition and the respective depth levels (1, 2, 3, 4, 5, and 6 km below the upper thermal boundary condition). According to our results, the geothermal gradient decreases with increasing depth and shows varying lateral trends and values for these two different margins. We compare the 3-D geological structural models and the geothermal gradient variations for both thermal models and show how radiogenic heat production, sediment insulating effect, and thermal lithosphere–asthenosphere boundary (LAB) depth influence the shallow thermal field pattern. The results indicate an ongoing process of oceanic mantle cooling at the young Norwegian margin compared with the old SW African passive margin that seems to be thermally equilibrated in the present day.